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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 非线性方程组的解法
61/1返回列表
查看: 2168  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

zouji1983

金虫 (正式写手)

[求助] 非线性方程组的解法

如果用fortran解下面两个方程组：
1. 5+3.2ln(x)+ln(1-x)-2x+5x^2+6x^3-3y=0;
2. 0.3x+0.7y=0.9
方程中既有ln，又有3次幂。
应该怎么编程呢？我是用别人的程序，现在需要加入这个。请高手指点，多谢！
回复此楼

» 猜你喜欢

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

zouji1983

金虫 (正式写手)
或者，如果用matlab的solve函数，怎么转化成fortran呢？
赞一下 回复此楼
2楼2012-11-21 21:25:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jingjingdou

银虫 (小有名气)

【答案】应助回帖

将上述两个方程化为：
y=(5+3.2ln(x)+ln(1-x)-2x+5x^2+6x^3)/3;
x=(0.9-0.7y)/0.3
然后将(x0,y0)赋初值采用上述式子迭代计算，得到的本次结果(xn,yn)作为下一次结果的初值。循环迭代至前后两次的结果满足精度要求，可以取误差的二范数或者无穷范数等作为判断。
该方法称为迭代法，具体的fortran算法，可以参考gauss-赛德尔迭代算法。网上有很多这样的程序，也可以参见徐士良的《fortran常用算法程序集》这本书。
以上仅介绍计算的主要过程，愿对你有帮助。
不足之处，请见谅……
赞一下(1人) 回复此楼
顺一点，一点就好！
3楼2012-12-22 09:03:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

11jxli1

禁虫 (著名写手)
本帖内容被屏蔽
4楼2012-12-22 21:22:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

615780236

禁虫 (知名作家)
本帖内容被屏蔽
5楼2012-12-24 18:36:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

573245006

至尊木虫 (正式写手)

【答案】应助回帖

解题思路：先对上述方程组第2式进行变形，把1式中y用x代替
最终得：x^22.4=exp(5x-35x^2-42x^3-8)/((1-x))^7
用迭代法VB程序如下：
Dim x, x1 As Single
x = 0.1
For i = 1 To 1000
x1 = ((Exp(5 * x - 35 * x ^ 2 - 42 * x ^ 3 - 8)) / ((1 - x) ^ 7)) ^ (1 / 22.4)
x = x1
Next i
Print x
解得：x=0.3031055
则 y=1.155811929(x带入方程组2式手算得)
赞一下 回复此楼
6楼2012-12-26 17:52:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zouji1983 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 生物学调剂招人！！！ +3 山海天岚 2026-03-17 4/200 2026-03-19 21:34 by 怎么释怀
[考研] 一志愿武汉理工材料工程专硕调剂 +5 Doleres 2026-03-19 5/250 2026-03-19 20:14 by 制度的
[考研] 一志愿天津大学化学工艺专业（081702）315分求调剂 +11 yangfz 2026-03-17 11/550 2026-03-19 15:06 by houyaoxu
[考研] 材料与化工求调剂 +7 为学666 2026-03-16 7/350 2026-03-19 14:48 by 尽舜尧1
[考研] 304求调剂 +3 曼殊2266 2026-03-18 3/150 2026-03-19 14:42 by peike
[考研] 材料080500调剂求收留 +4 一颗meteor 2026-03-13 4/200 2026-03-19 10:32 by 30660438
[考研] 332求调剂 +3 ydfyh 2026-03-17 3/150 2026-03-19 10:14 by 功夫疯狂
[考研] 材料专业求调剂 +5 hanamiko 2026-03-18 5/250 2026-03-18 20:19 by 星空星月
[考研] 一志愿武理材料305分求调剂 +5 想上岸的鲤鱼 2026-03-18 6/300 2026-03-18 17:53 by 无际的草原
[考研] 化学工程321分求调剂 +15 大米饭！ 2026-03-15 18/900 2026-03-18 14:52 by haxia
[考研] 材料与化工一志愿南昌大学327求调剂推荐 +8 Ncdx123456 2026-03-13 9/450 2026-03-18 14:40 by haxia
[考研] 331求调剂（0703有机化学 +7 ZY-05 2026-03-13 8/400 2026-03-18 14:13 by 007_lilei
[考研] 280求调剂 +6 咕噜晓晓 2026-03-18 7/350 2026-03-18 11:25 by 无际的草原
[考博] 环境领域全国重点实验室招收博士1-2名 +3 QGZDSYS 2026-03-13 5/250 2026-03-18 11:13 by QGZDSYS
[考研] 293求调剂 +11 zjl的号 2026-03-16 16/800 2026-03-18 08:10 by zhukairuo
[考研] 一志愿南京大学，080500材料科学与工程，调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 070300化学学硕求调剂 +6 太想进步了0608 2026-03-16 6/300 2026-03-16 16:13 by kykm678
[考研] 22408总分284求调剂 +3 InAspic 2026-03-13 3/150 2026-03-15 11:10 by zhq0425
[考研] 材料与化工 323 英一+数二+物化，一志愿：哈工大 本人本科双一流 +4 自由的_飞翔 2026-03-13 5/250 2026-03-14 19:39 by hmn_wj
[硕博家园] 085600 260分求调剂 +3 天空还下雨么 2026-03-13 5/250 2026-03-13 18:46 by 天空还下雨么
信息提示
关闭
请填处理意见
关闭