小波分解与重构代码

1、load leleccum;s = leleccum(1:3920);% 用 db1 小波函数对信号进行三尺度小波分解C,L=wavedec(s,2,'db1');figure(1);plot(s); title('leleccum 原始信号 ');% 提取尺度 1 的低频系数cA1 = appcoef(C,L,'db1',1); %用小波分解框架C.L计算1层低频系数的近似值，小波基 为 db1% 提取尺度 2 的低频系数cA2 = appcoef(C,L,'db1',2);figure(2);subplot(2,1,1);pl

2、ot(cA1);title(' 尺度 1 的低频系数 ');subplot(2,1,2);plot(cA2);title(' 尺度 2 的低频系数 ');% 提取尺度 1 的高频系数cD1 = detcoef(C,L,1); %用小波分解框架C.L计算1层高频系数的近似值，小波基为db1% 提取尺度 2 的高频系数cD2 = detcoef(C,L,2);figure(3);subplot(2,1,1);plot(cD1);title(' 尺度 1 的高频系数 ');subplot(2,1,2);plot(cD2);title(' 尺度

3、2 的高频系数 ');我给你大概标注了一下，但是你的程序有问题，% 小波图像压缩 - RGB 图像clear all;close all;% 读取图像im = input(' 输入图像 ');% 输入图像名称，要加分号X=imread(im);% 输入要分解的小波层数和小波n=input(' 输入要分解的小波层数 ');%输入所要分解的层数wname = input(' 输入小波名称 ');%输入小波名称，也要加分号x = double(X);NbColors = 255;map = gray(NbColors);x = uint8(x)

4、;%把 RGB 图像转换成灰度图% x = double(X);% xrgb = 0.2990*x(:,:,1) + 0.5870*x(:,:,2) + 0.1140*x(:,:,3);% colors = 255;% x = wcodemat(xrgb,colors);% map = pink(colors);% x = uint8(x);% 对图像 x 进行 n 维小波分解x=imread('D:a.jpg ');map=x;n=3wname= 'sym5 'c,s = wavedec2(x,n,wname);% 使用默认参数选择各层不同的阈值alpha =

5、 1.5; m = 2.7*prod(s(1,:);thr,nkeep = wdcbm2(c,s,alpha,m)% 使用上面的阈值和硬阈值处理进行图像压缩xd,cxd,sxd,perf0,perfl2 = wdencmp('lvd',c,s,wname,n,thr,'h');disp('压缩效率');disp(perf0);% 重构(下面这个地方有问题，你这里是原始图像小波变换后进行重构， xd 才是小波 阀值压缩后重构的图像，cxd, sxd,是c, s经过阀值处理后得到的小波分解结构，也就 是说xd=waverec2(cxd,sxd,wna

6、me);这个 wdencmp函数不需要另外进行重构，你下面那 些关于重构的都没用，而下面压缩后的图像才是重构后的图像，)R = waverec2(c,s,wname);rc = uint8(R);% 显示原始图像和压缩图像subplot(221), image(x);colormap(map);title(' 原始图像 ')subplot(222), image(xd);colormap(map);title(' 压缩后的图像 ')% 显示结果xlab1 = ' 图像压缩后保留能量百分比 ',num2str(perfl2);xlab2 = 

7、9;小波阀值压缩后置零系数百分比 ',num2str(perf0), ' %' xlabel(xlab1 xlab2);subplot(223), image(rc);colormap(map);title(' 重构图像 ');%计算图像大小disp('原始图像');imwrite(x,'original.tif);% 将图像 x 保存为 original.tif，下同 imfinfo('original.tif')% 显示图片 original.tif 详细信息，下同 disp('压缩后的图像');imwrite(xd,'compressed.tif');imfinfo('comp