连续时间LTI系统分析
连续时间LTI系统分析
连续时间L T I系统分析(总8页) -本页仅作为预览文档封面，使用时请删除本页-实验三 连续时间LTI 系统分析一、实验目的（一）掌握使用Matlab 进行连续系统时域分析的方法1、学会使用符号法求解连续系统的零输入响应和零状态响应2、学会使用数值法求解连续系统的零状态响应3、学会求解连续系统的冲激响应和阶跃响应（二）掌握使用Matlab 进行连续时间LTI 系统的频率特性及频域分析方法1、学会运用MATLAB 分析连续系统的频率特性2、学会运用MATLAB 进行连续系统的频域分析（三）掌握使用Matlab 进行连续时间LTI 系统s 域分析的方法1、学会运用MATLAB 求拉普拉斯变换（LT ）2、学会运用MATLAB 求拉普拉斯反变换（ILT ）3、学会在MATLAB 环境下进行连续时间LTI 系统s 域分析二、实验条件装有MATLAB 的电脑三、实验内容（一）熟悉三部分相关内容原理（二）完成作业1、已知某系统的微分方程如下：)(3)()(2)(3)(t e t e t r t r t r +'=+'+''其中，)(t e 为激励，)(t r 为响应。
（1） 用MATLAB 命令求出并画出2)0(,1)0(),()(3='==---r r t u e t e t 时系统的零状态响应和零输入响应（零状态响应分别使用符号法和数值法求解，零输入响应只使用符号法求解）；符号法求解零输入响应： >> eq='D2y+3*Dy+2*y=0';>> cond='y(0)=1,Dy(0)=2';>> yzi=dsolve(eq,cond);>> yzi=simplify(yzi)yzi =符号法求解零状态响应：exp(-2*t)*(4*exp(t) - 3)eq1='D2y+3*Dy+2*y=Dx+3*x';eq2='x=exp(-3*t)*heaviside(t)';cond='y=0,Dy=0';yzs=dsolve(eq1,eq2,cond);yzs=simplify(yzs)yzs =(exp(-2*t)*(exp(t) - 1)*(sign(t) + 1))/2图像如下：代码：subplot(211)ezplot(yzi,[0,8]);grid ontitle('ÁãÊäÈëÏìÓ¦')subplot(212)ezplot(yzs,[0,8]);grid ontitle('Áã״̬ÏìÓ¦')数值计算法：t=0::10;sys=tf([1,3],[1,3,2]);f=exp(-3*t).*uCT(t);y=lsim(sys,f,t);plot(t,y),grid on ;axis([0 10 ]);title('ÊýÖµ¼ÆËã·¨µÄÁã״̬ÏìÓ¦')（2）使用MATLAB命令求出并画出系统的冲激响应和阶跃响应（数值法）；用卷积积分法求系统的零状态响应并与（1）中结果进行比较；系统的冲激响应和阶跃响应（数值法）：代码：t=0::10;sys=tf([1,3],[1,3,2]);h=impulse(sys,t);g=step(sys,t);subplot(211)plot(t,h),grid on;axis([0 10 ]);title('³å¼¤ÏìÓ¦')subplot(212)plot(t,g),grid on;axis([0 10 ]);title('½×Ô¾ÏìÓ¦'卷积积分法求系统的零状态响应：Ctsconv函数的定义：function[f,t]=ctsconv(f1,f2,t1,t2,dt)f=conv(f1,f2);f=f*dt;ts=min(t1)+min(t2);te=max(t1)+max(t2);t=ts:dt:te;subplot(221)plot(t1,f1);grid onaxis([min(t1),max(t1),min(f1)-abs(min(f1)*,max(f1)+abs(max(f1)*])title('f1(t)');xlabel('t')subplot(222)plot(t2,f2);grid onaxis([min(t2),max(t2),min(f2)-abs(min(f2)*,max(f2)+abs(max(f2)*])title('f2(t)');xlabel('t')subplot(212)plot(t,f);grid onaxis([min(t),max(t),min(f)-abs(min(f)*,max(f)+abs(max(f)*])title('f(t)=f1(t)*f2(t)');xlabel('t')求系统的零状态响应代码：dt=;t1=0:dt:10;f1=exp(-3*t1).*uCT(t1);t2=t1;sys=tf([1,3],[1,3,2]);f2=impulse(sys,t2);[t,f]=ctsconv(f1,f2,t1,t2,dt)如图，根据两图相比较，两种方法做出的零状态响应大体相同。
MATLAB与信号实验——连续LTI系统的时域分析
MATLAB与信号实验——连续LTI系统的时域分析连续LTI系统的时域分析是信号与系统学中的重要课题。
MATLAB作为一种强大的科学计算软件，提供了丰富的工具和函数来进行信号与系统的分析。
下面将介绍MATLAB在连续LTI系统时域分析中的应用。
首先，我们需要了解连续LTI系统的基本概念。
一个连续域线性时不变系统（LTI系统）可以由它的冲激响应完全描述。
冲激响应是系统对单位冲激信号的响应。
在MATLAB中，可以使用impulse函数来生成单位冲激信号。
假设我们有一个连续LTI系统的冲激响应h(t)，我们可以使用conv 函数来计算系统对任意输入信号x(t)的响应y(t)。
conv函数实现了卷积运算，可以将输入信号与冲激响应进行卷积运算得到输出信号。
例如，我们假设一个连续LTI系统的冲激响应为h(t) = exp(-t)u(t)，其中u(t)是单位阶跃函数。
我们可以使用以下代码生成输入信号x(t)和计算输出信号y(t)：```matlabt=-10:0.1:10;%时间范围x = sin(t); % 输入信号h = exp(-t).*heaviside(t); % 冲激响应y = conv(x, h, 'same'); % 计算输出信号```这段代码首先定义了时间范围t，然后定义了输入信号x(t)和冲激响应h(t)。
接下来，使用conv函数计算输入信号和冲激响应的卷积，设置参数’same’表示输出信号与输入信号长度相同。
最后，得到了输出信号y(t)。
在得到输出信号后，我们可以使用MATLAB的绘图功能来可视化结果。
例如，使用以下代码可以绘制输入信号和输出信号的图像：```matlabfigure;plot(t, x, 'b', 'LineWidth', 2); % 绘制输入信号hold on;plot(t, y, 'r', 'LineWidth', 2); % 绘制输出信号xlabel('时间');ylabel('幅度');legend('输入信号', '输出信号');```除了卷积运算外，MATLAB还提供了许多其他函数来进行连续LTI系统的时域分析。
第五章1-连续LTI系统频域分析
连续时间LTI系统的频域分析 离散时间LTI系统的频域分析 信号的幅度调制和解调
时域分析的要点是，以冲激函数为基本信号，
任意输入信号可分解为一系列冲激函数；而系统零 状态响应yzs(t) = x(t)*h(t)。 由单位冲激函数δ (t)所引起的零状态响应称为单位 冲激响应，简称冲激响应，记为h(t)。
解： 利用H(j)与h(t)的关系
H ( j) F[h(t)] 1 1 j 1 j 2
1
( j)2 3( j) 2
只有当连续系统是稳定的LTI系统时，才存在H(j), 且可以由h(t)计算出H(j)。
电路系统的频率响应：
分析电路系统的频率响应，主要有两种方法。
H ( j) Yzs ( j)
( j) 3
X ( j) ( j)2 3( j) 2
在实际应用中， 只有当连续系统是稳定的LTI系统时，
才存在H(j),且频响函数才有意义。
例 已知某LTI系统的冲激响应为
h(t) = (e-t-e-2t) u(t)，求系统的频率响应H(j)。
vR (t) RiR (t)
VR ( jw) R IR ( jw)
ZR
VR ( IR(
jw) jw)
R
vL
(t)
L
diL (t) dt
VL ( jw) jwLIL ( jw)
ZL
VL ( jw) IL ( jw)
jwL
iC
(t)
C
d
vC (t) dt
IC ( jw) jwCVC ( jw)
例 已知某LTI系统的动态方程为 y"(t) + 3y'(t) + 2y(t) = x(t)，
[精品]连续时间LTI系统的频率特性及频域分析
![[精品]连续时间LTI系统的频率特性及频域分析](https://img.taocdn.com/s3/m/2bdc77054a73f242336c1eb91a37f111f1850d9f.png)
[精品]连续时间LTI系统的频率特性及频域分析连续时间LTI系统（Linear Time-Invariant System）是指可用于描述各种物理和工程系统运动规律的动态系统。
它们由一对连续时变系统（如模型、结构和控制）和一对线性运算符构成，其具有因变量（响应）和自变量（输入）之间的线性关联性、时间不变性、结构连续的性质，并且在响应上呈现出定义的平稳性，因而它们在描述众多系统运动规律中被广泛应用。
对于连续时间LTI系统的频域特性的研究，则涉及这些系统的相位特性、幅频特性、切趾特性等。
同时，也要探讨系统中不同频率分量的传输特性，因为有不同频率分量的信号既可以幅频分析也可以相位分析，可以衡量系统不同频率下的相应响应。
由于连续时间LTI系统在有限频率通道内传播信号时发生了部分信号丢失，因此我们引入了频域分析得到系统频响阻抗。
这样一来，它就可以用来测量系统频带上的增益，系统的模态表现，以及系统的传播属性和可控特性。
在频域分析过程中，由于信号可以被分解为离散频率分量，所以对于单个频率分量来说，有关连续时间LTI系统的分析可以比较容易地完成。
一般情况下，每一个频率分量的传播特性由一个线性系数连接，称之为频响函数，可以衡量一个系统的频率响应情况。
总的来说，对于连续时间LTI系统，研究其频率特性及频域分析具有重要的意义。
他可以提供一个系统的相位特性、幅频特性、切趾特性等详细的分析，而且由于信号可以分解为离散频率分量，因此可以很容易地实现频域分析，并衡量一个系统的频率响应情况。
此外，还可以利用频域分析来测量系统的增益，模态表现，以及系统的传播属性和可控特性，进而提高系统的性能，实现性能的优化。
实验三连续时间LTI系统的时域分析实验报告
实验三连续时间LTI系统的时域分析实验报告一、实验目的通过实验三的设计和实现，达到如下目的：1、了解连续时间LTI（线性时不变）系统的性质和概念；2、在时域内对连续时间LTI系统进行分析和研究；3、通过实验的设计和实现，了解连续时间LTI系统的传递函数、共轭-对称性质、单位冲激响应等重要性质。
二、实验原理在常见的线性连续时间系统中，我们知道采用差分方程的形式可以很好地表示出该系统的性质和特点。
但是，在本实验中，我们可以采用微分方程的形式来进行相关的研究。
设系统的输入为 x(t)，输出为 y(t)，系统的微分方程为：其中，a0、a1、…、an、b0、b1、…、bm为系统的系数，diff^n(x(t))和diff^m(y(t))分别是输入信号和输出信号对时间t的n阶和m阶导数，也可以记为x^(n)(t)和y^(m)(t)。
系统的单位冲激响应函数 h(t)=dy/dx| x(t)=δ(t),则有：其中，h^(i)(t)表示h(t)的第i阶导数定义系统的传递函数为：H(s)=Y(s)/X(s)在时域内，系统的输出y(t)可以表示为：其中，Laplace^-1[·]函数表示Laplace逆变换，即进行s域到t域的转化。
三、实验步骤1、在Simulink中，构建连续时间LTI系统模型，其中系统的微分方程为：y(t)=0.1*x(t)-y(t)+10*dx/dt2、对系统进行单位冲激响应测试，绘制出系统的单位冲激响应函数h(t)；4、在S函数中实现系统单位冲激响应函数h(t)的微分方程，并使用ODE45框图绘制出系统单位冲激响应函数h(t)在t=0~10s之间的图像；6、利用数据记录栏，记录系统在不同的参数下的变化曲线、阶跃响应函数u(t)和单位冲激响应函数h(t)的变化规律。
四、实验数据分析1、单位冲激响应测试那么，当输入信号为单位冲激函数δ(t)时，根据系统的微分方程，可以得知输出信号的形式为：即单位冲激响应函数h(t)为一个包含了单位冲激函数δ(t)在内的导数项序列。
MATLAB与信号实验——连续LTI系统的时域分析
MATLAB与信号实验-——-连续LTI系统的时域分析在信号处理中，MATLAB是一个强大的工具，它提供了许多功能，使我们能够模拟和分析各种信号系统。
对于连续LTI系统，时域分析是一个重要的方法，它允许我们直接观察系统的输入和输出信号之间的关系。
下面是一个关于连续LTI系统的时域分析的实验。
一、实验目的本实验的目的是验证连续LTI系统的时域响应，通过使用MATLAB模拟系统，我们可以观察到不同的输入信号产生的输出信号，从而了解系统的特性。
二、实验步骤1.定义系统：首先，我们需要定义我们的连续LTI系统。
这可以通过使用MATLAB中的lti函数来完成。
我们需要提供系统的传递函数，它描述了系统的输入和输出之间的关系。
2.设置输入信号：为了观察系统的行为，我们需要设置一个合适的输入信号。
在MATLAB中，我们可以使用square函数来生成一个方波信号，该信号具有固定的频率和幅度。
3.模拟系统：使用MATLAB的lsim函数，我们可以模拟我们的连续LTI系统。
这个函数将输入信号和系统的传递函数作为参数，然后计算出系统的输出信号。
4.分析结果：我们可以使用MATLAB的图形功能来观察输入和输出信号。
这可以帮助我们理解系统的行为，并验证我们的模型是否正确。
三、实验结果与分析在实验中，我们使用了不同的输入信号（如方波、正弦波等）来测试我们的连续LTI系统。
对于每种输入信号，我们都观察了系统的输出信号，并记录了结果。
通过对比不同的输入和输出信号，我们可以得出以下结论：1.对于方波输入，系统的输出信号是带有延迟的方波，这表明系统对突变信号的响应是瞬时的。
2.对于正弦波输入，系统的输出信号是与输入信号同频同相位的正弦波，这表明系统对正弦波的响应是具有稳定性的。
这些结果验证了连续LTI系统的基本特性：即对于单位阶跃函数（突变信号）的输入，系统的响应是瞬时的；而对于周期性输入（如正弦波），系统的响应具有稳定性。
这些结果与我们在理论上学到的知识相符，从而验证了我们的模型是正确的。
实验三 线性时不变(LTI)连续系统的时域分析
执行结果
实验任务 1：LTI系统的微分方程y''(t)+2y'(t)+y(t)=f'(t)+2f(t)，激励f (t)=e-2tε(t)，
（1） 利用 impulse 函数获得冲激响应； （2） 利用 lsim 函数求取零状态响应； （3） 用卷积分析法计算其零状态响应； 要求：在一个图形窗口里以 3 个子图形式绘制冲激响应和两种方法得到的零状 态响应的波形。 （4）改变系统的 a 系数矩阵，观察冲激响应和零状态响应时域波形的变化情况。 建议 a 系数向量分别如下取值讨论。
用。
一、实验目的
1. 掌握系统时域分析常用函数的使用方法； 2. 理解系统特征根对系统时域特性的影响；
二、实验原理及内容 2.1 连续系统的时域分析 2.1.1 连续系统时域分析的几个常用函数
设 LTI 连续时间系统的微分方程为
a 2 y''(t)+a 1 y'(t)+a 0 y(t)=b 2 f''(t)+b 1 f'(t)+b 0 f(t)
将积分变量离散化即将用n替代d用替代只要时域取样间隔足够小上式可近似为再把观察响应时刻离散化即将t用k替换只要足够小通常将fn简记为fnhkn简记为hknyk简记为yk这样上式便可表示为因此两个连续信号ft和ht的卷积yt可用ft和ht的取样信号f均取整数
实验三 线性时不变（LTI）连续系统的时域分析
2.1.2 连续信号卷积的近似计算
连续系统的零状态响应 y(t)可通过输入信号 f(t)与系统冲激响应 h(t)的卷积求 得；但是计算机只能处理数字信号，不能直接处理模拟信号，因此，卷积积分不 能直接用计算机计算。为了解决这个问题，可以将连续信号用取样信号来近似表 示，利用卷积和近似求得卷积积分。下面就连续信号卷积积分的近似计算进行简 单推导。
信号与系统连续时间LTI系统的几种响应求解方法及例题
谢谢您的聆听
THANKS
优点
能够直接得到系统在任意 时刻的响应值。
缺点
计算量大，需要逐个时间 点进行计算。
拉普拉斯变换法
定义
拉普拉斯变换法是一种将时域函 数转换为复频域函数的数学工具。
01
描述ห้องสมุดไป่ตู้
02 通过拉普拉斯变换，将系统的微 分方程转化为代数方程，然后求 解得到系统在复频域的响应。
优点
能够方便地求解高阶微分方程， 适用于具有复杂特性的系统。 03
拉普拉斯变换法
能够求解系统的零状态响应，但需要 已知系统传递函数，且变换过程可能 较为复杂。
05
结论
总结
本文介绍了求解连续时间LTI系统响应的几种方法，包括时域法和频域法。 通过具体实例，展示了这些方法在求解系统响应中的应用和优势。
时域法通过建立和求解微分方程来获取系统输出，具有直观和物理意义 明确的优点。而频域法则通过分析系统函数的频域特性来求解响应，具
信号与系统连续时间LTI系统的 几种响应求解方法及例
CONTENTS
• 引言 • 几种响应求解方法 • 例题解析 • 方法比较与选择 • 结论
01
引言
背景介绍
01
信号与系统是电子工程和通信工 程的重要基础学科，主要研究信 号和系统在时域和频域的行为和 特性。
02
在信号与系统中，线性时不变 （LTI）系统是最基本、最重要的 系统之一，其响应求解是研究的重 要内容。
LTI系统的基本概念
LTI系统是指系统的输出仅与输入和系统 的状态有关，而与时间无关。
LTI系统具有线性、时不变和因果性等基 本特性。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验三连续时间LTI系统分析一、实验目的（一）掌握使用Matlab进行连续系统时域分析的方法1、学会使用符号法求解连续系统的零输入响应和零状态响应2、学会使用数值法求解连续系统的零状态响应3、学会求解连续系统的冲激响应和阶跃响应（二）掌握使用Matlab进行连续时间LTI系统的频率特性及频域分析方法1、学会运用MATLAB分析连续系统的频率特性2、学会运用MATLAB进行连续系统的频域分析（三）掌握使用Matlab进行连续时间LTI系统s域分析的方法1、学会运用MATLAB求拉普拉斯变换（LT）2、学会运用MATLAB求拉普拉斯反变换（ILT）3、学会在MATLAB环境下进行连续时间LTI系统s域分析二、实验条件装有MATLAB的电脑三、实验内容（一）熟悉三部分相关内容原理（二）完成作业1、已知某系统的微分方程如下：)(3)()(2)(3)(t e t e t r t r t r +'=+'+''其中，)(t e 为激励，)(t r 为响应。
（1） 用MATLAB 命令求出并画出2)0(,1)0(),()(3='==---r r t u e t e t 时系统的零状态响应和零输入响应（零状态响应分别使用符号法和数值法求解，零输入响应只使用符号法求解）；符号法求解零输入响应： >> eq='D2y+3*Dy+2*y=0';>> cond='y(0)=1,Dy(0)=2';>> yzi=dsolve(eq,cond);>> yzi=simplify(yzi)yzi =符号法求解零状态响应：exp(-2*t)*(4*exp(t) - 3)eq1='D2y+3*Dy+2*y=Dx+3*x';eq2='x=exp(-3*t)*heaviside(t)';cond='y(-0.001)=0,Dy(-0.001)=0';yzs=dsolve(eq1,eq2,cond);yzs=simplify(yzs)yzs =(exp(-2*t)*(exp(t) - 1)*(sign(t) + 1))/2图像如下：代码：subplot(211)ezplot(yzi,[0,8]);grid ontitle('ÁãÊäÈëÏìÓ¦')subplot(212)ezplot(yzs,[0,8]);grid ontitle('Áã״̬ÏìÓ¦')数值计算法：t=0:0.01:10;sys=tf([1,3],[1,3,2]);f=exp(-3*t).*uCT(t);y=lsim(sys,f,t);plot(t,y),grid on;axis([0 10 -0.001 0.3]);title('ÊýÖµ¼ÆËã·¨µÄÁã״̬ÏìÓ¦')（2）使用MATLAB命令求出并画出系统的冲激响应和阶跃响应（数值法）；用卷积积分法求系统的零状态响应并与（1）中结果进行比较；系统的冲激响应和阶跃响应（数值法）：代码：t=0:0.01:10;sys=tf([1,3],[1,3,2]);h=impulse(sys,t);g=step(sys,t);subplot(211)plot(t,h),grid on;axis([0 10 -0.01 1.1]);title('³å¼¤ÏìÓ¦')subplot(212)plot(t,g),grid on;axis([0 10 -0.01 1.6]);title('½×Ô¾ÏìÓ¦'卷积积分法求系统的零状态响应：Ctsconv函数的定义：function[f,t]=ctsconv(f1,f2,t1,t2,dt)f=conv(f1,f2);f=f*dt;ts=min(t1)+min(t2);te=max(t1)+max(t2);t=ts:dt:te;subplot(221)plot(t1,f1);grid onaxis([min(t1),max(t1),min(f1)-abs(min(f1)*0.2),max(f1)+abs(max(f1)*0.2)])title('f1(t)');xlabel('t')subplot(222)plot(t2,f2);grid onaxis([min(t2),max(t2),min(f2)-abs(min(f2)*0.2),max(f2)+abs(max(f2)*0.2)])title('f2(t)');xlabel('t')subplot(212)plot(t,f);grid onaxis([min(t),max(t),min(f)-abs(min(f)*0.2),max(f)+abs(max(f)*0.2)])title('f(t)=f1(t)*f2(t)');xlabel('t')求系统的零状态响应代码：dt=0.01;t1=0:dt:10;f1=exp(-3*t1).*uCT(t1);t2=t1;sys=tf([1,3],[1,3,2]);f2=impulse(sys,t2);[t,f]=ctsconv(f1,f2,t1,t2,dt)如图，根据两图相比较，两种方法做出的零状态响应大体相同。
（3） 若已知条件同（1），借助MATLAB 符号数学工具箱实现拉普拉斯正反变换的方法求出并画出2)0(,1)0(),()(3='==---r r t u e t e t 时系统的零状态响应和零输入响应，并与（1）的结果进行比较。
普拉斯正反变换的方法求出系统的零状态响应和零输入响应：代码：syms t sRzis=(s+5)/(s^2+3*s+2);rzi=ilaplace(Rzis)rzi =4*exp(-t) - 3*exp(-2*t)et=exp(-3*t)*heaviside(t);es=laplace(et);Rzss=((3+s)*es)/(s^2+3*s+2);rzs=ilaplace(Rzss)rzs =exp(-t) - exp(-2*t)根据图像，同样也能看出拉普拉斯变换法得出的结果相同。
2、已知某RC 网络如下，()r t ()e t +-+-RC（1） 求出该网络的频域系统函数()H j ω；H （jw ）=a/(a+jw) 其中a=1/RC（2） 使用MATLAB 命令画出1RC=时系统的幅频特性和相频特性；代码：w=-3*pi:0.01:3*pi;b=[0,1];a=[1,1];h=freqs(b,a,w);subplot(211)plot(w,abs(h)),grid onaxis([-10 10 0 1.1]);title('H(w)µÄ·ùƵÌØÐÔ')subplot(212)plot(w,angle(h)),grid ontitle('H(W)µÄÏàƵÌØÐÔ')（3） 若1RC =，且激励信号()sin sin(3)e t t t =+，使用频域分析法求解()r t ，分别画出()e t 和()r t 波形，讨论经传输是否引起失真。
代码：t=0:0.1:20;w1=1;w2=3;H1=1/(1+1i*w1);H2=1/(1+1i*w2);f=sin(t)+sin(3*t);y=abs(H1)*sin(w1*t+angle(H1))+abs(H2)*sin(w2*t+angle(H2));subplot(2,1,1);plot(t,f);grid onylabel('f(t)'),xlabel('Time(s)')title('ÊäÈëÐźŵIJ¨ÐÎ')subplot(2,1,2);plot(t,y);grid onylabel('y(t)'),xlabel('Time(sec)')title('ÎÈ̬ÏìÓ¦µÄ²¨ÐÎ')如图，两组波形进行比较可以明显看出，二者不成线性关系，所以此传输系统失真。
3、已知某系统框图如下，∑+-))((121-+s s K)(1s V )2s V（1） 写出下图所示系统的s 域系统函数()H s ；H （s ）=1/(S 2+S-2+K)（2） 使用MATLAB 命令分别用两种方式画出90,1,2,,34K =时该系统的零极点分布图，并由图讨论K 从0增长时，该系统的稳定性变化情况。
代码：b1=[0 1];a1=[1 1 -2];sys1=tf(b1,a1);subplot(321)pzmap(sys1)axis([-2 2 -2 2])b2=[0 1];a2=[1 1 -1];sys1=tf(b2,a2);subplot(322)pzmap(sys1)axis([-2 2 -2 2])b3=[0 1];a3=[1 1 0];sys1=tf(b3,a3);subplot(323)pzmap(sys1)axis([-2 2 -2 2])b4=[0 1];a4=[1 1 0.25];sys1=tf(b4,a4);subplot(324)pzmap(sys1)axis([-2 2 -2 2])b4=[0 1];a4=[1 1 1];sys1=tf(b4,a4);subplot(325)pzmap(sys1)axis([-2 2 -2 2])根据图像，很明显的可以看出，随着K的逐渐增大，系统逐渐稳定。