解抛物型微分方程c语言程序,偏微分方程的matlab代码_偏微分方程matlab程序,matlab偏微分代码...

matlab偏微分方程含源代码，大家可以下载运行，里面的都能用

Pie Eais ae h t lat BLep

话k气飞“目基回

古班比比,一生热件身的些懂

1身定

回量1

古典显式格式稳定情况

于

吉氏,一生售挥导方栏的解里

点丁

空真

ecin.com

2、古典隐式格式求解抛物型偏微分方程(一维热传导方程)

function [U x t=PDEParabolicClassicalImplicit(uX, u T, phi, psil, psi2, M, N, C)

%古典隐式格式求解抛物型偏微分方程

"o[u x tPDEParabolicClassicallmplicit(ux, uT phi, psil psi, M, N,C)

%方程:ut=C*uxx0←=x

%初值条件:u(x.0)=phi(x

%边值条件;uOt)=pi.u(ux,t)=psi2(t

%输出参数:U-解矩阵,第一行表示初值,第一列和最后一列表示边值,第二行表示第2层

x-空间变量

t时间变量

%输入参数:ux空间变量x的取值上限

uT-时间变量t的取值上限

phi-初值条件,定义为内联函数

psil边值条件,定义为内联函数

psi2-边值条件,定义为内联函数

M沿x轴的等分区间数

N-沿t轴的等分区间数

C系数,默认情況下C=1

%应用举例

%uX=1;uT=0.2M=50N=50:C=1:

ophi inline( sin(pi*x)): psil -inline(O): psi inline(0)

oU x t-PDEParabolicClassicalImplicit(uX, uT, phi, psi l, psi, M, N,C)

%设置参数C的默认值

if margin

C=1:

ch

%计算步长

dx=uXM"x的步长

dt=uIN%的步长

x(O: M)*dx,

t(O: N)dt

rC* du dx/dx:%步长比

docin. com

Diag= zeros(1,M-1)矩阵的对角线元素

Low=