高手，请问如何用有限差分法求解抛物线型的偏微分方程，用matlab，能告诉我具体的编程程序了，万分感谢了~~~急 ？X/？t=？/？z(Deff？X/？z)；0<；z With the following conditions： Initial：t=0；0<；z；X=X0 Boundary：t>；0；z=0；？X/？z=0 t>；0；z=L；X=Xeq 这是我要求的偏微分方程，谢谢谢谢了
急求！！！ 大学数学，用matlab解决问题，题目是一维抛物型偏微分方程差分解法 显式前向欧拉法源程序：function[u，x，t]=EF_Euler(A，xf，T，it0，bx0，bxf，M，N)解方程 A u_xx=u_t，0，0初值：u(x，0)=it0(x)边界条件：u(0，t)=bx0(t)，u(xf，t)=bxf(t)M：x 轴的等分段数N：t 轴的等分段数dx=xf/M；x=[0：M]*dx；dt=T/N；t=[0：N]'*dt；for i=1：M+1u(i，1)=it0(x(i))；endfor j=1：N+1u([1 M+1]，j)=[bx0(t(j))；bxf(t(j))]；endr=A*dt/dx/dx，r1=1-2*r；if(r>；0.5)disp('r>；0.5，unstability')；endfor j=1：Nfor i=2：Mu(i，j+1)=r*(u(i+1，j)+u(i-1，j))+r1*u(i，j)；(9.2.3)endendu=u'；在MATLAB中编写脚本文件：A=0.5；方程系数it0=inline('sin(pi*x)'，'x')；初始条件bx0=inline('0')；bxf=inline('0')；边界条件xf=2；M=80；T=0.1；N=100；[u1，x，t]=EF_Euler(A，xf，T，it0，bx0，bxf，M，N)；figure(1)，clf，mesh(u1)xlabel('x')ylabel('t')zlabel('U')title('r>；0.5')M=50；[u1，x，t]=EF_Euler(A，xf，T，it0，bx0，bxf，M，N)；figure(2)，clf，mesh(u1)xlabel('x')ylabel('t')zlabel('U')title('r)隐式后向欧拉法源程序：function[u，x，t]=IB_Euler(A，xf，T，it0，bx0，bxf，M，N)解方程 A1 u_xx=u_t，0，0初值：u(x，0)=it0(x)边界条件：u(0，t)=bx0(t)，u(xf，t)=bxf(t)M：x 轴的。
如何用matlab解二维的非线性偏微分方程组， 其中每个方程是抛物线型的 如何用matlab解二维的非线性偏微分方程组，其中每个方程是抛物线型的 MATLAB提供了两种方法解决PDE问题：一是pdepe()函数，它可以求解一般的PDEs，据用较大的通用性，但只。
求解二维抛物线型偏微分方程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解抛物型偏微分方程的程序 1，不一定有效果，因为pdetool具体编程是不知道的，如果解决小问题两者的结果一样说明不了什麽问题，尤其对于偏微分方程。2有限元的边界必须固定，从数理方程上讲静态有限元问题就是边值问题，如果边界变化的话，初始一下别的专业有限元软件，比如anasys，adima等。
