应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » C/C++ » 隐式欧拉法求解一阶常微分方程
汕头大学海洋科学接受调剂
51/1返回列表
查看: 4266  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

微尘、梦想

木虫 (知名作家)

[求助] 隐式欧拉法求解一阶常微分方程

算法：


MATLAB程序：
CODE:
Format long;
N=(b-a)/h;
y=zeros(N+1,1);
y(1)=y0;
x=a:h:b;
var=findsym(f);
for i=2:N+1
   fx=Funval(f,var(1),x(i));
   gx=y(i-1)+h*fx-varvec(2);
   y(i)=NewtonRoot(gx,-10,10,eps);%用牛顿法得出下步的迭代值
end
format short;

其中，Funval是求函数f在点(var(1),x(i))处的值

PS：显示法求解常微分方程的算法都用C语言实现，但隐式法不知道怎么实现，主要是不知道如何求后一个点的导数值，希望高人给指点一下，呵呵……
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
任风云变幻，我笑对人生！
1楼 2012-02-23 12:20:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sudo

木虫 (正式写手)
余泽成(金币+1): 鼓励讨论！ 2012-02-25 04:34:54
引用回帖:
7楼: Originally posted by 微尘、梦想 at 2012-02-23 16:50:00:
嗯，我也就是想把它搞明白它到底是如何工作的，解常微分方程的显示解法，我都搞明白了，从最简单的欧拉法到四阶龙格-库塔法，不过我在网上找到一篇关于算法分析的文章，结论是隐式算法的稳定性普遍比显式算法高， ...

嗯，隐式欧拉的确是有它的意义的。。。不过你好像误会了我的意思我说的是估计y[n+1]那里不用牛顿法而用显式欧拉法来做会方便点
赞一下 回复此楼
8楼2012-02-23 17:19:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

sudo

木虫 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1
余泽成(金币+1): 谢谢参与应助！ 2012-02-25 04:34:01
微尘、梦想(金币+30): ★★★★★最佳答案 谢谢，sudo教会了我很多东西呀，呵呵…… 2012-02-27 13:57:13
隐式欧拉法迭代式右边f(x[n+1], y[n+1])里面的y[n+1]，其实也只是用估计的方法得出来的

一般可以用欧拉法（或者牛顿法和其他别的什么方法）得到它的估计值，所以这个过程实际上是：
CODE:
yt[n+1] = y[n] + h * f(x[n], y[n]);
y[n+1] = y[n] + h * f(x[n+1], yt[n+1]);

梯形法什么的类似，迭代式中f里面的y[n+1]，无论怎样都是需要某种方法先估计出来的
赞一下 回复此楼
2楼2012-02-23 13:25:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sudo

木虫 (正式写手)
恩，如果不要求通用的方法，只需要针对特定函数的话，也可以先化简迭代式，或者可能可以直接导出y[n+1]的公式什么的~
赞一下(1人) 回复此楼
3楼2012-02-23 13:29:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

微尘、梦想

木虫 (知名作家)
引用回帖:
3楼: Originally posted by sudo at 2012-02-23 13:29:11:
恩，如果不要求通用的方法，只需要针对特定函数的话，也可以先化简迭代式，或者可能可以直接导出y[n+1]的公式什么的~

那个MATLAB程序就是用牛顿法得到的下一个值，不过，不知道牛顿法的算法是什么呀，我只知道用牛顿法解方程的公式：
x[n+1]=x[n]-f(x[n])/f`(x[n])
这是一种切线法，从一端向一个方向逼近方程的根，可是跟这个联系不到一块呀
赞一下 回复此楼
任风云变幻，我笑对人生！
4楼2012-02-23 13:38:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 调剂求收留 +32 果然有我 2026-04-10 33/1650 2026-04-14 08:49 by 木木mumu～
[考研] 310求调剂 +15 666真好 2026-04-11 17/850 2026-04-13 22:27 by pies112
[考研] 302求调剂 +10 易！? 2026-04-13 10/500 2026-04-13 19:04 by lbsjt
[考研] 339求调剂 +4 hanwudada 2026-04-12 4/200 2026-04-13 12:03 by 蓝云思雨
[考研] 0854调剂 +10 长弓傲 2026-04-11 11/550 2026-04-13 10:38 by wp06
[考研] 一志愿浙大生物325分求调剂 +9 zysheng 2026-04-12 9/450 2026-04-12 22:31 by yuyin1233
[考研] 22408 352分求调剂 +5 努力的夏末 2026-04-09 5/250 2026-04-12 19:17 by wj165256
[考研] 材料工程日语考生求调剂 +7 0856?调剂 2026-04-10 7/350 2026-04-11 21:33 by 蓝云思雨
[考研] 303求调剂 +14 SereinQ 2026-04-10 15/750 2026-04-11 20:43 by 蓝云思雨
[考研] 085410-273求调剂 +6 X1999 2026-04-10 6/300 2026-04-11 10:32 by Delta2012
[考研] 288求调剂 +15 代fish 2026-04-09 16/800 2026-04-11 10:26 by wwj2530616
[考研] 346，工科求调剂 +3 moser233 2026-04-09 3/150 2026-04-11 10:04 by zhq0425
[考研] 085506-求调剂-285分 +3 雷欧飞踢 2026-04-08 3/150 2026-04-11 08:37 by zhq0425
[考研] 材料类284调剂 +40 想换手机不想解� 2026-04-08 48/2400 2026-04-10 23:28 by 314126402
[考研] 初试261 +3 Asht少 2026-04-10 6/300 2026-04-10 16:38 by Asht少
[考研] 一志愿京区985，085401电子信息，本科电子信息 +3 阳光开朗的男孩 2026-04-10 3/150 2026-04-10 16:29 by sophia_93
[考研] 085800 能源动力求调剂 +6 阿biu啊啊啊啊啊 2026-04-10 6/300 2026-04-10 15:03 by hemengdong
[考研] 085404，334分，求调剂 +5 sunjie8888 2026-04-08 8/400 2026-04-09 07:26 by sunjie8888
[考研] 软件工程求调剂22软工296分求调剂，接受跨调 +4 yangchen2017 2026-04-08 5/250 2026-04-08 21:56 by 土木硕士招生
[考研] 259求调剂 +5 就爱吃土豆呀呀 2026-04-07 5/250 2026-04-07 22:40 by JourneyLucky
信息提示
关闭
请填处理意见
关闭