function y=myidwt2(LL,HL,LH,HH);
% 函數(shù) MYIDWT2() 對(duì)輸入的子矩陣序列進(jìn)行逆小波變換，重構(gòu)出矩陣 y
% 輸入?yún)?shù)：LL,HL,LH,HH —— 是四個(gè)大小均為 r*c 維的矩陣
% 輸出參數(shù)：y —— 是一個(gè)大小為 2r*2c 維的矩陣

lpr=[1 1];hpr=[1 -1];           % 默認(rèn)的低通、高通濾波器
tmp_mat=[LL,HL;LH,HH];          % 將輸入的四個(gè)矩陣組合為一個(gè)矩陣
[row,col]=size(tmp_mat);        % 求出組合矩陣的行列數(shù)

for k=1:col                     % 首先對(duì)組合矩陣tmp_mat的每一列，分開成上下兩半
    ca1=tmp_mat(1:row/2,k);     % 分開的兩部分分別作為平均系數(shù)序列ca1、細(xì)節(jié)系數(shù)序列cd1
    cd1=tmp_mat(row/2+1:row,k);
    tmp1=myidwt(ca1,cd1,lpr,hpr);   % 重構(gòu)序列
    yt(:,k)=tmp1;                % 將重構(gòu)序列存入待輸出矩陣 yt 的相應(yīng)列，此時(shí) y=[L|H]
end

for j=1:row                     % 將輸出矩陣 y 的每一行，分開成左右兩半
    ca2=yt(j,1:col/2);           % 分開的兩部分分別作為平均系數(shù)序列ca2、細(xì)節(jié)系數(shù)序列cd2
    cd2=yt(j,col/2+1:col);
    tmp2=myidwt(ca2,cd2,lpr,hpr);   % 重構(gòu)序列
    yt(j,:)=tmp2;                % 將重構(gòu)序列存入待輸出矩陣 yt 的相應(yīng)行，得到最終的輸出矩陣 y=yt
end
y=yt;