生成原始信号
为了方便起见,生成一个包含三个频率的复信号,分别是5Hz、10Hz、15Hz,初相位分别为
采样率取100,采样点取512个,代码如下。
clc;clear;
Fs =100;%采样率
N=512;%序列长度
T = 1/Fs;%采样间隔
t = 0:T:(N-1)*T;%时间序列
s1 = cos(2*pi*5.*t+pi/4)+cos(2*pi*10.*t+pi/8)+cos(2*pi*15.*t+pi/4);%信号实部
s2 = sin(2*pi*5.*t+pi/4)+sin(2*pi*10.*t+pi/8)+sin(2*pi*15.*t+pi/4);%信号虚部
ss = complex(s1,s2);%合成复信号
使用FFT的方法
对信号做fft,并生成相位谱和幅度谱:
y = fft(ss);
f = Fs.*((0:n-1)-(n/2))./(n);
ys = fftshift(y);
plot(f,abs(ys))