24小时热门版块排行榜    

查看: 2488  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 时间戳又变了8-15 +12 archvillain 2026-08-15 23/1150 2026-08-15 20:47 by foyo404
[基金申请] filecode=后面第一个是大写字母 +6 wangze12014 2026-08-14 7/350 2026-08-15 20:12 by gltch
[基金申请] filecode +8 cratir 2026-08-14 12/600 2026-08-15 18:08 by zyfgau
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +12 Tide man 2026-08-10 13/650 2026-08-15 16:34 by 氺木
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +7 Ldrop2023 2026-08-13 7/350 2026-08-15 15:57 by mzhh2000
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +23 Lanmanbaby 2026-08-09 37/1850 2026-08-14 13:40 by Equinoxhua
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
信息提示
请填处理意见