24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2019  |  回复: 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 的主题更新
信息提示
请填处理意见