文档介绍：实验三连续时间系统的时域分析
[实验目的]
1、掌握卷积计算方法及其性质。
2、运用课堂上学到的理论知识，从RC、RL一阶电路的响应中正确区分零输入响应、零状态响应、冲激响应和阶跃响应。
[实验原理]
描述线性非时变连续时间系统的[0,1]
hold off
title(单位阶跃信号)
axis([t1,t2,-,])
jieyue(-2,5,0)
程序执行后绘出波形如图2所示：
图2
二. 连续时间信号卷积的MATLAB实现
图3
解：首先自编通用函数sconv()实现连续时间信号的卷积，该程序在计算出卷积积分的数值近似的同时，还绘出f(t)的时域波形图。需要注意的是程序中如何构造f(t)的对应时间向量k。另外，程序在绘制f(t)波形图时采用的是plot命令而不是stem命令。 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)=*h(3)
set(gca,‘position’，h) ％
title(‘f(t)=f1(t)*f2(t)')
xlabel(‘t’)
ylabel(‘f(t)’)