信号与系统实验78.docx

信号与系统实验78.docx

《信号与系统实验78.docx》由会员分享，可在线阅读，更多相关《信号与系统实验78.docx（9页珍藏版）》请在冰豆网上搜索。

信号与系统实验78

信号与系统

上

机

仿

真

实

验

报

告

实验七、实验八

班级：

08021103

学号：

2011301886

姓名：

石盼

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

一、实验目的

（1）掌握离散时间信号时域运算的基本实现方法。

（2）熟悉相关函数的调用格式及运算。

（3）掌握离散信号的基本运算。

（4）掌握信号的分解，会将任意离散信号分解为单位脉冲信号的信号组合。

二、实验原理

离散时间信号的运算包括信号的相加和相乘。

信号的时域变换包括信号的平移、反折、倒相以及尺度变换。

三、涉及的MATLAB函数

Fliplr函数

功能：

实现矩阵行元素的左右翻转。

调用格式：

B=fliplr（A）:

其中A指要翻转的矩阵

四、实验内容与方法

1.验证性实验

1）序列的加法

MATLAB程序：

clearall;

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（k1）,min（k2）-1,max（max（k1）,max（k2））+1,min（f）+0.5]）;

结果截图如下

2）序列的翻转

MATLAB程序：

clearall;

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]）;

结果截图

3）序列的倒相

MATLAB程序：

clearall;

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程序：

clearall;

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=

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）的图形

1.f1（k）=f（k-2）的程序为：

clearall;

x1=[2,3,1,2,3,4,3,1];k1=-3:

4;k=k1;f=x1;

stem（k,f,'filled'）;axis（[min（k）-1,max（k）+1,min（f）-0.5,max（f）+0.5]）;

f1（k）的图形

2.f2（k）=f（-k）的程序为：

clearall;

x1=[2,3,1,2,3,4,3,1];k1=-3:

4;

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]）;

f2（k）的图形

3.f3（k）=f（k-1）U（k）的程序为：

clearall;

x1=[2,3,1,2,3,4,3,1];

k1_1=-3:

4;k0=1;k1=k1_1+k0;

N=6;k2=0:

（N-1）;

x2=ones（1,N）;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（k）-1,max（k）+1,min（f）-0.5,max（f）+0.5]）;

f3（k）的图形

4.f4（k）=f（-k+2）的程序为：

clearall;

x1=[2,3,1,2,3,4,3,1];k1=-3:

4;

k0=2;k2=k1+k0;k=-fliplr（k2）;f=fliplr（x1）;

stem（k,f,'filled'）;

axis（[min（k）-1,max（k）+1,min（f）-0.5,max（f）+0.5]）;

f4（k）的图形

5.f5（k）=f（k+1）的程序为：

clearall;

x1=[2,3,1,2,3,4,3,1];k1=-3:

4;

k0=-1;k=k1+k0;f=x1;

stem（k,f,'filled'）;axis（[min（k）-1,max（k）+1,min（f）-0.5,max（f）+0.5]）;

f5（k）的图形

6.f6（k）=f（k-2）U（k）的程序为：

clearall;

x1=[2,3,1,2,3,4,3,1];k1_1=-3:

4;k0=2;k1=k1_1+k0;

N=7;k2=0:

（N-1）;

x2=ones（1,N）;

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（k）-1,max（k）+1,min（f）-0.5,max（f）+0.5]）;

f6（k）的图形

7.f7（k）=f（k+2）U（k）的程序为：

clearall;

x1=[2,3,1,2,3,4,3,1];

k1_1=-3:

4;k0=-2;k1=k1_1+k0;

N=3;k2=0:

（N-1）;

x2=ones（1,N）;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（k）-1,max（k）+1,min（f）-0.5,max（f）+0.5]）;

f7（k）的图形