该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
x=imread('F:/flower.bmp');
wname='db5'
[Lo_D,Hi_D] = wfilters(wname);
lf=length(Lo_D);
% 画出原始图像
imshow(x);title('Original Image');
% 标明图像大小
[r,c]=size(x);
xlabel(['Size : ',num2str(r),'*',num2str(c)]);
% 将矩阵x的数据格式转换为适合数值处理的double格式
xd=double(x);
[rx,cx]=size(x);
[o2sa,f1sa,rsx]=sizcoef([rx,cx],lf,N); a=[o2sa,f1sa,rsx]
coef=[];
scf=[rx,cx];
for i=1:N
[cA,cV,cH,cD]=mydwt2(xd,Lo_D,Hi_D); % 第 i 级小波分解
xd=cA; % 将第 i 级分解得到的低频系数矩阵作为第 i+1 级分解的源矩阵
outmp={cV;cH;cD}; % 将第 i 级分解得到的高频系数矩阵cV,cH,cD存入细胞矩阵 outmp
scf=[size(cV);scf]; % 将各级分解系数矩阵的大小存入矩阵 scf
coef=[outmp;coef]; % 将细胞矩阵 outmp 存入输出矩阵 coef,coef将由空矩阵变为细胞矩阵
end
% 迭代结束后&