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