Курсовая работа: Проектирование системы автоматического управления
x1=0;
i=0;
As=0.0125*(x.^3)+0.3*(x.^2)+4.886*x+61.72;
x=x-(As.^2)./((0.0125*((x+As).^3)+0.3*((x+As).^2)+4.886*(x+As)+61.72)+As);
As=0.0125*(x.^3)+0.3*(x.^2)+4.886*x+61.72;
A(1)=x;
i=i+1;
while abs(x-x1)>e
x1=x;
x=x-(As.^2)./((0.0125*((x+As).^3)+0.3*((x+As).^2)+4.886*(x+As)+61.72)+As);
As=0.0125*(x.^3)+0.3*(x.^2)+4.886*x+61.72;
A(i+1)=x;
i=i+1;
end
plot(1,A(1));
hold on
for n=1:i
plot(n,A(n),'b-o')
end
grid on
xlabel('iteraciya')
ylabel('roots')
disp('ответ');
disp(x);
disp('число итераций');
disp(i);
Программа 2.
Решение дифференциального уравнения численным способом.
clc
a2=24;