2022年实验三连续时间LTI系统的频域分析实用 .pdf

实验三连续时间 LTI 系统的频域分析一、实验目的1、掌握系统频率响应特性的概念及其物理意义；2、掌握系统频率响应特性的计算方法和特性曲线的绘制方法，理解具有不同频率响应特性的滤波器对信号的滤波作用；3、学习和掌握幅度特性、相位特性以及群延时的物理意义；4、掌握用MATLAB 语言进行系统频响特性分析的方法。基本要求：掌握LTI 连续和离散时间系统的频域数学模型和频域数学模型的MATLAB 描述方法，深刻理LTI系统的频率响应特性的物理意义，理解滤波和滤波器的概念，掌握利用 MATLAB 计算和绘制LTI 系统频率响应特性曲线中的编程。二、实验原理及方法1 连续时间 LTI 系统的频率响应所谓 频率特性 ，也称为 频率响应特性，简称 频率响应（ Frequency response） ，是指系统在 正弦信号 激励下的稳态响应随频率 变化的情况，包括响应的幅度随频率的变化情况和响应的相位随频率的变化情况两个方面。上图中x(t)、y(t)分别为系统的时域激励信号和响应信号，h(t) 是系统的单位冲激响。它们三者之间的关系为：)(*)()(thtxty，由傅里叶变换的时域卷积定理可得到：)()()(jHjXjY3.1 或者：)()()(jXjYjH 3.2 )(jH为系统的频域数学模型，它实际上就是系统的单位冲激响应h(t)的傅里叶变换。即dtethjHtj)()( 3.3 由于 H(j)实际上是系统单位冲激响应h(t)的傅里叶变换，如果h(t) 是收敛的，或者说是 绝对可积（ Absolutly integrabel）的话，那么H(j) 一定存在，而且H(j)通常是复数 ，因此，也可以表示成复数的不同表达形式。在研究系统的频率响应时，更多的是把连续时间 LTI 系统的时域及频域分析图系统LTI)(th)( jH)(ty)( jX)( jY)(tx名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 9 页 - - - - - - - - - 它表示成极坐标形式：)()()(jejHjH 3.4 上式中，)j (H称为 幅度频率相应（Magnitude response） ，反映信号经过系统之后，信号各频率分量的幅度发生变化的情况，)(称为 相位特性（ Phase response ） ，反映信号经过系统后，信号各频率分量在相位上发生变换的情况。)( jH和)(都是频率的函数。对于一个系统，其频率响应为H(j) ，其幅度响应和相位响应分别为)( jH和)(，如果作用于系统的信号为tjetx0)(，则其响应信号为tjejHty0)()(0tjjeejH00)(0)()(000)(tjejH3.5 若输入信号为正弦信号，即x(t) = sin(0t) ，则系统响应为)(sin(|)(|)sin()()(00000tjHtjHty3.6 可见，系统对某一频率分量的影响表现为两个方面，一是信号的幅度要被)( jH加权，二是信号的相位要被)(移相。由于)( jH和)(都是频率的函数，所以，系统对不同频率的频率分量造成的幅度和相位上的影响是不同的。三、实验内容及步骤Q3-1 修改程序Program3_1，并以 Q3_1存盘，使之能够能够接受键盘方式输入的微分方程系数向量。并利用该程序计算并绘制由微分方程Eq.3.1 、Eq.3.2和 Eq.3.3描述的系统的幅度响应特性、相位响应特性、频率响应的实部和频率响应的虚部曲线图。抄写程序 Q3_1 如下：%Program Q3_1 %This program is used to compute and draw the the plots of the frequency response %Of a continuous-time system clc,clear,close all; b = input(Type in the right coefficient vector of differential equation b); %The coefficient vector of the right side of the differential equation a = input(Type in the left coefficient vector of differential equation a); %The coefficient vector of the left side of the differential equation 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 9 页 - - - - - - - - - H,w= freqs(b,a);%Compute the frequency response H Hm= abs(H);%Compute the magnitude response Hm phai= angle(H);%Compute the phase response phai Hr= real(H);%Compute the real part of the frequency response Hi= imag(H);%Compute the imaginary part of the frequency response subplot(221) plot(w,Hm), grid on, title(Magnitude response), xlabel(Frequency in rad/sec) subplot(223) plot(w,phai), grid on, title(Phase response), xlabel(Frequency in rad/sec) subplot(222) plot(w,Hr), grid on, title(Real part of frequency response), xlabel(Frequency in rad/sec) subplot(224) plot(w,Hi), grid on, title(Imaginary part of frequency response), xlabel(Frequency in rad/sec) %Eq.3.1 b=1 0;a=1 1 25; %Eq.3.2 b=1 -1;a=1 1; %Eq.3.3 b=262;a=1 10 48 148 306 401 262; 执行程序 Q3_1 ，绘制的系统1 的频率响应特性曲线如下：从系统 1 的幅度频率响应曲线看，系统1 是低通、高通、全通、带通还是 带阻滤波器？答：带通滤波器。执行程序 Q3_1 ，绘制的系统2 的频率响应特性曲线如下：名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 9 页 - - - - - - - - - 从系统 2 的幅度频率响应曲线看，系统2 低通、高通、全通、带通还是 带阻滤波器？答：低通滤波器执行程序 Q3_1 ，绘制的系统3 的频率响应特性曲线如下：从系统 3 的幅度频率响应曲线看，系统3 是低通、高通、全通、带通还是 带阻 滤波器？答：带阻滤波器Q3-2编写程序Q3_2 ，使之能够能够接受键盘方式输入的输入信号x(t)的数学表达式，系统微分方程的系数向量，计算输入信号的幅度频谱，系统的幅度频率响应，系统输出信号 y(t)的幅度频谱，系统的单位冲激响应h(t) ，并按照下面的图Q3-2 的布局，绘制出各个信号的时域和频域图形。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 9 页 - - - - - - - - - 图 Q3-2 你编写的程序Q3_2抄写如下：% Program Q3_2 % This Program is used to compute and draw the plots of the frequency response % of a continuous-time system clc,clear,close all T =0.01;dw =0.1; t=0:T:40; w=-4*pi:dw:4*pi; b=input(b=); a = input(a=); x = input(Type in the expression of the input signal x(t); X=x*exp(-j*t*w)*T; X1=abs(X); y=lsim(b,a,x,t); Y=y*exp(-j*t*w)*T; Y1=abs(Y); h=impulse(b,a,40); H,w= freqs(b,a); Hm= abs(H); subplot(324) plot(w,Hm),axis(0 4*pi -0 1.2);grid on,title(Magnitude response of the system) subplot(321) plot(t,x),axis(0 40 -3 3);grid on, title(Input sihnal x(t) 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 9 页 - - - - - - - - - subplot(323) impulse(b,a,40),axis(0 40 -0.2 1);grid on,title(Impulse response h(t) of the system) subplot(325) lsim (b,a,x,t),axis(0 40 -1 1);grid on,title(Output signal y(t) w=-4*pi:dw:4*pi; subplot(322) plot(w,X1),axis(-4*pi 4*pi 0 20);grid on,title(Magnitude response of input signal x(t) subplot(326) plot(w,Y1),axis(-4*pi 4*pi 0 20);grid on,title(Magnitude response of output signal y(t);xlabel(Frequency in rad/sec) %Eq.3.3 b=262;a=1 10 48 148 306 401 262; % x(t)=sin(t)+sin(8*t); % xlabel(Time t); 执行程序 Q3_2，输入信号x(t) = sin(t) + sin(8t) ，输入由Eq.3.3 描述的系统。得到的图形如下：执行程序 Q3_2 所得到的图形请手工绘制出信号x(t) = sin(t) + sin(8t) 的幅度频谱图如下：名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 6 页，共 9 页 - - - - - - - - - 你手工绘制的信号x(t) = sin(t) + sin(8t) 的幅度频谱图与执行程序Q3_2 得到的x(t) = sin(t) + sin(8t) 的幅度频谱图是否相同？如不同，是何原因造成的？答： ：不相同。手工绘制的幅度频谱图是无法做到那么密集的取值，所以看起来比较离散，而用程序绘制的图比较密集。Q3-3 编写程序Q3_3，能够接受从键盘输入的系统微分方程系数向量，并分别绘制所给三个系统的群延时曲线图。抄写程序 Q3_3如下：% Program Q3_3 % This Program is used to compute and draw the plots of the frequency response % of a continuous-time system clc,clear,close all; b = input(b=); a = input(a=); H,w= freqs(b,a); phai= angle(H); phi=-diff(phai); subplot(211) plot(w,phai), grid on, title(Phase response), xlabel(Frequency in rad/sec) w(1,:)=; subplot(212) plot(w,phi), grid on, title(Group delay), xlabel(Frequency in rad/sec) %Eq.3.1 b=1 0;a=1 1 25; %Eq.3.2 b=1 -1;a=1 1; %Eq.3.3 b=262;a=1 10 48 148 306 401 262; 系统 Eq.3.1 的群延时曲线图名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 7 页，共 9 页 - - - - - - - - - 系统 Eq.3.2 的群延时曲线图系统 Eq.3.3 的群延时曲线图根据上面的群延时曲线图可以看出，对系统Eq.3.1，当频率为5 弧度 /秒时，群延时为秒，当频率为10 弧度 /秒时，群延时为 0.2 秒，如何解释这两个群延时时间？根据上面的群延时曲线图，说明这两个系统是否会造成对信号的相位失真？为什么？名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 8 页，共 9 页 - - - - - - - - - 从系统 Eq.3.3 的群延时曲线图中可以看出，当信号的频率为1 弧度 /秒时，系统Eq.3.3 对这一频率的信号的延时是秒。所以，执行程序Q3_2 时，当作用于系统Eq.3.3 的输入信号为 x(t) = sin(t) + sin(8t) 时，其输出信号y(t)的数学表达式为：本实验完成时间：2010 年 6 月 6 日名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 9 页，共 9 页 - - - - - - - - -