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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » C/C++ » 隐式欧拉法求解一阶常微分方程
汕头大学海洋科学接受调剂
81/1返回列表
查看: 4274  |  回复: 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
余泽成(金币+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的回帖

sudo

木虫 (正式写手)
余泽成(金币+2): 谢谢参与应助！ 2012-02-25 04:34:19
引用回帖:
4楼: Originally posted by 微尘、梦想 at 2012-02-23 13:38:29:
那个MATLAB程序就是用牛顿法得到的下一个值，不过，不知道牛顿法的算法是什么呀，我只知道用牛顿法解方程的公式：
x[n+1]=x[n]-f(x[n])/f`(x[n])
这是一种切线法，从一端向一个方向逼近方程的根，可是跟这个联 ...

牛顿法是用来做什么的呢？一般是用来求非线性方程的根的：

对于f(x)=0
迭代式x[n+1] = x[n] - f(x[n])/f`(x[n])
迭代的结果就是f(x)=0的数值解

=============================================

对于y[n+1] = y[n] + h * f(x[n+1], y[n+1])
CODE:
   fx=Funval(f,var(1),x(i)); %里面的x(i)就相当于x[n+1]，var(1)是一个符号变量（其实是等下要求解的y[n+1]），把它们代入到f中之后，得到一个符号表达式
   gx=y(i-1)+h*fx-varvec(2); %其实这里就是方程g(x)=0左边的g(x)，就是y[n] + h*fx - y[n+1] = 0这个方程，里面y[n]是上一步知道的值，fx是含有变量y[n+1]的式子，外面的y[n+1]也视为变量
   y(i)=NewtonRoot(gx,-10,10,eps); %这就是用牛顿法来解非线性方程了

也就是说，这段MATLAB程序里面涉及了符号的运算
赞一下 回复此楼
5楼2012-02-23 15:34:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sudo

木虫 (正式写手)
余泽成(金币+1): 鼓励交流！ 2012-02-25 04:34:33
如果它的方法要通用的话，在NewtonRoot还必须对gx进行一次符号求导（用diff函数），所以这个方法其实比较蛋疼

还不如直接用前向欧拉给估计了...
赞一下 回复此楼
6楼2012-02-23 15:41:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

微尘、梦想

木虫 (知名作家)
引用回帖:
6楼: Originally posted by sudo at 2012-02-23 15:41:38:
如果它的方法要通用的话，在NewtonRoot还必须对gx进行一次符号求导（用diff函数），所以这个方法其实比较蛋疼

还不如直接用前向欧拉给估计了...

嗯，我也就是想把它搞明白它到底是如何工作的，解常微分方程的显示解法，我都搞明白了，从最简单的欧拉法到四阶龙格-库塔法，不过我在网上找到一篇关于算法分析的文章，结论是隐式算法的稳定性普遍比显式算法高，且只有隐式欧拉方法才是绝对稳定的。
也没仔细看，反正用四阶龙格-库塔法的精度已经非常高了，呵呵……
这是文章的地址：
http://wenku.baidu.com/view/8bcb52492e3f5727a4e96207.html
赞一下 回复此楼
任风云变幻，我笑对人生！
7楼2012-02-23 16:50:00
已阅   回复此楼   关注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的回帖
相关版块跳转 我要订阅楼主 微尘、梦想 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085801电气专硕272求调剂 +16 电气李 2026-04-13 18/900 2026-04-15 09:17 by heqizheng
[考研] 通信工程求调剂！！！ +4 zlb770521 2026-04-14 4/200 2026-04-14 18:19 by lbsjt
[考研] 302分求调剂 +10 凡语祈愿 2026-04-08 11/550 2026-04-14 16:50 by jiangguiquan11
[考研] 恳请有学校收留 +3 柯淮然 2026-04-12 3/150 2026-04-14 16:25 by 逆水乘风
[考研] 材料专业344求调剂 +17 hualkop 2026-04-10 22/1100 2026-04-14 16:21 by sxdj2
[考研] 药学305求调剂 +10 玛卡巴卡boom 2026-04-10 10/500 2026-04-14 15:55 by zs92450
[考研] 297求调剂 +23 ORCHID1 2026-04-10 26/1300 2026-04-14 13:52 by 陈皮皮
[考研] 本科西工大 324求调剂 +5 wysyjs25 2026-04-10 5/250 2026-04-13 23:08 by pies112
[考研] 293求调剂 +16 我爱高数高数爱� 2026-04-12 18/900 2026-04-13 21:47 by 学员JpLReM
[考研] 0854调剂 +10 长弓傲 2026-04-11 11/550 2026-04-13 10:38 by wp06
[考研] 求调剂 +16 张番茄不炒蛋 2026-04-10 17/850 2026-04-12 13:58 by 熬夜成！
[找工作] 山东高校教师考核超级无底线，员工过不下去啦 +4 qut2026 2026-04-09 9/450 2026-04-12 00:54 by qut2026
[考研] 0859，337求调剂 +4 研s. 2026-04-10 4/200 2026-04-11 11:34 by caotw2020
[考研] 0854调剂 +5 音像店听花鼓戏 2026-04-10 5/250 2026-04-11 10:49 by qingpingzhu
[考研] 309求调剂 +14 wdhw 2026-04-10 15/750 2026-04-10 21:06 by zhouxiaoyu
[考研] 吉大计算机技术331分，英语六级，求调剂 +3 峰峰021116 2026-04-09 3/150 2026-04-10 20:01 by chemisry
[考研] 298求调剂 +13 钉叮咚冬瓜 2026-04-09 13/650 2026-04-10 15:49 by jiajinhpu
[考研] 江苏大学 工科调剂 捡漏 +3 Evan_Liu 2026-04-09 5/250 2026-04-10 10:22 by Evan_Liu
[考研] 求机械专硕297第二批调剂 +5 拾柒12。 2026-04-08 5/250 2026-04-09 16:43 by 允当适度
[考研] 331求调剂 +5 luoxin0706. 2026-04-08 5/250 2026-04-08 22:15 by zhouyuwinner
信息提示
关闭
请填处理意见
关闭