24小时热门版块排行榜    

查看: 1063  |  回复: 3

莫selina

新虫 (初入文坛)

[求助] 这个方程怎样实现差分?(偶是新人,目前只有两个金币,不要嫌少啊)

这样一个方程怎样实现差分?
在n点t时刻,方程如图:

未命名.jpg
http://kuai.xunlei.com/d/GXKYUGYYBUXX?p=130497
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangww2011

木虫 (著名写手)

这个方程应该可以直接求出解析解的吧 不需要求查分的啊

即使求差分,你的v_{n}和v_{n-1}是什么意思?是不同时刻的v吗还是同一时刻有几个的v?
2楼2012-09-03 20:52:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

莫selina

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by wangww2011 at 2012-09-03 20:52:01
这个方程应该可以直接求出解析解的吧 不需要求查分的啊

即使求差分,你的v_{n}和v_{n-1}是什么意思?是不同时刻的v吗还是同一时刻有几个的v?

我做的车辆的模拟,n是指第n辆车,n-1 指前面那一辆车。
3楼2012-09-04 14:32:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangww2011

木虫 (著名写手)

【答案】应助回帖

★ ★
csgt0: 金币+2, 应助指数+1, 多谢应助 2012-09-04 22:49:28
引用回帖:
3楼: Originally posted by 莫selina at 2012-09-04 14:32:35
我做的车辆的模拟,n是指第n辆车,n-1 指前面那一辆车。...

其实查分本身非常容易,就是因为太容易,所以有可能求解你这个方程是不稳定的,所以也不是很建议你直接用,不过你是不是应该提供更多的信息?

你给出T为常数,其实这个常数是n的函数吧?譬如第一辆车T应该为零的,否则整个方程组没法开始。从第二辆车开始,加速度可以用你给的式子计算,但是第一辆车的的加速度就不能了,所以在你的模拟中,第一辆车的速度其实应该是你提前假定好的(当然可以变化,但你是不是应该提前给出个函数呢),然后也作为方程的初始条件对吧?

你把加速度写成dv/dt是应该的,也是求解微分方程组的第一步。

至于你的这个方程组,假设有n辆车,你需要求解的是一个有n-1个方程的方程组,本身很容易求解,如果你懂点计算机的话,你可以用matlab的ode45,或者在mathematica里直接敲出式子来给出初始条件就可以直接得出结果了。或者如果你用python,直接用scipy.integrate.odeint就可以了,很简单的,而且貌似前段时间有本一个家伙写的python数值计算的中文教程,你拿来读一读(当字典来查就行了),你这个问题也就几十行的代码吧,加上漂亮的绘图,应该也不长的。
4楼2012-09-04 17:35:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 莫selina 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 让我中一个面上吧! +5 大萍1987 2026-08-20 5/250 2026-08-20 20:46 by RRRKKK
[基金申请] 言之凿凿,明天官网会通知放榜? +15 医学老男孩 2026-08-20 17/850 2026-08-20 20:35 by CaoLG
[基金申请] 时间戳又变了 +12 wuchongjun 2026-08-20 17/850 2026-08-20 20:05 by spzhao
[基金申请] 今天系统多次维护,明天很可能放榜! +10 zju2000 2026-08-16 11/550 2026-08-20 20:02 by cl479861084
[基金申请] 放榜前的不淡定 40+3 snowwithsea 2026-08-19 12/600 2026-08-20 15:56 by lfy8008
[基金申请] 今天基金会出结果吗?20260819 +14 kkkl_v 2026-08-19 15/750 2026-08-20 08:21 by sunzitan
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
信息提示
请填处理意见