function eigvector_pixel=eigvector_wavelet(pixel )
[c,s]=wavedec2(pixel,3,'db3');%多尺度二维小波变换
%CA3=appcoef2(c,s,'db3',3);%提取低频系数
CA1=appcoef2(c,s,'db3',1);%提取低频系数
CH3=detcoef2('h',c,s,3);%提取3层高频系数,水平方向
CV3=detcoef2('v',c,s,3);%提取3层高频系数,垂直方向
CD3=detcoef2('d',c,s,3);%提取3层高频系数,对角方向
CH2=detcoef2('h',c,s,2);%提取2层高频系数,水平方向
CV2=detcoef2('v',c,s,2);%提取2层高频系数,垂直方向
CD2=detcoef2('d',c,s,2);%提取2层高频系数,对角方向
CH1=detcoef2('h',c,s,1);%提取1层高频系数,水平方向
CV1=detcoef2('v',c,s,1);%提取1层高频系数,垂直方向
CD1=detcoef2('d',c,s,1);%提取1层高频系数,对角方向
%提取各子带的能量
eigvector_pixel(1)= sum(sum(CA1.^2));
eigvector_pixel(2)= sum(sum(CH3.^2));
eigvector_pixel(3)= sum(sum(CV3.^2));
eigvector_pixel(4)= sum(sum(CD3.^2));
eigvector_pixel(5)= sum(sum(CH2.^2));
eigvector_pixel(6)= sum(sum(CV2.^2));
eigvector_pixel(7)= sum(sum(CD2.^2));
eigvector_pixel(8)= sum(sum(CH1.^2));
eigvector_pixel(9)= sum(sum(CV1.^2));
eigvector_pixel(10)= sum(sum(CD1.^2));
wavelet小波变换MATLAB程序
需积分: 33 156 浏览量
2015-12-09
21:25:01
上传
评论 2
收藏 1KB ZIP 举报
xiaojiaoyagong
- 粉丝: 0
- 资源: 3