《典型序列频谱分析》由会员分享，可在线阅读，更多相关《典型序列频谱分析（10页珍藏版）》请在金锄头文库上搜索。
1、第1章设计任务与要求111对于三种典型序列单位采样序列、实指数序列、矩形序列11.2自行设计一个周期序列1第2章原理及过程21设计原理2第3章设计内容41.1单位采样序列4时域波形41.2傅里叶变换41.3幅度谱及相位谱51.4频移61.5时移72.1时域图形72.2傅里叶变换82.3幅度谱与相位谱92.4频移102.5时移103.1时域图形错误!未定义书签。3.2傅里叶变换错误!未定义书签。3.3幅度谱与相位谱错误!未定义书签。3.4时移错误!未定义书签。3.5频移错误!未定义书签。4.1幅度特性曲线114.4周期序列的DFS124.5傅里叶变换12第4章心得与体会13参考文献14第1章设计
2、任务与要求11对于三种典型序列单位采样序列、实指数序列、矩形序列要求：（1）画出以上序列的时域波形图；（2）求出以上序列的傅里叶变换；（3）画出以上序列的幅度谱及相位谱，并对相关结果予以理论分析；（4）对以上序列分别进行时移，画出时移后序列的频谱图，验证傅里叶变换的时移性质；（5）对以上序列的频谱分别进行频移，求出频移后频谱所对应的序列，并画出序列的时域波形图，验证傅里叶变换的频移性质。1.2自行设计一个周期序列要求：（1）画出周期序列的时域波形图；（2）求周期序列的DFS,并画出幅度特性曲线；（3）求周期序列的FT，并画出幅频特性曲线；（4）比较DFS和FT的结果，从中可以得出什么结论。第2
3、章原理及过程1设计原理运用Matlab及相关软件，对三种典型序列进行频谱分析，得到并验证相应的结果。MATLAB是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学争论、工程设计以及必需进行有效数值计算的众多科学领域供应了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。1MATLAB和Mathematica、Maple并称为三大数学软件。它在数
4、学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算接其他编程语言的程序等，主要应用于工程计算、把握设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。MATLAB的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相像，故用MATLAB来解算问题要比用C，FORTRAN等语言完成相同的事情简捷得多，并且MATLAB也吸取了像Maple等软件的优点，使MATLAB成为一个强大的数学软件。在新的版本中也加入了对C，FORTRAN,C+，JAVA的支持。可以直接调用，用户也可以将自己编写的有用程序导入到MATLAB函数库中便利自
5、己以后调用，此外许多的MATLAB爱好者都编写了一些经典的程序，用户可以直接进行下载就可以用。采纳傅里叶变换的方法对振动的信号进行分解，并按频率挨次开放，使其成为频率的函数，进而在频率域中对信号进行争论和处理的一种过程，称为频谱分析。将信号在时间域中的波形转变为频率域的频谱，进而可以对信号的信息作定量解释。对信号进行频谱分析，是对其进行傅里叶变换，得到其振幅谱与相位谱。分析软件主要为Matlab。对于信号来说，分模拟信号与数字信号。进行频谱分析时，对于模拟信号来说，首先对其进行抽样，使其离散化，然后采纳离散傅里叶变换（DFT）或者快速傅里叶变换（FFT），然后对其幅度（ABS）和相位（ANGL
6、E）的图像进行分析，而对于数字信号来说，则可直接进行离散傅里叶变换或快速傅里叶变换第3章设计内容1.1单位采样序列1.1.1时域波形程序：n=-5:10;y=zeros(1,5),1,zeros(1,10);stem(n,y)axis(-5,10,0,2);title(单位取样序列)图11.2傅里叶变换程序：clearclc%输入信号n=-5:10;%时域信号的时间范围x=zeros(l,5),l,zeros(l,10);%时域信号xw=0:le-2:2;%想要观看的频率范围%预定义y=w;a=w;j=sqrt(-l);%先定义变量维度，提高运算速度%计算频点fori=1:length(w)f
7、=trapz(n,x.*exp(-j*w(i)*n);y(i)=abs(f);a(i)=angle(f);end%输出subplot(3,1,1),plot(n,x)subplot(3,1,2),plot(w,y)subplot(3,1,3),plot(w,a)图21.3幅度谱及相位谱n=1:50;%定义序列的长度是50x=zeros(l,50);%留意：MATLAB中数组下标从1开头x(1)=1;closeall;subplot(3,1,1);stem(x);title(单位冲击信号序列);k=-25:25;X=x*(exp(-j*pi/12.5).(n*k);magX=abs(X);%绘制
8、x(n)的幅度谱subplot(3,1,2);stem(magX);title(单位冲击信号的幅度谱);angX二angle(X);%绘制x(n)的相位谱subplot(3,1,3);stem(angX);title(单位冲击信号的相位谱)图31.4频移程序：fs=256;N=256;%采样频率和数据点数n=0:N-1;t=n/fs;%时间序列f=20;%定义信号中心频率(hz)x=zeros(1,5),1,zeros(1,10);%信号%程序区y1=fft(x,N);%对信号进行快速Fourier变换y2=fftshift(y1);mag2=abs(y2);f1=n*fs/N;%频率序列f2
9、二n*fs/N-fs/2;%shift后频率序列subplot(4,1,1),plot(f2,mag2,r);%绘出随频率变化的振幅图41.5时移程序N=50;n=0:N-1;M=20;xn=1zeros(1,N);nm=mod(n-M),N);xm=xn(nm+1);subplot(3,1,1),stem(xm);title(|X(k)|);k=-25:25;X=x*(exp(-j*pi/25).(nm)*k);magX=abs(X);%绘制x(m)的幅度谱subplot(3,l,2);stem(magX);title(矩形序列的幅度谱);angX二angle(X);%绘制x(m)的相位谱s
10、ubplot(3,1,3);stem(angX);title(矩形序列相位谱)图52.1时域图形程序：n=-5:10;y=zeros(1,5),ones(1,5),zeros(1,6);plot(n,y)stem(n,y)axis(-5,10,0,2);title(矩形序列)；拒形信号序列图62.2傅里叶变换程序：clearclc%输入信号n=-5:10;%时域信号的时间范围x=zeros(l,5),l,zeros(l,10);%时域信号xw=0:le-2:2;%想要观看的频率范围%预定义y=w;a=w;j=sqrt(-l);%先定义变量维度，提高运算速度%计算频点fori=1:length(
11、w)f=trapz(n,x.*exp(-j*w(i)*n);y(i)=abs(f);a(i)=angle(f);end%输出subplot(3,1,1),plot(n,x)subplot(3,1,2),plot(w,y)subplot(3,1,3),plot(w,a)图72.3幅度谱与相位谱程序：n=1:50;x=sign(sign(10-n)+1);closeall;subplot(3,l,l);stem(x);title(矩形信号序列)；k=-25:25;X=x*(exp(-j*pi/25).(n*k);magX=abs(X);%绘制x(n)的幅度谱subplot(3,l,2);stem(
12、magX);title(矩形序列的幅度谱);angX二angle(X);%绘制x(n)的相位谱subplot(3,1,3);stem(angX);title(矩形序列相位谱)矩形序列II1nr-1|I0.5-10*!*0510*15*025*0*540矩形序列的幅度诸TP?J91111(小竹和詢WeJJJJJ_*_*_*051015202630354045矩形序列的相位谱IlliIllijo-jj11沁iPiyf1111110S1015202630354045图82.4频移程序：fs=256;N=256;%采样频率和数据点数n=0:N-1;t=n/fs;%时间序列f=20;%定义信号中心频率(
13、hz)x=zeros(1,5),ones(1,5),zeros(1,6);%信号%程序区y1=fft(x,N);%对信号进行快速Fourier变换y2=fftshift(y1);mag2=abs(y2);f1=n*fs/N;%频率序列f2二n*fs/N-fs/2;%shift后频率序列2.5时移程序：N=50;n=0:N-1;M=20;xn=sign(sign(10-n)+1);nm=mod(n-M),N);xm=xn(nm+1);subplot(3,1,1),stem(xm);title(|X(k)|);k=-25:25;X=x*(exp(-j*pi/25).(nm)*k);magX=abs
14、(X);%绘制x(m)的幅度谱subplot(3,l,2);stem(magX);title(矩形序列的幅度谱);angX二angle(X);%绘制x(m)的相位谱subplot(3,1,3);stem(angX);title(矩形序列相位谱)图104.1幅度特性曲线程序：n=1:50;%定义序列的长度是50x=sin(pi*n);%留意：MATLAB中数组下标从1开头closeall;subplot(3,1,1);stem(x);title(单位冲击信号序列);k=-25:25;X=x*(exp(-j*pi/12.5).(n*k);magX=abs(X);%绘制x(n)的幅度谱subplot
15、(3,1,2);stem(magX);title(单位冲击信号的幅度谱);angX二angle(X);%绘制x(n)的相位谱subplot(3,1,3);stem(angX);title(单位冲击信号的相位谱)x|0-14单位冲击信号库列2IIIIIIIIIIII1Ipl一一甲一甲-9c诉宀早d4T4亠況:心:冬彗护O-忙J占也“也111111116id12III|i+i061015202630354D4650”关单位冲击信号的幅度谱1IIIIII0102030405060单位冲击信号的相位谱11111。MF旷1北屈必11111_5111110102030405D60图1544周期序列的DFS程序：n=0:9;k=0:1:9;x=sin(5*pi*n);WN=exp(-j*2*pi/9);nk=n*k;WNnk=WN.nk;Xk=x*WNnk;xlabel(Xk)4.5傅里叶变换程序：fs=100;%设定采
