24小时热门版块排行榜    

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

莫selina

新虫 (初入文坛)

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

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

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

wangww2011

木虫 (著名写手)

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

即使求差分,你的v_{n}和v_{n-1}是什么意思?是不同时刻的v吗还是同一时刻有几个的v?
2楼2012-09-03 20:52:01
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2026国自然函评费到账 +11 羊腰板 2026-08-21 11/550 2026-08-22 22:03 by lsbin3733
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 3/150 2026-08-22 22:03 by OEbVnUOu01ol
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-22 21:42 by OEbVnUOu01ol
[基金申请] 范进中举一文的中心思想 +3 zjuhero 2026-08-22 4/200 2026-08-22 20:02 by mycsru
[基金申请] 什么时候开奖? +7 CrisMessi 2026-08-18 8/400 2026-08-22 18:29 by 淀粉搬运工
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-21 4/200 2026-08-22 16:52 by sunzitan
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 建议基金发布提前给出明确的时间点 +10 kulium 2026-08-21 13/650 2026-08-21 21:58 by alongwaytogo
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 看来今天不会放榜了? +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 紫杉醇
[基金申请] 让我中一个面上吧! +11 大萍1987 2026-08-20 13/650 2026-08-21 14:48 by 20120902066
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
信息提示
请填处理意见