实验四 离散时间信号的时域基本运算.pptVIP

实验四 离散时间信号的时域基本运算 实验目的 通过本实验的练习，了解离散时间信号时域运算的基本实现方法。 了解函数fliplr的调用格式及作用。 通过该实验，掌握离散时间信号时域基本运算的原理及编程思想。 实验内容 验证性实验 采用fliplr函数编程，实现对离散序列信号的反折运算； 利用向量的平移运算，实现对离散信号的平移运算； 实现对离散序列信号的平移、反折运算的综合。 程序设计实验。 预备知识 向量和矩阵的左右反折运算用fliplr函数来实现，调用格式如下 y=fliplr（x） 例如： 结果变为： 实验要求 在计算机中输入程序，验证实验结果，并将实验结果存入指定存储区域。 对于程序设计实验，要求通过对验证性实验的练习，自行编制完整的实验程序，实现对信号的模拟，并得出实验结果。 在实验报告中写出完整的自编程序，并给出实验结果。 实验步骤 一）验证性实验 离散序列的时域运算和变换不能用符号运算来实现，而必须用向量的表示方法，因此参加运算的两序列向量必须有相同的维数。 1、反折 向量的反折用fliplr来实现，实现序列反折的子程序如下所示： function[f,k]=lsfz(f1,k1) f=fliplr(f1) k=-fliplr(k1) stem(k,f,filled) axis([min(k)-1,max(k)+1,min(f)-0.5,max(f)+0.5]) 例：绘出离散序列的波形及其反折后的波形。 k=-3:3; f=2.^k subplot(1,2,1) stem(k,f) axis([-4,4,-0.5,8.5]) subplot(1,2,2) lsfz(f,k) 2、平移 离散序列的平移可认为是将离散序列的时间序号向量平移，而序列样值不变。 子函数程序如下： function[f,k]=lsyw(ff,kk,k0) k=kk+k0 f=ff stem(k,f,filled) axis([min(k)-1,max(k)+1,min(f)-0.5,max(f)+0.5]) %ff,kk分别为平移前离散序列样值及其对应的时间向量序列，k0为平移量。f,k分别为平移后的序列样值和时间向量 例：绘出序列、f(k)及f(k-2)的波形 k=-4:4 f=k.^2 subplot(2,1,1) stem(k,f,filled) axis([-5,5,-0.5,16.5]) subplot(2,1,2) lsyw(f,k,2) 形如 序列移位的实现 已知序列 分别绘出下列信号的波形。 k=-4:3 f=[3 3 1 2 3 3 3 3] subplot(2,3,1) stem(k,f) title(f(k)) f1=f subplot(2,3,2) [f2,k]=lsanpw(f1,k,2) % 移位运算f(k-2) title(f1(k)) subplot(2,3,3) f2 k m=0 % 求f(k-2)因果信号 n=length(k) for i=1:n if(k(i)0) m=m+1 end end m for j=1:m f2(j)=0 end f2 subplot(2,3,3) stem(k,f2) title(f2(k)) 3.反折、平移运算的综合实现 已知序列， 分别绘出下列信号的图形。 k=-3:4 f=[2 3 1 2 3 4 3 1] subplot(2,3,1) stem(k,f) title(f(k)) f1=f subplot(2,3,2) lsanpw(f1,k,2) % 移位运算f(k-2) title(f1(k)) subplot(2,3,3) f2=f lsanfz(f2,k) % 反折运算 title(f2(k)) subplot(2,3,4) f3=f [f4,k]=lsanpw(f3,k,-2) % f(k+2) lsanfz(f4,k) % f(-k+2) title(f4(k)) subplot(2,3,5) k=-3:4 % 重新设定输入序列，以防内部变量混淆 f=[2 3 1 2 3 4 3 1] f5=f [f6,k]=lsanpw(f5,k,2) % f(k-2) title(f(k-2)) subplot(2,3,6) m=0 % 求f(k-2)因果信号 n=length(k) for i=1:n if(k(i)0) m=m+