24小时热门版块排行榜    

查看: 2697  |  回复: 4

有机物

至尊木虫 (著名写手)


[交流] 【求助】同样一个问题用matlab和fortran编程计算,结果相差比较大,不知道为什么

求解一个偏微分方程,过程一样,分别用matlab和fortran编程,不知道为什么,结果相差很大。有没有虫友碰到过类似的问题啊,纠结中。
回复此楼

» 猜你喜欢

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

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijinfeng042

木虫 (小有名气)


★ ★
有机物(金币+2):谢谢参与
wuming524(金币+1):谢谢应助 2011-01-08 01:04:59
matlab对初值依赖很大 有时候甚至于没有结果...尝试maple对比一下
如果是自己写的算法 那就应该可以每一步都可以检查啊
2楼2011-01-06 10:49:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★
有机物(金币+2):谢谢参与
wuming524(金币+1):谢谢应助 2011-01-08 01:05:09
引用回帖:
Originally posted by 有机物 at 2011-01-05 22:43:42:
求解一个偏微分方程,过程一样,分别用matlab和fortran编程,不知道为什么,结果相差很大。有没有虫友碰到过类似的问题啊,纠结中。

MATLAB只能求解一些比较简单的偏微方程,对复杂一些的PDE,需要自己写程序实现。
不懂FORTRAN,但是用c写过简单的计算程序,用的有限差分进行的离散,对规则的求解域,还比较容易实现。

不知道你的PDE是什么样的,建议用COMSOL计算一下作为比较。

如果算法都正确,计算结果应该只有精度误差的。
3楼2011-01-06 22:39:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

有机物

至尊木虫 (著名写手)


引用回帖:
Originally posted by xiegangmai at 2011-01-06 22:39:37:


MATLAB只能求解一些比较简单的偏微方程,对复杂一些的PDE,需要自己写程序实现。
不懂FORTRAN,但是用c写过简单的计算程序,用的有限差分进行的离散,对规则的求解域,还比较容易实现。

不知道你的PDE是什 ...

因为我的PDE的边界条件比较特殊,不能直接用MATLAB里面的PDE工具解。我也是用有限差分法自己写程序算,算法都一样,只是用MATLAB和FORTRAN来计算,结果相差很大。仔细检查了每一步的计算,实在是找不到哪里有问题,郁闷啊。
4楼2011-01-07 10:52:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

信彼南山

木虫 (著名写手)


★ ★
有机物(金币+2):谢谢参与
xiegangmai(金币+1):谢谢参与 2011-01-09 22:36:23
把每一步的值都记录一下,看看偏差出现在什么地方了
如果是我的话,我选择相信matlab。

看看你的误差产生原因是什么,是累积误差还是截断误差
或者尝试修改步长试试看,步长与两个计算的结果成什么样的关系
5楼2011-01-09 22:34:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 有机物 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 面上函评意见出来了,像什么等级? 20+3 Tsingking1 2026-08-27 14/700 2026-09-01 08:29 by 尼古拉斯小虫
[基金申请] 国社科又开始会评了,不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 21:30 by xiangy672
[基金申请] 基金未中,这种答复是模板吗? +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[文学芳草园] 梦想 +5 myrtle 2026-08-26 5/250 2026-08-31 14:40 by hahaboy
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
请填处理意见