24小时热门版块排行榜    

Znn3bq.jpeg
北京石油化工学院2026年研究生招生接收调剂公告
查看: 865  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学308分求调剂 +13 你好明天你好 2026-04-07 15/750 2026-04-07 19:41 by wenqingshen
[考研] 材料求调剂 +18 一样YWY 2026-04-05 18/900 2026-04-07 15:49 by dxlg
[考研] 328求调剂 +4 ghhh88888 2026-04-06 5/250 2026-04-07 14:45 by ghhh88888
[考研] 362求调剂 +7 我要考大 2026-04-06 11/550 2026-04-07 09:32 by 我要考大
[考研] 277工科求调剂 +10 1915668 2026-04-04 11/550 2026-04-06 23:53 by 心媛意北
[考研] 287分求调剂 有专利国奖一志愿哈工大085406 +6 白易辰 2026-04-06 7/350 2026-04-06 22:46 by 875465
[考研] 生物学求调剂 一志愿沪9,326分 +6 刘墨墨 2026-04-06 6/300 2026-04-06 19:36 by lijunpoly
[考研] 071000生物学调剂 +7 拉提桃 2026-04-06 7/350 2026-04-06 18:55 by 52305043001
[考研] 362求调剂一志愿中国石油大学 +4 我要考大 2026-04-06 6/300 2026-04-06 14:11 by 无际的草原
[考研] 315求调剂 +5 &123456789 2026-04-05 5/250 2026-04-05 19:55 by nepu_uu
[考研] 08专硕275调剂 +5 AaAa7420 2026-04-05 5/250 2026-04-05 18:01 by jkddd
[考研] 本科211 分数293请求调剂 +4 莲菜就是藕吧 2026-04-01 4/200 2026-04-04 22:32 by hemengdong
[考研] 材料专硕322分 +13 哈哈哈吼吼吼哈 2026-04-01 13/650 2026-04-03 16:08 by 哦哦123
[考研] 求调剂 +3 usbdndj 2026-04-03 3/150 2026-04-03 14:10 by dxiaoxin
[考研] 专硕085601求调剂 +7 suyifei 2026-04-03 8/400 2026-04-03 14:00 by 欣喜777
[考研] 求材料调剂 一志愿南昌大学 328分 +5 yyy..... 2026-04-03 5/250 2026-04-03 13:46 by 百灵童888
[考研] 专硕 351 086100 也是考的材科基 本科也是材料 +8 202451007219 2026-04-02 8/400 2026-04-03 09:50 by 蓝云思雨
[考研] 一志愿山东大学,085600,344 +7 魏子per 2026-04-02 8/400 2026-04-02 21:12 by 百灵童888
[考研] 301求调剂 +8 axibli 2026-04-01 8/400 2026-04-01 09:51 by 我的船我的海
[考研] 求调剂:一志愿:南京大学 专业:0705 总分320 ,本科985,四六级已过 +3 lfy760306 2026-03-31 3/150 2026-04-01 01:57 by Creta
信息提示
请填处理意见