信号处理实验报告Word文档格式.docx

信号处理实验报告Word文档格式.docx

《信号处理实验报告Word文档格式.docx》由会员分享，可在线阅读，更多相关《信号处理实验报告Word文档格式.docx（10页珍藏版）》请在冰豆网上搜索。

z变换和零极点的求法

实验目的：

1、掌握z变换和反z变换的应用

2、掌握用局部分式法求z变换

3、掌握零极点的求法

实验步骤与内容：

1、求Z变换。

2.求z反变换。

3.

4.

1.将上述的程序进展验证，并进展思考，掌握其编写的方法，将程序和结果复制到下面：

1、

symszaw0n

x1=a^n;

X1=ztrans（x1）

x2=n;

X2=ztrans（x2）

x3=（n*（n-1））/2;

X3=ztrans（x3）

x4=exp（j*w0*n）;

X4=ztrans（x4）

x5=1/n*（n-1）;

X5=ztrans（x5）

X1=-z/（a-z）

X2=z/（z-1）^2

X3=（z^2+z）/（2*（z-1）^3）-z/（2*（z-1）^2）

X4=z/（z-exp（i*w0））

X5=z/（z-1）-ztrans（1/n,n,z）

2、

symsnza

X1=z/（z-1）;

x1=iztrans（X1）

X2=a*z/（a-z）^2;

x2=iztrans（X2）

X3=z/（z-1）^3;

x3=iztrans（X3）

X4=（1-z^-n）/（1-z^-1）;

x4=iztrans（X4）

x1=1

x2=piecewise（[a<

>

0,a*（a^n/a-kroneckerDelta（n,0）/a）+a^2*（kroneckerDelta（n,0）/a^2+（a^n*（n-1））/a^2）]）

x3=n+binomial（n-1,2）-1

x4=iztrans（1/（z^n*（1/z-1））,z,n）+1

3、

b=[1,0,0];

a=[1,-1.5,0.5];

[rpc]=residuez（b,a）

N=20;

n=0:

N-1;

x=r

（1）*p

（1）.^n+r

（2）*p

（2）.^n

stem（n,x）

结果

4、

b=[0.2,0.1,0.3,0.1,0.2]

a=[1,-1.1,1.5,-0.7,0.3]

rz=roots（b）

rp=roots（a）

subplot（1,2,1）,zplane（b,a）;

subplot（1,2,2）,impz（b,a,20）;

title（'

Ï

µ

Í

³

Ä

å

¼

¤

ì

Ó

¦

'

）;

xlabel（'

n'

）,ylabel（'

h（n）'

）

2.编写程序实现

1.求z变换：

〔1〕

〔2〕

代码：

x1=n*a^n;

x2=sin（w0*n）;

结果：

X1=z/（a*（z/a-1）^2）

X2=（z*sin（w0））/（z^2-2*cos（w0）*z+1）

2.求z反变换：

〔2〕

X1=z/（z-a）;

X2=z/（z-a）^2;

x1=piecewise（[a<

0,a*（a^n/a-kroneckerDelta（n,0）/a）+kroneckerDelta（n,0）]）

0,a*（kroneckerDelta（n,0）/a^2+（a^n*（n-1））/a^2）+a^n/a-kroneckerDelta（n,0）/a]）

3.用局部分式法求解以下系统函数的z反变换，写出x〔n〕的表示式，用图形表示出来，取前10个点作图。

c=[10,20,0,0]

d=[1,8,19,12]

[r,p,c]=residuez（c,d）

N1=10;

N1-1;

x=r

（1）*p

（1）.^n+r

（2）*p

（2）.^n+r（3）*p（3）.^n

4.离散时间系统函数分为：

求该系统的零极点分布图，并判断系统的因果稳定性。

b=[4,-1.6,-1.6,4]

a=[1,0.4,0.35,-0.4]