通信11002班16号李俊良实验一 序列的基本运算.docx

通信11002班16号李俊良实验一 序列的基本运算.docx

《通信11002班16号李俊良实验一 序列的基本运算.docx》由会员分享，可在线阅读，更多相关《通信11002班16号李俊良实验一 序列的基本运算.docx（12页珍藏版）》请在冰豆网上搜索。

通信11002班16号李俊良实验一序列的基本运算

一、实验目的

1．在理论学习的基础上，通过本实验，加深对序列基本运算的理解；

2．熟悉利用matlab软件做简单的仿真实验；

例1用MATLAB下列序列的离散卷积：

{-201–13}，{120-1}

解MATLAB程序如下：

a=[-201-13];

b=[120-1];

c=conv（a,b）;

M=length（c）-1;

n=0:

1:

M;

stem（n,c）;

xlabel（'n'）;ylabel（'幅度'）;

运行结果如下：

卷积结果的图形，求得的结果存放在数组c中为：

{-2-413151-3}。

分析：

对着例题代码画图soeasy

第二题：

N=41;

a=[0.8-0.440.360.22];

b=[10.7-0.45-0.6];

x=[1zeros（1,N-1）];

k=0:

1:

N-1;

y=filter（a,b,x）;

stem（k,y）

xlabel（'n'）;ylabel（'幅度'）

运行结果如下：

N=41;

a=[0.8-0.440.360.22];

b=[10.7-0.45-0.6];

x=[1zeros（1,N-22）];!

N-22

k=0:

1:

N-22;

y=filter（a,b,x）;

stem（k,y）

xlabel（'n'）;ylabel（'幅度'）

N=41;

a=[0.8-0.440.360.22];

b=[10.7-0.45-0.6];

x=[1zeros（1,N+22）];!

N+22

k=0:

1:

N+22;

y=filter（a,b,x）;

stem（k,y）

xlabel（'n'）;ylabel（'幅度'）

分析：

N的值变大了就越密集了，N的值变小了就越稀疏了！

第三题：

0到pai

分析：

画图比较简单。

0到2pai

k=256;

num=[0.8-0.440.360.02];

den=[10.7-0.45-0.6];

w=0:

2*pi/k:

2*pi;

h=freqz（num,den,w）;

subplot（2,2,1）;

plot（w/pi,real（h））;grid

title（'实部'）

xlabel（'\omega/\2*pi'）;ylabel（'幅度'）

subplot（2,2,2）;

plot（w/2*pi,imag（h））;grid

title（'虚部'）

xlabel（'\omega/\2*pi'）;ylabel（'Amplitude'）

subplot（2,2,3）;

plot（w/2*pi,abs（h））;grid

title（'幅度谱'）

xlabel（'\omega/\2*pi'）;ylabel（'幅值'）

subplot（2,2,4）;

plot（w/pi,angle（h））;grid

title（'相位谱'）

xlabel（'\omega/\2*pi'）;ylabel（'弧度'）

分析：

改变pai的值就可以了！

第五题

i=0:

1:

49;

x=[sin（2*pi*i/50）];

h1=[zeros（1,10）,1,zeros（1,20）];

h2=[zeros（1,20）,1,zeros（1,10）];

h3=h1+h2;

y1=conv（x,h1）;

y2=conv（x,h2）;

y3=conv（x,h3）;

M=length（y1）-1;

n=0:

1:

M;

m1=length（h1）-1;

n1=0:

1:

m1;

m2=length（h2）-1;

n2=0:

1:

m2;

m3=length（h3）-1;

n3=0:

1:

m3;

!

stem（x,y1,x,y2）;

!

xlabel（'x'）;ylabel（'幅度'）;

figure,

subplot（3,1,1）;

stem（n,y1）;

title（'卷积y1（n）'）

xlabel（'n'）;ylabel（'幅度'）

subplot（3,1,2）;

stem（n,y2）;

title（'卷积y2（n）'）

xlabel（'n'）;ylabel（'幅度'）

subplot（3,1,3）;

stem（n,y3）;

title（'卷积y3（n）'）

xlabel（'n'）;ylabel（'幅度'）

figure,

subplot（3,1,1）;

stem（n1,h1）;

title（'h1（n）'）

xlabel（'n'）;ylabel（'幅度'）

subplot（3,1,2）;

stem（n2,h2）;

title（'h2（n）'）

xlabel（'n'）;ylabel（'幅度'）

subplot（3,1,3）;

stem（n3,h3）;

title（'h3（n）'）

xlabel（'n'）;ylabel（'幅度'）

figure,

stem（i,x）;

title（'x（n）'）

xlabel（'n'）;ylabel（'幅度'）

分析：

画图需要时间比较长！

线性相关不会！