基于matlab的偏微分方程差分解法(更新版)
【正文】 figure(2)。 U(i,1)=sin(pi*h*(i1))。 % x方向上的步長k=b/(m1)。如果差分格式有解，且當網(wǎng)格無限變小時其解收斂于原微分方程定解問題的解，則差分格式的解就作為原問題的近似解（數(shù)值解）?；贛ATLAB的偏微分方程差分解法學院：核工程與地球物理學院專業(yè)：勘查技術(shù)與工程班級：1120203學號：姓名:2014/6/11在科學技術(shù)各領域中，有很多問題都可以歸結(jié)為偏微分方程問題。它的基本思想是：先對求解區(qū)域作網(wǎng)格剖分，將自變量的連續(xù)變化區(qū)域用有限離散點（網(wǎng)格點）集代替；將問題中出現(xiàn)的連續(xù)變量的函數(shù)用定義在網(wǎng)格點上離散變量的函數(shù)代替；通過用網(wǎng)格點上函數(shù)的差商代替導數(shù)，將含連續(xù)變量的偏微分方程定解問題化成只含有限個未知數(shù)的代數(shù)方程組（稱為差分格式）。Alembert方法function U=hyperbolic(a,b,c,n,m)% a為x的取值范圍% b為t的取值范圍% c為系數(shù)% n為x方向上的節(jié)點數(shù)% m為t上的節(jié)點數(shù)h=a/(n1)。for i=2:(n1)。 surf(U)。2 拋物線方程中熱傳導方程的有限差分解法。U=zeros(n,m)。 surf(U)。這里使用零法線導數(shù)條件：對于熱傳導而言，這表示邊是熱絕緣的而且經(jīng)過邊的熱通量為零，從而得到：得出點的Laplace差分方程為：利用差分方程：在（4）上式使用（5）式這個近似值，這結(jié)果為：這個公式講函數(shù)聯(lián)系起來。ave=(a*(20+180)+b*(80+0))/(2*a+2*b)。U(n,1)=(U(n1,1)+U(n,2))/2。 U(i,j)=U(i,j)+relx。figure(2)。3 附錄雙曲型及拋物線方程MATLAB語言提供的pdepe()函數(shù)程序。,39。)xlabel(39。% 目標PDE函數(shù)function [c,f,s]=pdefun (x,t,u,du)%p1414c=1。pb=ub。11
環(huán)評公示相關推薦
鄂ICP備17016276號-1