matlab实验报告

1、 实 验 报 告 实验课程： MATLAB综合实验 学生姓名： 学 号： 专业班级： 2012年 5 月 20 日 实验四 模拟滤波器频率特性测试一、实验目的· 1、掌握低通无源滤波器的设计；· 2、学会将无源低通滤波器向带通、高通滤波器的转换；· 3、了解常用有源低通滤波器、高通滤器、带通滤波器、带阻滤波器的结构与特性；二、预备知识 · 1、 学习“模拟滤波器的逼近”；· 2、 系统函数的展开方法；· 3、低通滤波器的结构与转换方法； 预习报告中回答以下问题：· 1、实际中常用的滤波器电路类型有哪些，有何特点？·

2、 2、有源滤波器、无源滤波器的概念，优缺点和各自的应用场合？· 3、绘出低通、带通、带阻、高通四种滤波器的理想频响曲线及实际频响曲线，两者有何根本区别，产生原因？三、实验原理模拟滤波器根据其通带的特征可分为：· （1）低通滤波器：允许低频信号通过，将高频信号衰减；· （2）高通滤波器：允许高频信号通过，将低频信号衰减；· （3）带通滤波器：允许一定频带范围内的信号通过，将此频带外的信号衰减；· （4）带阻滤波器：阻止某一频带范围内的信号通过，而允许此频带以外的信号衰减；各种滤波器的频响特性图： 图2一1低通滤波器 图2一2高通滤波器 图2一3带

3、通滤波器 图2一4带阻滤波器· 在这四类滤波器中，又以低通滤波器最为典型，其它几种类型的滤波器均可从它转化而来。1、系统的频率响应特性是指系统在正弦信号激励下系统的稳态响应随激励信号频率变化的情况。用矢量形式表示： 其中：H(j)为幅频特性，表示输出信号与输入信号的幅度比随输入信号频率的变化关系；()为相频特性，表示输出信号与输入信号的相位差随输入信号频率的变化关系。2、H(j)可根据系统函数H(s)求得：H(j)= H(s)s=j因此，对于给定的电路可根椐S域模型先求出系统函数H(s)，再求H(j)，然后讨论系统的频响特性。3、频响特性的测量可分别测量幅频特性和相频特性，幅频特性的

4、测试采用改变激励信号的频率逐点测出响应的幅度，然后用描图法描出幅频特性曲线；相频特性的测量方法亦可改变激励信号的频率用双踪示波器逐点测出输出信号与输入信号的延时，根椐下面的公式推算出相位差 当响应超前激励时为 正，当响应落后激励时为负。四、实验原理图B函数发生=生生生器CH1 示波器RRR/2CC2CINPUTAIN1IN2OUT1OUT2GNDGND 图2一5实验电路图中：R=38k，C=3900pF，红色框内为实验板上的电路。五、实验前预习内容： 1、写出原理图中高、低通及并联后滤波器网络的电压移函数，计算截止频率，并画出幅频特性及相频特性曲线； 2、测试频率特性时，测试点频率应如何选取。

5、六、实验内容及步骤： 将信号源CH1的信号波形调为正弦波，信号的幅度调为Vpp=10V 。 1、RC高通滤波器的频响特性的测量： 将信号源的输出端(A)接实验板的IN1端，滤波后的信号OUT1接示波器的输入(B) 。根据被测电路的参数及系统的频特性，将输入信号的频率从低到高逐次改变十 次以上(幅度保持Vipp=10v) ， 逐个测量输出信号的峰峰值大小(Vopp)及输出信号与输入信号的相位差 ，并将测量数据填入表一：表一Vi（V）1010101010101010101010101010f(Hz)4505486247208229259731023117512761420156020004000V

6、o(v)4.04.65.46.06.46.87.07.27.67.88.08.28.258.3 得到幅频特性为：2RC低通滤波器的频响特性的测量： 将信号源的输出(A)接实验板的IN2，滤波后的输出信号OUT2接示波器的输入(B) 。根据被测电路的参数及系统的幅频特性，将输入信号的频率从低到高逐次改变十 次以上(幅度保持Vipp=10v) ， 逐个测量输出信号的峰峰值大小(Vopp) 及()，并将测量数据填入表二：表二Vi(V)10101010101010101010f(Hz)2050105200286392504602700796Vo(v)9.39.39.39.298.68.27.87.47

7、.0()0-8-15-24-27-30-34-38-44-47幅频特性：相频特性：3双TRC带阻滤波器的频响特性的测量：将实验板上的两输入端IN1与IN2短接，输出端OUT1与OUT2短接；并将信号源的输出 (A)接实验板输入(IN1 )或(IN2 )，滤波后的输出OUT1或OUT2接示波器的输入(B) 。根据被测电路的参数及系统的幅频特性，将输入信号的频率从低到高逐次改变二十 次以上(幅度保持Vipp=10v) ， 逐个测量输出信号的峰峰值大小(Vopp)及() ，并将测量数据填入表三：Vi(V) 1010101010101010101010101010f(Hz)30037847955263

8、2790850938112012541440192322672517Vo(v)5.64.23.22.41.61.00.70.61.42.02.63.64.45.0幅频特性：七、实验仪器：函数发生器一台 ，双踪示波器一台，实验板一块 八、误差分析：(1)由于示波器的稳定性和读数的偏差，导致实验有误差。 (2)实验中 测了大量的数据，由于要求不多于15组，故省去了一些数据。注：在matlab里有些图是用插值命令interp1画的。九、实验小结：通过模拟滤波器频率特性测试实验，我初步了解了有源低通滤波器、高通滤器、带通滤波器、带阻滤波器的结构与特性，同时对仪器的使用有了较熟练的把握。为今后进一步学习

9、奠定了基础。实验五 连续时间系统的模拟一、 目的学习根据给定的连续系统的传输函数，用基本运算单元组成模拟装置。二、 原理1 线性系统的模拟系统的模拟就是用基本运算单元组成的模拟装置来模拟实际的系统。这些实际的系统可以是电的或非电的物理量系统，也可以是社会、经济和军事等非物理量系统。模拟装置可以与实际系统的内容完全不同，但是两者之间的微分方程完全相同，输入输出关系即传输函数也完全相同。模拟装置的激励和响应是电物理量，而实际系统的激励和响应不一定是电物理量，但它们之间的关系是一一对应的。所以，可以通过对模拟装置的研究来分析实际系统，最终达到在一定条件下确定最佳参数的目的。对于那些用数学手段较难处理

10、的高阶系统来说，系统模拟就更为有效。2 传输函数的模拟若已知实际系统的传输函数为： （1）分子、分母同乘以得： （2）式中和分别代表分子、分母的负幂次方多项式。因此： （3）令： （4）则 （5） (6) (7)根据式（6）可以画出如图1所示的模拟框图。在该图的基础上考虑式（7）就可以画出如图2所示系统模拟框图。在连接模拟电路时，用积分器，、及、均用标量乘法器，负号可用倒相器，求和用加法器。值得注意的问题是，积分运算单元有积分时间常数，即积分运算单元的实际传递函数为，所示标量乘法器的标量应分别乘以。同理，应分别乘以。此外，本实验采用的积分器是反相积分器，即传递函数为，所以还应分别乘以，同理，也

11、应分别乘。对于图3(a)所示的电路，其电压传输函数为： （8）如值等于积分器的时间常数，则可以用图3(b)所示的模拟装置来模拟，该装置只用了一个加法器和一个积分时间常数为的反相积分器。附：用信号流图法，有整理成梅森(Mason)公式形式，得： （9）由Mason公式的含义，可画出此系统的信号流图如图4所示，其中和可以用加法器实现，可以用积分器实现，常数及可以用标量乘法器实现。因此，根据此信号的流图可画出图2所示的模拟系统的方框图。图31 模拟框图图32 系统模拟框图图33 一阶RC电路模拟(a) 一阶RC电路；(b) 模拟电路图34 系统信号流图图35 RC低通电路图36 运算单元连接方式，其

12、中该连接方式中的四个运放可采用LM324实现。LM324芯片的管脚如图7所示。图37 LM324芯片的管脚图三、 实验内容图35 所示的 RC低通电路的系统函数： 其中：用基本运算单元模拟图35 所示的 RC低通电路的传输特性。在运算单元连接方式中，反相积分器的时间常数，与图3-5中的RC值一致。实验时分别测量RC电路及其模拟装置的幅频特性，并比较两者是否一致。四、 实验仪器1 GDS-806C数字存储示波器；2 GPD-3303直流电源；3 EE1640C系列函数信号发生器/计数器；4 LM324芯片、相应的电阻、电容和面包板。五、 思考题如反相积分器的积分时间常数与RC电路中的RC值不相等

13、，应如何处理？答：如反相积分器的积分时间常数与RC电路中的RC值不相等，应该通过调整R或C的值使其值与积分时间常数相等。六、 实验数据：F/Hz7899123140172220292328351370410494573U/v4.964.804.644.524.324.043.683.543.323.2032.722.366797979401198130318122.081.801.641.321.160.84RC低通电路：(Vp-p=5V)幅频特性：相频特性：连续时间系统模拟：得到数据如下所示：F/Hz8095128147185227298335360376405487585U/v4.934.

14、824.604.484.313.983.583.483.293.183.12.802.246758059641204135017562.101.751.571.281.080.95 实验六 序列的傅里叶变换1 实验原理：序列的离散时间傅里叶变换（DTFT）和逆变换定义为：DTFTx(n)=X()=IDTFTx()=x(n)=和傅里叶变换、拉氏变换及z变换不同，MATLAB没有直接提供计算DTFT的符号函数方法。但可先用ztrans计算出z变换，然后将z=代入，即可得到DTFT变换。若序列x(n)是有限长，即当且仅当0nN时，x(n)0,就可以用数字方法计算DTFT中-,区间内N个抽样点的值，X

15、(k)=X().事实上，从实序列到这些频域抽样点点的过程被赋予了一个专用名称离散傅里叶变换（DFT）。 DFTx(n)=X(k)= IDFTX(k)=x(n)=为了实现DFT,在MATLAB中可以用下述公式描述 X=Wx X=X注：这里W矩阵是通过n和k从0到N-1变换得到的N*N的矩阵。上诉公式准确的实现了离散傅里叶变换，故有W=NI连续傅里叶变换的求解方法可调用fft函数代替矩阵乘法近似计算连续时间的傅里叶变换。因而有 = 2 实验内容下面是实现矩形脉冲波形和频谱的程序：clcT=2; %定义信号抽样区间，长度N=100; %定义时域和频域抽样点数OMG=100*pi; %定义频域抽样区间

16、t,omg,FT,IFT=prefourier(-T/2,T/2,N,-OMG/2,OMG/2,N);f=0*t; %初始经时域波形f(t>-1/2&t<1/2)=1;F=FT*f; %用矩阵做傅里叶变换fs=IFT*F%用矩阵做傅里叶逆变换f1=f.*exp(-j*omg(1)*t); %生成辅助时域函数f1F1=T*exp(-j*omg(1)*t(1)/N*fft(f1); %用fft函数得到F1F_fft=F1.*exp(j*omg*t(1); %从F1恢复到频谱F_fftf1_ifft=OMG*exp(-j*omg(1)*t(1)/2/pi*ifft(F1); %利

17、用ifft函数做逆变换f_ifft=f1_ifft.*exp(j*omg(1)*t); %由得到的时域辅助信号恢复f_ifftsubplot(2,2,1); %绘制输出波形和频谱plot(t,f,t,fs,'r-.');legend('f','fs');xlabel('t');ylabel('f(t)');subplot(2,2,2)plot(omg,F,'r');xlabel('omega');ylabel('F(w)');subplot(2,2,3);plot(t,f,t,f_ifft,'g-.');legend('f','f_ifft');xlabel('t');ylabel('