信号的采样与重建

通过用MATLAB对f(t)= 5sin(2*pi*30*t)+2sin(2*pi*60*t)+0.5sin(2*pi*90*t)进行设计仿真，让我们通过试验论证理论的正确性，同时学会使用并掌握MATLAB软件的使用，进一步熟悉掌握连续时间信号的傅立叶变换、采样定理等。

1.建立源信号：f(t)= 5sin(2*pi*30*t)+2sin(2*pi*60*t)+0.5sin(2*pi*90*t)，对f（t）进行采样，其结果显示如图1所示：

t=-1:pi/100:1;

x1=5*sin(2*pi*30*t);x2=2*sin(2*pi*60*t);

x3=0.5*sin(2*pi*90*t);f=x1+x2+x3;

subplot(221),plot(t,x1);subplot(222),plot(t,x1);

subplot(223),plot(t,x1);subplot(224),plot(t,f);

2、采样：用120Hz的频率对f(t)进行采样，其采样图如图（2）所示；用240Hz的频率对f(t)进行采样，其采样图如图（3）所示：

fs1=120;t1=-1:1/fs1:1;

 

f1=5*sin(2*pi*30*t1)+2*sin(2*pi*60*t1)+0.5*sin(2*pi*90*t1);

figure(1);plot(t1,f1);

axis([-0.1 0.1 -8 8]);

hold off;

fs2=240;t2=-1:1/fs2:1;

f2=5*sin(2*pi*30*t2)+2*sin(2*pi*60*t2)+0.5*sin(2*pi*90*t2);

figure(2);plot(t2,f2);

axis([-0.1 0.1 -8 8]);

hold off;

3、将二个采样信号进行快速离散傅里叶变换(FFT),观察频谱图,指出是否产生频谱混迭现象. 用120Hz的频率对f(t)进行采样，其采样后快速傅立叶变换频谱图图（4）所示；用240Hz的频率对f(t)进行采样，其采样后快速傅立叶变换频谱图图（5）所示：

f1=30;f2=60;f3=90;fs=120;

N=120;W=2*pi*5;k=0:N-1;w=k*W/N;

t=0:1/fs:0.1;

x1=5*sin(2*pi*f1*t)+2*sin(2*pi*f2*t)+0.5*sin(2*pi*f3*t);

xf1=fft(x1,N);xf1=abs(xf1);w1=120*k/N

figure(1);plot(w1,xf1);

f1=30;f2=60;f3=90;fs=240;

N=240;W=2*pi*5;

k=0:N-1;w=k*W/N;t=0:1/fs:0.1;

x2=5*sin(2*pi*f1*t)+2*sin(2*pi*f2*t)+0.5*sin(2*pi*f3*t);

xf2=fft(x2,N);xf2=abs(xf2);

w2=240*k/N

figure(2);plot(w2,xf2);

4、因为用信号fs=120HZ进行采样时，fs<2f3,其采样频率太小，所以采样信号重建无法复原，

其重建如图（6）所示。当fs=240Hz时，fs>2f3,可以将采样信号重建，对f(t)进行采样的信号

重建如图（7）所示；

Wm=180*pi;Wc=Wm;fs1=120;Ws=2*pi*fs1;n=-500:500;nTs1=n/fs1;

f1=5*sin(2*pi*30*nTs1)+2*sin(2*pi*60*nTs1)+0.5*sin(2*pi*90*nTs1);

Dt=pi/100;t1=-1:Dt:1;

fa1=f1/fs1*Wc/pi*sinc((Wc/pi)*(ones(length(nTs1),1)*t1-nTs1'*ones(1,length(t1))));

t2=-0.1:1/fs1:0.1;

f2=5*sin(2*pi*30*t2)+2*sin(2*pi*60*t2)+0.5*sin(2*pi*90*t2);

figure(1);plot(t1,fa1);

Wm=180*pi;Wc=Wm;fs2=240;Ws=2*pi*fs2;n=-500:500;nTs2=n/fs2;

f1=5*sin(2*pi*30*nTs2)+2*sin(2*pi*60*nTs2)+0.5*sin(2*pi*90*nTs2);

Dt=pi/100;t1=-1:Dt:1;

fa2=f1/fs2*Wc/pi*sinc((Wc/pi)*(ones(length(nTs2),1)*t1-nTs2'*ones(1,length(t1))));

t2=-0.1:1/fs2:0.1;

f2=5*sin(2*pi*30*t2)+2*sin(2*pi*60*t2)+0.5*sin(2*pi*90*t2);

figure(2);plot(t1,fa2);

grid;

5、f(t)用120Hz的频率采样进行快速傅立叶变换后向右平移100Hz如图（8）所示；将f(t)用

240Hz的频率采样进行快速傅立叶变换后向右平移100Hz如图（9）所示：

f1=30;f2=60;f3=90;fs1=120;t=0:1/fs1:0.1;n=500;k=0:n-1;

x1=5*sin(2*pi*f1*t)+2*sin(2*pi*f2*t)+0.5*sin(2*pi*f3*t);

y1=exp(i*100*pi*t);y2=x1.*y1;

y3=fft(y2,n);y3=abs(y3);y4=fft(x1,n);y4=abs(y4);w=300*k/n;

figure(1);plot(w,y3,w,y4);

f1=30;f2=60;f3=90;fs2=240;t=0:1/fs2:0.1;n=500;k=0:n-1;

x1=5*sin(2*pi*f1*t)+2*sin(2*pi*f2*t)+0.5*sin(2*pi*f3*t);

y1=exp(i*100*pi*t);y2=x1.*y1;

y3=fft(y2,n);y3=abs(y3);y4=fft(x1,n);y4=abs(y4);w=300*k/n;

figure(2);plot(w,y3,w,y4);

 

1、信号f(t)= 5sin(2*pi*30*t)+2sin(2*pi*60*t)+0.5sin(2*pi*90*t)。 对f(t) 进行采样，其

采样图如图（1）所示：

 

2、用120Hz的频率对f(t)进行采样，其采样图如图（2）所示；用240Hz的频率对f(t)进行

采样，其采样图如图（3）所示:

图2  f=120HZ时,f(t)的采样信号

图3  f=240HZ时,f(t)的采样信号

3、将二个采样信号进行快速离散傅里叶变换(FFT),观察频谱图,指出是否产生频谱混迭现象。 用120Hz的频率对f(t)进行采样，其采样后快速傅立叶变换频谱图图（4）所示；用240Hz的频率对f(t)进行采样，其采样后快速傅立叶变换频谱图图（5）所示：

图4  f=120HZ时,采样后快速傅立叶变换频谱图

图5  f=240HZ时,采样后快速傅立叶变换频谱图

4、当fs=120HZ时，fs<2f3,其采样频率太小，所以采样信号重建无法复原，其重建如图（6）所

示。当fs=240Hz时，fs>2f3,可以将采样信号重建，对f(t)进行采样的信号重建如图（7）所示：

图6  fs=240Hz时，对f(t)进行采样的信号的重建图谱

图7  fs=240Hz时，对f(t)进行采样的信号的重建图谱

5、f(t)用120Hz的频率采样进行快速傅立叶变换后向右平移100Hz如图（7）所示；将f(t)用

240Hz的频率采样进行快速傅立叶变换后向右平移100Hz如图（8）所示：

图8  fs=120Hz的频率采样快速傅立叶变换后右平移100Hz频谱图

图9  fs=240Hz的频率采样快速傅立叶变换后右平移100Hz频谱图