http://blog.csdn.net/compression/article/details/9039619
看了幾篇文章,關(guān)于信號(hào)的稀疏逼近,gabor原子庫一般用于語音信號(hào)的表示和去噪,由于其良好的時(shí)頻逼近性能,以及原子庫的過完備性,可以用于表示一維信號(hào)。
其參數(shù)有4個(gè),(s,u,v,w)分別表示尺度拉伸、位移中心、頻率變化、相位變化。將參數(shù)離散化,對(duì)于32位長(zhǎng)的信號(hào),用200個(gè)點(diǎn)來逼近得到10804個(gè)原子。matlab代碼如下
a=2;
delta_u=1/2;
delta_v=pi;
delta_w=pi/6;
N=32;
M=log2(N);
k=1;
T=2e-6;
Fs=100e6;
Ts=1/Fs;
t=0:Ts:T-Ts;
L=length(t);
for j=1:M
s(j)=a^j;
P=N*(2^(-j+1));
Q=2^(j+1);
for p=1:(P+1)
u(p)=(p-1)*(a^j)*delta_u;
for q=1:(Q+1)
v(q)=(q-1)*(a^j)*delta_v;
for i=1:13
w(i)=(i-1)*delta_w;
grt(k,1:L)=1/sqrt(s(j))*exp(-pi*((t-u(p))/s(j)).^2).*cos(v(q)*t+w(i));
k=k+1;
end
end
end
end
% 一共得到10803個(gè)原子,選取了10個(gè)原子顯示
figure;
for tu=1:10
subplot(5,5,tu)
plot(t,grt((tu-1)*500+1,:));
end
圖片如下
聯(lián)系客服