信号实验二离散信号的频谱分析.doc

实验二 离散信号的频谱分析 一、[实验目的] （1）加深对采样定理的理解和掌握，以及对信号恢复的必要性； （2）掌握对连续信号在时域的采样与重构的方法 （3）理解和加深傅里叶变换的概念及其性质。 （4）离散时间傅里叶变换(DTFT)的计算和基本性质。 （5）离散傅里叶变换(DFT)的计算和基本性质。 二、[实验内容] 1.实验原理验证 （一）．采样定理及采样后信号的频谱 对Sa(t)的采样后信号的频谱 （二）．信号重建 对cos(t)的采样与重建信号cos(t) cos(t)重建信号与原信号的比较及误差 ．离散时间信号的傅立叶变换及频谱分析 （1））离散时间傅里叶变换的概念及其性质。 有限长序列x(n)={1,2,3,4,5} 离散傅里叶变换的概念及其性质 x(n)=sin(n*pi/8)+sin(n*pi/4)，N=16的序列傅里叶变换。 2. 选取信号 f(t)= cos(t)作为被采样信号（最高频率为f=8Hz），取理想低通的截止 频率wc=1/2*ws。实现对信号f(t)= cos(t)的采样及由该采样信号的恢复重建，按 要求完成以下内容: (1) 分别令采样角频率ws=1.5*wm 及ws=3*wm，给出在欠采样及过采样条件下冲 激取样后信号的频谱，从而观察频谱的混叠现象。 答：实验程序如下 clc,clear dt=0.01; t=0:dt:1; f=8; %信号频率 wm=2*pi*f; %信号角频率 ft=cos(wm*t); %时域信号 %bs=1.5; %采样角频率，欠采样 bs=3; %采样角频率，大于两倍采样 ws=bs*wm; Ts=2*pi/ws; %采样时间间隔 wc=1/2*ws; %理想低通截止频率 nTs=0:Ts:1; Tf=0.01; nTf=-10:Tf:10; f_nTs=cos(wm*nTs); %时域采样信号 Fs=funexer4_1(f_nTs,nTs,Ts,nTf); figure(1); plot(nTf,Fs); title(cos(t)的3倍采样信号频谱); xlabel(ω); ylabel(F(jw)); grid on %//////////////////1.5倍采样 figure(2) bs=1.5; %采样角频率，大于两倍采样 ws=bs*wm; Ts=2*pi/ws; %采样时间间隔 wc=1/2*ws; %理想低通截止频率 nTs=0:Ts:1; Tf=0.01; nTf=-10:Tf:10; Fs=funexer4_1(f_nTs,nTs,Ts,nTf); plot(nTf,Fs); title(cos(t)的1.5倍采样信号频谱); xlabel(ω); ylabel(F(jw)); grid on (2) 若采样角频率取为ws=3*wm，欲使输出信号与输入信号一致为cos(t)，试根据采样信号恢复信号的误差，确定理想低通滤波器H ( jw)的截止角频率Wc的取值范围应为多大？ 答：截止频率wc应满足：?wmwc≤ws/2。? (3) 以f(t)= cos(t)为被采样信号，以矩形脉冲（门宽t =0.5）作为取样信号，取样 周期ws=4*wm，画出采样后信号的奈奎斯特采样频谱图。 clc,clear all wm=1; %信号宽度 bs=4;% 4倍采样角频率 ws=bs*wm; %采样频率 Ts=2*pi/ws; %采样时间间隔 nTs=0:Ts:10 dw=0.1; w=-3*ws:dw:3*ws; dt=0.01 t=-8:dt:8; f=sinc(t/pi);%Sa(t)函数 for n=1:length(t) b(n)=0; for k=1:length(nTs) b(n)=b(n)+rectpuls(t(n)+1*(k-1),0.5)+rectpuls(t(n)-1*k,0.5)；end end xlabel(t); ylabel(P(t)); grid on subplot(224) plot(t,y); title(Sa(t)函数); xlabel(t); y=f.*b;%采样后的时域信号 yf1=fun4_1(b,t,dt,w);%计算矩形脉冲的频谱 yf2=fun4_1(f,t,dt,w);%计算原函数的频谱 yf3=fun4_1(y,t,dt,w);%计算采样后信号的频谱 figure(1) subplot(221) plot(w,yf1); title(矩形脉