基于matlab語(yǔ)言求偏微分方程畢業(yè)設(shè)計(jì)(doc畢業(yè)設(shè)計(jì)論文)(編輯修改稿)
2024-07-15 14:48 本頁(yè)面
【文章內(nèi)容簡(jiǎn)介】 為1，演化時(shí)間為130。動(dòng)畫(huà)演示了波包向勢(shì)壘行進(jìn)，在勢(shì)壘上發(fā)生透射和反射，最后分為透射波和反射波二部分向不同的方向運(yùn)動(dòng)。編寫(xiě)MATLAB程序如下：clear all %清除NPTS=220。 %設(shè)定格點(diǎn)范圍sigma=10。 %設(shè)置奇異值頻域圖k0=。 %設(shè)置平均波數(shù)x0=40。 %設(shè)置中心位置time=130。 %設(shè)定時(shí)間v(NPTS)=0。 %設(shè)定初值v(105:115)=+。 %替換A=diag(2+2i+v)+diag(ones(NPTS1,1),1)+diag(ones(NPTS1,1),1)。 %建立對(duì)角陣運(yùn)算PHI0(1)=0。PHI0(NPTS)=0。 %設(shè)定初值for x=2:NPTS1。 %循環(huán)PHI0(x)=exp(*x)*exp((xx0)^2*log10(2)/sigma^2)。 End %結(jié)束PHI(:,1)=PHI0.39。 %for y=2:timeCHI(:,y)=4i*(A\PHI(:,y1))。PHI(:,y)=CHI(:,y)PHI(:,y1)。endmo=moviein(time)。for j=1:timeplot([1:NPTS],abs(PHI(:,j)).^2/norm(PHI(:,j)).^2,[1:NPTS],v/2)。%輸出mo(:,j)=getframe。 endmovie(mo) %播放圖NPTS=220。 %設(shè)定格點(diǎn)范圍通過(guò)上述程序得到運(yùn)行結(jié)果如圖2所示 圖2 某一時(shí)刻透射波與反射波運(yùn)動(dòng)圖像差分法與松弛法解橢圓型方程利用二階導(dǎo)數(shù)的中心差分公式得到二維拉普拉斯方程 的差分公式為 （21） 記得顯式 （22）穩(wěn)定條件是. 問(wèn)題3： 差分法解平面溫度場(chǎng),設(shè)正方形的一邊溫度為10度,其余各邊的溫度為零度,求穩(wěn)定的溫度場(chǎng)。按照上面的差分公式，任意設(shè)定內(nèi)部各點(diǎn)的溫度為零，就停止計(jì)算。編寫(xiě)MATLAB程序如下：clear all %清除u=zeros(100,100)。 %產(chǎn)生100行100列的全零矩陣 u(100,:)=10。 %取十個(gè)數(shù)uold=u+1。 %定義表達(dá)式 unew=u。 %把u賦給unewfor k=1:500 %循環(huán)語(yǔ)句if max(max(abs(uuold)))= %循環(huán)語(yǔ)句unew(2:99,2:99)=*(u(3:100,2:99)+u(1:98,2:99)+...u(2:99,3:100)+u(2:99,1:98))。 %疊加語(yǔ)句uold=u。 %賦值 u=unew。 %賦值End %結(jié)束endsurf(u) %繪制三維圖形運(yùn)行上述程序，結(jié)果如圖3所示圖3 差分法解平面溫度運(yùn)行結(jié)果 如果在差分公式中，隨時(shí)將上一步算得的格點(diǎn)上的值替代舊值，并且每次算出的新值也替換成新值與舊值的“組合”，則得到下列松弛法的計(jì)算公式,其中。這個(gè)公式可以用變分原理去證明。 （23）問(wèn)題4： 松弛法計(jì)算平面溫度場(chǎng),已知定解問(wèn)題如下 (24)為了
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1