信号与系统之 基础演练.docx

1、信号与系统之 基础演练实验一 连续系统的时域分析一、实验目的（1）掌握零状态响应和冲激响应的计算。（2）了解卷积的计算。二、实验原理连续时间线性时不变系统（LTI），可以用线性常系数微分方程描述。系统分析即求解一定输入信号所产生的输出响应，并根据响应判断系统特性。系统响应可分为由初始状态产生的零输入响应和由输入信号产生的零状态响应。当输入信号为冲激信号、阶跃信号时，产生的零状态响应分别为冲激响应和阶跃响应。（1） 直接求解建立系统函数模型tf、零状态响应lsim、冲激响应impulse、阶跃响应step。这些函数属于Control System Toolbox。系统方程的系数将作为建立系统函数

2、的参数。零输入响应较少用到，可利用roots来求。（2） 卷积求解系统的零状态响应，等于其冲激响应与输入激励信号的卷积。采用数值计算，用离散序列的卷积和conv近似这实际是用宽度为1的矩形脉冲之和代替积分，因此最后的结果需要乘上实际的数据间隔时间。这是一个近似结果，间隔时间越窄，结果越接近积分值。利用卷积求响应主要意义体现在理论上，实际做数值计算时使用现成函数更简便，也更准确。此外，利用符号计算可以求得响应的表达式。三、验证性实验已知系统。（1） e(t)=u(t)时，求零状态响应程序：b=3 9; a=1 6 8;sys=tf(b,a); % Create transfer function

3、 modelt=0:0.1:10;e=ones(1,length(t);r=lsim(sys,e,t); % Simulate LTI model response%figure(Position,100,100,400,200);plot(t,r); 问题：本例所求响应是原方程的特解还是齐次解？方程的特征根有哪些？为什么响应最终趋于常数？总响应中有没有暂态分量？（2） 求冲激响应h(t)程序：b=3 9; a=1 6 8;sys=tf(b,a); % Specify transfer functionst=0:0.1:10;h=impulse(sys,t); % Compute the im

4、pulse response of LTI modelsplot(t,h);问题：本例所求响应是原方程的什么解？为什么响应最终趋于0？系统的因果性、稳定性如何？（3） e(t)=u(t)时，用卷积法求零状态响应程序：clear all;b=3 9; a=1 6 8;sys=tf(b,a); % Specify transfer functionsdt=0.1; tf=10;t=0:dt:tf;h=impulse(sys,t); % Compute the impulse response of LTI modelse=ones(1,length(t);r=conv(h,e); % Convol

5、ution multiplicationr=r*dt;tr=0:dt:2*tf;plot(tr,r);axis(0,tf,0,1.5);问题：该例结果与例（1）是否一致？原因何在？四、设计性实验已知系统，当e(t)=cos(t)u(t)时，求系统的零状态响应。将得到的响应与激励作对比分析。（直流分量，幅度、相位、频率等）还可以尝试改变激励信号的幅度、频率，看响应有何变化。五、实验要求1运行验证性实验，观察记录结果，回答问题。2完成设计性实验，记录程序和结果，并分析。六、补充：符号计算（1）已知系统，求零输入响应。r=dsolve(D2r+6*Dr+8*r=0, r(0)=1,Dr(0)=0,

6、t)r = 2/exp(2*t) - 1/exp(4*t)（2）冲激信号dirac(x)抽样syms xint(dirac(x-5)*sin(x),-inf, inf)ans = sin(5)（3）阶跃信号heaviside(t)求导syms x;diff(heaviside(x),x)ans = dirac(x)实验二 离散系统的分析一、实验目的（1）熟悉零状态响应和冲激响应的计算。（2）熟悉z变换。（3）熟悉系统的频率响应。二、实验原理（1）离散线性时不变系统（LTI），可以用线性常系数差分方程描述。当输入信号为冲激信号时，产生的零状态响应为冲激响应。（2）离散LTI系统，在z域可以用系统

7、函数H(z)描述，其实质是冲激响应h(n)的z变换。（3） 如果系统函数H(z)的收敛域包含单位圆，则令，得到系统的频率响应。三、验证性实验（1）已知系统y(n)-y(n-1)+0.8y(n-2)=x(n)，当时，求零状态响应。程序：b=1; a=1,-1,0.8; % 差分方程系数n=0:15; % 输入序列长度xn=0.5.n; % 输入序列yn=filter(b,a,xn); % 零状态响应figure(Position,100,100,400,200);stem(n,yn,filled);grid on问题：求方程的特征根，分析其零状态响应（齐次解+特解）的特征。（2）已知下式，求其部

8、分分式展开程序：clear;p=0.6,0.3;a=poly(p); %分母多项式系数b=2.5,-0.9; %分子多项式系数r,p,k=residue(b,a)运行结果：r = 2.0000 0.5000p = 0.6000 0.3000k = 则问题：根据极点分布判断该系统的稳定性。（3） 已知系统函数画出零、极点图；求冲激响应；求幅频响应和相频响应。b=0,1,2,1; a=1,0.5,-0.005,0.3;figure(Position,100,100,400,300);subplot(2,2,1),zplane(b,a); % 画零极点title(Zero-Pole);hn n=im

9、pz(b,a,16); % 求冲激响应subplot(2,2,2),stem(n,hn,filled);title(h(n);grid on;h w=freqz(b,a,16); % 求频率响应subplot(2,2,3),plot(w,abs(h);title(Amplitude);grid on;subplot(2,2,4),plot(w,angle(h);title(Phase);grid on;问题：该系统是何种滤波器？对信号延迟还是超前？四、设计性实验已知系统函数画出零、极点图，求冲激响应，并判断系统稳定性；求幅频响应，并判断系统属于就何种滤波器。 求相频响应，分析系统对信号的相移。

10、五、实验要求1运行验证性实验，观察记录结果，回答问题。2完成设计性实验，在实验报告上记录程序和结果。实验三 信号的抽样和重构实验目的（1）熟悉抽样信号及其频谱。（2）掌握抽样定理。（3）了解理想低通滤波器。一、实验原理1抽样信号抽样信号相当于连续信号与周期性的冲击序列相乘。 在Matlab中可以很方便的用不同的时间间隔实现对连续信号不同频率的抽样。抽样信号的频谱等于原始信号的频谱与冲击序列的频谱的卷积。 抽样信号的频谱是对原始信号的频谱的周期性延拓，周期大小为抽样品率，其中每一个周期都复制了原始信号的频谱。2抽样定理一个带宽为wm的带限信号f(t)，可唯一地由它的均匀取样信号fs(nTs)确定

11、，其中，取样间隔Ts1/的分量所具有的能量占有很少的比重，所以一般定义fm=1/为门信号的截止频率。其中的为门信号在时域的宽度。在本例中选取fm=0.5，临界采样频率为fs=2fm=1，过采样频率为fs1（为了保证精度，可以将其值提高到该值的50倍），欠采样频率为fs2fm=1t=-4:Ts:4;f=rectpuls(t,2);% 宽度为2的门信号w1=2*pi*10; % 频谱范围-20*pi 20*piN=1000; % 计算出2*1000+1个频率点k=0:N;wk=k*w1/N;F=f*exp(-j*t*wk)*Ts; % 计算Fourier变换F=abs(F); % 计算频谱的幅度w

12、k=-fliplr(wk),wk(2:1001);F=fliplr(F),F(2:1001); % 补充对应负频率的频谱subplot(2,1,1); plot(t,f);xlabel(t); ylabel(f(t);title(f(t)=u(t+1)-u(t-1);subplot(2,1,2); plot(wk,F);xlabel(w); ylabel(F(jw);title(f(t)的幅度谱);由于抽样周期很小，抽样频率（100Hz）远远大于信号带宽，结果相当于连续信号。请回答：门信号的频谱第一个过零点是多大的频率？2绘制宽度为2的门信号G2(t)=u(t+1)-u(t-1)的的抽样信号的

13、图形和频谱。MATLAB程序：Ts=0.25;%采样周期=0.25,fs=42fm=1t=-4:Ts:4;f=rectpuls(t,2);% 宽度为2的门信号w1=2*pi*10; % 频谱范围-20*pi 20*piN=1000; % 计算出2*1000+1个频率点k=0:N;wk=k*w1/N;F=f*exp(-j*t*wk); % 计算Fourier变换F=abs(F); % 计算频谱的幅度wk=-fliplr(wk),wk(2:1001);F=fliplr(F),F(2:1001); % 补充对应负频率的频谱subplot(2,1,1); stem(t/Ts,f); % 采样信号的离散

14、时间显示xlabel(n); ylabel(f(n); title(f(n);subplot(2,1,2); plot(wk,F); % 显示采样信号的连续的幅度谱xlabel(w); ylabel(F(jw); title(f(n)的频谱);请回答：抽样序列f(n)的脉冲宽度是多少，为什么不是2？抽样信号的频谱周期是多少？3求宽度为2的门信号G2(t)=u(t+1)-u(t-1)的的抽样信号，利用低通滤波器重构门信号，并计算误差。MATLAB程序：%计算抽样信号Ts=0.25; % 采样周期t=-4:Ts:4; %时间区间-4，4f=rectpuls(t,2);%抽样信号%计算抽样信号的频谱

15、w1=2*pi*10;N=1000;k=0:N;wk=k*w1/N;F=f*exp(-j*t*wk); % 利用数值计算连续信号的Fourier变换wk=-fliplr(wk),wk(2:1001);F=fliplr(F),F(2:1001);%抽样信号的频谱%计算理想低通滤波器Tw=w1/N; % 频率采样间隔w=-2*pi*10:Tw:2*pi*10;H=Ts*rectpuls(w,2.*pi/Ts); % 理想低通滤波器频率特性% 重构原始信号Fa=F.*H; % 采样信号通过滤波器后的频谱Dt=0.01;t1=-4:Dt:4;%时间区间-4，4fa=Tw/(2*pi)*(Fa*exp(

16、j*wk*t1);% Fourier逆变换%绘重构信号subplot(2,1,1) plot(t1,fa);xlabel(t); ylabel(fa(t);title(f(t)的重构信号);%绘误差图err=fa-rectpuls(t1,2);%误差subplot(2,1,2) plot(t1,err); xlabel(t);ylabel(err(t);title(f(t)的重构误差信号);三、设计性实验求宽度为2的门信号G2(t)=u(t+1)-u(t-1)的的抽样信号，利用低通滤波器重构门信号，并计算误差。比较抽样周期分别为0.1、0.5、1、1.5时的不同结果，分析其原因。四、实验要求1运行验证性实验，观察记录结果。2完成设计性实验，在实验报告上记录程序和结果。