(改进)上机实验7离散时间信号的时域分析.docVIP

上机实验7 离散时间信号的时域分析 实验目的 掌握离散时间信号时域运算的基本实现方法。 熟悉相关函数的调用格式及运算。 掌握离散信号的基本运算。 掌握信号的分解，会将任意离散信号分解为单位脉冲信号的信号组合。 实验原理 离散时间信号的运算包括信号的相加和相乘。信号的时域变换包括信号的平移、反折、倒相以及尺度变换。 涉及的MATLAB函数 Fliplr函数 功能：实现矩阵行元素的左右翻转。 调用格式： B=fliplr(A):其中A指要翻转的矩阵 四 、实验内容与方法 1. 验证性实验 序列的加法 MATLAB程序： clear; x1=-2:2;k1=-2:2; x2=[1 -1 1];k2=-1:1; k=min([k1,k2]):max([k1,k2]); f1=zeros(1,length(k)); f2=zeros(1,length(k)); f1(find((k=min(k1)k=max(k1))==1))=x1; f2(find((k=min(k2)k=max(k2))==1))=x2; f=f1+f2;subplot(2,3,1); stem(k,f,filled); axis([min(min(k1),min(k2))-1,max(max(k1),max(k2))+1,min(f)-1,max(f)+1]); 两个序列的加法如下所示： 2）序列的乘法 clear all; x1=-2:2; k1=-2:2; x2=[1,-1,1]; k2=-1:1;k=min([k1,k2]):max([k1,k2]); f1=zeros(1,length(k));f2=zeros(1,length(k)); f1(find((k=min(k1))(k=max(k1))==1))=x1; f2(find((k=min(k2))(k=max(k2))==1))=x2; f=f1*f2;stem(k,f,filled); axis([min(min(k1),min(k2))-1,max(max(k1),max(k2))+1,min(f)-0.5,max(f)+0.5); 两个序列的乘法如下所示： 3）序列的翻转 MATLAB程序： clear all; x1=-2:2; k1=-2:2; k=-fliplr(k1); f=fliplr(x1); stem(k,f,filled); axis([min(k)-1,max(k)+1,min(f)-0.5,max(f)+0.5]); 序列的翻转如下所示: 4）序列的倒相 MATLAB程序： clear all; x1=-2:2; k1=-2:2; k=k1; f=-x1; stem(k,f,filled); axis([min(k)-1,max(k)+1,min(f)-0.5,max(f)+0.5]); 序列倒相图如下： 4）序列的平移 MATLAB程序： clear all; x1=-2:2; k1=-2:2; k0=2; k=k1+k0;f=x1; stem(k,f,filled); axis([min(k)-1,max(k)+1,min(f)-0.5,max(f)+0.5]); 序列的平移图如下所示： 2. 程序设计实验 已知序列f(k)={2,3,1,2,3,4,3,1},对应的k值为-3=k=4,分别绘出下列信号的图形： f1(k)=f(k-2),f2(k)=f(-k),f3(k)-f(k-1)U(k),f4(k)=f(-k+2),f5(k)=f(k+1),f6(k)=f(k-2)U(k),f7(k)=f(k+2)U(k) f(k)图形 f1(k)=f(k-2)的程序为： k1=k1+2; subplot(2,4,2); stem(k1,x1,filled); axis([min(k1)-1,max(k1)+1,min(x1)-1,max(x1)+1]); f1(k)的图形 f2(k)=f(-k) 的程序为： k1=-3:4;k1=-k1; subplot(2,4,3); stem(k1,x1,filled); axis([min(k1)-1,max(k1)+1,min(x1)-1,max(x1)+1]); f2(k)的图形 f3(k)=f(k-1)U(k) 的程序为： k1=-3:4;k1=k1+1; k1(find(k10))=0; for temp=1:length(k1) if k1(temp)==0 x1(temp)=x1(3); temp=temp+1; end if temp==length(k1)+1 break; end end subplot(2,4,4); stem(k1,x1,filled); axis([min(k1)-1