function?u_fine=interpolate（u_coarse）
%Copyright@Peng?Yj?NWPU2016.3.3
%由粗網格延拓到細網格，采用線性插值算子
%將下一層近似解的修正量插值到上一層

%輸入項：u_coarse（粗網格列向量）
%輸出項：u_fine（密網格列向量）

%示例如下：
%下一層共9個點，首尾函數值已知，去除，共7個點
%上一層共17個點，首尾函數值已知，去除，共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）;?????%粗網格向量的長度，則細網格向量的長度為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）;??%偶數項
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;??%奇數項（除首尾項外）
end
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件??????186114??2016-03-04?19:27??多重網格法求解微分方程-matlab\Readme.pdf
?????文件?????????594??2016-03-03?15:21??多重網格法求解微分方程-matlab\SOR.m
?????文件????????1277??2016-03-04?19:20??多重網格法求解微分方程-matlab\interpolate.m
?????文件????????7052??2017-05-12?10:39??多重網格法求解微分方程-matlab\main.m
?????文件????????2234??2017-05-15?15:20??多重網格法求解微分方程-matlab\mytestV.m
?????文件????????1080??2016-03-04?19:21??多重網格法求解微分方程-matlab\restrict.m
?????目錄???????????0??2017-09-25?10:43??多重網格法求解微分方程-matlab\