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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 非线性方程组的解法
61/1返回列表
查看: 2146  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 伙伴们，祝我生日快乐吧 +15 myrtle 2026-03-10 24/1200 2026-03-15 21:16 by 苏州_逗号
[考研] 0703化学调剂，求各位老师收留 +7 秋有木北 2026-03-14 7/350 2026-03-15 17:30 by 小物理化学
[考研] 274求调剂 +4 时间点 2026-03-13 4/200 2026-03-15 15:29 by Rambo13
[考研] 255求调剂 +3 李嘉慧， 2026-03-12 4/200 2026-03-14 16:58 by 有只狸奴
[考研] 337一志愿华南理工材料求调剂（有希望2吗？） +3 mysdl 2026-03-09 3/150 2026-03-14 02:53 by JourneyLucky
[考研] 288求调剂 +14 王晓阳- 2026-03-09 19/950 2026-03-14 02:05 by JourneyLucky
[考研] 一志愿天津大学，英一数二305分求调剂，四六级已过 +8 小小番的茄 2026-03-09 8/400 2026-03-14 01:53 by JourneyLucky
[考研] 环境调剂 +6 晓看天暮看云 2026-03-09 6/300 2026-03-14 01:16 by JourneyLucky
[考研] 一志愿安徽大学材料工程专硕313分，求调剂的学校 +8 Yu先生 2026-03-10 10/500 2026-03-14 01:04 by JourneyLucky
[基金申请] 有必要更换申报口吗 20+3 fannyamoy 2026-03-11 3/150 2026-03-14 00:52 by zhanghaozhu
[考研] 307求调剂 +7 超级伊昂大王 2026-03-10 7/350 2026-03-14 00:49 by JourneyLucky
[考研] 一志愿华中农业大学071010，总分三百二,求调剂 +3 困困困困坤坤 2026-03-10 3/150 2026-03-14 00:35 by JourneyLucky
[考研] 341求调剂 +3 番茄头--- 2026-03-10 3/150 2026-03-13 23:07 by JourneyLucky
[考研] 四川大学085601材料工程专硕 初试294求调剂 +4 祝我们好在冬天 2026-03-11 4/200 2026-03-13 21:39 by peike
[考研] 085600材料与化工 309分请求调剂 +7 dtdxzxx 2026-03-12 8/400 2026-03-13 14:43 by jxchenghu
[考研] 一志愿河海大学085900土木水利专硕279求调剂不挑专业 +4 SunWwWwWw 2026-03-10 8/400 2026-03-13 02:23 by SunWwWwWw
[考博] 2026年博士申请 +3 QwQwQW10 2026-03-11 3/150 2026-03-12 17:58 by gxch43
[考博] 读博申请 +5 感dd 2026-03-10 7/350 2026-03-11 17:02 by QGZDSYS
[考研] 298求调剂 +3 Vv呀！ 2026-03-10 3/150 2026-03-10 22:40 by 剑诗杜康
[考研] 一志愿山东大学，总分327，英语二79，有论文，有竞赛，已过四六级 +3 木木目目1 2026-03-09 3/150 2026-03-09 19:52 by yuningshan
信息提示
关闭
请填处理意见
关闭