信号与系统实验10Word下载.docx

信号与系统实验10Word下载.docx

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

MatlabR2016

三、实验内容及分析

（a）

代码如下：

a=[132];

b=[1-10];

draw（b,a）;

title（'

H（z）µ

Ä

Á

ã

¼

«

µ

Í

'

）;

draw.m:

functiondraw（num,den）

p=roots（den）;

q=roots（num）;

p=p'

;

q=q'

x=max（abs（[pq]））;

x=x+1;

y=x;

plot（[-xx],[00],'

:

holdon;

plot（[00],[-yy],'

plot（real（p）,imag（p）,'

*'

plot（real（q）,imag（q）,'

o'

xlabel（'

Ê

²

¿

ylabel（'

Ð

é

axis（[-xx-yy]）;

（b）

a=[110.5];

b=[1];

Â

Ë

¨

Æ

÷

（c）

a=[1-1.250.75-0.125];

b=[10.5];

b=[100];

a=[1-0.90.81];

zs=roots（b）;

ps=roots（a）;

plot（real（zs）,imag（zs）,'

plot（real（ps）,imag（ps）,'

x'

gridon;

axis（[-44-44]）;

H1（z）µ

[b,a]=ellip（4,.2,40,[.41.47]）;

[H,w]=freqz（b,a,4096）;

plot（w/pi,20*log10（abs（H）））;

axis（[01-8010]）;

Hµ

¶

Ô

ý

·

ù

Ö

axis（[01-16020]）;

[b,a]=ellip（4,.2,40,[.41.47]）;

M=max（abs（[ba]））;

a16=quant（a,16,M）;

b16=quant（b,16,M）;

[H,w]=freqz（b16,a16,4096）;

quant.m:

functionqc=quant（x,N,M）

[mm,nn]=size（x）;

qc=zeros（mm,nn）;

levels=2^（N-1）;

maxlevel=2^N-1;

fork=1:

mm,

tmp=fix（（x（k,:

）+M）./（M/levels））

q=zeros（1,nn）;

q（tmp<

=maxlevel）=tmp（tmp<

=maxlevel）;

q（tmp>

maxlevel）=maxlevel*ones（1,length（tmp（tmp>

maxlevel）））;

0）=zeros（tmp（tmp<

0））;

q=（q-levels）*M/levels;

qc（k,:

）=q;

end

（d）

（e）

draw（a16,b16）

a1p=[1/8];

b1p=[110.51/8];

w=[-10:

0.25:

10];

[H,w]=freqz（b1p,a1p,w）;

plot（w,abs（H））;

Ì

Ø

figure;

draw（b1p,a1p）;

H1p（s）µ

四、实验总结

掌握了如何使用fft、freqz、compose、dpzplot简单函数的使用，学会了使用Matlab实现Z变换。

教师评语：

签名：

日期：

成绩：