24小时热门版块排行榜    

查看: 940  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 怎么查啊 +4 huang1991js 2026-08-26 4/200 2026-08-26 11:26 by panjy.cn
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 能否退出参与的面上项目解除限项 +22 koalala 2026-08-24 25/1250 2026-08-26 11:14 by dxcharlary
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 出来了 +8 trojank 2026-08-26 8/400 2026-08-26 09:00 by gxc
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 有没有大神帮我看看基金代码 28+4 1234567wang 2026-08-24 10/500 2026-08-25 19:15 by lfy8008
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 19/950 2026-08-25 17:57 by chick875
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 只有每年这种时候来逛逛小木虫 +25 yaoyewhu2008 2026-08-20 27/1350 2026-08-25 08:01 by Equinoxhua
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
请填处理意见