Question

如何用matlab求解微分方程组 我 要 求的是让模型=0时的解（即平衡点）求各位大神帮帮忙

%%min f(x)=(2*x1^2-2*x1*x2+2x2^2+x1+x2-3)/(3*x1+4)
%%s.t. 3*x1+x2-3=0
function fyp=li1(t,x)
N=2;
e0=2;
%%f

f1=(6*x(1)^2+16*x(1)-6*x(2)^2-11*x(2)+13)/(3*x(1)+4)^2;
f2=(-2*x(1)+4*x(2)+1)/(3*x(1)+4);
f=[f1;f2];

A=[3,1];
b=3;
B1=A'*inv(A*A')*A;
I=[1,0;
0,1];
B=(I-B1);
C=A'*A;
D=3*A';
e=1/((t+1)^(1/3));
fyp=-B*(e*f)-C*x+D;
我要问的是 我的模型是-B*(e*f)-C*x+D 我要令-B*(e*f)-C*x+D=0 求出模型=0时候的所有解

Answer 1

通过这个子函数可以求出fyp的所有值。然后在主程序中添加一个选择fyp中为0的元素，另外需要考虑求解精度，fyp中的符合条件的元素未必能完全等于0，误差范围在10e-6，差不多就行了。
在主程序中的代码大致如此
n=length(fyp);
p=1;
for i=1:n
if abs(fyp(i))<10e-6;
Answer(p)=fyp(i);%选择符合条件的值
p=p+1;
end

追问

运行不出来

追答

如果可以的话，你把你的全部代码贴出来看看。