24小时热门版块排行榜    

查看: 859  |  回复: 6

利维坦6021

新虫 (初入文坛)

[求助] 【求助】使用fortran95牛顿法求根

最近被一个超越方程困扰了很久,因为对fortran不太熟悉,因此求教各位如何使用fortran语言求这个方程的根呢?
公示在图中,其中p为常数。
在网上搜到的程序包有点看不明白,希望各位不吝赐教

【求助】使用fortran95牛顿法求根
捕获.PNG
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

scoot

新虫 (初入文坛)

2楼2017-09-15 18:08:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijf517

木虫 (小有名气)

这个方程这么简单,还搜什么程序包啊。自己写代码,十行内搞定。标准的牛顿迭代法。选个合适的初值x0,然后x_(k+1)=x_k-f‘_k。随便找本数值分析的书都会讲得很清楚。

发自小木虫Android客户端
3楼2017-09-16 06:25:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ylsxz2012

铜虫 (初入文坛)

引用回帖:
3楼: Originally posted by lijf517 at 2017-09-16 06:25:12
这个方程这么简单,还搜什么程序包啊。自己写代码,十行内搞定。标准的牛顿迭代法。选个合适的初值x0,然后x_(k+1)=x_k-f‘_k。随便找本数值分析的书都会讲得很清楚。
...

的确是最简单的。一般微积分书上都有。给您更正一个小小的笔误。
x_(k+1)=x_k-f‘_k*dx,
这里dx是步长。
4楼2017-09-16 08:57:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ylsxz2012

铜虫 (初入文坛)

引用回帖:
3楼: Originally posted by lijf517 at 2017-09-16 06:25:12
这个方程这么简单,还搜什么程序包啊。自己写代码,十行内搞定。标准的牛顿迭代法。选个合适的初值x0,然后x_(k+1)=x_k-f‘_k。随便找本数值分析的书都会讲得很清楚。
...

我刚才应该是dt不是dx,你用x做因变量了,没有仔细看。呵呵。
x_(k+1)=x_k-x‘_k*dt
dt步长。
5楼2017-09-16 09:01:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijf517

木虫 (小有名气)

引用回帖:
5楼: Originally posted by ylsxz2012 at 2017-09-16 09:01:54
我刚才应该是dt不是dx,你用x做因变量了,没有仔细看。呵呵。
x_(k+1)=x_k-x‘_k*dt
dt步长。...

谢谢指正,不过也不是dt或dx,这个迭代过程与步长无关。应该为x_(k+1)=x_k-f(x_k)/f’(x_k)。其中,f(x)=0为带求解的方程。

发自小木虫Android客户端
6楼2017-09-16 15:01:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ylsxz2012

铜虫 (初入文坛)

引用回帖:
6楼: Originally posted by lijf517 at 2017-09-16 15:01:16
谢谢指正,不过也不是dt或dx,这个迭代过程与步长无关。应该为x_(k+1)=x_k-f(x_k)/f’(x_k)。其中,f(x)=0为带求解的方程。
...

嗯,牛顿法是求根。我那个有步长的是求函数值,不是牛顿法。呵呵,我犯了个方法错误,虽然两个公式都用一个公式推导。谢谢指正,也免得误导其他初学者。
7楼2017-09-17 01:39:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 利维坦6021 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 341求调剂 +5 捣蛋猪猪 2026-03-11 7/350 2026-03-17 19:09 by 捣蛋猪猪
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
[考研] 286求调剂 +3 lemonzzn 2026-03-16 5/250 2026-03-16 20:43 by lemonzzn
[考研] 一志愿211 0703方向310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-16 16:44 by houyaoxu
[考研] 0703化学调剂 +6 妮妮ninicgb 2026-03-15 9/450 2026-03-16 16:40 by houyaoxu
[考研] 304求调剂 +3 曼殊2266 2026-03-14 3/150 2026-03-16 16:39 by houyaoxu
[考研] 求老师收留调剂 +4 jiang姜66 2026-03-14 5/250 2026-03-15 20:11 by Winj1e
[考研] 328求调剂 +3 5201314Lsy! 2026-03-13 6/300 2026-03-14 15:31 by hyswxzs
[考研] 一志愿哈工大材料324分求调剂 +5 闫旭东 2026-03-14 5/250 2026-03-14 14:53 by 木瓜膏
[考研] 297求调剂 +4 学海漂泊 2026-03-13 4/200 2026-03-14 11:51 by 热情沙漠
[考研] 0703,333分求调剂 一志愿郑州大学-物理化学 +3 李魔女斗篷 2026-03-11 3/150 2026-03-13 22:24 by JourneyLucky
[考研] 336求调剂 +6 Iuruoh 2026-03-11 6/300 2026-03-13 22:06 by JourneyLucky
[考研] 0856材料与化工301求调剂 +5 奕束光 2026-03-13 5/250 2026-03-13 22:00 by 星空星月
[考研] 求调剂 +7 18880831720 2026-03-11 7/350 2026-03-13 16:10 by JourneyLucky
[考研] 工科278分求调剂 +5 周慢热啊 2026-03-12 7/350 2026-03-13 15:49 by JourneyLucky
[考研] 一志愿211化学学硕310分求调剂 +8 努力奋斗112 2026-03-12 9/450 2026-03-13 15:41 by JourneyLucky
[考研] 0856化学工程280分求调剂 +4 shenzxsn 2026-03-11 4/200 2026-03-13 11:55 by ymwdoctor
[考博] 福州大学杨黄浩课题组招收2026年专业学位博士研究生,2026.03.20截止 +3 Xiangyu_ou 2026-03-12 3/150 2026-03-13 09:36 by duanwu655
[考博] 26读博 +4 Rui135246 2026-03-12 10/500 2026-03-13 07:15 by gaobiao
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-03-12 4/200 2026-03-12 19:33 by 求调剂zz
信息提示
请填处理意见