基于matlab的偏微分方程差分解法-免費閱讀
【正文】 % pb=ub。% f=du。)zlabel(39。,39。sol=pdepe(m,pdefun,pdeic,pdebc,x,t)。 h=。endU=flipud(U39。t=0。U(1:n,1)=20。對所有網(wǎng)格點應(yīng)用上式（8）直到。 通過MATLAB語言提供了pdepe()函數(shù)，可以直接求解一般偏微分方程(組)，它的調(diào)用格式為：sol=pdepe(m,pdefun,pdeic,pdebc,x,t)，具體程序見附錄得出的結(jié)果為： 結(jié)果對比通過編寫MATLAB的差分方程程序求取結(jié)果和MATLAB自帶函數(shù)求取結(jié)果進行對比，發(fā)現(xiàn)這兩種方法求得到的結(jié)果是非常相似的，差距不大證明程序編寫是成功的。% U(I,1)=4*(i1)*h4*((i1)*h)^2。題：其中初始條件為其中邊界條件為：解：第一步，分析并帶入（3）并編寫MATLAB求解程序如下：function U=forwdif(c1,c2,a,b,c,n,m)clch=a/(n1)。b=。 for i=2:(n1)。r22=r^2/2。167。隨著計算機技術(shù)的發(fā)展，采用數(shù)值計算方法，可以得到其數(shù)值解。近些年來，求解偏微分方程的數(shù)值方法取得進展，特別是有限差分區(qū)域分解算法，此類算法的特點是在內(nèi)邊界處設(shè)計不同于整體的格式， 將全局的隱式計算化為局部的分段隱式計算。1 雙曲型方程中波動方程的有限差分解法。s1=1r^2。 U(i,j)=s2*U(i,j1)+r2*(U(i+1,j1)+U(i1,j1))U(i,j2)。c=4n=11m=11執(zhí)行hyperbolic(1,4,11,11)；第三步：得出結(jié)果并畫圖入下167。k=b/(m1)。endfor j=2:mfor i=2:n1 U(I,j)=s*U(I,j1)+r*(U(i1,j1)+U(i+1,j1))。167。其中： 橢圓型方程中波動方程例題的差分解法結(jié)果及程序。U(1:n,m)=180。while((err
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1