用matlab求非线性规划和常微分方程组

Question

11．求下列非线性规划，初值为（0.5，0.5,0.5）

12.求常微分方程组的解。
范围[0 10]

Answer 1

11

%fun1.m文件
 function f=fun1(x)
f=sum(x.^2)+8;
end

%fun2.m文件
function [g,h]=fun2(x)
g=[-x(1)^2+x(2)-x(3)^2,x(1)+x(2)^2+x(3)^2-20];
h=[x(2)+2*x(3)^2-3,-x(1)-x(2)^2+2];
end

命令行输入

clear
x0=[0.5,0.5,0.5];    %初值
lb=[0,0,0];
ub=[inf,inf,inf];
A=[];b=[];Aeq=[];beq=[];
[x,fval]=fmincon('fun1',x0,A,b,Aeq,beq,lb,ub,'fun2');

22

[x,y]=ode45(@cal,[0 1],[0;1;1]);
plot(x,y(:,1));
 
function res=cal(~,y)
    res = zeros(3,1);
    res(1)=y(2)*y(3);
    res(2)=-y(1)*y(3);
    res(3)=-0.5*y(1)*y(2);
end