數(shù)值方法期末作業(yè)有限體積WENO+北大高彬


高同學講得很細也很清楚,報告時的自信從容、有條不紊、對所做格式熟悉到信手拈來??紤]到班里的很多同學對偏微分方程數(shù)值數(shù)值解不熟悉,建議高同學能夠把自己代碼的注釋進一步的完善,使得其他同學以后再遇到這樣的問題的時候能夠方便地參考高同學的代碼。
---下面是我看了高同學的代碼之后拋磚引玉,增加的一些注釋。 ----

在北太天元上運行得到數(shù)值結(jié)果如上圖所示。
另外,可以做一些小改可以使得計算得更快,
把 WENO函數(shù) 修改的可以處理向量參數(shù),
實際山就是把 ^2 改成 .^2, 把 * 改成 .*
然后把下面的for循環(huán)
%for i=1:N %重構(gòu)
%uL(i+1)=WENO(uL2(i),uL1(i),u(i),uR1(i),uR2(i));
%uR(i+1)=WENO(uR3(i),uR2(i),uR1(i),u(i),uL1(i));
%end
改成
uL(2:end)=WENO(uL2,uL1,u,uR1,uR2);
uR(2:end)=WENO(uR3,uR2,uR1,u,uL1);
標簽: