24小时热门版块排行榜    

查看: 708  |  回复: 8
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

dxyhn1979

银虫 (小有名气)

[交流] 【求助】差分法解微分方程

请问:附件中的微分方程可以用中心差分法来解么?能否提供以下解得思路啊,高手们?要是不考虑步长对稳定性的影响,用显式还是隐式的差分格式?
谢了啊!
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

senlia

木虫 (小有名气)

★ ★ ★
nono2009(金币+2,VIP+0):谢谢专家!好young的特聘专家:) 9-23 21:12
dxyhn1979(金币+1,VIP+0):说的很详细了,我试试。感谢指导了。 9-23 21:16
引用回帖:
Originally posted by dxyhn1979 at 2009-9-23 18:30:
"线性方程组左端的系数矩阵"中的“系数矩阵”指的是什么?

你把方程进行差分离散 考虑所有的初、边值条件 最后不是可以得到一个大致形式为A*v=f线性方程组么?(可能比这个复杂些 如果原方程有非线性项之类的话)其中v和f可能和时间层有关 而A和时间层无关 这个A就是系数矩阵啊 你预先把A放到方右端去 变成 v(t_k)=B*f(t_k) 其中B=A^(-1) 那么每一步(每一时间层)只要计算出f(t_k) 再进行“矩阵和向量的乘法”就可以了 这样的运行速度不会慢多少的
书到用时方恨少, 天涯何处无芳草.
7楼2009-09-23 21:06:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

dxyhn1979

银虫 (小有名气)

以上方程中,Cm,  d  和  roua  都是常数,V是电压,i(V)是电流,是电压的函数
2楼2009-09-22 23:31:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liang_2004

木虫 (正式写手)

始作俑者

★ ★ ★ ★
nono2009(金币+1,VIP+0):谢谢交流!欢迎常来。 9-23 06:08
wuguocheng(金币+2,VIP+0): 回答的不错,鼓励一下新虫 9-23 07:34
dxyhn1979(金币+1,VIP+0):谢谢回答啊 9-23 18:26
电压对时间的偏导用欧拉公式离散,对x的二阶偏导可以用最近的四点离散,为了更均匀的话,可以用九点格式
我认为这样比较好做
如果你想用中心差分提高误差阶的话,也可以。至于隐式和显式来说,隐式更稳定,但计算速度慢,显式计算速度快,有时候不稳定
走着,看看再说
3楼2009-09-22 23:43:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxyhn1979

银虫 (小有名气)

隐式是否一定稳定呢?显式要达到稳定对时间步长和空间步长有什么要求呢?比如理论上二者要满足什么样的关系?谢谢回答
4楼2009-09-23 10:18:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见