北理工信号与系统实验.docVIP

实验3 信号的频域分析 实验目的 深入理解信号频谱的概念，掌握信号的频域分析方法。 观察典型周期信号和非周期信号的频谱，掌握其频谱特性。 实验原理 连续周期信号的频谱分析 如果周期信号满足Dirichlet条件，就可以展开为傅里叶级数的形式，即 ， 式中，表示基波周期， 为基波周期，表示任一个基波周期内的积分。上述两式定义为周期信号复指数形式的傅里叶级数，系数Ck称为傅里叶级数。周期信号的级数还可以由三角函数的线性组合来表示，即 第一式中同频率正弦项和余弦项可以合并，从而得到三角函数形式的傅里叶级数，即 可见，任何满足Dirichlet条件的周期信号都可以一组谐波关系的复指数函数或三角函数的叠加。一般来说表示为傅里叶时需要无限多项才能完全逼近原信号，但是在实际应用中采用有限项，所选越多就越逼近原信号。 对于非周期连续时间信号，信号的傅里叶变换和傅里叶逆变换定义为 , 把信号的时域特性和频域特性联系起来，确定了非周期信和频谱之间的关系。 MATLAB可以方便的求取非周期连续时间信号的傅里叶变换。1）符号运算法 的符号数学工具箱提供了直接求解傅里叶变换反变换的函数，fourier函数和ifourier函数，基本调用格式为， 默认的时域变量为t，频域变量为。 2数值积分法 符号运算的方法外，我们还可以利用MATLAB的quad函数，采用数值积分的方法来进行连续信号的频谱分析。quad是一个用来计算数值积分的函数利用函数可以计算非周期连续时间信号的频谱。q的一般调用格式为：y = quad(fun,a,b)，y = quad(fun,a,b,TOL,TRACE,p1,p2,…) fun是指定被积函数，可以采用inline命令来创建，通过传递函数句柄的形式来制定，、b表示定积分的下限，TOL表示允许的相对或绝对积分误差，TRACE表示以被积函数的点绘图形式来跟踪该函数的返回值，如果TOL和TRACE为空矩阵，则使用缺省值，“1，p2，…”被积函数除时间t之外所需的其他额外输入参数。 3数值近似法 可以近似计算 当t)为时限信号，且足够小，则上式可以演变为 而式中求和部分又可以表示成一个行向量和一列向量的乘积 可以很方便地用ATLAB实现。 基波周期为N的周期序列(n)可以用N个成谐波关系的复指数序列的加权和表示，即 序列表示成上式的形式，称为离散傅里叶级数，而Ck则称为离散傅里叶系数。离散系数Ck可由下式确定。 系数Ck也n)的频谱系数，而且可以证明Ck是以N为的离散频率序列。这说明了离散时间函数对应于频域为周期的离散频率。 用周期与傅里叶系数Ck的来表示周期离散时间信号的频谱即 可以利用MATLAB提供的函数fft用来计算，调用格式为 返回X(k)一个周期内的值，x表示x(n)的样本值。序列(n)可以表示成一组复指数序列的连续和 第二式称为n)的离散时间傅里叶变换，两式确立了非周期离散时间信号n)及其离散时间傅里叶变换关系。频率的函数，为频谱函数，且周期的连续频率函数，其周期为。可见非周期离散时间函数对应于频域中是一个连续的周期频率函数。 可以很方便地利用t)是周期矩阵脉冲信号 解：,, 该信号的傅里叶级数为： 利用MATLAB绘出由前N次谐波合成的信号波形，观察随着N的变化合成信号波形的变化规律 解：设t0= t=-1.5:0.01:1.5; N=input(N=);