| 查看: 1150 | 回复: 3 | ||
[求助]
这个方程怎样实现差分？（偶是新人，目前只有两个金币，不要嫌少啊）
|
|
这样一个方程怎样实现差分？ 在n点t时刻，方程如图： 未命名.jpg http://kuai.xunlei.com/d/GXKYUGYYBUXX?p=130497 |
回复此楼» 猜你喜欢
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有186人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有120人回复
-
深圳大学2027级光电信息工程/物理学研究生推免
已经有70人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有5人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有257人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 这样一个微分方程怎么实现差分？？？？？？
已经有7人回复
- 这个方程该怎么采用有限差分法求解呢？
已经有4人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 求助一个差分方程的MATLAB实现
已经有3人回复
- 【求助】高维PDE的求解
已经有11人回复
- 【求助】多重网格预处理子如何构造？
已经有5人回复
wangww2011
木虫 (著名写手)
- 应助: 11 (小学生)
- 金币: 4023.1
- 散金: 2709
- 红花: 18
- 沙发: 1
- 帖子: 1915
- 在线: 1537.1小时
- 虫号: 772953
- 注册: 2009-05-17
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2012-09-03 20:52:01
3楼2012-09-04 14:32:35
wangww2011
木虫 (著名写手)
- 应助: 11 (小学生)
- 金币: 4023.1
- 散金: 2709
- 红花: 18
- 沙发: 1
- 帖子: 1915
- 在线: 1537.1小时
- 虫号: 772953
- 注册: 2009-05-17
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★ ★
csgt0: 金币+2, 应助指数+1, 多谢应助 2012-09-04 22:49:28
csgt0: 金币+2, 应助指数+1, 多谢应助 2012-09-04 22:49:28
|
其实查分本身非常容易，就是因为太容易，所以有可能求解你这个方程是不稳定的，所以也不是很建议你直接用，不过你是不是应该提供更多的信息？ 你给出T为常数，其实这个常数是n的函数吧？譬如第一辆车T应该为零的，否则整个方程组没法开始。从第二辆车开始，加速度可以用你给的式子计算，但是第一辆车的的加速度就不能了，所以在你的模拟中，第一辆车的速度其实应该是你提前假定好的（当然可以变化，但你是不是应该提前给出个函数呢），然后也作为方程的初始条件对吧？ 你把加速度写成dv/dt是应该的，也是求解微分方程组的第一步。 至于你的这个方程组，假设有n辆车，你需要求解的是一个有n-1个方程的方程组，本身很容易求解，如果你懂点计算机的话，你可以用matlab的ode45，或者在mathematica里直接敲出式子来给出初始条件就可以直接得出结果了。或者如果你用python，直接用scipy.integrate.odeint就可以了，很简单的，而且貌似前段时间有本一个家伙写的python数值计算的中文教程，你拿来读一读（当字典来查就行了），你这个问题也就几十行的代码吧，加上漂亮的绘图，应该也不长的。 |
4楼2012-09-04 17:35:00
