【matlab】信号频谱分析之研究采样率对信号混叠的影响

  具体内容：利用FFT，分析连续周期方波信号的频谱：通过改变采样率研究频谱的混叠；

%b混叠
clc;
clear;
close all;
fs1=400;
fs2=4000;
TT1=0.2;
TT2=0.6;
Freq=5;
ts1=1/fs1;
ts2=1/fs2;

t1=0:ts1:(TT1-ts1);
t2=0:ts2:(TT2-ts2);
N=length(t1);
M=length(t2);

x1=(square(2*pi*Freq*t1,50)+1)/2;
x2=(square(2*pi*Freq*t2,50)+1)/2;
y1=fft(x1)/N;
y2=fft(x2)/M;
Z1=abs(y1);
Z2=abs(y2);

k1=-N/2:N/2-1;
k2=-M/2:M/2-1;
f1=k1*fs1/N;
f2=k2*fs2/M;

figure(1)

subplot(2,1,1);
plot(t1,x1);
axis([0,0.6,-0.5,1.5]);
hold on;
plot(t2,x2);
axis([0,0.6,-0.5,1.5]);

subplot(2,1,2);
stem(f1,fftshift(Z1),'b');
xlim([-100,100]);
hold on;
stem(f2,fftshift(Z2),'r');
xlim([-100,100]);
%%%%%%%%%%%%%%%%%%%%%%
k1=-40:-1;
k2=1:40;
k=[k1,0,k2];
f3=k*Freq;

ak1=(sin(k1*pi/2)./k1/pi).*exp(-1i*k1*2*pi*5*0.04875);
ak2=(sin(k2*pi/2)./k2/pi).*exp(-1i*k2*2*pi*5*0.04875);

ak=[ak1,0.5,ak2];
akabs=abs(ak);

figure(2);
stem(f3,akabs);
hold on;
stem(f1,fftshift(Z1),'r');

 

图1：两采样率下的原信号时域图像

图2：两采样率下的原信号频域图像

 理论幅度谱与实际谱的叠加