29

第2题： fc=1000; Fs=10000;

Wn=fc/(Fs/2); %归一化 n=4; Rs=30;

[b,a]=cheby2(n,Rs,Wn,'high','s'); sys=tf(b,a) freqs(b,a); 结果：

s^4 - 1.42e-017 s^3 + 0.04 s^2 - 2.938e-018 s + 0.0002 ------------------------------------------------------ s^4 + 0.6442 s^3 + 0.2475 s^2 + 0.05581 s + 0.006325

第3题： fc1=1000; fc2=3000; Fs=10000;

Wn=[fc1,fc2]/(Fs/2); %归一化 n=6;

[b,a]=butter(n/2,Wn,'bandpass'); sys=tf(b,a)

freqz(b,a,512); 结果：

Transfer function:

0.09853 - 0.2956 z^-2 + 0.2956 z^-4 - 0.09853 z^-6

《Matlab语言及其在电子信息科学中的应用》实验指导书

30

--------------------------------------------------------- 1 - 1.366 z^-1 + 1.245 z^-2 - 0.8777 z^-3 + 0.6537 z^-4 - 0.2256 z^-5 + 0.0563 z^-6

第4题：

[x,Fs,Bit]=wavread('Test2.wav'); Wn=500/(Fs/2); %归一化 n=6;

[b,a]=butter(n,Wn); y=filter(b,a,x);

subplot(211); plot(x) title('原始语音信号'); subplot(212); plot(y)

title('500Hz以下子带信号'); 结果：

《Matlab语言及其在电子信息科学中的应用》实验指导书

31

第5题： n=65;

Wn=[0.45,0.65,0.85]; b=fir1(n,Wn,'DC-1'); a=1;

sys=tf(b,a)

freqz(b,a,512); 结果：

《Matlab语言及其在电子信息科学中的应用》实验指导书

32

实验十：数字图像处理初步

一、实验目的

1．掌握图像文件的读取与写入。

2．掌握图像文件的显示、翻转、缩放、叠加

3．掌握黑白和彩色图像在Matlab中矩阵存储格式。

二、实验内容（边做实验，边将生成结果和图形拷贝到Word文档中）

1．将Test3.jpg和Test4.jpg两个图像文件复制到C:\\Matlab6p5\\work\\目录下。 2．读取Test3.jpg到矩阵A中，观察A的维数，并显示该图像。

3．对A做一定的处理和变换，使之成为下面的图像，并存储到文件PicOut.bmp中。

4．读取Test4.jpg到矩阵X中，去掉红色图层，然后显示出来。

三、设计提示

1．首先计算A的维数，然后列数放大1倍，再考虑将矩阵左半边的图像数据镜像到右半边去。

2．对于彩色图像，矩阵的维数是：行数×列数×图层数（1-红，2-绿，3-蓝）。

四、实验报告要求

1．编写实现第二节实验内容所使用的M文件，并记录相应的生成结果。

2．书写实验报告时要结构合理，层次分明，在分析描述的时候，需要注意语言的流畅。

五、参考答案

《Matlab语言及其在电子信息科学中的应用》实验指导书