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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 计算数学 » matlab求解高阶方程组
41/1返回列表
查看: 1540  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

fy640323

新虫 (初入文坛)

[求助] 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就可以了
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-04-22 22:59:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
内容已删除
赞一下 回复此楼
2楼2013-04-23 08:53:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fy640323

新虫 (初入文坛)
引用回帖:
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

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
fy640323: 金币+10, ★★★很有帮助, 很感谢，用你说的方法，虽然方程解出来了，但是解和你给的有点差别，但是用参数估计把问题解决了，非常感谢， 2013-04-23 20:11:37
自己用1stOpt算吧，简单方便。
赞一下 回复此楼
4楼2013-04-23 11:52:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fy640323 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +7 3n8v2C8RimXI 2026-09-13 7/350 2026-09-16 15:12 by i3KBRD19ERkG
[找工作] 浙江师范大学是怎么坑我的 +7 recruit123 2026-09-09 7/350 2026-09-16 15:07 by 防腐剂科技部
[考博] 上海工程技术大学激光智能制造课题组｜2027级博士研究生招生公告 +7 水士口 2026-09-16 8/400 2026-09-16 13:47 by kkklens
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-16 08:51 by SjQobnC04j06
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 Aj1rhIDL5ixY 2026-09-14 4/200 2026-09-16 07:48 by i0IGWqQuchto
[考研] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-16 06:11 by LH5NkK5Ud4bh
[找工作] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-16 05:15 by LH5NkK5Ud4bh
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 6F5UbRU2I5hL 2026-09-14 5/250 2026-09-16 05:11 by LH5NkK5Ud4bh
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-16 04:47 by LH5NkK5Ud4bh
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-13 7/350 2026-09-16 03:26 by DgNGHc3h5tPl
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 3n8v2C8RimXI 2026-09-13 5/250 2026-09-16 03:02 by DgNGHc3h5tPl
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +7 QUjhNVAcOSff 2026-09-13 8/400 2026-09-16 02:23 by DgNGHc3h5tPl
[考研] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +6 QUjhNVAcOSff 2026-09-13 6/300 2026-09-16 02:02 by DgNGHc3h5tPl
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 20:26 by zNcEhTcH7ihq
[公派出国] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 19:59 by CgyNCDVNhGVg
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 23jxep3nCNZb 2026-09-14 4/200 2026-09-15 18:11 by CgyNCDVNhGVg
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:59 by CgyNCDVNhGVg
[公派出国] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +7 QUjhNVAcOSff 2026-09-13 9/450 2026-09-15 14:10 by VMDqgnlKimtZ
[教师之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:46 by C79jjtjAKjEn
[教师之家] 浙江师范大学是怎么坑我的 +10 recruit123 2026-09-09 15/750 2026-09-13 10:45 by seaskyy
信息提示
关闭
请填处理意见
关闭