function y=myidwt2(LL,HL,LH,HH);
% 函數 MYIDWT2() 對輸入的子矩陣序列進行逆小波變換，重構出矩陣 y
% 輸入參數：LL,HL,LH,HH —— 是四個大小均為 r*c 維的矩陣
% 輸出參數：y —— 是一個大小為 2r*2c 維的矩陣

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

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

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