24小时热门版块排行榜    

查看: 2735  |  回复: 9

ws89

金虫 (著名写手)

悠然书生

[求助] fortran程序不收敛如何调程序 已有2人参与

希望大家能够帮助我,我学的是工程力学,用Fortran程序来计算中空圆柱的力学,热力学性能,包括位移应变应力等。程序的实现过程就是在拉普拉斯域内求解出变量的表达式,然后在用经验公式近似反变换求出实数变量,然后再用origin画图。
     我尝试过以下操作:1,稍稍改变起点位置,比如空心圆柱半径从1到3,我改为1.01到3,
                                       或者直接把半径范围改到从5到15
                                  2,改变计算的精度,误差度从10的-10次改到10的-2或-3次,
但是起点或者第二三点仍旧有不收敛的情况,不能得出计算结果,而且精度改大的话,不一定满足计算可靠性,向大家跪拜了,程序高手们,力学精英们,你们来救救我吧
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

天佑一帆&soul.St.realand
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็

为啥用Fortran?
ส็็็็็็็็็็็็็็็็็็็็
2楼2014-07-25 22:36:59
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

ws89

金虫 (著名写手)

悠然书生

引用回帖:
2楼: Originally posted by cooooldog at 2014-07-25 22:36:59
为啥用Fortran?

我们的论文就是这么做的呀,程序调不出来的话只能干耗着,愁啊
天佑一帆&soul.St.realand
3楼2014-07-25 23:29:16
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็

引用回帖:
3楼: Originally posted by ws89 at 2014-07-25 23:29:16
我们的论文就是这么做的呀,程序调不出来的话只能干耗着,愁啊...

换个语言, 比如matlab 或python
否则算好还得用origin,多别扭
ส็็็็็็็็็็็็็็็็็็็็
4楼2014-07-26 06:57:36
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

ws89

金虫 (著名写手)

悠然书生

引用回帖:
4楼: Originally posted by cooooldog at 2014-07-26 06:57:36
换个语言, 比如matlab 或python
否则算好还得用origin,多别扭...

不会呀
天佑一帆&soul.St.realand
5楼2014-07-26 10:57:58
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

锐利的碎片

木虫 (正式写手)

star watcher

算法,程序实现bug都有可能。具体要看程序。

[ 发自手机版 http://muchong.com/3g ]
6楼2014-07-26 11:50:50
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

ws89

金虫 (著名写手)

悠然书生

引用回帖:
6楼: Originally posted by 锐利的碎片 at 2014-07-26 11:50:50
算法,程序实现bug都有可能。具体要看程序。

有没有一些切实可行的方法呀,可以让我操作尝试一下
天佑一帆&soul.St.realand
7楼2014-07-28 21:41:05
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

绿遍山原

铜虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
引用回帖:
7楼: Originally posted by ws89 at 2014-07-28 21:41:05
有没有一些切实可行的方法呀,可以让我操作尝试一下...

把你的程序分解成几个独立部分,分别测试,看看有没有错误。

[ 发自手机版 http://muchong.com/3g ]
要夢遊,不要催眠。
8楼2014-07-28 22:17:04
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

ws89

金虫 (著名写手)

悠然书生

引用回帖:
8楼: Originally posted by 绿遍山原 at 2014-07-28 22:17:04
把你的程序分解成几个独立部分,分别测试,看看有没有错误。
...

我们研究的是分数阶导数的问题,整数阶导数就收敛了,但是分数阶不收敛,不知道如何处理呀
天佑一帆&soul.St.realand
9楼2014-07-29 22:49:58
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

carey-only

铁虫 (初入文坛)

【答案】应助回帖

这种情况需要在反复检查方程的求解正确与否,即方程的离散方法是否对?并试着调初始值
10楼2014-09-29 13:17:42
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ws89 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 考博求助 10+3 社会实践对别 2026-11-08 3/150 2026-11-13 14:04 by 沐小虫252
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 gRND4NVAZ8VG 2026-11-11 5/250 2026-11-13 11:15 by Equinoxhua
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 FChKX8meu1ds 2026-11-09 10/500 2026-11-13 11:08 by Equinoxhua
[基金申请] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +6 C6CRn66RYYXS 2026-11-10 7/350 2026-11-13 10:55 by Equinoxhua
[考博] 博士申请 +4 z_l_l_666 2026-11-06 4/200 2026-11-13 06:19 by wlf1982
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 TgvkbZZl9tqB 2026-11-11 6/300 2026-11-12 12:42 by ECLTXnJnRSXf
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 cdqGD6ftHmbl 2026-11-11 5/250 2026-11-12 12:28 by ECLTXnJnRSXf
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 TgvkbZZl9tqB 2026-11-11 4/200 2026-11-12 12:28 by ECLTXnJnRSXf
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 Z7qfCsdIiD1T 2026-11-10 9/450 2026-11-12 11:22 by HajghKOMynmz
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 czAiolzR0IJC 2026-11-10 6/300 2026-11-12 11:03 by HajghKOMynmz
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +7 czAiolzR0IJC 2026-11-10 9/450 2026-11-12 11:00 by HajghKOMynmz
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 czAiolzR0IJC 2026-11-10 8/400 2026-11-12 11:00 by HajghKOMynmz
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 G2rhSvqvhNzl 2026-11-10 6/300 2026-11-12 10:49 by HajghKOMynmz
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 czAiolzR0IJC 2026-11-10 8/400 2026-11-12 10:49 by HajghKOMynmz
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 C6CRn66RYYXS 2026-11-10 6/300 2026-11-12 06:03 by 9fOTqj074tIg
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 FChKX8meu1ds 2026-11-09 10/500 2026-11-12 05:19 by 9fOTqj074tIg
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 zSbvVX9JFl9J 2026-11-09 8/400 2026-11-12 04:57 by 9fOTqj074tIg
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 Z7qfCsdIiD1T 2026-11-10 4/200 2026-11-12 02:23 by 9fOTqj074tIg
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 12/600 2026-11-09 06:32 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 5sJFDKp2kwuI 2026-11-06 9/450 2026-11-08 22:46 by l48hFd8Nwae2
信息提示
请填处理意见