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

登录注册

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

24小时热门版块排行榜

>基金申请 (8260)
>论坛更新日志 (5109)
>休闲灌水 (1471)
>文献求助 (392)
>导师招生 (372)
>虫友互识 (333)
>论文投稿 (324)
>有奖起名 (165)
>硕博家园 (161)
>考博 (125)
>健康生活 (89)
>博后之家 (47)
>考研 (45)
>留学DIY (39)
>绿色求助(高悬赏) (37)
>教师之家 (34)

返回列表

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就可以了

回复此楼

» 猜你喜欢

有多少人是今天查系统知道结果的？已经有16人回复
我不理解！已经有14人回复
国自然面上复盘~欢迎讨论已经有4人回复
系统查不到已经有10人回复
基金不中，共勉已经有3人回复
怎么查啊已经有5人回复
范进中举一文的中心思想已经有10人回复
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战已经有9人回复
2026国自然函评费到账已经有25人回复
2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？已经有28人回复

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

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 的主题更新

返回列表

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

[基金申请] 我不理解！ +9	Edward_pc 2026-08-26	14/700	2026-08-26 18:03 by hsding2026
[基金申请] 国自然面上复盘~欢迎讨论 +3	晴天加油 2026-08-26	4/200	2026-08-26 17:59 by Tide man
[基金申请] 系统查不到 +10	董八千 2026-08-26	10/500	2026-08-26 16:30 by Equinoxhua
[基金申请] 怎么查啊 +5	huang1991js 2026-08-26	5/250	2026-08-26 16:09 by Equinoxhua
[基金申请] 出来了 +9	trojank 2026-08-26	9/450	2026-08-26 14:25 by 宝贝虫子
[基金申请] 只有每年这种时候来逛逛小木虫 +26	yaoyewhu2008 2026-08-20	28/1400	2026-08-26 13:49 by zzuzxg
[基金申请] 哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识 +4	yuleib84 2026-08-26	5/250	2026-08-26 13:40 by yuleib84
[基金申请] 为什么国自然不能直接公布 +4	bjdxyxy 2026-08-26	4/200	2026-08-26 13:12 by qingmu1201
[基金申请] 国合可查了 +3	paperzjh 2026-08-26	3/150	2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了，明天出结果吗 +19	angus9576 2026-08-25	23/1150	2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗？ +10	chengyan1220 2026-08-24	20/1000	2026-08-26 08:57 by peasantsprig
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+10) +6	ziyangfang 2026-08-25	9/450	2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了！？ +6	chengyan1220 2026-08-23	6/300	2026-08-25 19:45 by zfd97
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +17	医学老男孩 2026-08-20	22/1100	2026-08-25 15:36 by 医学老男孩
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8	scalable 2026-08-24	8/400	2026-08-25 12:52 by jnhyjjm
[基金申请] 人气不行了 +11	fansofjerry 2026-08-21	11/550	2026-08-25 11:04 by 孤独的英雄6
[教师之家] 跳槽后在研项目怎么办？ +5	简单化xn 2026-08-22	10/500	2026-08-23 12:38 by 简单化xn
[基金申请] 看来今天不会放榜了？ +8	chengyan1220 2026-08-21	11/550	2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13	wuchongjun 2026-08-20	19/950	2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧？ +4	哈哈蛤？ 2026-08-21	4/200	2026-08-21 10:58 by Vivilian