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

1、上机实验7 离散时间信号的时域分析一、 实验目的(1) 掌握离散时间信号时域运算的基本实现方法。(2) 熟悉相关函数的调用格式及运算。(3) 掌握离散信号的基本运算。(4) 掌握信号的分解，会将任意离散信号分解为单位脉冲信号的信号组合。二、 实验原理离散时间信号的运算包括信号的相加和相乘。信号的时域变换包括信号的平移、反折、倒相以及尺度变换。三、 涉及的matlab函数fliplr函数功能：实现矩阵行元素的左右翻转。调用格式：b=fliplr(a):其中a指要翻转的矩阵四 、实验内容与方法1. 验证性实验1) 序列的加法matlab程序：clear;x1=-2:2;k1=-2:2;x2=1 -

2、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); 两个序列的加法如下所示

3、：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(

4、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,ma

5、x(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),f

6、6(k)=f(k-2)u(k),f7(k)=f(k+2)u(k) f(k)图形1. 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)的图形2. 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)的图形3.

7、 f3(k)=f(k-1)u(k) 的程序为：k1=-3:4;k1=k1+1;k1(find(k1<0)=0;for temp=1:length(k1)if k1(temp)=0 x1(temp)=x1(3); temp=temp+1;endif temp=length(k1)+1 break;endendsubplot(2,4,4);stem(k1,x1,'filled');axis(min(k1)-1,max(k1)+1,min(x1)-1,max(x1)+1);f3(k)的图形4. f4(k)=f(-k+2) 的程序为：clear all;x1=2 3 1 2 3

8、4 3 1;k1=-3:4;k1=-k1;k1=k1+2;subplot(2,4,5);stem(k1,x1,'filled');axis(min(k1)-1,max(k1)+1,min(x1)-1,max(x1)+1);f4(k)的图形5. f5(k)=f(k+1) 的程序为：k1=-3:4;k1=k1-1;subplot(2,4,6);stem(k1,x1,'filled');axis(min(k1)-1,max(k1)+1,min(x1)-1,max(x1)+1);f5(k)的图形6. f6(k)=f(k-2)u(k) 的程序为：k1=-3:4;k1=k1

9、+2;subplot(2,4,7);k1(find(k1<0)=0;for temp=1:length(k1)if k1(temp)=0 x1(temp)=x1(2); temp=temp+1;endif temp=length(k1)+1 break;endendstem(k1,x1,'filled');axis(min(k1)-1,max(k1)+1,min(x1)-1,max(x1)+1);f6(k)的图形7. f7(k)=f(k+2)u(k) 的程序为：k1=-3:4;k1=k1-2;subplot(2,4,8);k1(find(k1<0)=0;for temp=1:length(k1)if k1(temp)=0 x1(temp)=x1(6); temp=temp+1;endif temp=length(k1)+1 break;endendstem(k1,x1,'filled&#