ZKX's LAB

急求!!! 大学数学,用matlab解决问题,题目是一维抛物型偏微分方程差分解法 求一维抛物型偏微分方程组

2020-11-26知识6

急求!!! 大学数学,用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解决问题,题目是一维抛物型偏微分方程差分解法 求一维抛物型偏微分方程组

抛物型偏微分方程的反应扩散

你见过哪些有趣的偏微分方程组? 说说我一直在做的一个偏微分方程系统吧,该模型为生物领域的趋化性(chemotaxis)模型,也叫Keller-Segel模…

偏微分问题 椭圆型偏微分方程:二维平面稳定场方程,如稳定浓度分布,稳定温度分布,静电场方程,无旋稳恒电流场方程,无旋稳恒流动方程等抛物型偏微分方程:一维输运方程,如扩散方程,热传导方程等双曲型偏微分方程:一维波动方程,如弦振动方程,杆振动方程,电报方程等它们是分别描述二维平面稳定场,一维输运,一维波动问题的方程

帮忙求解以下偏微分方程(急) 好难的,没学。

抛物型偏微分方程的格林函数 基本解是点热源的影响函数。如果在t=0时刻在(ξ,η,ζ)处给定单位点热源,即u0(x,y,z,0)=δ(ξ,η,ζ)(δ是狄喇克函数),则当t>;0时由它引起的在全空间 R3的温度分布(即热。

随机阅读

qrcode
访问手机版