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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 计算数学 » matlab求解高阶方程组
41/1返回列表
查看: 1538  |  回复: 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一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-16 03:23 by DgNGHc3h5tPl
[找工作] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 3n8v2C8RimXI 2026-09-13 3/150 2026-09-16 02:38 by DgNGHc3h5tPl
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-16 02:14 by DgNGHc3h5tPl
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 20:59 by zNcEhTcH7ihq
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +6 两三岁ss 2026-09-14 6/300 2026-09-15 20:17 by 只吃肉葡萄
[考研] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 20:14 by zNcEhTcH7ihq
[公派出国] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 19:59 by CgyNCDVNhGVg
[硕博家园] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 40ms4Wlo4umh 2026-09-14 3/150 2026-09-15 18:59 by CgyNCDVNhGVg
[考研] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 17:23 by CgyNCDVNhGVg
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:35 by CgyNCDVNhGVg
[考博] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +7 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:23 by CgyNCDVNhGVg
[找工作] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 15:02 by CgyNCDVNhGVg
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 09:38 by VMDqgnlKimtZ
[教师之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 07:50 by C79jjtjAKjEn
[教师之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 07:02 by 5BDX0d0WFp7t
[博后之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:50 by 5BDX0d0WFp7t
[考研] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-14 22:25 by vZfe6xYu34yj
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 21:37 by vZfe6xYu34yj
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:25 by vZfe6xYu34yj
[考博] 西北工业大学材料学院电化学传感与催化课题组招收2027年推免研究生、博士生、博士后 +3 马方园 2026-09-10 4/200 2026-09-14 09:29 by Chem张zz
信息提示
关闭
请填处理意见
关闭