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

登录注册

应《网络安全法》要求，自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就可以了

回复此楼

» 猜你喜欢

球磨粉体时遇到了大的问题，请指教！已经有9人回复
江汉大学解明教授课题组招博士研究生/博士后已经有3人回复

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

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

返回列表