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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » C/C++ » 隐式欧拉法求解一阶常微分方程
51/1返回列表
查看: 4262  |  回复: 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

木虫 (正式写手)
余泽成(金币+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的回帖
查看全部 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 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +7 何气正 2026-04-13 8/400 2026-04-13 09:33 by 零零二
[考研] 291 求调剂 +30 化工2026届毕业� 2026-04-09 30/1500 2026-04-13 08:52 by lhj2009
[考研] 电子信息270求调剂 +18 terminal469 2026-04-07 18/900 2026-04-12 16:23 by ajpv风雷
[考研] 296求调剂 +14 汪！？！ 2026-04-10 16/800 2026-04-12 10:48 by zhouyuwinner
[考研] 070300化学279求调剂 +19 哈哈哈^_^ 2026-04-08 20/1000 2026-04-11 20:43 by stoner78
[考研] 调剂求助 +6 果然有我 2026-04-11 7/350 2026-04-11 16:22 by 明月此时有
[考研] 087100初试311求调剂 +4 任雅琴 2026-04-09 4/200 2026-04-11 10:33 by zhq0425
[考研] 一志愿985机械学硕380求调剂 +5 关关雎鸠10 2026-04-11 5/250 2026-04-11 10:10 by 知念。A
[考研] 一志愿东北大学控制工程085406数二英二385，求调剂 +8 Ezra_Zhang 2026-04-09 8/400 2026-04-11 09:15 by 猪会飞
[考研] 求调剂 +5 不会飞的鱼@ 2026-04-10 5/250 2026-04-10 19:07 by chemisry
[考研] 085404，285分求调剂 +12 薇薇考研 2026-04-07 14/700 2026-04-09 23:10 by parmtree
[考研] 材料调剂 +10 18815505510 2026-04-09 11/550 2026-04-09 17:07 by 544594351
[考研] 材料307分求大佬组收留 +17 Hll胡 2026-04-07 17/850 2026-04-09 10:53 by liuhuiying09
[考研] 求调剂 +13 柒luck 2026-04-07 13/650 2026-04-08 22:46 by 猪会飞
[考研] 生物学328分求调剂 +9 闪电kkl 2026-04-08 10/500 2026-04-08 21:42 by liuhuiying09
[考研] 材料考研求调剂总分280 +30 mkjlz1 2026-04-06 35/1750 2026-04-08 21:25 by cyh—315
[考研] 277求调剂 +4 考研调剂lxh 2026-04-06 6/300 2026-04-08 10:40 by 逆水乘风
[考研] 287求调剂 +6 Fnhc 2026-04-07 6/300 2026-04-08 10:05 by xingguangj
[考研] 22408 318分求调剂 +4 勤奋的小笼包 2026-04-06 6/300 2026-04-07 15:05 by 纸鹤555
[考研] 求调剂到材料 +5 程9915 2026-04-06 5/250 2026-04-06 15:21 by yulian1987
信息提示
关闭
请填处理意见
关闭