clear?all;
clc;
tic;
addpath（genpath（?‘D:\已完成仿真\FRSR\重建初值‘））;%加入路徑
addpath（genpath（?‘D:\已完成仿真\FRSR\產(chǎn)生待重建的低分辨率圖像‘））;%加入路徑
number=16;
scale=4;
sd=[7.5035e-04;0.0034;0.0023;0.0039;0.0042;1.6262e-04;0.0044;0.0031];
sp=[0.0147;0.0188;0.0166;0.0058;0.0259;0.0088;0.0133;0.0245?];
noisesd=[sd;sp];
delta=[1.2459?0.2544;1.0818?1.5652;3.2187?0.1709;...
??????????1.1261?0.8884;3.7900?4.4700;0.8561?1.4863;...
??????????1.3039?0.0585;4.3434?0.7524;1.6179?0.9584;...
?????????2.4087?4.1360;3.9513?1.0876;2.4135?1.9259;...
?????????3.3602?2.5350;4.9883?4.2239;3.2744?0.6202;...
?????????3.0191?1.4393];
?%phi=rand（number1）*5;
phi=[2.1312;4.0317;1.2249;2.5983;1.4762;0.4953;0.1746;...
??????????3.7099;0.1464;4.2639;1.4812;2.3967;3.7673;4.3981;...
??????????1.7816;1.8918];

y=cell（number1）;
f=cell（number1）;
%g=[0.00010.00020.00030.00040.00050.0030.0040.0050.0040.0030.00050.00040.00030.00020.0001];




x=imread（‘lena.tiff‘）;
%x=x（::3）;
[mn]=size（x）;
LR=imread（‘初值X0.bmp‘）;%讀入待重建的低分辨率圖像
%插值
s=imresize（LRscale‘bicubic‘）;%雙三次插值，效果對比
%f0=deconvreg（spsfm*n*0.002）;
f{1}=double（s）;
%f{1}=imresize（f0scale‘bicubic‘）;
%s=imresize（f0scale‘bicubic‘）;
s1=imresize（LRscale‘bilinear‘）;%雙線性插值
x=double（x）;

for?i=1:number
????imageName?=?[?‘LR1_‘?num2str（i）?‘.bmp‘];
???y{i}=double（imread（imageName））;
end%讀入16幅同一場景下的LR圖像
%Sc=generate_L（f{1}number）;
D_f=cell（number1）;
deD_f=cell（number1）;
sadd=zeros（size（x））;
iterative_number=15;%設(shè)定迭代次數(shù)
RMS=zeros（1iterative_number）;
PSNR=zeros（1iterative_number）;
sigma_d=12;sigma_r=0.25;
for?i=1:iterative_number
????i
Sc=generate_L（f{i}numberdeltaphisdspscale）;%在原有的退化模型下退化
delta=-delta;phi=-phi;
??for?k=1:number
??????n=k;
?????%?y{k}=double（y{k}）;
?????Sc{k}=255*Sc{k};
%size（Sc{k}）
%size（y{k}）
????D_f{k}=sign（Sc{k}-y{k}）;????
????deD_f{k}=huifu（D_f{k}nscalenoisesddeltaphi）;????
????sadd=sadd+deD_f{k};
??end
?sadd=sadd/number;
%frec=fswt（f{i}）;
frec=composite_cs（f{i}）;
?beta=80;
f{i+1}=f{i}-8*sadd-frec;%公式見文章中的迭代算法

???for?m1=1:m
?????for?n1=1:n
????????if?f{i+1}（m1n1）<0
?????????????f{i+1}（m1n1）=0;
????????elseif?f{i+1}（m1n1）>255
??????????????????f{i+1}（m1n1）=255;
????????end
??????end
???end%將像素值統(tǒng)一到0-255
???
[MN]=size（f{i+1}）;

%f{i+1}（M-1:）=f{i+1}（M-2:）;?
%f{i+1}（M:）=f{i+1}（M-2:）;?
??%f{i}?=?mybilateral（f{i}sigma_dsigma_r）;
??[PSNR（i）RMS（i）]=indicator（xf{i}）
???figure;imshow（f{i}/255[]）;
???xlabel（[‘第‘num2str（i）‘次迭代重建圖像?PSNR=‘num2str（PSNR（i））]）;
???delta=-delta;phi=-phi;
?
end
???


toc
figure
imshow（x/255[]）;xlabel（‘原始HR圖像‘）
%雙三次插值結(jié)果
[PSNR1RMS1]=indicator（xdouble（s））
figure
imshow（s[]）;xlabel（[‘雙三次插值?PSNR=‘num2str（PSNR1）]）
%雙線性插值結(jié)果
[PSNR2RMS2]=indicator（xdouble（s1））
figure
imshow（s1[]）;xlabel（[‘雙線性插值?PSNR=‘num2str（PSNR2）]）

figure
imshow（LR[]）;xlabel（‘輸入LR圖像‘）
%保存評價指標
dict_path=[‘PSNR+RMSE.mat‘];
save（dict_path?‘PSNR‘?‘RMS‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-12-28?09:49??像素上超分辨重構(gòu)\
?????文件?????????899??2012-04-09?09:20??像素上超分辨重構(gòu)\Affine_trans.asv
?????文件?????5315639??2011-12-26?16:58??像素上超分辨重構(gòu)\Fast?and?robust?multiframe?super?resolution.pdf
?????文件????????3120??2012-07-25?10:43??像素上超分辨重構(gòu)\chongjian.m
?????文件?????????405??2012-04-10?20:35??像素上超分辨重構(gòu)\composite_cs.m
?????文件????????2639??2012-04-09?08:36??像素上超分辨重構(gòu)\dic_training.m
?????文件????????1487??2012-02-21?22:10??像素上超分辨重構(gòu)\generate_L.m
?????文件????????1249??2012-07-25?10:24??像素上超分辨重構(gòu)\get_LR.m
?????文件?????????196??2011-08-31?17:53??像素上超分辨重構(gòu)\goback.m
?????文件?????????483??2012-04-10?15:17??像素上超分辨重構(gòu)\huifu.m
?????文件?????????310??2012-04-10?20:31??像素上超分辨重構(gòu)\indicator.m
?????文件????????2886??2012-02-07?09:28??像素上超分辨重構(gòu)\shift_fun1.m
?????目錄???????????0??2018-03-27?16:24??像素上超分辨重構(gòu)\產(chǎn)生待重建的低分辨率圖像\
?????文件??????298097??2011-07-26?21:38??像素上超分辨重構(gòu)\產(chǎn)生待重建的低分辨率圖像\1.gif
?????文件?????????730??2012-04-09?09:45??像素上超分辨重構(gòu)\產(chǎn)生待重建的低分辨率圖像\generate_L1.m
?????文件?????????427??2017-09-19?20:09??像素上超分辨重構(gòu)\產(chǎn)生待重建的低分辨率圖像\input_LR.m
?????文件??????249438??2011-07-26?21:38??像素上超分辨重構(gòu)\產(chǎn)生待重建的低分辨率圖像\lena.tiff
?????文件?????????345??2012-07-25?10:38??像素上超分辨重構(gòu)\產(chǎn)生待重建的低分辨率圖像\para.mat
?????文件???????17462??2012-07-25?10:38??像素上超分辨重構(gòu)\產(chǎn)生待重建的低分辨率圖像\初值X0.bmp
?????目錄???????????0??2018-03-27?16:24??像素上超分辨重構(gòu)\測試圖像\
?????文件??????298097??2011-07-26?21:38??像素上超分辨重構(gòu)\測試圖像\1.gif
?????文件??????298097??2011-07-26?21:38??像素上超分辨重構(gòu)\測試圖像\10.gif
?????文件??????298097??2011-07-26?21:38??像素上超分辨重構(gòu)\測試圖像\11.gif
?????文件??????298097??2011-07-26?21:38??像素上超分辨重構(gòu)\測試圖像\12.gif
?????文件??????298097??2011-07-26?21:38??像素上超分辨重構(gòu)\測試圖像\13.gif
?????文件??????298097??2011-07-26?21:38??像素上超分辨重構(gòu)\測試圖像\14.gif
?????文件??????298097??2011-07-26?21:38??像素上超分辨重構(gòu)\測試圖像\15.gif
?????文件??????257214??2011-07-26?21:38??像素上超分辨重構(gòu)\測試圖像\16.tiff
?????文件??????224348??2011-07-26?21:38??像素上超分辨重構(gòu)\測試圖像\17.tiff
?????文件??????398064??2011-07-26?21:38??像素上超分辨重構(gòu)\測試圖像\18.tiff
?????文件??????292638??2011-07-26?21:38??像素上超分辨重構(gòu)\測試圖像\19.tiff
............此處省略373個文件信息