Контрольный расчет к курсовой работе в среде Wxmaxima

fpprintprec:5; tn:10;tk:35;t1:22.5;a:12;b:12;N:11;dt:(tk-tn)/(N-1),numer;

// определение начальных условий для расчета Uvx(t)

Uvx(t):=if t<t1 then a*(t-tn) else a*(t1-tn)-b*(t-t1); // определение значений Uvx(t)

for t:tn thru tk step dt do display(Uvx(t)),numer; // вывод значений на экран морнитора

Uvx1:5;U1:20;Uvx2:25;U2:150; // определение начальных условий для преобразования графика передаточной характеристики в аналитическое выражение

Для определения Uvix(t) можно воспользоваться формулой 3 или использовать более простой матричный способ.

A: matrix([5,1], [25,1]); B: matrix([20], [150]); x:invert(A).B,numer;

Uvix(t):=if Uvx(t)<Uvx1 then U1 else if Uvx(t)<=Uvx2 then 6.5*Uvx(t)-12.5 else U2,numer;

for t:tn thru tk step dt do display(Uvix(t)),numer;

wxplot2d([Uvx(t)], [t,tn,tk],[gnuplot_preamble, "set grid;"])

wxplot2d([Uvix(t)], [t,tn,tk],[y,0,160], [gnuplot_preamble, "set grid;"])