小波变换与逆变换
2018-03-21
0 0 0
暂无评分
其他
如何获取积分?
%实例
% 产生一维信号
clc;clear
wname='db3';
randn('seed',531316785);
s=2+kron(ones(1,8),[1 -1])+((1:16).^2)/32+0.2*randn(1,16);
figure;
plot(s);title('(a) 原始图像');
% 利用db3小波进行二尺度分解
% 小波分解函数dwt
[cA1,cD1]=dwt(s,wname);
figure;plot(cA1);title('(b) 近似系数cA1');
figure;plot(cD1);title('(c) 细节系数cD1');
% 对一尺度上的近似系数再次进行小波分解
[cA2,cD2]=dwt(cA1,wname);
figure;plot(cA2);title('(d) 近似系数cA2');
figure;plot(cD2);title('(e) 细节系数cD2');
[Lo_R,Hi_R]=wfilters(wname,'r');
ss=idwt(cA2,cD2,Lo_R,Hi_R);
sss=idwt(ss,cD1,Lo_R,Hi_R);
% 计算机重构信号与原始信号的误差
err=norm(s-sss)
figure;plot(sss);title('(f) 重构的原始信号');
xlabel(['相对误差=',num2str(err)]);
% 产生一维信号
clc;clear
wname='db3';
randn('seed',531316785);
s=2+kron(ones(1,8),[1 -1])+((1:16).^2)/32+0.2*randn(1,16);
figure;
plot(s);title('(a) 原始图像');
% 利用db3小波进行二尺度分解
% 小波分解函数dwt
[cA1,cD1]=dwt(s,wname);
figure;plot(cA1);title('(b) 近似系数cA1');
figure;plot(cD1);title('(c) 细节系数cD1');
% 对一尺度上的近似系数再次进行小波分解
[cA2,cD2]=dwt(cA1,wname);
figure;plot(cA2);title('(d) 近似系数cA2');
figure;plot(cD2);title('(e) 细节系数cD2');
[Lo_R,Hi_R]=wfilters(wname,'r');
ss=idwt(cA2,cD2,Lo_R,Hi_R);
sss=idwt(ss,cD1,Lo_R,Hi_R);
% 计算机重构信号与原始信号的误差
err=norm(s-sss)
figure;plot(sss);title('(f) 重构的原始信号');
xlabel(['相对误差=',num2str(err)]);
相关源码推荐
使用Java开发Android AOA Android开放式附件
0
0
暂无评分
Golang AOA Android开放式附件HID控制
0
0
暂无评分
pclite
0
0
暂无评分
GMSK线性接收机
0
0
暂无评分
传世AFT中变打金版版本--以前很火爆的,收藏很久
0
0
暂无评分
暂无评论