信号实验报告2.doc

计算机信号与线性系统 实验报告 班 级： 学 号： 姓 名： 指导教师： 完成日期： 计算机科学与技术学院 实验二 信号的Fourier分析 一 实验目的 1 掌握周期信号的频谱-------Fourier级数的分析方法及物理意义； 2 深入理解信号频谱的概念，掌握典型信号的频谱以及Fourier变换的主要性质； 3通过阅读。修改并调试本实验系统所给周期信号频谱分析的源程序，加强Matlab的编程能力； 二 实验原理及方法 1 Fourier级数的理论告诉我们：任何周期信号只要满足Dirichlet条件就可以分解成许多指数分量之和或许多正弦，余弦分量之和直流分量及许多正弦，余弦分量之和； 2 一般来说，任意周期信号表示为Fourier级数时需要无限多项才能完全逼近原函数。但在实际应用中，我们经常采用有限项级数来代替无限级数； 3 通过选择不同的周期函数，及输入不同的周期值和谐波次数，可以得到四种不同得周期信号的频谱分析结果； 4 把上述理论推广到非周期信号中去，就可以导出Fourier变换； 三 实验内容及步骤 1 在上机之前，认真复习信号的Fourier分析有关内容，仔细阅读本实验原理和步骤。准备实验所需数据； 2 周期信号的频谱分析 （1）选择周期信号波形 （2）输入参数 （3）显示波形 （4）编程 （a）1．TTao/TT=1/10 % design2_1.m TT=10; % 周期信号周期 N=20; % 希望看到的谐波次数 syms TTa TT k; % 定义时间 TT、周期 TTa、下标 k 三个符号 y=subs(sym(Heaviside(TT+TT/10)-Heaviside(TT-TT/10)),TT,TT); % 锯齿脉冲基础波形，TTao/TT=1/10 A0=int(y,TT,-TT/2,TT/2)/TT; % 直流分量 A0 % 利用符号法得到各频率分量的复数 Fourier 系数 Ak 的符号表达式 Ak=int(y*exp(-2*i*pi*k*TT/TT),TT,-TT/2,TT/2)/TT; % 利用 symmul 命令得到 Ak*exp(2*j*k*pi*TT/TT) 的符号表达式 fk=sym(Ak)*sym(exp(2*i*k*pi*TT/TT)); % 计算 k 在区间 [-N,N] 内的 Fourier 系数，存入向量 a（长度为 2*N+1） for m=-N:-1 a(m+N+1)=double(sym(subs(Ak,k,m))); end a(N+1)=double(sym((A0))); for m=1:N a(m+N+1)=double(sym((subs(Ak,k,m)))); end % for % 利用 symmul 命令进行符号求和，得到综合信号 f f=symsum(fk,k,-N,-1)+A0+symsum(fk,k,1,N); % 图形显示结果 figure n=-N:N; as=abs(a)*2; % 由 Fourier 系数得到频谱幅度 subplot(3,1,1) ezplot(y,[-TT,TT]) ylabel(原函数); subplot(3,1,2) ezplot(f,[-TT,TT]) ylabel(合成函数); subplot(3,1,3) stem(n,as) ylabel(幅度频谱图); 2．TTao/TT=1/20 % design2_1.m TT=10; % 周期信号周期 N=20; % 希望看到的谐波次数 syms TTa TT k; % 定义时间 TT、周期 TTa、下标 k 三个符号 y=subs(sym(Heaviside(TT+TT/10)-Heaviside(TT-TT/10)),TT,TT); % 锯齿脉冲基础波形，TTao/TT=1/20 A0=int(y,TT,-TT/2,TT/2)/TT; % 直流分量 A0 % 利用符号法得到各频率分量的复数 Fourier 系数 Ak 的符号表达式 Ak=int(y*exp(-2*i*pi*k*TT/TT),TT,-TT/2,TT/2)/TT; % 利用 symmul 命令得到 Ak*exp(2*j*k*pi*TT/TT) 的符号表达式 fk=sym(Ak)*sym(exp(2*i*k*pi*TT/TT)); % 计算 k 在区间 [-N,N] 内的 Fourier 系数，存入向量 a（长度为 2*N+1） for m=