Добавил:
lightcorpz
nikolozzz15@gmail.com
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:fmsfunction3
.m function J=fmsfunction3(Q)
global u;
u = Q;
OPTIONS = odeset('RelTol', 100000);
[t,x] = ode23s('ode3',[0 4],[1 0 0],OPTIONS);
n=length(t);
for k = 1:n
U(k)=u(1)*exp(u(2)*t(k)).*(cos(u(3)*t(k)+u(4)));
end
J=x(n,3)+20*(x(n,1)^2);
figure(2)
plot(t,x(:,1),'--',t,U)
legend('x(t)','U(t)')
grid on;
end
global u;
u = Q;
OPTIONS = odeset('RelTol', 100000);
[t,x] = ode23s('ode3',[0 4],[1 0 0],OPTIONS);
n=length(t);
for k = 1:n
U(k)=u(1)*exp(u(2)*t(k)).*(cos(u(3)*t(k)+u(4)));
end
J=x(n,3)+20*(x(n,1)^2);
figure(2)
plot(t,x(:,1),'--',t,U)
legend('x(t)','U(t)')
grid on;
end
Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]