/* [wxMaxima batch file version 1] [ DO NOT EDIT BY HAND! ]*/ /* [ Created with wxMaxima version 16.04.1 ] */ /* [wxMaxima: title start ] Dimmer [wxMaxima: title end ] */ /* [wxMaxima: comment start ] Normalized voltage [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ m: sin(x); /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] RMS [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ Mrms: sqrt(1 / %pi * integrate(m^2, x, alpha, %pi)); /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ Mrms: ratsimp(Mrms); /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] First Harmonic [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ M1s: 1 / %pi * (integrate(m * sin(x), x, alpha, %pi) + integrate(m * sin(x), x, %pi + alpha, 2 * %pi)); /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ M1s: ratsimp(M1s); /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ ev(M1s, alpha = 0); /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ M1c: 1 / %pi * (integrate(m * cos(x), x, alpha, %pi) + integrate(m * cos(x), x, %pi + alpha, 2 * %pi)); /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ ev(M1c, alpha = 0); /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ M1: sqrt(M1s^2 + M1c^2); /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ M1: ratsimp(M1); /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ M1: trigsimp(M1); /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ M1rms: M1 / sqrt(2); /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] Power [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ P: Vm^2 / R * Mrms^2; /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ ev(P, alpha = 0); /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ wxplot2d(P * R / Vm^2, [alpha, 0, %pi-1e-3], [box, false]); /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] Apparent Power [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ S: Vm / sqrt(2) * Vm / R * Mrms; /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ wxplot2d([P * R / Vm^2, S * R / Vm^2], [alpha, 0, %pi-1e-3], [box, false], [legend, false]); /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] Power Factor [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ PF: P / S; /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ wxplot2d(PF, [alpha, 0, %pi-1e-3], [box, false], [legend, false]); /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] Displacement Power Factor [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ DPF: M1s / M1; /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ wxplot2d(DPF, [alpha, 0, %pi-1e-3], [box, false], [legend, false]); /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] verify PF-DPF relation [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ wxplot2d([PF, M1rms / Mrms * DPF], [alpha, 0, %pi-1e-3], [box, false], [legend, false]); /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] THD [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ THD: sqrt(Mrms^2 / (M1^2 / 2) - 1) * 100; /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ wxplot2d(THD, [alpha, 0, %pi-1e-3], [y, 0, 300], [box, false], [legend, false]); /* [wxMaxima: input end ] */ /* Maxima can't load/batch files which end with a comment! */ "Created with wxMaxima"$