| 查看: 3668 | 回复: 1 | |||
[交流]
【求助】用matlab解含第二类边界条件的偏微分方程 已有1人参与
|
|
一个简单的抛物线型偏微分方程，求u（x，t） u_t 表示对t的一阶导数，u_x为对x的一阶导数，u_xx表示对x的二阶导数 方程 u_t =u_xx （1） 边界条件 x=0,t>0，u_x=0； （2） x=1,t>0, u(1,t)=0; （3） 初始条件 t=0, u(x,0)=1; （4） 我用有限差分法，主要是第二类边界条件（2）式怎么在matlab里面利用 我大致解了一下，但觉得结果有问题，而且我认为问题是出在（2）式上。真心求教这个问题怎么解，谢谢虫友！ 附：我写的程序-------------------------------------------------------------------------------------- function [u,x,t]=model2(xf,T) % solve u_xx=u_t for 0<=x<=xf,0<=t<=T % Initial condition:u(x,0)=1 % boundary conditon:u(1,t)=0;u_x(0,t)=0 dx=xf/50;x=[0:50]'*dx; dt=T/10000;t=[0:10000]*dt; for i=1:51,u(i,1)=1;end % Initial condition:u(x,0)=1 for n=1:10001,u(51,n)=0;end % boundary conditon:u(1,t)=0 r=dt/dx/dx;r1=1-2*r; for k=1:10000 u(1,k+1)=2*r*u(2,k)+r1*u(1,k);end % boundary conditon:u_x(0,t)=0 ? for k=1:10000 for i=2:50 u(i,k+1)=r1*u(i,k)+r*(u(i-1,k)+u(i+1,k));% solve u_xx=u_t for 0<=x<=xf,0<=t<=T end end -------------------------------------------------------------------------------------- 命令窗口输入： xf=1; T=1; [u,x,t]=model2(xf,T); figure(1),clf,mesh(t,x,u); ============================================== [ Last edited by zju-hjm on 2009-12-17 at 16:28 ] |
回复此楼» 猜你喜欢
各位大佬，求一份最新Wien2k版本作为学习用途
已经有0人回复
-
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后
已经有523人回复
-
物理学I论文润色/翻译怎么收费?
已经有297人回复
-
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有130人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 椭圆形偏微分方程标准化后该如何解？
已经有6人回复
- Matlab求解偏微分方程组
已经有22人回复
- matlab的函数调用问题
已经有3人回复
- matlab数值求解边界条件微分方程组
已经有7人回复
- 一个偏微分方程组（包括边界条件）如何求解得到自振频率和振型？
已经有4人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
- 【求助】使用MATLAB怎么实现拟合力场参数的程序？【已解决】
已经有15人回复
- 【求助】偏微分方程的基本解
已经有5人回复
- 【求助】偏微分方程解法
已经有19人回复
2楼2010-05-05 16:34:41
