-
大小: 166KB文件類型: .rar金幣: 2下載: 0 次發(fā)布日期: 2022-09-07
- 語言: Matlab
- 標(biāo)簽:
資源簡介
以某偏微分方程為例,應(yīng)用多重網(wǎng)格法進(jìn)行了求解,編寫了matlab程序,注釋詳細(xì),并對計算結(jié)果進(jìn)行了分析。

代碼片段和文件信息
function?u_fine=interpolate(u_coarse)
%Copyright@Peng?Yj?NWPU2016.3.3
%由粗網(wǎng)格延拓到細(xì)網(wǎng)格,采用線性插值算子
%將下一層近似解的修正量插值到上一層
%輸入項:u_coarse(粗網(wǎng)格列向量)
%輸出項:u_fine(密網(wǎng)格列向量)
%示例如下:
%下一層共9個點,首尾函數(shù)值已知,去除,共7個點
%上一層共17個點,首尾函數(shù)值已知,去除,共15個點
%下一層首尾的修正量為0
%???1????1____1____1____1____1____1____1____1____1____1____1____1____1____1____1????1
%???????/?\???|???/?\???|???/?\???|???/?\???|???/?\???|???/?\???|???/?\???|???/?\
%????0?/???\??|??/???\??|??/???\??|1?/???\?1|??/???\??|??/???\??|??/???\??|??/???\
%?????/?????\?|?/?????\?|?/?????\?|?/?????\?|?/?????\?|?/?????\?|?/?????\?|?/????0\
%????/??1/2??\|/???????\|/???????\|/??1/2??\|/???????\|/???????\|/???????\|/???????\
%???1?????????1_________1_________1_________1_________1_________1_________1?????????1??????????
%
u_coarse_size=size(u_coarse1);?????%粗網(wǎng)格向量的長度,則細(xì)網(wǎng)格向量的長度為2*u_coarse_size+1
u_fine(11)=u_coarse(11)/2;??%首項
u_fine(2*u_coarse_size+11)=u_coarse(u_coarse_size1)/2;??%尾項
for?i=1:u_coarse_size
????u_fine(2*i1)=u_coarse(i1);??%偶數(shù)項
end
for?i=3:2:2*u_coarse_size-1
????u_fine(i1)=u_coarse((i-1)/21)/2+u_coarse((i+1)/21)/2;??%奇數(shù)項(除首尾項外)
end
end
?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件???????1277??2016-03-04?19:20??多重網(wǎng)格法實例及matlab程序\interpolate.m
?????文件???????7041??2016-03-04?19:10??多重網(wǎng)格法實例及matlab程序\main.m
?????文件?????186114??2016-03-04?19:27??多重網(wǎng)格法實例及matlab程序\Readme.pdf
?????文件???????1080??2016-03-04?19:21??多重網(wǎng)格法實例及matlab程序\restrict.m
?????文件????????594??2016-03-03?15:21??多重網(wǎng)格法實例及matlab程序\SOR.m
?????目錄??????????0??2016-03-04?19:28??多重網(wǎng)格法實例及matlab程序
-----------?---------??----------?-----??----
???????????????196106????????????????????6
評論
共有 條評論