function WIk=CalculateWIk(I,Wk,width,height,step);
%計算圖像I經過小波變換Wk后的結果圖像；步長step表示，每隔step個點計算一次，即計算寬度為step的離散網格上的小波變換；
WIkAll=conv2(I,Wk,'same');
for x1=1:width/step
    for x2=1:height/step
        WIk(x1,x2)=abs(WIkAll(x1*step,x2*step));
    end
end
%Id = im2double(I);
%for x1=1:width/step
%    for x2=1:height/step
%        WIk(x1,x2)=0;
%        for x3=1:x1*step-1
%            for x4=1:x2*step-1
%                WIk(x1,x2)=WIk(x1,x2)+Wk(x1*step-x3,x2*step-x4)*Id(x3,x4);
%            end
%        end
%    end
%end