实验信号的MATLAB表示及系统的时域分析要点.pptVIP

信号的MATLAB 表示及系统的时域分析 实验目的 实验原理一 实验原理二 连续时间信号的卷积 Matlab的库函数conv()只能计算离散序列f1与f2的卷积和，构造函数sconv()实现两连续函数的卷积积分。 举例1 function[f,k]=sconv(f1,f2,k1,k2,p) ％计算连续信号卷积积分f(t)=f1(t)*f2(t) ％f: 卷积积分f(t)对应的非零样值向量 ％k: f(t)的对应时间向量 ％f1: f1(t)非零样值向量 ％f2: f2(t)非零样值向量 ％k1: f1(t)的对应时间向量 ％k2: f2(t)的对应时间向量 ％p: 取样时间间隔 续： f=conv(f1,f2); 计算序列f1和f2的卷积和f f=f*p; k0=k1(1)+k2(1); 计算序列f非零样值的起点位置 k3=length(f1)+length(f2)-2 ％计算卷积和f非零样值的宽度 k=k0:p:k3*p ％确定卷积和f非零样值的时间量 subplot(2,2,1) plot(k1,f1) ％在子图1绘出f1(t)的时域波形图 title(‘f1(t)) xlabel(‘t) ylabel(‘f1(t)) 续： subplot(2,2,2) plot(k2,f2) ％在子图2绘出f2(t)的时域波形图 title(‘f2(t)) xlabel(‘t) ylabel(‘f2(t)) subplot(2,2,3) plot(k,f) ％画卷积f(t)的时域波形图 h=get(gca,‘position’) h(3)=2.5*h(3) set(gca,‘position’，h) ％将第三个子图的横坐标范围扩为原来的2.5倍 title(‘f(t)=f1(t)*f2(t)) xlabel(‘t’) ylabel(‘f(t)’) 举例1 已知两个连续时间信号如下图所示，试用MATLAB求f(t)=f1(t)*f2(t)，并绘出f(t)的时域波形图。 MATLAB命令如下： p=0.01; k1=0:p:2; f1=0.5*k1; k2=k1; f2=f1; [f,k]=sconv(f1,f2,k1,k2,p) 仿真结果 实验原理三 连续时间系统零状态响应的MATLAB实现 Matlab库函数中的lsim()能对微分方程描述的LIT连续时间系统响应进行仿真。 函数lsim()的调用格式： lsim(b,a,x,t) 该调用格式中，a和b是由描述系统的微分方程左边和右边系数构成的两个行向量；t为输入信号时间范围的向量，x为输入信号在向量t定义的时间点上的抽样值 举例2 解：MATLAB命令如下： a=[1 2 1]; b=[1 2]; p=0.5; %定义取样时间间隔 t=0:p:5; %定义时间范围 x=exp(-2*t); %定义输入信号 lsim(b,a,x,t); %对系统输出信号进行仿真 hold on p=0.3; t=0:p:5; x=exp(-2*t); lsim(b,a,x,t); p=0.01; t=0:p:5; x=exp(-2*t); lsim(b,a,x,t); hold off legend(p=0.5, p=0.3,p=0.01); 系统零状态响应的仿真波形图如图2.4所示：系统零状态响应的仿真波形图如下图所示： 仿真结果 实验原理四 连续时间系统冲激响应的MATLAB实现 MATLAB为用户提供了专门用于连续系统冲激响应并绘制其时域波形的impulse函数 函数impulse()的调用格式： impulse(b,a) 连续时间系统阶跃响应的MATLAB实现 MATLAB为用户提供了专门用于连续系统