tiaopin .m

调频信号的时-频分析 针对时变信号有很好的效果

x=fmlin(256,0,0.5);
% plot(real(x));
N=length(x);
K=3;
ny=0:K*N-1;
y=x(mod(ny,N)+1);
%  plot(real(y));
% axis([0 1050 -3 3]);
% ft=fftshift(abs(fft(y)).^2);
% plot((-512:511)/1024,ft);
% xlabel('Normalized frequency');
% ylabel('Squared modulus');
h=window(@kaiser,21);
t=1:768;
% [tfr,f]=tfrpwv(y,t,1024,h);
tfrpwv(y,t,768,h);
imagesc(abs(tfr(1:768,':')))
axis('xy');
% tfrwv(y,t,1024,h);
% htl(tfr,512,256,1);
% nn=0:N-1;
% I=cos(nn*2*pi*0.25);
% Q=sin(nn*2*pi*0.25);
% II=x.*I;
% QQ=x.*Q;
% window=chebwin(51,40);
% [b,a]=fir1(50,2*0.05,window);
% II=filter(b,a,II);
% QQ=filter(b,a,QQ);
% plot(Q);