matlab抛物偏微分方程,用matlab解抛物型方程 求解二维抛物线型偏微分方程matlab程序...

求解抛物线型偏微分方程matlab程序

82064956f575c906587622efe6967ed2.png

用matlab求解抛物型方程,急啊!!用最简隐格式(向后差分格式)求解抛物型方程 用matlab求解抛物型方程,急啊!用最简隐格式(向后差分格式)求解抛物型方程 要用matlab求解,但是不能用里面的求微分方程的工具来求解,就是自己编程序,要有图示的啊,。

ffdf4834ce7622cd2556db8783db40e8.png

跪求MATLAB解抛物型偏微分方程的程序 1,不一定有效果,因为pdetool具体编程是不知道的,如果解决小问题两者的结果一样说明不了什麽问题,尤其对于偏微分方程。2有限元的边界必须固定,从数理方程上讲静态有限元问题就是边值问题,如果边界变化的话,初始一下别的专业有限元软件,比如anasys,adima等。

66202f9399115f51ba5f9c3c083fda75.png

如何用Matlab解偏微分方程组,该方程组由两个抛物型偏微分方程组成

求解二维抛物线型偏微分方程matlab程序 function[u,x,y,t]=TDE(A,D,T,ixy0,bxyt,Mx,My,N)解方程 u_t=c(u_xx+u_yy)for D(1)(2),D(3)(4),0初值:u(x,y,0)=ixy0(x,y)边界条件:u(x,y,t)=bxyt(x,y,t)for(x,y)cBMx/My:x轴和y轴的等分段数N:t 轴的等分段数dx=(D(2)-D(1))/Mx;x=D(1)+[0:Mx]*dx;dy=(D(4)-D(3))/My;y=D(3)+[0:My]'*dy;dt=T/N;t=[0:N]*dt;初始化ufor i=1:Mx+1for j=1:My+1u(i,j)=ixy0(x(i),y(j));endendrx=A*dt/(dx*dx);rx1=1+2*rx;rx2=1-2*rx;ry=A*dt/(dy*dy);ry1=1+2*ry;ry2=1-2*ry;for i=1:Mx-1%(11.2.21a)P(i,i)=ry1;if i>;1P(i-1,i)=-ry;P(i,i-1)=-ry;endendfor j=1:My-1%(11.2.21b)Q(j,j)=rx1;if j>;1Q(j-1,j)=-rx;Q(j,j-1)=-rx;endendfor k=1:Nu_1=u;t=k*dt;for i=1:Mx+1%边界条件u(i,1)=feval(bxyt,x(i),y(1),t);u(i,My+1)=feval(bxyt,x(i),y(My+1),t);endfor j=1:My+1u(1,j)=feval(bxyt,x(1),y(j),t);u(Mx+1,j)=feval(bxyt,x(Mx+1),y(j),t);endif mod(k,2)=0for i=2:Mxj=2:My;bx=[ry*u(i,1)zeros(1,Mx-3)ry*u(i,My+1)]+rx*(u_1(i-1,j)+u_1(i+1,j))+rx2*u_1(i,j);u(i,j)=linsolve(P,bx');(11.2.21a)endelsefor j=2:Myi=2:Mx;by=[rx*u(1,j);zeros(My-3,1);rx*u。

请求各位大虾,帮忙指点一下,如何用MATLAB解一维抛物型方程?谢谢! 可以用命令pdetool进入pde工具箱来用数值方法求解

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值