ZKX's LAB

抛物型偏微分方程的格林函数 一维抛物型微分方程

2020-09-27知识10

急求!!! 大学数学,用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求解,但是不能用里面的求微分方程的工具来求解,就是自己编程序,要有图示的啊,。

#偏微分方程#微分

随机阅读

qrcode
访问手机版