第3章32 频率周期测量和频谱分析
- 格式:ppt
- 大小:2.19 MB
- 文档页数:36
信号与系统实验报告实验三周期信号的频谱分析学院专业班级学号指导教师实验报告评分:_______实验三 周期信号的频谱分析一、实验目的1、掌握连续时间周期信号的傅里叶级数的物理意义和分析方法;2、观察截短傅里叶级数而产生的“Gibbs 现象”,了解其特点以及产生的原因;3、掌握各种典型的连续时间非周期信号的频谱特征。
二、实验容实验前,必须首先阅读本实验原理,读懂所给出的全部例程序。
实验开始时,先在计算机上运行这些例程序,观察所得到的信号的波形图。
并结合例程序应该完成的工作,进一步分析程序中各个语句的作用,从而真正理解这些程序。
实验前,一定要针对下面的实验项目做好相应的实验准备工作,包括事先编写好相应的实验程序等事项。
Q3-1 编写程序Q3_1,绘制下面的信号的波形图:-+-=)5cos(51)3cos(31)cos()(000t t t t x ωωω∑∞==10)cos()2sin(1n t n n nωπ其中,ω0 = 0.5π,要求将一个图形窗口分割成四个子图,分别绘制cos(ω0t)、cos(3ω0t)、cos(5ω0t) 和x(t) 的波形图,给图形加title ,网格线和x 坐标标签,并且程序能够接受从键盘输入的和式中的项数。
抄写程序Q3_1如下: clear,%Clear all variablesclose all,%Close all figure windowsdt = 0.00001; %Specify the step of time variable t = -2:dt:4; %Specify the interval of time w0=0.5*pi; x1=cos(w0.*t); x2=cos(3*w0.*t); x3=cos(5*w0.*t);N=input('Type in the number of the harmonic components N='); x=0; for q=1:N;x=x+(sin(q*(pi/2)).*cos(q*w0*t))/q; endsubplot(221)plot(t,x1)%Plot x1axis([-2 4 -2 2]);grid on,title('signal cos(w0.*t)')subplot(222)plot(t,x2)%Plot x2axis([-2 4 -2 2]);grid on,title('signal cos(3*w0.*t))')subplot(223)plot(t,x3)%Plot x3axis([-2 4 -2 2])grid on,title('signal cos(5*w0.*t))')执行程序Q3_1所得到的图形如下:Q3-2给程序Program3_1增加适当的语句,并以Q3_2存盘,使之能够计算例题1中的周期方波信号的傅里叶级数的系数,并绘制出信号的幅度谱和相位谱的谱线图。
实验报告一、实验目的和要求谱分析即求信号的频谱。
本实验采用DFT/FFT技术对周期性信号进行谱分析。
通过实验,了解用X(k)近似地表示频谱X(ejω)带来的栅栏效应、混叠现象和频谱泄漏,了解如何正确地选择参数(抽样间隔T、抽样点数N)。
二、实验内容和步骤2-1 选用最简单的周期信号:单频正弦信号、频率f=50赫兹,进行谱分析。
2-2 谱分析参数可以从下表中任选一组(也可自定)。
对各组参数时的序列,计算:一个正弦周期是否对应整数个抽样间隔?观察区间是否对应整数个正弦周期?2-3 对以上几个正弦序列,依次进行以下过程。
2-3-1 观察并记录一个正弦序列的图形(时域)、频谱(幅度谱、频谱实部、频谱虚部)形状、幅度谱的第一个峰的坐标(U,V)。
2-3-2 分析抽样间隔T、截断长度N(抽样个数)对谱分析结果的影响;2-3-3 思考X(k)与X(e jω)的关系;2-3-4 讨论用X(k)近似表示X(ejω)时的栅栏效应、混叠现象、频谱泄漏。
三、主要仪器设备MATLAB编程。
四、操作方法和实验步骤(参见“二、实验内容和步骤”)五、实验数据记录和处理clc;clf;clear;%清除缓存%第一组数据的MATLAB程序(之后几组只需要将参数改变即可) T=0.000625;length=32;n=0:length-1;t=0:0.0001:31;%原序列和采样序列xn=sin(2*pi*50*n*T);xt=sin(2*pi*50*t);%画第一幅图(原序列和采样序列)figure(1);subplot(2,1,1);plot(t,xt);xlabel('t');ylabel('xt');axis([0,0.2,-1.1,1.1]);title('原序列时域');subplot(2,1,2);stem(n,xn ,'filled');xlabel('n');ylabel('xn');axis([0,length,-1.1,1.1]);title('采样后序列时域');%画第二幅图(采样序列实部、虚部、模和相角)figure(2);subplot(2,2,1);stem(n,real(xn) ,'filled');xlabel('n');ylabel('real(xn)');axis([0,length,-1.1,1.1]);title('采样序列的实部');subplot(2,2,2);stem(n,imag(xn) ,'filled');xlabel('n');ylabel('imag(xn)');axis([0,length,-1.1,1.1]);title('采样序列的虚部');subplot(2,2,3);stem(n,abs(xn) ,'filled');xlabel('n');ylabel('abs(xn)');axis([0,length,-1.1,1.1]);title('采样序列的模');subplot(2,2,4);stem(n,angle(xn) ,'filled');xlabel('n');ylabel('angle(xn)');axis([0,length,-(pi+0.5),pi+0.5]);title('采样序列的相角');%计算DFTDFT=fft(xn,length);%画第三幅图(DFT的幅度、实部和虚部)figure(3);subplot(3,1,1);stem(n,abs(DFT) ,'filled');xlabel('k');%DFT后的频域变量为kylabel('abs(DFT)');title('DFT 幅度谱');subplot(3,1,2);stem(n,real(DFT) ,'filled');xlabel('k');ylabel('real(DFT)');title('DFT的实部');subplot(3,1,3);stem(n,imag(DFT) ,'filled');xlabel('k');ylabel('imag(DFT)');title('DFT的虚部');六、实验结果与分析实验结果:第一组数据:实验名称:DFT/FFT的应用之一 确定性信号谱分析姓名:张清学号:3110103952 P.4第二组数据:第三组数据:第四组数据:第五组数据:第六组数据:6-1 实验前预习有关概念,并根据上列参数来推测相应频谱的形状、谱峰所在频率(U)和谱峰的数值(V)、混叠现象和频谱泄漏的有无。