[原創(chuàng)]偏微分方程數(shù)值解法的MATLAB源碼【更新完畢】

說明：由于偏微分的程序都比較長，比其他的算法稍復雜一些，所以另開一貼，專門上傳偏微分的程序謝謝大家的支持！

其他的數(shù)值算法見：

..//Announce/Announce.asp?BoardID=209&id=8245004

1、古典顯式格式求解拋物型偏微分方程(一維熱傳導方程)

function [U x t]=PDEParabolicClassicalExplicit(uX,uT,phi,psi1,psi2,M,N,C)

%古典顯式格式求解拋物型偏微分方程

%[U x t]=PDEParabolicClassicalExplicit(uX,uT,phi,psi1,psi2,M,N,C)

%

%方程：u_t=C*u_xx 0 <= x <= uX,0 <= t <= uT

%初值條件：u(x,0)=phi(x)

%邊值條件：u(0,t)=psi1(t), u(uX,t)=psi2(t)

%

%輸出參數(shù)：U -解矩陣，第一行表示初值，第一列和最后一列表示邊值，第二行表示第2層……

% x -空間變量

% t -時間變量

%輸入?yún)?shù)：uX -空間變量x的取值上限

% uT -時間變量t的取值上限

% phi -初值條件，定義為內(nèi)聯(lián)函數(shù)

% psi1 -邊值條件，定義為內(nèi)聯(lián)函數(shù)

% psi2 -邊值條件，定義為內(nèi)聯(lián)函數(shù)

% M -沿x軸的等分區(qū)間數(shù)

% N -沿t軸的等分區(qū)間數(shù)

% C -系數(shù)，默認情況下C=1

%

%應用舉例：

%uX=1;uT=0.2;M=15;N=100;C=1;

%phi=inline('sin(pi*x)');psi1=inline('0');psi2=inline('0');

%[U x t]=PDEParabolicClassicalExplicit(uX,uT,phi,psi1,psi2,M,N,C);

%設置參數(shù)C的默認值

if nargin==7

C=1;

end

%計算步長

dx=uX/M;%x的步長

dt=uT/N;%t的步長

總結(jié)

以上是生活随笔為你收集整理的偏微分方程数值解的matlab程序,偏微分方程数值解法的MATLAB源码的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。