ZKX's LAB

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

2020-10-09知识12

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

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

偏微分方程的分类 二阶偏微分方程的一般形式为A*Uxx+2*B*Uxy+C*Uyy+D*Ux+E*Uy+F*U=0其特征方程为A*(dy)^2-2*B*dx*dy+C*(dx)^2=0若在某域内B^2-A*C0则在此域内称为双曲形方程其实主要是按特征方程的曲线类型分的注:Uxx表示U对x求二阶.

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

抛物型偏微分方程的定解问题 为了确定一个具体的热传导过程,除了列出方程(1)以外,还必须知道物体Ω的初始温度(初始条件)和在它的边界嬠Ω上所受到的外界的影响(边界条件)。初始条件:边界条件,最通常的形式有三类。第一边界条件(或称狄利克雷条件):即表面温度为已知函数。第二边界条件(或称诺伊曼条件):式中n是Ω的外法向,即通过表面的热量已知。第三边界条件(或称罗宾条件):式中α≥0;即物体表面给定热交换条件。除了以上三类边界条件外还可以在边界嬠Ω上给定其他形式的边界条件,如斜微商条件、混合边界条件等。方程(1)连同初始条件(2)以及边界条件(3)、(4)、(5)中的任意一个一起构成了一个定解问题,根据边界条件的不同形式,分别称为第一、二、三边值问题,统称为热传导方程的初边值问题或混合问题。若Ω呏R3,则由方程(1)和初始条件(2)构成的定解问题称为热传导方程的初值问题或柯西问题。

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

抛物型偏微分方程的介绍 简称抛物型方程,一类重要的偏微分方程。热传导方程是最简单的一种抛物型方程。热传导方程 研究热传导过程的一个简单数学模型。根据热量守恒定律和傅里叶热传导实验定律导致热传导方程

抛物型偏微分方程数值解怎么给出第三类边界条件

抛物型偏微分方程数值解怎么给出第三类边界条件 沿外法线的导数与边界内外函数值之差成正比dy/dn=k(y-f)

急求!!! 大学数学,用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 轴的。

#偏微分方程#热传导#边界条件

随机阅读

qrcode
访问手机版