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

%%minf(x)=(2*x1^2-2*x1*x2+2x2^2+x1+x2-3)/(3*x1+4)%%s.t.3*x1+x2-3=0functionfyp=li1(t,x... %%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时候的所有解
展开
 我来答
chenbin6699
2014-09-01
知道答主
回答量:12
采纳率:0%
帮助的人:16.2万
展开全部
通过这个子函数可以求出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
追问
运行不出来
追答
如果可以的话,你把你的全部代码贴出来看看。
富港检测技术(东莞)有限公司_
2024-04-02 广告
正弦振动多用于找出产品设计或包装设计的脆弱点。看在哪一个具体频率点响应最大(共振点);正弦振动在任一瞬间只包含一种频率的振动,而随机振动在任一瞬间包含频谱范围内的各种频率的振动。由于随机振动包含频谱内所有的频率,所以样品上的共振点会同时激发... 点击进入详情页
本回答由富港检测技术(东莞)有限公司_提供
推荐律师服务: 若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询

为你推荐:

下载百度知道APP,抢鲜体验
使用百度知道APP,立即抢鲜体验。你的手机镜头里或许有别人想知道的答案。
扫描二维码下载
×

类别

我们会通过消息、邮箱等方式尽快将举报结果通知您。

说明

0/200

提交
取消

辅 助

模 式