找回密码
 注册
Simdroid-非首页
查看: 128|回复: 3

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

[复制链接]
发表于 2010-8-25 14:15:12 | 显示全部楼层 |阅读模式 来自 山东烟台
我想解这个方程组,有点棘手,借这个问题好像只能使用有限差分法
z代表长度段,总长为mt代表时间段。

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

想把方程组进行离散化,然后用网格法进行zt方向上的求解。可是问题是,我查找了关于解偏微分方程的书,大部分是解
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

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

×
发表于 2010-8-25 14:29:39 | 显示全部楼层 来自 黑龙江哈尔滨
Simdroid开发平台
n是什么?
回复 不支持

使用道具 举报

发表于 2010-8-26 16:34:30 | 显示全部楼层 来自 陕西西安
你不用化成0啊,用差分法做
你去看计算方法或者去弄一本matlab数值解方面的书,边界是0其实比较好弄
回复 不支持

使用道具 举报

 楼主| 发表于 2010-9-15 14:32:53 | 显示全部楼层 来自 山东烟台
你不用化成0啊,用差分法做
你去看计算方法或者去弄一本matlab数值解方面的书,边界是0其实比较好弄
wayne190214 发表于 2010-8-26 16:34

是啊,现在的关键是等式右边的u1*u2+u2*u3
,怎么处理,我发的帖说明的是现在的计算方法的书讲的全是等式右边=0的情况,我这种情况怎么做?望指点
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Simapps系列直播

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-10-6 19:20 , Processed in 0.042579 second(s), 13 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表