信号与系统实验课件—卷积和与卷积的计算.pptVIP

实验十六 卷积和/卷积的计算 一、实验目的 1. 熟悉Matlab的卷积和函数conv( )的调用方法； 2. 会计算任意两个离散时间序列的卷积和； 3. 掌握连续时间信号卷积积分的数值计算原理和计算方法。 二、实验原理 （一）离散时间序列的卷积和 n -1 0 1 0 1 2 -2 -1 0 1 2 0 1 n 例16-1 已知序列f1(n)和f2(n)如图16-1所示， 求f(n)=f1(n)*f2(n)，并画出f(n)的波形。 f1=[1 2 1]; n1=[-1 0 1] f2=ones(1,5); %1*5的矩阵 n2=[-2:2] [f,n]=dconv(f1,f2,n1,n2)，%求离散时间序列卷积和函数 stem(n,f) 运行结果: f =1 3 4 4 4 3 1 n =-3 -2 -1 0 1 2 3 需要调用dconv( )函数：求离散时间序列卷积和函数 function [f, n]=dconv(f1, f2, n1, n2) %f1、f2 ：序列的非零样值向量 %n1、n2 ：序列的对应时间序号向量 %f：卷积和序列的非零样值向量 %n：卷积和序列的对应时间序号向量 f=conv(f1, f2); %内置函数，求f1、f2的卷积和 N0=n1(1)+n2(1); %该矩阵的第一个元素（左边界） L=length(n1)+length(n2)-1; %行矩阵n1、n2的长度 n=N0: N0+L-1; 例16-2 已知描述某离散时间系统的差分方程为 a=[1 -2/3]; b=[1]; n1=0:10; %h(n)的时间序号向量 h=impz(b, a, n1); %求离散系统单位样值响应 x=ones(1, 5); n2=0: 4; %激励x(n)的样值向量及对应的时间序号向量 [y,n]=dconv(h,x,n1,n2) %y(n)=h(n)*x(n)求卷积和 subplot(3,1,1), stem(n1, h), hold on,title(h(n)) subplot(3,1,2), stem(n2, x),title(x(n)) subplot(3,1,3), stem(n,y,.k), hold on,title(y(n)) 试画出该系统的单位样值响应h(n)和在x(n)=U(n)-U(n-5)激励下的零状态响应y(n)的波形。 在t=nΔ时刻的取值 （二）连续时间信号的卷积 例16-3 已知f1(t)=U(t+1)-U(t-1)，f2(t)=2t[U(t)-U(t-1)]，画出f(t)=f1(t)*f2(t)的时间波形。 p=0.01; n1=-1:p:1; f1=ones(1,length(n1)); %f1(t) 的时间及样值向量 n2=0:p:1; f2=2*n2; % f2(t)的时间及样值向量 [f, n]=sconv(f1, f2, n1, n2, p); %f(t)=f1(t)*f2(t) 求连续时间系统卷积积分 subplot(3,1,1), plot(n1, f1);title(f1(t)) subplot(3,1,2), plot(n2, f2);title(f2(t)) subplot(3,1,3), plot(n, f)；title(f(t)) 求连续时间信号卷积积分的通用函数sconv( )。 function [f, n]=sconv(f1, f2, n1, n2, p) %计算连续时间信号卷积积分f(t)=f1(t)*f2(t) %f1、f2：f1(t)、f2(t)的非零样值向量 %n1、n2：f1(t) 、f2(t)的对应时间向量 %f：卷积结果f(t)的非零样值向量 %n：f(t)的对应时间向量; %p：取样时间间隔 f=conv(f1, f2); f=f*p; %计算卷积结果f(t)=f1(t)*f2(t)对应的样值向量 N0=n1(1)+n2(1); %左边界 L=length(n1)+length(n2)-1; %行矩阵n1、n2的长度 n=N0:p:N0+(L-1)*p; %生成f(t)的时间向量 1．描述离散时间系统的差分方程为y(n)+y(n-1)=2x(n-1)+x(n-2)，试画出系统的单位样值响应h(n)和单位阶跃响应g(n)的波形。 a=[1,1]; b=[0 2 1];n1=0:10; %的时间序号向量 h=impz(b,a,n1); %求解离散系统单位样值响应 n2=0:9; %激励的时间序号向量 x=ones(1,length(n2));