24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

自己的歌

银虫 (初入文坛)

应助: 0 (幼儿园)
金币: 747.5
红花: 1
帖子: 30
在线: 18小时
虫号: 1211560
注册: 2011-02-24
性别: GG
专业: 工程热物理相关交叉领域

[求助] 求助：非线性方程组的求解（郁闷中）已有10人参与

各位大侠好，本人工科出身，数学功底实在一般。最近自编一个计算程序，涉及到一个非线性方程组的求解，无奈解法不理想，很多时候不收敛，求大侠指导一二。

方程形式：方程组中的每个方程的形式都是这样的，x + A = f1(x) + f2(x) ,其中f(x)的形式为 f(x) = (x+ B) / ln (x + C)
其中A，B，C为常数。

解法：采用不动点迭代法，即假设一组初值，带入方程的右边，从而得到一组新的值。如误差值大于允许误差，采用加权因子的方式获得新的迭代值，加权因子从0.1到0.625已尝试过多个。

问题：有时迭代过程中变量计算值超过边界条件。例如x的允许范围为8<= x <=12 , 迭代过程中x会超过12或小于8，会导致计算出错。因此我限定如果x超过12，则等于12；小于8，则等于8。但是没有效果，最后x值总是一直超过12就是小于8。费解。

Help：各位大侠对于这个不动点解法，有没有优化的方法，或者有没有其它可行的解法？谢谢先！

回复此楼

» 猜你喜欢

【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
物理学I论文润色/翻译怎么收费? 已经有152人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有23人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复

1楼 2013-12-22 21:45:03

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

专家经验: +518
应助: 942 (博后)
贵宾: 1.275
金币: 3430
散金: 58785
红花: 532
沙发: 11
帖子: 24215
在线: 2601.8小时
虫号: 2139575
注册: 2012-11-21
专业: 光学信息获取与处理
管辖: 数学

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
自己的歌: 金币+5, ★★★很有帮助 2013-12-23 15:28:51
fegg7502: 金币+2, 3ks 2013-12-26 09:05:20

用1stop, 代码如下，各个constant 需要先赋值。
Parameters  x;
Constant  a1= , a2= , a3= , a4= , a5= , a6= ;
Constant  b1= , b2= , b3= , b4= , b5= , b6= ;
Constant  c1= , c2= , c3= , c4= , c5= , c6= ;
Constant  d1= , d2= , d3= , d4= , d5= , d6= ;
Constant  e1= , e2= , e3= , e4= , e5= , e6= ;
Constant  f1= , f2= , f3= , f4= , f5= , f6= ;
Function z1=(a3-x)/ln(a4-x)+ (a5-x)/ln(a6-x)- a1*(x-a2);
Function  z2=(b3-x)/ln(b4-x)+ (b5-x)/ln(b6-x)- b1*(x-b2);
Function  z3=(c3-x)/ln(c4-x)+ (c5-x)/ln(c6-x)- c1*(x-c2);
Function  z4=(d3-x)/ln(d4-x)+ (d5-x)/ln(d6-x)- d1*(x-d2);
Function  z5=(e3-x)/ln(e4-x)+ (e5-x)/ln(e6-x)- e1*(x-e2);
Function  z6=(f3-x)/ln(f4-x)+ (f5-x)/ln(f6-x)- f1*(x-f2);
MinFunction z1^2+z2^2+z3^2+ z4^2+z5^2+z6^2;