24小时热门版块排行榜    

查看: 1068  |  回复: 3
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

莫selina

新虫 (初入文坛)

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

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

未命名.jpg
http://kuai.xunlei.com/d/GXKYUGYYBUXX?p=130497
回复此楼
已阅   回复此楼   关注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的回帖
查看全部 4 个回答

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 6/300 2026-08-23 14:45 by KM5EcsNQRBPn
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +3 process2012 2026-08-23 4/200 2026-08-23 12:54 by huixian257
[基金申请] 让我中一个面上吧! +12 大萍1987 2026-08-20 14/700 2026-08-23 10:30 by wrm
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 07:41 by OEbVnUOu01ol
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 10/500 2026-08-23 07:17 by OEbVnUOu01ol
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 04:04 by OEbVnUOu01ol
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 03:52 by OEbVnUOu01ol
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 7/350 2026-08-23 03:43 by OEbVnUOu01ol
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 DpSrDtM079iu 2026-08-22 3/150 2026-08-23 00:55 by OEbVnUOu01ol
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 建议基金发布提前给出明确的时间点 +10 kulium 2026-08-21 13/650 2026-08-21 21:58 by alongwaytogo
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 今天系统多次维护,明天很可能放榜! +10 zju2000 2026-08-16 11/550 2026-08-20 20:02 by cl479861084
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
信息提示
请填处理意见