版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

fy640323

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 21
帖子: 2
在线: 1.2小时
虫号: 2425480
注册: 2013-04-19
专业: 能源化工

[求助] matlab求解高阶方程组

方程组如下：

CODE:

F=[x(1)*(1/(x(4)*x(5)*(x(2)+1)*exp(x(5)*1)))*((x(3)+x(4)*(exp(x(5)*1)))^(x(2)+1))-0.029785;

x(1)*(1/(x(4)*x(5)*(x(2)+1)*exp(x(5)*2)))*((x(3)+x(4)*(exp(x(5)*2)))^(x(2)+1))-0.06402;

x(1)*(1/(x(4)*x(5)*(x(2)+1)*exp(x(5)*3)))*((x(3)+x(4)*(exp(x(5)*3)))^(x(2)+1))-0.090489;

x(1)*(1/(x(4)*x(5)*(x(2)+1)*exp(x(5)*4)))*((x(3)+x(4)*(exp(x(5)*4)))^(x(2)+1))-0.111586;

x(1)*(1/(x(4)*x(5)*(x(2)+1)*exp(x(5)*5)))*((x(3)+x(4)*(exp(x(5)*5)))^(x(2)+1))-0.135672];

首先用sove函数解了一下：

CODE:

[k,n,A,B,y]=solve( ...,

'k*(1/(A*B*(n+1)*exp(B*1)))*((y+A*(exp(B*1)))^(n+1))==0.029785', ...,'k*(1/(A*B*(n+1)*exp(B*2)))*((y+A*(exp(B*2)))^(n+1))==0.06402', ...,'k*(1/(A*B*(n+1)*exp(B*3)))*((y+A*(exp(B*3)))^(n+1))==0.090489', ...,'k*(1/(A*B*(n+1)*exp(B*4)))*((y+A*(exp(B*4)))^(n+1))==0.111586', ...,'k*(1/(A*B*(n+1)*exp(B*5)))*((y+A*(exp(B*5)))^(n+1))==0.135672');

有解，但是带回方程验证发现每一个等式的左边都无限于零，和右边是不等的，而且更改方程组右边的常数对方程组的解没有影响，这显然不对，这是为什么？

然后又尝试用fsolve：（其中的x1 x2 x3 x4 x5分别于k n y A B对应）

CODE:

function F=test(x)

F=[(x(1)/910/(0.434*10^(-10)))*(1/(x(4)*x(5)*(x(2)+1)*exp(x(5)*1)))*((x(3)+x(4)*(exp(x(5)*1)))^(x(2)+1))-0.029785;

(x(1)/910/(0.434*10^(-10)))*(1/(x(4)*x(5)*(x(2)+1)*exp(x(5)*2)))*((x(3)+x(4)*(exp(x(5)*2)))^(x(2)+1))-0.06402;

(x(1)/910/(0.434*10^(-10)))*(1/(x(4)*x(5)*(x(2)+1)*exp(x(5)*3)))*((x(3)+x(4)*(exp(x(5)*3)))^(x(2)+1))-0.090489;

(x(1)/910/(0.434*10^(-10)))*(1/(x(4)*x(5)*(x(2)+1)*exp(x(5)*4)))*((x(3)+x(4)*(exp(x(5)*4)))^(x(2)+1))-0.111586;

(x(1)/910/(0.434*10^(-10)))*(1/(x(4)*x(5)*(x(2)+1)*exp(x(5)*5)))*((x(3)+x(4)*(exp(x(5)*5)))^(x(2)+1))-0.135672];

复制代码命令是：

CODE:

op=optimset('MaxIter',5000)

[x,fval]=fsolve(@test,[0.02,0.02,0.02,0.02,0.02],op);

复制代码更改初值，不是提示迭代次数不够，就是不收敛，由于是新模型，各参数范围不清楚

请各位帮忙看看，怎么解，本人学化工的，需要用实验数据检验模型是否正确，现在是用五组数据五个方程解出其中的系数，谢谢了！或者大家还有什么其他的方法，能保证方程左边的精度在0.001就可以了

回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

MATLAB求解非线性方程组已经有5人回复
Matlab求解二阶偏微分方程组，希望能给出相应的求解代码已经有10人回复
请问这个4阶的非线性偏微分方程组在PDE中怎么解？等待高手解答。已经有17人回复
求Matlab高手解决线性方程组的迭代求解问题已经有17人回复
matlab求解非线性方程组，错误提示怎么解决已经有5人回复
微分方程解，解的稳定性及相关分析已经有7人回复
matlab求解偏微分方程组已经有3人回复
求助 MATLAB解方程组-fslove 已经有7人回复
求助matlab---fsolve解非线性方程组已经有6人回复
matlab求解非线性方程组，求助！已经有6人回复
matlab有限差分法求解复杂的偏微分方程组已经有19人回复
matlab微分方程组参数拟合，以周为单位求解，汇总后以年为单位进行数值比较已经有4人回复
Matlab求解偏微分方程组已经有22人回复
matlab求解非线性方程组已经有16人回复
matlab数值求解边界条件微分方程组已经有7人回复
求matlab求解一个方程组的问题已经有6人回复
求高人指点用matlab求解非线性方程组，解决了追加100金币；已经有11人回复
matlab的fsove 命令求解非线性方程组已经有6人回复
MATLAB非线性方程组求解~~~反应釜模拟已经有9人回复
【求助】MATLAB 解方程组【已解决】已经有7人回复
【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题已经有21人回复
【求助】matlab怎么求解偏微分方程组啊，先谢谢了已经有13人回复
【求助】求教matlab解非线性方程组已经有9人回复

1楼 2013-04-22 22:59:15

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

dingd

铁杆木虫 (职业作家)

应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.7小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

感谢参与，应助指数 +1

内容已删除

赞一下

回复此楼

2楼2013-04-23 08:53:52

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

fy640323

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 21
帖子: 2
在线: 1.2小时
虫号: 2425480
注册: 2013-04-19
专业: 能源化工

引用回帖:

2楼: Originally posted by dingd at 2013-04-23 08:53:52
不唯一：
1：
x1: 6.63491531705323E-10
x2: 0.170494561773157
x3: -5172.65986018108
x4: 3722.18334080518
x5: 0.98435674669011

2：
x1: 1.41600788120601E-9
x2: 0.170494570410365
x3: -60.62931 ...

这个解差不多，请问是怎么做的呢？可以把m文件给我发一下吗？谢谢了，毕业设计急着用

赞一下

回复此楼

3楼2013-04-23 09:50:41

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

dingd

铁杆木虫 (职业作家)

应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.7小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
fy640323: 金币+10, ★★★很有帮助, 很感谢，用你说的方法，虽然方程解出来了，但是解和你给的有点差别，但是用参数估计把问题解决了，非常感谢， 2013-04-23 20:11:37

自己用1stOpt算吧，简单方便。

赞一下

回复此楼

4楼2013-04-23 11:52:17

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

相关版块跳转我要订阅楼主 fy640323 的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[基金申请] 生命口面上会评结束了吗大家有知道结果吗？现在还没有结果正常不？ +4	碧水00 2026-07-16	6/300	2026-07-16 08:11 by gltch
[基金申请] 生生命学部会评 +4	wuke100666 2026-07-15	4/200	2026-07-16 00:12 by 啦啦小妮子
[有机交流] 伯胺上甲基求助 100+4	002tmumu 2026-07-14	4/200	2026-07-15 22:01 by 学员oN9JQD
[基金申请] H口，面上，时间戳7月7日。 +3	plumage5 2026-07-15	3/150	2026-07-15 21:06 by 6543yes
[基金申请] 小木虫没落了，除了祈祷帖子，几乎看不到有价值的帖子 +13	botar 2026-07-14	14/700	2026-07-15 19:28 by 寒流32
[基金申请] 会评 +5	无名者登山 2026-07-15	5/250	2026-07-15 18:15 by 杠就是你对
[基金申请] E口会评 +11	布布和一二 2026-07-13	13/650	2026-07-15 16:34 by Kittylucky
[基金申请] E口青C会评结束了吗？ 40+3	小小书虫c 2026-07-14	9/450	2026-07-15 16:04 by 安静的知了
[基金申请] 生命口C13面上会评时间 +6	luckinging 2026-07-14	9/450	2026-07-15 14:39 by 雨冰共舞
[教师之家] 内心匮乏 +13	水冰月月野兔 2026-07-12	13/650	2026-07-15 14:21 by jiduquegai
[文学芳草园] 你可曾有过这种感受？ +3	Jeanya 2026-07-09	6/300	2026-07-15 11:07 by mbtwt
[论文投稿] 农业机械学报投稿周期 +3	G0DLIN 2026-07-14	7/350	2026-07-15 09:26 by xs74101122
[基金申请] 国自然面上D口祈祷 +4	monkeynana 2026-07-14	4/200	2026-07-15 09:14 by 白水煮
[论文投稿] 跨出版社商投稿 20+4	倪好520 2026-07-13	4/200	2026-07-15 08:35 by bobvan
[基金申请] 没收到消息，再次陪跑 +3	天空星辰fly 2026-07-14	3/150	2026-07-14 19:14 by 6543yes
[基金申请] 不要再数国自然申请书的 filecode 的分隔符个数了 +13	sadpencil 2026-07-12	20/1000	2026-07-14 13:12 by 杠就是你对
[基金申请] 明天E口面上会评 +8	布布和一二 2026-07-12	9/450	2026-07-13 22:59 by QTSorange
[基金申请] 祈祷青基必中 50+6	hadengry 2026-07-09	16/800	2026-07-12 20:31 by 登山虫虫
[考博] 27届辽宁大学应届毕业生申博 +3	可乐微风大海 2026-07-09	3/150	2026-07-10 21:23 by 贪玩的研究僧
[基金申请] 关于如何从代码看上不上会 +17	且听虎啸 2026-07-09	23/1150	2026-07-10 19:51 by sdfapple719