《信号分析与处理》实验报告
四、思考： 1、代数运算符号*和.*的区别是？ *是矩阵相乘，是矩阵A行元素与B的列元素相乘的和 .*是数组相乘，表示数组A和数组B中的对应元素相乘 实 验 内 容 实验三 连续时间信号的卷积 一、实验内容 1、已知两连续时间信号如下图所示，绘制信号f1(t)、f2(t)及卷积结果f(t)的波形；设时间变化步长dt分别取为0.5、0.1、0.01，当dt取多少时，程序的计算结果就是连续时间卷积的较好近似？ 《信号分析与处理》实验报告
2、、计算信号f1?t??e?atu?t??a?1?和f2?t??sintu?t?的卷积f(t)，f1(t)、f2(t)的时间范围取为0~10，步长值取为0.1。绘制三个信号的波形。 二、实验方法与步骤 1、绘制信号f1(t)、f2(t)及卷积结果f(t)的波形，当dt取0.01时程序的计算结果就是连续时间卷积的较好近似 程序代码如下： clear all close all clc dt=0.01 t1=0:dt:2; t2=-1:dt:1; f1=0.5*t1; f2=0.5*(t2+1); y=dt*conv(f1,f2); %计算卷积 t0=t1(1)+t2(1); %计算卷积结果的非零样值的起点位置 t3=length(y); %计算卷积结果的非零样值的宽度 ty=t0:dt:(t0+(t3-1)*dt); %确定卷积结果的非零样值的时间向量 subplot(3,1,1) plot(t1,f1) %绘制信号f1(t)的时域波形 title('f1') xlabel('t1') axis([-1.5,3,-0.2,1.1]) grid on subplot(3,1,2) plot(t2,f2) %绘制信号f2(t)的时域波形 《信号分析与处理》实验报告
title('f2') xlabel('t2') axis([-1.5,3,-0.2,1.1]) grid on subplot(3,1,3) plot(ty,y) %绘制两信号卷积结果y(t)的时域波形 title('y') xlabel('ty') axis([-1.5,3.2,-0.2,0.7]) grid on 2、计算信号f1?t??e?atu?t??a?1?和f2?t??sintu?t?的卷积f(t) 程序代码如下： clear all close all clc dt=0.1 t1=0:dt:10; t2=0:dt:10; f1=exp(-1*t1).*heaviside(t1); f2=sin(t2).*heaviside(t2); y=dt*conv(f1,f2); %计算卷积 t0=t1(1)+t2(1); %计算卷积结果的非零样值的起点位置 t3=length(y); %计算卷积结果的非零样值的宽度 ty=t0:dt:(t0+(t3-1)*dt); %确定卷积结果的非零样值的时间向量 subplot(3,1,1) plot(t1,f1) %绘制信号f1(t)的时域波形 title('f1') 《信号分析与处理》实验报告
xlabel('t1') axis([-0.2,10.2,-0.2,1]) grid on subplot(3,1,2) plot(t2,f2) %绘制信号f2(t)的时域波形 title('f2') xlabel('t2') axis([-0.2,10.2,-1.2,1.2]) grid on subplot(3,1,3) plot(ty,y) %绘制两信号卷积结果y(t)的时域波形 title('y') xlabel('ty') axis([-0.2,20.2,-1.2,1.2]) grid on 三、实验数据与结果分析 1.
