碰撞求解

无水1324 · 发表于 2007-7-21 11:15

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

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

x

function dydt = vv(t,y)
dydt = [y(2);-0.1*y(2)-y(1)+cos(2.8*(t+0.019))];

上面这个方程，假设给他加一个控制条件。
当位移x=0.5时候产生碰撞，然后反向。但是此时速度变为碰撞前的0.8倍。
请问用ode45求解的时候怎么控制？

octopussheng · 发表于 2007-7-21 12:22

加一个判断语句，类似于分段线性函数的写法，不知道可不可行！

无水1324 · 发表于 2007-7-21 13:55

这样不行

octopussheng · 发表于 2007-7-21 14:39

那考虑在求解程序中添加判断语句，

想了下，还真不好实现！难哦！

是不是考虑做一个循环，在循环中加入一个判断语句，我再想想！

hohoo · 发表于 2007-7-21 17:05

我这样考虑给一个位移的误差比如 0.5+(1e－5)和 0.5-(1e-5), 如果位移进入这个范围内，就反向当然位移的误差越小越好计算时间就长了:@)

无水1324 · 发表于 2007-7-21 17:36

原帖由 hohoo 于 2007-7-21 17:05 发表

 登录/注册后可看大图

我这样考虑给一个位移的误差比如 0.5+(1e－5)和 0.5-(1e-5), 如果位移进入这个范围内，就反向当然位移的误差越小越好计算时间就长了:@)

想法是这样的，但是直接用ode45求解的时候是控制不住的

中原 · 发表于 2007-7-21 20:52

用fortran或c编个程序就ok，龙特库塔求解模块很好找，然后几条判断语句ok

octopussheng · 发表于 2007-7-22 08:36

确实修改Runge－Kutta法计算程序比我们在这里想着怎么修改函数去让RK识别要好很多，呵呵！

在下面贴两个程序模块，一个是Fotran的，一个是C的，省得无水去找了！

直接作为模块使用，我测试过，没有问题的！

无水1324 · 发表于 2007-7-22 09:21

原帖由 octopussheng 于 2007-7-22 08:36 发表

 登录/注册后可看大图

确实修改Runge－Kutta法计算程序比我们在这里想着怎么修改函数去让RK识别要好很多，呵呵！

在下面贴两个程序模块，一个是Fotran的，一个是C的，省得无水去找了！

直接作为模块使用，我测试过，没有问题的！

你说的是已经试着解决上面那个问题了？

octopussheng · 发表于 2007-7-22 09:38

不是，我的意思是这两个模块我用别的方程算过，模块本身是没有问题的

但是里面还需要作些修改，像中原说的，加一些判断语句！

我这两天在算些东西，有些忙不过头，所以就不能帮着修改啦！

无水加油哦，呵呵！

无水1324 · 发表于 2007-7-23 16:32

已经想到一个比较笨的方法解决。就是直接自己编写一个定步长RK 法。如：
http://forum.vibunion.com/forum/ ... p%3Bfilter%3Ddigest

octopussheng · 发表于 2007-7-23 19:51

确实，这个因为是在计算过程中进行判断的，只能根据你的要求去修改RK方法了

xiejinling · 发表于 2008-1-15 18:19

可以用解析法来做吧！但是我现在也还没有调出来！！

无水1324 · 发表于 2008-1-15 19:55

但是我这个解析解比较困难啊

beel · 发表于 2008-12-12 22:35

请问碰撞需要哪些数学方面的知识

账号		自动登录	找回密码
密码			我要加入

[近似分析] 碰撞求解

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

回复 #2 octopussheng 的帖子

评分

回复 13楼的帖子

[近似分析] 碰撞求解

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

回复 #2 octopussheng 的帖子

评分

回复 13楼 的帖子

回复 13楼的帖子