% 定义带宽和中心频率
fb = 2; fc = 1;
% 设置有效支撑和网格
lb = -4; ub = 4; n = 1000;
% 计算复Morlet小波cmor2-1
[psi,x] = cmorwavf(lb,ub,n,fb,fc);
% 画出复Morlet小波
figure(1)
plot(x,real(psi),x,imag(psi),'m-.')
grid on;
title('复Morlet小波cmor2-1')
xlabel('时间t/s')
ylabel('幅值A/dB')
legend('实部','虚部');
%%%%%%进行FFT变换并做频谱图%%%%%%%
N=10000;
fs=1000;%设定采样频率
y=fft(real(psi),N);%进行fft变换
mag=abs(y);%求幅值
m=length(y);
f=(0:m/2-1)'*fs/m;%进行对应的频率转换
figure(2);
plot(f,mag(1:m/2));%做频谱图
axis([0,16,0,100]);
xlabel('频率(Hz)');
ylabel('幅值');
grid on;
这个你试一下 |