信号系统实验报告.doc

实 验 报 告 实验课程： MATLAB综合实验 学生姓名： 胡 慧 敏 学 号： 6100209056 专业班级： 电 子 093 班 2011年 5 月 24 日目录软件实验实验一、Simulink仿真（1）实验二、Simulink仿真（2）实验三、傅里叶变换实验四、拉普拉斯变换，连续时间系统的s域分析实验五、傅里叶变换利用于通信系统硬件实验实验一、周期信号的频谱测试实验二、模拟滤波器频率特性测试实验三、连续时间系统的模拟实验一、Smulink仿真（1）一、实验目的1、熟悉matlab 编程语言；2、掌握用simulink 进行系统仿真和输出波形的功能；3、进一步掌握matlab 的一些其他功能。二、实验任务1、组建一个函数，并用它求x= 的值；2、试利用simulink 对下列微分方程进行系统仿真，并得到输出波形，比较在 e(t)=u(t)、和e(t)=sin(20t)时的输出情况。 1、 用simulink 对微分方程进行系统仿真，求。3、 实验内容与结果1）function函数的调用要求：组建一个函数，并用它求x= 的值；内容：在编写MATLAB程序的时候，稍微拓宽了函数的适用范围，循环相加次数不限于100次，而变为可调。编写.m文件程序如下： function x=add(N,k) % N为循环次数，k为n的幂 x=0; for n=1:N x=x+nk; end在command window中调用该函数运行的结果如下： x= add(100,2)+ add(100,1)+ add(100,-1) x = 3.4341e+0052）Simulink求以下函数（1） 进行仿真并比较在 e(t)=u(t)、和e(t)=sin(20t)时的输出情况。e(t)=u(t)时仿真模型输e(t)=u(t)时出输出波形e(t)=sin(20t)时的仿真模型e(t)=sin(20t)时的输出波形（2）对微分方程 进行系统仿真，求。解：设则当 e(t)=(t)时得系统的相应的方框图如下图所示：仿真示波器波形如图： 实验二、Simulink仿真（2）一、实验目的.1，通过实验，熟悉Simulink基本用法。2，观察信号的波形特点及产生方法。二、实验仪器1、计算机2、MATLAB 软件三、实验内容已知某系统的二阶微分方程r(t)+r(t)+r(t)=e(t)+e(t)，分别用两种方法计算其冲激响应和阶跃响应，对比理论结果进行验证。使用Simulink创建微分方程：当e(t)=(t)时，得系统框图：simulink画图如下：r(t）=(-6)*exp(-2*t)+10*exp(-3*t); 再用matlab程序： t=-1:0.1:10; y=(-6)*exp(-2*t)+10*exp(-3*t); a=plot(t,y); set(a,linewidth,2); set(a,color,red);grid on axis(-1,10,-0.5,4) 可见计算结果与simulink画的响应图一致当e(t)=u(t)时，得系统框图： simulink画图如下：再用matlab画图得下图： 当e(t)=sin(t)时，得系统框图：经用simulink画图得下图：2、 使用传递函数模块求解二阶微分方程i+7i+10i=e+6e+4e的方程解，其中u(t)是脉冲信号。需要使用Simulink求解g(t)。得系统框图：用线性系统lsim仿真的下列程序：a=1,7,10;b=1,6,4;sys=tf(b,a);t=0:0.01:3; figure;subplot(2,2,1);step(sys); subplot(2,2,2); x_step=zeros(size(t);x_step(t0)=1;x_step(t=0)=1/2; lsim(sys,x_step,t);subplot(2,2,3);h1,t1=impulse(sys,t); plot(t1,h1,k);title(Impulse Response);xlabel(Time(sec);ylabel(Amplitude); subplot(2,2,4); x_delta=zeros(size(t);x_delta(t=0)=100;y1,t=lsim(sys,x_delta,t); y2=y1-x_delta;plot(t,y2,k);title(Impulse Response);xlabel(Time(sec);ylabel(Amplitude);故得下图结果：五、实验小结此实验中实用的仿真软件simulink使用简便但是高能功能强大，是信号系统实验的重要工具。必须重视该软件包的学习。实验三、傅里叶变换一、 实验目的将图形化的方式显示傅里叶变换的性质，进一步理解傅里叶变换二、 实验仪器计算机MATLAB软件三、 实验内容2.请编写函数F=fsana(t,f,N),计算周期信号f的前N个指数形式的傅立叶级数系数，t表示f对应的抽样时间（均为一个周期）；再编写函数f=fssyn（F，t），由傅立叶级数系数F合成抽样时间t对应的函数。设计信号验证这两个是否正确。3利用fsana和fssyn计算习题1中x1（t）的钱10个傅立叶级数系数Fi，0i10,用这些系数合成周期为0.5的锯齿波y（t），0t0&ttau)=1R=1;C=0.3;Alpha=1/R/C;H=aplpha./(aplpha+j*omg);V1=FT*V1;V2=V1*HH=itf*hH=ift*v2运行结果：2、 编写程序为：figure,subplot(2,2,1), hold on, box on;plot(t,h,k-);set(gca,FontSize,16,XLim,-20 20,YLim,-0.1 1);xlabel(t);ylabel(h(t);subplot(2,2,2), hold on, box on;plot(omg,real(H),k-);set(gca,FontSize,16,XLim,-10*pi 10*pi,YLim,0 1);xlabel(omega);ylabel(realH(omega);subplot(2,2,3), hold on, box on;plot(t,real(h1),k-);set(gca,FontSize,16,XLim,-20 20,YLim,-0.1 1);xlabel(t);ylabel(h_1(t);subplot(2,2,4), hold on, box on;plot(omg,imag(H1),k-);set(gca,FontSize,16,XLim,-10*pi 10*pi,YLim,-0.5 0.5);xlabel(omega);ylabel(imag(H_1omega);五、实验总结： 通过这几次实验学习了MATLAB软件，可见计算机技术在人们的学习生活中发挥了很重要的作用，作为当代大学生应该认真女里学习计算机技术。在对信号的操作过程中，我加深了对信号与系统这门课程的认识和了解，也产生了兴趣。实验一、周期信号的频谱测试一、实验目的1、掌握周期信号频谱的测试方法；2、了解典型信号频谱的特点，建立典型信号的波形与频谱之间的关系。二、实验原理及方法1、信号的频谱可分为幅度谱、相位谱和功率谱，分别是将信号的基波和各次谐波的振幅、相位和功率按频率的高低依次排列而成的图形。2、连续周期信号的频谱具有离散性、谐波性、收敛性。例如正弦波、周期矩形脉冲、三角波的幅度谱分别如图1-1、1-2、1-3所示：图1-1(a) 正弦波信号图1-1(b) 相应的幅度谱图1-2(a) 周期矩形脉冲 图1-2(b) 相应的幅度谱 图1-3（a） 三角波 图1-3（b）相应的幅度谱因此，信号的频谱测试方法可用频谱分析仪直接测量亦可用逐点选频测量法进行测量。本实验使用GDS-806C型号的数字存储示波器直接测试幅度谱。用示波器直接测试，就是将其与EE1460C函数信号发生器连好。分别输入相应频率和幅度的正弦波,三角波和矩形波，此时示波器将显示按频率由低到高的各输入信号的谐波分量。GDS-806C数字存储示波器测频谱的方法，就是将MATH键按下，F1键选择FFT(快速傅立叶转换)功能可以将一个时域信号转换成频率构成，显示器出现一条红颜色的频谱扫描线。当示波器输入了不同信号的波形时就显示它们相应的频谱, 参数的测量由调试水平（即频率）与垂直（即增益）游标获取，从而得到输入信号的频谱图。三、实验设备GDS-806C数字存储示波器和EE1640函数信号发生器/计数器.四、实验前预习内容：1、计算重复频率为500HZ的方波，三角波的频谱，并画出频谱图；2、计算重复频率为500HZ，脉冲宽度分别为0.4ms和1ms的对称矩形脉冲的频谱，并画出频谱图。五、实验内容及步骤按图1-4示连接电路；1、 测试正弦波，频率为10HZ的幅度频谱 图1-4 实验原理图 将信号源、示波器、按图1-4连接好；信号源CH1的输出波形调为正弦波，输出频率自选，输出信号幅度自选 ，并记录幅度与频率的参数.测出前五次谐波分量.将其数据填入表一。表一：正弦波前五次谐波的幅度谱(HZ)100200300400500(V)0.3944.79*10（-5）1.908*10（-5）4.37*10（-6）4.37*10（-6）2、测试三角波的幅度频谱在实验步骤1的基础上将信号源CH1的输出波形调为三角波(T) ，频率自选,幅度自选.并记录幅度和周期的参数.测出前五次谐波分量。将测量数据填入表二。表二：三角波前五次谐波的幅度谱(HZ)100300500700900(dB)0.2733.315*10（-3）4.37*10（-4）1.204*10（-4）6.927*10（-5）3、测试周期矩形方波脉冲的幅度频谱（2） 将信号源的输出线接“脉冲”输出端 ，信号频率,幅度和脉宽自选，测出信号的前5次谐波分量，填入表三.表三：周期矩形方波脉冲前五次谐波的幅度谱(HZ)100300500700900(dB)0.575.254*10（-2）2.092*10（-2）1.203*10（-2）5.761*10（-2）六、数据分析与处理2、 由表1数据画出正弦波的频谱图，如图1-5示 图1-5 正弦波频谱图3、 由表2数据画出三角波的频谱图，如图1-6示 图1-6 三角波频谱图4、 由表3数据画出方波的频谱图，如图1-7示 图1-7 方波频谱图七、误差分析 1、对于信号周期的采集不是很精确； 2、应用FFT测试时，信号点数有限； 3、测试谐波分量时，利用频谱图形找点记录数据存在较大的偶然误差。八、实验总结做实验之前，必须先了解线路的解法，如果不清楚，在做实验室才去摸索，这将极大的浪费时间，使事倍功半。做实验时，一定要亲力亲为，务必要将每个步骤、每个细节弄清楚，搞明白，试验后，还要复习、思考，这样，印象才深刻，记得才牢固。实验二 模拟滤波器频率特性测试一、实验目的3、 1、掌握低通无源滤波器的设计；4、 2、学会将无源低通滤波器向带通、高通滤波器的转换；5、 3、了解常用有源低通滤波器、高通滤器、带通滤波器、带阻滤波器的结构与特性；二、预备知识 3、 1、 学习“模拟滤波器的逼近”；4、 2、 系统函数的展开方法；5、 3、低通滤波器的结构与转换方法； 三、实验原理模拟滤波器根据其通带的特征可分为：2、 （1）低通滤波器：允许低频信号通过，将高频信号衰减；3、 （2）高通滤波器：允许高频信号通过，将低频信号衰减；4、 （3）带通滤波器：允许一定频带范围内的信号通过，将此频带外的信号衰减；5、 （4）带阻滤波器：阻止某一频带范围内的信号通过，而允许此频带以外的信号衰减；各种滤波器的频响特性图： 图2一1低通滤波器 图2一2高通滤波器 图2一3带通滤波器 图2一4带阻滤波器 在这四类滤波器中，又以低通滤波器最为典型，其它几种类型的滤波器均可从它转化而来。1、系统的频率响应特性是指系统在正弦信号激励下系统的稳态响应随激励信号频率变化的情况。用矢量形式表示： 其中：H(j)为幅频特性，表示输出信号与输入信号的幅度比随输入信号频率的变化关系；()为相频特性，表示输出信号与输入信号的相位差随输入信号频率的变化关系。2、H(j)可根据系统函数H(s)求得：H(j)= H(s)s=j因此，对于给定的电路可根椐S域模型先求出系统函数H(s)，再求H(j)，然后讨论系统的频响特性。3、频响特性的测量可分别测量幅频特性和相频特性，幅频特性的测试采用改变激励信号的频率逐点测出响应的幅度，然后用描图法描出幅频特性曲线；相频特性的测量方法亦可改变激励信号的频率用双踪示波器逐点测出输出信号与输入信号的延时，根椐下面的公式推算出相位差 当响应超前激励时为 正，当响应落后激励时为负。四、实验原理图B函数发生=生生生器CH1 示波器RRR/2CC2CINPUTAIN1IN2OUT1OUT2GNDGND 图2一5实验电路图中：R=38k，C=3900pF，红色框内为实验板上的电路。五、实验内容及步骤： 将信号源CH1的信号波形调为正弦波，信号的幅度调为Vpp=10V 。 1、RC高通滤波器的频响特性的测量： 将信号源的输出端(A)接实验板的IN1端，滤波后的信号OUT1接示波器的输入(B) 。根据被测电路的参数及系统的频特性，将输入信号的频率从低到高逐次改变十 次以上(幅度保持Vipp=10v) ， 逐个测量输出信号的峰峰值大小(Vopp)及输出信号与输入信号的相位差 ，并将测量数据填入表一：Vi(V)101010101010101010101010f(Hz)101002005001K10K50K100K200K500K1M10MVo(v)1.152.144.56.688.728.72T()32ms2.2ms1.1ms320us100us300us0.031250.45450.90910.0031250.010.0033.2RC低通滤波器的频响特性的测量： 将信号源的输出(A)接实验板的IN2，滤波后的输出信号OUT2接示波器的输入(B) 。根据被测电路的参数及系统的幅频特性，将输入信号的频率从低到高逐次改变十 次以上(幅度保持Vipp=10v) ， 逐个测量输出信号的峰峰值大小(Vopp) 及()，并将测量数据填入表二：Vi(V)101010101010101010101010f(Hz)101002005001K10K50K100K200K500K1M10MVo(v)8.69.28.87.525.52760mv138mv53.6mv12.1mv8.4mv6mv0()0280us200us1901402863.61.920ns000.0035710.0050.005260.007140.0357140.166670.277780.5263150.053双TRC带阻滤波器的频响特性的测量：将实验板上的两输入端IN1与IN2短接，输出端OUT1与OUT2短接；并将信号源的输出 (A)接实验板输入(IN1 )或(IN2 )，滤波后的输出OUT1或OUT2接示波器的输入(B) 。根据被测电路的参数及系统的幅频特性，将输入信号的频率从低到高逐次改变二十 次以上(幅度保持Vipp=10v) ， 逐个测量输出信号的峰峰值大小(Vopp)及() ，并将测量数据填入表三：Vi(V)101010101010101010101010f(Hz)101002005001K10K50K100K200K500K1M10MVo(v)8.47.65.443.64.48.88.68.48.69.09.07.0()0700us40018016091.6600ns3601206000.0014280.00250.0055560.006250.1111110.6250.0016670.0027780.0083330.01667七、实验仪器：函数发生器一台 ，双踪示波器一台，实验板一块八、实验小结 本次试验中，我学到了很多东西，加强了我的动手能力，并且在实验过程中遇到一些困难，培养了我的独立思考能力。实验中我往往受益匪浅，独立思考，最终解决问题，培养理论联系实际的作风，增强创新意识。实验三 连续时间系统的模拟一、 实验目的学习根据给定的连续系统的传输函数，用基本运算单元组成模拟装置。二、 实验原理1 线性系统的模拟系统的模拟就是用基本运算单元组成的模拟装置来模拟实际的系统。这些实际的系统可以是电的或非电的物理量系统，也可以是社会、经济和军事等非物理量系统。模拟装置可以与实际系统的内容完全不同，但是两者之间的微分方程完全相同，输入输出关系即传输函数也完全相同。模拟装置的激励和响应是电物理量，而实际系统的激励和响应不一定是电物理量，但它们之间的关系是一一对应的。所以，可以通过对模拟装置的研究来分析实际系统，最终达到在一定条件下确定最佳参数的目的。对于那些用数学手段较难处理的高阶系统来说，系统模拟就更为有效。2 传输函数的模拟若已知实际系统的传输函数为： （1）分子、分母同乘以得： （2）式中和分别代表分子、分母的负幂次方多项式。因此： （3）令： （4）则： （5） (6) (7)根据式（6）可以画出如图1所示的模拟框图。在该图的基础上考虑式（7）就可以画出如图2所示系统模拟框图。在连接模拟电路时，用积分器，、及、均用标量乘法器，负号可用倒相器，求和用加法器。值得注意的问题是，积分运算单元有积分时间常数，即积分运算单元的实际传递函数为，所示标量乘法器的标量应分别乘以。同理，应分别乘以。此外，本实验采用的积分器是反相积分器，即传递函数为，所以还应分别乘以，同理，也应分别乘。对于图3(a)所示的电路，其电压传输函数为： （8）如值等于积分器的时间常数，则可以用图3