24小时热门版块排行榜

返回列表

【悬赏金币】回答本帖问题，作者Claire3741将赠送您 5 个金币

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

Claire3741

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 58.6
帖子: 28
在线: 4.2小时
虫号: 23386947
注册: 2020-08-18
专业: 电力电子学

[求助] 求解高次方程组（5次，matlab&1stopt）已有1人参与

我需要解一个五元方程组，首先使用matlab中的fsolve，我所设置的迭代次数为10000次，但是每次跑到400次就停了，而且还会报“solver stopped prematurely.fsolve stopped because it exceeded the iteration limit, options.MaxIterations=400(the default value).”我不知道原因。后面考虑到fsolve初值选取比较难，考虑用1stopt去解，仍然解不出来，我想着没有精确解，有个看得过去的近似解也行啊。但是始终没解。有没有人能帮忙看看为什么前面两种方法解不出来，或者是提供一些解答的新思路啊，谢谢大家了！
1、matlab代码
x = sym('x', [1,2,3,4,5],'positive');
% 生成符号变量向量
w1=sqrt((x(1)+x(3)*x(4))./(2*x(2)*x(4)*x(5)));
z=((2*x(2)*x(4)+x(3)*x(4)*x(5))./(2*(x(1)+x(3)*x(4))))*w1;
w2=sqrt(1-z^2)*w1;
a=sqrt((1-2*x(5)*z*w1+x(5)^2*w1^2)./(1-z^2));
f1=atan((w2*x(5))./(1-z*w1*x(5)))-atan((sqrt(1-z^2)./(-z)));
% 参数组成的表达式
eq1=-0.5*x(4)*a./(x(1)+x(3)*x(4))+0.7959;
eq2=z*w1-1.275;
eq3=w2-2.284;
eq4=f1+0.4204;
eq5=-0.5*x(4)./(x(1)+x(3)*x(4))+0.3504;
f = matlabFunction([eq1;eq2;eq3;eq4;eq5], 'vars', {[x(1) x(2) x(3) x(4) x(5)]});
x0 = [0.5 2.5 1 0.05 0.2];
options=optimset('Display','iter','MaxFunEvals',10000,'TolFun',1e-10);
[x,fval,exitflag,output] = fsolve(f,x0,options);
2、1stopt
const delta=-0.5;
parameter x1[0,5],x2[0,5],x3[0,5],x4[0,1],x5[0,5];
conststr
w1=sqrt((x1+x3*x4)/(2*x2*x4*x5));
z=((2*x2*x4+x3*x4*x5)/(2*(x1+x3*x4)))*w1;
w2=sqrt(1-z^2)*w1;
a=sqrt((1-2*x5*z*w1+x5^2*w1^2)/(1-z^2));
f1=atan((w2*x5)/(1-z*w1*x5))-atan(sqrt(1-z^2)./(-z));
//参数组成的表达式
A=delta*x4*a/(x1+x3*x4)+0.7959;
B=z*w1-1.275;
C=w2-2.284;
E=f1+0.4204;
F=delta*x4/(x1+x3*x4)+0.3504；

matlab出错图

回复此楼

» 猜你喜欢

【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
物理学I论文润色/翻译怎么收费? 已经有58人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有23人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有0人回复

1楼 2020-08-22 20:49:39

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

Claire3741

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 58.6
帖子: 28
在线: 4.2小时
虫号: 23386947
注册: 2020-08-18
专业: 电力电子学

引用回帖:

5楼: Originally posted by 独孤神宇 at 2020-08-25 08:29:27
5个参数，需要正版软件才能计算。

Constant delta=-0.5;
Parameter x1=,x2=,x3=,x4=,x5=;
ConstStr
w1=sqrt((x1+x3*x4)/(2*x2*x4*x5)),
z=((2*x2*x4+x3*x4*x5)/(2*(x1+x3*x4)))*w1,
w2=sqrt(1-z^2)*w1,
a ...

谢谢你的回复，我又重做了实验，用fsolve还是没办法解出来，我没1stopt正版软件，能不能麻烦你帮忙跑一下啊。我这个其实是6个参数，5个方程，我想弄两组试看看能不能接出来。
1、6个参数，5个方程
Constant delta=-0.225；
Parameter x1=[0,5],x2=[0,5],x3=[0,5],x4=[0,5],x5=[0,5],x6=[0,5];
//生成变量K,F,T,H,D,R
ConstStr
w1=sqrt((x1+x5*x6)/(2*x3*x4*x6)),
z=((2*x6*x4+(x5*x6+x1*x2)*x3)/(2*(x1+x5*x6)))*w1,
w2=sqrt(1-z^2)*w1,
aa=sqrt((1-2*x3*z*w1+x3^2*w1^2)./(1-z^2)),
f1=atan((w2*x3)./(1-z*w1*x3))-atan((sqrt(1-z^2)./(-z)));
//参数组成的表达式
Function
delta*x6*aa/(x1+x5*x6)+0.02436=0;
z*w1-1.135=0;
w2-2.3=0;
f1+0.5595=0;
delta*x6./(x1+x5*x6)+0.01119=0;
2、5个参数，5个方程
Constant delta=-0.225；
Parameter x1=[0,5],x2=[0,5],x3=[0,5],x4=[0,5],x5=[0,5];
//生成变量K,F,T,H,D
ConstStr
w1=sqrt((x1+x5*0.05)/(2*x3*x4*0.05)),
z=((2*0.05*x4+(x5*0.05+x1*x2)*x3)/(2*(x1+x5*0.05)))*w1,
w2=sqrt(1-z^2)*w1,
aa=sqrt((1-2*x3*z*w1+x3^2*w1^2)./(1-z^2)),
f1=atan((w2*x3)./(1-z*w1*x3))-atan((sqrt(1-z^2)./(-z)));
//参数组成的表达式
Function
delta*0.05*aa/(x1+x5*0.05)+0.02436=0;
z*w1-1.135=0;
w2-2.3=0;
f1+0.5595=0;
delta*0.05./(x1+x5*0.05)+0.01119=0;
十分感谢！