信号与系统实验三实验报告.docxVIP

实验三利用DFT分析连续信号频谱 一、实验目的 应用离散傅里叶变换(DFT),分析模拟信号龙⑴的频谱。深刻理解利用DFT分析模拟信号频 谱的原理，分析过程屮出现的现象及解决方法。 二、实验原理 连续周期信号相对于离散周期信号，连续非周期信号相对于离散非周期信号，都可以通过时 域抽样定理建立相互关系。因此，在离散信号的DFT分析方法基础上，增加时域抽样的步 骤，就可以实现连续信号的DFT分析。 三、实验内容 利用FFT分析信号*⑴=eF⑴的频谱。 (1)确定DFT计算的各参数(抽样间隔，截短长度，频谱分辨率等)； 答：选取fm=25Hz为近似的最高频率，则抽样间隔T」/(2佥)=0.02s 选取Tp=10s分析，则截短点数为N =叩丁 = 5oo 釆用矩形窗，确定频域抽样点数为512点。 fsam=50;Tp=10; N=600; T=l/fsam; t=O:T:Tp; x=cxp(-2*t); X=T*fft(x,N); subplot(2,l ,1 );plot(t,x); xlabel(f);title(咐域波形 J; w=(-N/2:N/2-1)*(2 *pi/N)* fsam; y=l./(j*w+2); subplot(2,l ,2);plot(w,abs(flftshift(X)),w,abs(y)/r-.); titled 幅度谱);xlabel(w); legend(理论值丁计算值;0); axis([-10,10,0,l-4]) 当fsam为50HZ时 时域波形 1 幅度谱 1 时域波形 (2)比较理论值与计算值，分析误差原因，提出改善误差的措施。 2.分析周期信号的E) = cos(10M) + 2sin(l畑)频谱时，如果分析长度不为 整周期，利用fft函数计算并绘出其频谱，总结对周期信号进行频谱 分析时，如何选取信号的分析长度。 信号基频心％信号周期 信号最高次谐频为9? = 18兀，所以収W ? 2x9 +1 = 19 如果分析长度不为整周期，比如取N=10,则函数代码为：TO=1;N=1O; T=TO/N; %周期 TO、FFT的点数N、抽样间隔T t=O:T:TO; x=cos(2*pi*5*t)+2*sin(2*pi*9*t); %周期信号 Xm=fft(x,N)/N; %利用FFT计算其频谱 H-(N-l)/2:(N-l)/2)/N/T; %若 N 为偶数 f^l/T/N*(-N/2:(N/2-l)); stem(f,abs(fTtshift(Xm))); %画出幅度谱 xlabclff (Hz),);ylabcl(,magnitudc,); titled幅度谱J; 截 时域波形幅度谱1 4 时域波形 幅度谱 1 4 0.8 0.6 04 02 °5―I~0 °5 f(Hz) 频谱不对称，故应当周期来取，信号基频?=2兀，信号周期％ = 信号最高次谐频为9? = 18兀 所以取N 12x9 + 1 = 19 按 照 N=19 , 取 得 的 函 数 的 截 图 为 时域波形 0.80.60.40.200.5I 0.8 0.6 0.4 0.2 0 0.5 I 2.5 幅度谱 1 —— 1 1—— 1 1 1 1 1 心—— 0.8 - 一 0.6 - 一 $ P 0.4 - 一 0.2 - 一 0 ——6—?g——0———e—— ——~G—— ?10 ?8 -6 -4 -2 0 2 4 6 8 10 f(Hz) 照 N=38 , 取 得 的 函 数 的 截 图 为 幅度谱 1 ~1 《1 1 T f1 1 0.9 - - 0.8 - - 0.7 - - 0.6 - - 0.5 - - 0.4 ■ ■ 0.3 - - 0.2 ■ - 0.1 ■ - ■GO? -15 -10 -5 0 5 10 15 20 f(Hz) 3.假设一实际测得的一段信号的长度为0.4s,其表达式为 x(t)=cos(2pif11)+0.75 cos(2pif2t) 其中f1=100Hz, f2=110HZo当利用FFT近似分析该信号的 频谱时，需要对信号进行时域抽样。试确定一合适抽样频率 人，”，利用DFT分析信号x(t)的频谱。 若在信号截短时使用Hamming窗，由实验确定能够分辨最 小谱峰间隔6和信号长度几的关系。若采用不同参数的 Kaiser窗，重新确定能够分辨最小谱峰间隔6和信号长度乙 的关系。 答：最高次谐频为 pw=220HZ 1R w0=10pi = p=22 = N=46 选取抽样频率为fasm=2pw=440HZ; fasm=440;T0=0.4; N=46; T=l/fasm; % 周期 TO、FFT 的点数 N、抽样间隔 T LOTTO; x=cos(2*pi* 100*t)+