linshaojie 发表于 2010-8-25 14:15:12

解这个非线性方程组,有点棘手

我想解这个方程组,有点棘手,借这个问题好像只能使用有限差分法
z代表长度段,总长为m。t代表时间段。

n(du1/dt)+(du1/dt)=u1*u2+u2*u3
n(du2/dt)+(du2/dt)=u1*u3+u1*u2
n(du3/dt)+(du3/dt)=u1*u3+u1*u2

设初始值为
u1(0,t)= u1(m,t)=0
u2(0,t)= u2(m,t)=0
u3(0,t)= u3(m,t)=0

想把方程组进行离散化,然后用网格法进行z、t方向上的求解。可是问题是,我查找了关于解偏微分方程的书,大部分是解
n(du1/dt)+(du1/dt)=0,也就是后面没有u1*u2+u2*u3这一项。
我给出的方程组是双曲型方程组,如何才能把我给的方程组转化为求解n(du1/dt)+(du1/dt)=0,或者说,
怎么求解n(du1/dt)+(du1/dt)=u1*u2+u2*u3?

messenger 发表于 2010-8-25 14:29:39

n是什么?

wayne190214 发表于 2010-8-26 16:34:30

你不用化成0啊,用差分法做
你去看计算方法或者去弄一本matlab数值解方面的书,边界是0其实比较好弄

linshaojie 发表于 2010-9-15 14:32:53

你不用化成0啊,用差分法做
你去看计算方法或者去弄一本matlab数值解方面的书,边界是0其实比较好弄
wayne190214 发表于 2010-8-26 16:34 http://forum.simwe.com/images/common/back.gif
是啊,现在的关键是等式右边的u1*u2+u2*u3
,怎么处理,我发的帖说明的是现在的计算方法的书讲的全是等式右边=0的情况,我这种情况怎么做?望指点
页: [1]
查看完整版本: 解这个非线性方程组,有点棘手