通信试验指导书(09级).docVIP

实验四 离散系统分析 一、实验目的 １．熟悉离散时间系统的频域分析方法。 ２．掌握离散时间系统频域分析的MATLAB实现方法。 二、实验内容 １．三阶归一化的Butterworth低通滤波器的频率响应为 试画出系统的幅度响应和相位响应。 ２．已知RC电路如图所示，系统的输入电压信号为f(t),输出信号为电阻两端的电压y(t)。当RC=0.04，f(t)=cos5t+cos100t,-t。试求该系统的响应y(t)。 三、实验预备知识 利用MATLAB分析系统的频率特性，当系统的频率响应H（jw）是jw的有理式时，有： MATLAB信号处理工具箱提供的freqs函数可直接计算函数系统的频率响应。 格式：H=freqs(b,a,w) 说明： （１）b是上式中分子多项式的系数； a是上式中分母多项式的系数。 （２）w为需计算的H（jw）的抽样点。（数组w中最少需包含两个w的抽样点）。 四、实验步骤 １．编写程序。 ２．调试程序。 ３．写出程序运行结果。 五、思考题 １．总结连续时间系统的频响特性分析方法。 ２．总结离散时间系统的频响特性分析方法。 实验五 信号﹑系统及系统响应 一、实验目的 熟悉连续周期、非周期信号的频域分析方法及MATLAB编程实现方法。 掌握离散周期、非周期信号的频域分析方法及MATLAB编程实现方法。 二、实验内容 １．试用MATLAB计算如图所示周期矩形波序列的DFS系数。 ２．求如图所示周期矩形脉冲信号的Fourier级数表达式，并用MATLAB画出由前几项Fourier级数Fourier级数系数得出的信号近似波形（N值由键盘输入）。 ３．试画出时，幅度频谱。 4．试画出如图所示矩形序列的频谱。 三、实验预备知识 １．离散周期信号傅立叶级数DFS分析 若设定DFS和1DFS的求和范围为0到N-1， (1) 　　　(2) 则MATLAB提供的函数F=fft(f),可用来计算（1）式定义的N个DFS系数。 N由上式中序列f长度确定。F给出的是DFS系数。MATLAB提供的函数f=ifft(F)由DFS系数F(m)按（2）式计算出时域信号f[k]。 quad8是计算数值积分的函数 1：y=quad8(F,a,b) 格式2：y=quad8(F,a,b,[　],[　],p) 说明： （１）F是一个字符串，它表示被积分函数的文件名，通常使用函数调用。 a,b分别表示定积分的下限和上限。 quad8的返回值是用自适应Simpson算法得出来的积分值P表示被积分函数中的一个参数function　y=sf1(t,w); y=(t=-1t=1).*(1-abs(t)).*exp(-j*w*t); Program for k=1:N F(k)=quad8(‘sf1’,-1,1,[　],[　],w(k)); End ３．离散非周期信号的傅里叶变换 当序列的DTFT可写成的有理多项式时，MATLAB Signal ProcessinToolbox中的fregz函数可用来计算DTFT的值。另外MATLAB提供的abs,angle,real,imad等基本函数可用来计算DTFT的幅度，相位，实部，虚部。设DTFT的有理多项式为： 则freqz的调用形式为：h=frez(b,a,w)　　　　　　　　　　　(2) 说明： （１）b和a分别为（1）式中分子多项式和分母多项式函数的向量 （２）W为抽样的频率点 在以（2）式形式调用freqz函数时，W中至少要有2个频率点返回的值h就是DTFT在抽样点W上的值H的值一般是复数 一般来说是实变量的复值函数，可用实部和虚部将其表示为：=+，其中和分别是的实部和虚部；也可用幅度和相位将表示为：=||，其中||和分别为序列幅度谱和相位谱。 四、实验步骤 １．编写程序。 ２．调试程序。 ３．写出程序运行结果。 五、思考题 １．总结周期连续、离散信号的傅里叶级数分析方法。 ２．总结非周期连续、离散信号的傅里叶变换分析方法。 实验六 FFT对信号进行频谱分析 一．实验目的 1．在理论学习的基础上，通过本次实验，加深对FFT的理解，熟悉FFT算法及其程序的编写。 2．熟悉应用FFT对典型信号进行频谱分析的方法。 3．了解应用FFT进行信号频谱分析过程中可能出现的问题，以便在实际中正确应用FFT。 二．实验内容 (一) 编制实验用主程序及相应子程序 1．编制信号产生子程序及本实验的频谱分析主程序。实验中需要用到的基本信号包括： （1）高斯序列： （2）衰减正旋序列： ? （3）三角波序列 （4）反三角序列 1．观察高斯序列的时域和频域特性 2．观察衰减正旋序列的时域和幅频特性 3．观测三角波序列和反三角波序列的时域和幅频特性 三、试验预备知识 一个连续信号的频谱可以用它的傅立叶变换表示为