实验3信号的频域分析（综合型实验）一、实验目的1）深入理解信号频谱的概念，掌握信号的频域分析方法。2）观察典型周期信号和非周期信号的频谱，掌握其频谱特性。二、实验原理与方法1.连续周期信号的频谱分析如果周期信号满足Dirichlet条件，就可展开为傅里叶级数的形式，即0(t)jktkkxce（1）0001(t)ejktkTcxdtT（2）其中0T表示基波周期，002/T为基波频率，0(...)T表示任一个基波周期内的积分。上面两式为周期信号复指数形式的傅里叶级数，系数kc成为x(t)的傅里叶系数。周期信号的傅里叶级数还可由三角函数的线性组合来表示，即00011(t)cossinkkkkxaaktbkt（3）其中000000000122(t),(t)cosk,(t)sinkkkTTTaxdtaxtdtbxtdtTTT（4）（3）式中同频率的正弦、余弦项合并可以得到三角函数形式的傅里叶级数，即001(t)cos(kt)kkkxAA（5）其中2200,,arctankkkkkkbAaAaba（6）任何满足Dirichlet条件的周期信号都可以表示成一组谐波关系的复指数函数或三角函数的叠加。周期信号表示为傅里叶级数时需要无限多项才能完全逼近原信号，但在实际应用中常采用有限项级数代替，所选级数项越多就越接近原信号。2.连续非周期信号的频谱分析对于非周期连续时间信号，信号的傅里叶变换和傅里叶逆变换定义为()(t)ejtXxdt（7）1(t)()e2jtxXd（8）以上两式把信号的时频特性联系起来，确立了非周期信号(t)x和频谱()X之间的关系。利用MATLAB可以方便地求出非周期连续时间信号的傅里叶变换，几种常见方法如下：1）符号运算法MATLAB的符号数学工具箱提供了直接求解傅里叶变换和反变换的函数，fourier函数和ifourier函数，基本调用格式为(x)Xfourierx(X)ifourier默认的时域变量为t，频域变量为。例：求2|t|(t)ex的傅里叶变换，代码及运行结果如下：>>symst>>x=exp(-2*abs(t));>>X=fourier(x)X=4/(w^2+4)所以傅里叶变换结果为24()4X也可利用int函数直接根据式（7）求傅里叶变换。2）数值积分法除了采用符号运算的方法外，还可以采用MATLAB的quad函数，采用数值积分的方法来进行连续信号的频谱分析。quad函数是一个用来计算数值积分的函数。利用quad函数可以计算非周期连续时间信号的频谱。Quad函数的一般调用格式为：y=quad(fun,a,b)y=quad(fun,a,b,TOL,TRACE,p1,p2,...)其中fun指定被积函数，可以用inline命令来创建，也可通过传递函数句柄的形式来指定，a、b表示定积分的上下限，TOL表示允许的相对或绝对积分误差，TRACE表示以被积函数的点绘图形式来跟踪该函数的返回值，如果TOL和TRACE为空矩阵，则使用缺省值，“p1，p2，...”表示被积函数除时间t之外所需的其他额外输入参数。3）数值近似法还可以利用MATLAB的数值计算的方法近似计算连续时间傅里叶变换。傅里叶变换()X可以由（9）式近似计算0()(t)elim(k)ejtjkkXxdtx（9）当x(t)为时限信号，且足够小，则（9）可以演变成()(k)ebjkkaXx（10）（10）中求和部分又可以表示成一个行向量和一个列向量的乘积(1)ee(k)e[x(),x((a1)),...,x(b)]...ejajabjkkajbxa（11）上式可以很方便地利用MATLAB实现。3.离散周期信号的分频域分析基波周期为N的周期序列(n)x可以用N个成谐波关系的复指数序列的加权和表示，即(2/N)n(n)jkkkNxce（12）这里k=表示求和仅需包括一个周期内的N项，周期序列在一个周期内的求和与起点无关。将周期序列表示成式（12）的形式，称为离散傅里叶级数，而系数kc则称为傅里叶系数。离散傅里叶系数kc可由（13）式确定。(2/N)n1(n)ejkkkNcxN（13）傅里叶系数kc也称为(n)x的频谱系数，而且可以证明kc是以N为周期的离散频率序列。这说明了周期的离散时间函数对应于频域为周期的离散频率。这里，我们用周期N与傅里叶系数kc的乘积来表示周期离散时间信号的频谱，即(2/N)n(k)N(n)ejkkkNXcx（14）X(k)可以利用MATLAB提供的函数fft用来计算，调用格式为X=fft(x)该函数返回X(k)一个周期内的值，其中x表示x(n)一个周期内的样本值。4.离散非周期信号的频域分析非周期序列x(n)可以表示成一组复指数序列的连续和21(n)(e)2jjnxXed（15）其中(e)(n)ejjnnXx（16）式（16）成为x(n)的离散时间傅里叶变换，式（15）和（16）确立了非周期离散时间信号x(n)及其离散时间傅...
