24小时热门版块排行榜

返回列表

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

我是xingxing

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 782.6
散金: 17
红花: 2
帖子: 142
在线: 26.5小时
虫号: 4435625
注册: 2016-02-25
性别: MM
专业: 电化学

[求助] 10参数非线性方程组求助已有2人参与

a+d+2*e+2*f+2*g+2*p=0.05
b+d+2*e+2*g+h+m+n+p=0.025
d/(a*b*c)=10^9.55
e/((a^2)*(b^2))=10^14.43
f*c/((a^2)*(b^2))=10^10.85
g*(c^2)/((a^2)*(b^2))=10^5.87
p*c/(a^2*b)=10^4.92
m*c/h=10^(-2.92)
n*c/m=10^(-4.34)
b*c/n=10^(-5.69)

C值从3变到9 间隔为0.1 。每取一个c值，求出对应的其他值。‘

有大神可以通过编程帮忙解决吗？不胜感激。’

[ Last edited by 独孤神宇 on 2018-7-19 at 15:23 ]

回复此楼

» 猜你喜欢

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

1楼 2018-07-18 11:36:19

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

我是xingxing

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 782.6
散金: 17
红花: 2
帖子: 142
在线: 26.5小时
虫号: 4435625
注册: 2016-02-25
性别: MM
专业: 电化学

引用回帖:

2楼: Originally posted by 独孤神宇 at 2018-07-18 22:21:32
初始值太难试了，可以考虑用 1stOpt ，这还是只是一个参数的 c=3的情况。

function example_18_7_18
c=3;
x0=; %% 设定初始值
options=optimset('MaxFunEvals',20000,'MaxIter',2000);
=fsolve(@f,x0,o ...

这个求出来全部正值吗？

发自小木虫Android客户端

赞一下

回复此楼

3楼2018-07-18 23:46:30

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

查看全部 6 个回答

独孤神宇

版主 (知名作家)

应助: 490 (硕士)
贵宾: 0.008
金币: 31014.8
散金: 802
红花: 122
沙发: 1
帖子: 5600
在线: 855.5小时
虫号: 3522474
注册: 2014-11-06
性别: GG
专业: 机械动力学
管辖: 计算模拟

【答案】应助回帖

感谢参与，应助指数 +1

初始值太难试了，可以考虑用 1stOpt ，这还是只是一个参数的 c=3的情况。

function example_18_7_18
c=3;
x0=[1;1;1;1;1;1;1;0.1;0.1;0.1]; %% 设定初始值
options=optimset('MaxFunEvals',20000,'MaxIter',2000);
[x,fval,exitflag]=fsolve(@f,x0,options)
function F=f(x) %%% 待求方程组
F1=x(1)+x(3)+2*x(4)+2*x(5)+2*x(6)+2*x(10)-0.05;
F2=x(2)+x(3)+2*x(4)+2*x(6)+x(7)+x(8)+x(9)+x(10)-0.025;
F3=x(3)/(x(1)*x(2)*c)-10^9.55;
F4=x(4)/((x(1)^2)*(x(2)^2))-10^14.43;
F5=x(5)*c/((x(1)^2)*(x(2)^2))-10^10.85;
F6=x(6)*(c^2)/((x(1)^2)*(x(2)^2))-10^5.87;
F7=x(10)*c/(x(1)^2*x(2))-10^4.92;
F8=x(8)*c/x(7)-10^(-2.92);
F9=x(9)*c/x(8)-10^(-4.34);
F10=x(2)*c/x(9)-10^(-5.69);
F=[F1;F2;F3;F4;F5;F6;F7;F8;F9;F10];
end
end

赞一下

回复此楼

数值计算