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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 非线性方程组的系数的数量级差很大，结果发散
51/1返回列表
查看: 1588  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

RK1224

新虫 (初入文坛)

[求助] 非线性方程组的系数的数量级差很大，结果发散 已有3人参与

小生最近做点东西，模型简化为非线性方程组，输入一组参数后，发现方程的系数数量级相差很大，用牛顿迭代结果发散。请问各位大神有遇到过这种情况吗？有什么好的方法吗？
简化后的方程代入一组参数如下：
(4.1018e-07)*x1^4+38716e*x1-38708*x2-(9.2333e-08)*x3^4-297.7296)=0;

38708*x1-38943*x2+1174.2e=0;

-23679*x2+47828e*x3-24149e*x4=0;

(9.2333e-08)*x1^4+23914e*x2-47834*x3-(9.2333e-08)*x3^4+23914*x4 -5.3514*x5+4613.8=0;

5.3514*x3-83.6514*x5+1566=0;
回复此楼

» 收录本帖的淘帖专辑推荐

matlab编程绘图

» 猜你喜欢
1楼 2015-08-10 19:50:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

引用回帖:
6楼: Originally posted by RK1224 at 2015-08-11 13:32:07
谢谢你的回复，确实是我没写清楚，我改了一下如下：
(4.1018e-07)*x1^4+38716*x1-38708*x2-(9.2333e-08)*x3^4-297.7296=0;

38708*x1-38943*x2+1174.2=0;

-23679*x2+47828*x3-24149*x4=0;

(9.2333e-08)*x ...

MATLAB代码：
CODE:
function solve_eq
clear all;clc
format long

x0=[0 0 0 0 0];

options=optimset('TolFun',1e-9);
[x fval]=fsolve(@fun,x0,options);
[x;fval]'

function f=fun(u)



x1=u(1);
x2=u(2);
x3=u(3);
x4=u(4);
x5=u(5);

f(1)=(4.1018e-07)*x1^4+38716*x1-38708*x2-(9.2333e-08)*x3^4-297.7296;

f(2)=38708*x1-38943*x2+1174.2;

f(3)=-23679*x2+47828*x3-24149*x4;

f(4)=(9.2333e-08)*x1^4+23914*x2-47834*x3-(9.2333e-08)*x3^4+23914*x4 -5.3514*x5+4613.8;

f(5)=5.3514*x3-83.6514*x5+1566;

计算结果：
           x1~x5                         方程残差
   6.063569995329909     -0.000000000031662
   6.057131381229750     -0.000000000011596
  15.543130816422794     -0.000000000116415
  24.844508953237408      0.000000000042746
  19.714882360020336  -0.000000000000227
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
8楼2015-08-11 14:34:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

hsdwlxhx

金虫 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1
5维迭代，你也太bug了

[ 发自小木虫客户端 ]
赞一下 回复此楼
BringLightonChip
2楼2015-08-10 20:09:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

RK1224

新虫 (初入文坛)
引用回帖:
2楼: Originally posted by hsdwlxhx at 2015-08-10 20:09:31
5维迭代，你也太bug了

我是菜鸟。。。希望你指点一下
赞一下 回复此楼
3楼2015-08-10 20:43:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hsdwlxhx

金虫 (正式写手)

【答案】应助回帖

引用回帖:
3楼: Originally posted by RK1224 at 2015-08-10 20:43:33
我是菜鸟。。。希望你指点一下...

指点谈不上，你的方程可以简化，先用数值方法跑一跑，看看有没有解，大概在什么区间，然后再用迭代计算，matlab有fsolve函数可以直接用

[ 发自小木虫客户端 ]
赞一下 回复此楼
BringLightonChip
4楼2015-08-10 22:02:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 生物学308分求调剂（一志愿华东师大）做过分子实验 +9 相信必会光芒万� 2026-04-07 10/500 2026-04-13 10:20 by 可淡不可忘
[考研] 271求调剂 +23 2261744733 2026-04-11 26/1300 2026-04-13 09:59 by szhize
[考研] 085801电气专硕272求调剂 +9 电气李 2026-04-13 10/500 2026-04-13 09:31 by fqwang
[考研] 327求调剂 +14 Xxjc1107. 2026-04-06 14/700 2026-04-13 08:31 by Delta2012
[考研] 313求调剂 +5 余摆摆是大哥 2026-04-11 5/250 2026-04-12 10:34 by 猪会飞
[考研] 305求调剂 +7 玛卡巴卡boom 2026-04-11 7/350 2026-04-12 07:35 by zhouxiaoyu
[考研] 347求调剂 +4 mhyqyy 2026-04-06 4/200 2026-04-12 02:27 by 秋豆菜芽
[考研] 303求调剂 +14 SereinQ 2026-04-10 15/750 2026-04-11 20:43 by 蓝云思雨
[考研] 求调剂 +3 胃痉挛累了 2026-04-11 5/250 2026-04-11 14:13 by luhong1990
[考研] 300分求调剂 （085501机械专硕，本科扬大） +8 xu@841019 2026-04-11 8/400 2026-04-11 10:46 by qingpingzhu
[考研] 087100初试311求调剂 +4 任雅琴 2026-04-09 4/200 2026-04-11 10:33 by zhq0425
[考研] 一志愿985机械学硕380求调剂 +5 关关雎鸠10 2026-04-11 5/250 2026-04-11 10:10 by 知念。A
[考研] 346，工科求调剂 +3 moser233 2026-04-09 3/150 2026-04-11 10:04 by zhq0425
[考研] 311求调剂 +13 xyp想读书 2026-04-10 14/700 2026-04-11 09:41 by 猪会飞
[考研] 本科211 工科085400 280分求调剂 可跨专业 +11 LZH（等待调剂中 2026-04-10 11/550 2026-04-11 08:39 by zhq0425
[考研] 计算机类求调剂，22408-274分 +7 上岸de小虫 2026-04-09 8/400 2026-04-10 19:56 by fxue1114
[考研] 284求调剂 +9 让我上岸吧阿西 2026-04-09 11/550 2026-04-10 19:18 by 靖jing
[考研] 考研调剂 +13 冰冰，，， 2026-04-07 13/650 2026-04-09 17:01 by Lilly_Li
[考研] 求调剂 +3 猪肉墩粉条cc 2026-04-08 4/200 2026-04-09 10:05 by 猪肉墩粉条cc
[考研] 材料调剂 +14 一样YWY 2026-04-06 14/700 2026-04-08 23:00 by 猪会飞
信息提示
关闭
请填处理意见
关闭