matlab中使用pspectrum心得以及用fft函数画频谱图

1、基础使用pspectrum(y,fs)其中y是信号，也就是一维向量，fs是采样频率2、画出chirppspectrum(y,t,'spectrogram','TimeResolution',10e-4, ...'OverlapPercent',99,'Leakage',0.85,'FrequencyLimits',[20e3 50e3])其中 'FrequencyLimits',[20e3 5

youzjuer

11152人浏览 · 2022-04-17 13:59:33

youzjuer · 2022-04-17 13:59:33 发布

1、基础使用

pspectrum(y,fs)

其中y是信号，也就是一维向量，fs是采样频率

2、画出chirp

pspectrum(y,t,'spectrogram','TimeResolution',10e-4, ...
    'OverlapPercent',99,'Leakage',0.85,'FrequencyLimits',[20e3 50e3])

其中 'FrequencyLimits',[20e3 50e3] ,表示频率限制的范围

其中10e-4决定中间黄色的宽度

3、原始频谱

fs = 50e3;
f = 5e3;
t = [0:1/fs:1-1/fs];
y = sin(2*pi*f*t);
L= length(y)
f_sin = fft(y)
f0 = (0:L-1)*fs/L;
plot(f0,abs(f_sin));

画出的的谱线会关于fs/2对称。

剩下的部分可以参考:

matlab绘制频谱图代码_Matlab的FFT绘制频谱图_我是牛哥哥的博客-CSDN博客

4、实信号和复信号的频谱

都是单边谱

fs = 100e4;
fc = 25e3;
t = [0:1/fs:1-1/fs];
y_sin = sin(2*pi*fc*t);
y_exp = exp(1j*2*pi*fc*t);
subplot(2,1,1)
pspectrum(y_exp,fs);
subplot(2,1,2)
pspectrum(y_sin,fs)

当引入频偏，或者叫频谱搬移的时候，则会出现两根谱线：

fs = 100e4;
fc = 25e3;
t = [0:1/fs:1-1/fs];
cfo = 5e3;

y_sin = sin(2*pi*fc*t);
y_exp = exp(1j*2*pi*fc*t);

y_sin_cfo = y_sin.*exp(j*2*pi*cfo*t);
y_exp_cfo = y_sin.*exp(j*2*pi*cfo*t);

subplot(2,1,1)
pspectrum(y_exp_cfo,fs);
subplot(2,1,2)
pspectrum(y_sin_cfo,fs)
% cfo = 10e3;

其中负的谱线是：-fc + cfo

如果使用fft直接计算，则实际上还是两根谱线，也是实际中的样子

fs = 100e4;
fc = 25e3;
t = [0:1/fs:1-1/fs];
cfo = 5e3;

y_sin = sin(2*pi*fc*t);
y_exp = exp(1j*2*pi*fc*t);

L= length(y_sin)
f = (0:L-1)*fs/L-fs/2;
f1 = (0:L-1)*fs/L;

y_sin_cfo = y_sin.*exp(j*2*pi*cfo*t);
y_exp_cfo = y_sin.*exp(j*2*pi*cfo*t);

subplot(2,1,1)
% pspectrum(y_exp_cfo,fs);
plot(f,abs(fftshift(fft(y_sin))))
subplot(2,1,2)
% pspectrum(y_sin_cfo,fs);
plot(f1,abs(fft(y_sin)));