声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1439|回复: 8

在求解常微分方程的过程中,可以调用两个时刻的函数值吗?很急

[复制链接]
发表于 2006-3-22 16:16 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
在求解常微分方程时,需要调用计算时刻的前两个时刻的函数值作为一个判定条件,在求解微分方程的过程中,应该如何实现?谢谢大家
回复
分享到:

使用道具 举报

发表于 2006-3-22 19:32 | 显示全部楼层
怎么没有人回答<BR>我也有同样的问题<BR>顶上希望有人帮忙<BR>谢谢
发表于 2006-3-22 21:28 | 显示全部楼层

回复:(yuqing)在求解常微分方程的过程中,可以调用...

不懂你说的意思
 楼主| 发表于 2006-3-23 09:25 | 显示全部楼层

回复:(happy)回复:(yuqing)在求解常微分方程的过...

意思大体是这样:在关于x的常微分方程的求解过程中,是根据所取的tapan取了很多点来算得吗,在k步计算时,如何调用k-1 ,k-2 时刻的x的值<BR><BR>不知这样你明白了吗?
发表于 2006-3-23 09:43 | 显示全部楼层

回复:(yuqing)回复:(happy)回复:(yuqing)在求...

<DIV class=quote><B>以下是引用<I>yuqing</I>在2006-3-23 9:25:02的发言:</B><BR>意思大体是这样:在关于x的常微分方程的求解过程中,是根据所取的tapan取了很多点来算得吗,在k步计算时,如何调用k-1 ,k-2 时刻的x的值<BR><BR>不知这样你明白了吗?</DIV>
<br>没仔细想过,直接好像很难实现,需要该程序<BR>不过可以做定步长计算,相关定步长的m代码可以在matlab Central站点找到
 楼主| 发表于 2006-3-23 22:44 | 显示全部楼层
多谢
发表于 2006-3-24 10:28 | 显示全部楼层
我也遇到类似的问题,我现在解决的办法是:用ode45分段处理,把上一段的最后一个值作为下一段初始值,并且做条件判断,不过计算量会增加很多。
发表于 2006-3-24 10:37 | 显示全部楼层
在调用ode45的时候,如果把时间设定为[0 T/n T],就可以返回n=1,2,...时刻的值了。
发表于 2006-3-30 07:49 | 显示全部楼层
n=length(t);这样可以算出点数,再一个for(i=1:m),m&lt;=n,由需要取吧。
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-26 09:52 , Processed in 0.068606 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表