Курсовая работа: Исследование метода продолжения решения по параметру для нелинейных САУ

if (m < 4)

dt=dt*2;

end;

t=dt+t;

it=it+1;

end;

disp ('it='); %количество итераций

disp (it);

disp ('t='); %время выполнения итерационного процесса

disp (t);

pause;

xout - конечное значение x;

dxout - конечное значение вектора ошибки.

m - номер итерации.

Prog. m

trace=1;

dt=0.1;

x0=0;

edop=0.1;

[xout,dxout,m] = mpr ('funf','funj',x0,dt,edop,trace);

plot (m,xout); %График значений x

pause;

plot (m,dxout); %График ошибки

pause;

Funf. m

function [f] =funf (x)

f= [0.0001*exp (30*x) +x-6];

end

Funj. m