matlab大作业Matlab在数字信号处理及图像处理中的应用

<>

(1)产生简单信号 程序如下: close all clear all

t1=-5:0.01:5;%定义t的取值范围 y1=sinc(t1);%调用sinc函数

subplot(221);plot(t1,y1);grid; title('sinc函数');

y2=rectpuls(t1,2);%调用门函数 subplot(222);plot(t1,y2);grid; title('门函数');

axis（[-5,5,-0.5,1.5]）;%定义坐标轴的范围 y3=tripuls(t1,3);%调用三角脉冲函数 subplot(223);plot(t1,y3);grid; title（'三角脉形冲'）; axis（[-5,5,-0.5,1.5]）;

（2）求模型的零输入响应

????x1?=?????0.7814x2????0.5572?0.781401?0?u1?x1??????+???? ???x2??02??u2?x10????Y=[1.9691 6.4493]??,初始条件为x(0)=??

?1??x2?

程序如下

a=[-0.5572 -0.7814;0.7814 0]; c=[1.969 6.4493]; x0=[1 ; 0]%初始状态?? x0 =

1 0

sys=ss(a,[],c,[]);%状态方程?? initial(sys,x0)

(3）滤波器设计

1.设计一个butterworth数字低通滤波器，抽样频率为1000HZ，滤波器的3db截至频率为40HZ，阻带截止频率为１５０HZ，阻带最小衰减为６０ｄｂ．绘制该滤波器的频率响应曲线。

程序如下：

Wp=40/500;%通带截止频率

Ws=150/500;%阻带截止频率

[n,Wn]=buttord(Wp,Ws,3,60)%求butterworth滤波器的阶数和3db带宽的截止频率 [b,a]=butter(n,Wn);%调用butter函数设计模拟滤波器 freqz(b,a,512,1000);%画频率响应图 title('巴特沃斯滤波器')

２.设计一个ｃｈｅｂｙｓｈｅ Ｉ数字低通滤波器，抽样频率为1000HZ，滤波器的3db截至频率为40HZ，阻带截止频率为１５０HZ，阻带最小衰减为６０ｄｂ．绘制该滤波器的频率响应曲线。 程序如下： Wp=40/500; Ws=150/500;

Rp=3;Rs=60;%通带和阻带的波纹系数 [n,Wn]=cheb1ord(Wp,Ws,Rp,Rs) [b,a]=cheby1(n,Rp,Wn); freqz(b,a,512,1000); title('切比雪夫滤波器')

(4)图像处理

1.实现两图像相减 程序如下：

I=imread('pepper.jpg');%用imread函数读入图片 K=imread('fruit.jpg');

Z=imsubtract(I,K);%用imsubtract函数实现两图片相减的功能 imshow(Z)%显示相减后的图片

2.实现辆图像相加 程序如下：

I=imread('pepper.jpg');%用imread函数读入图片 J=imread('fruit.jpg');

K=imadd(I,J,'uint8');%用imadd函数实现两图片相加的功能

imshow(K)

title('两图像相加')

3.图像切割 程序如下：

I=imread('pepper.jpg');

i=imcrop(I,[40,30,50,60]);%用imcrop函数实现切割图片的功能

subplot(121); imshow(I);

title('辣椒原图'); subplot(122); imshow(i);

title('切割后的图');

4.图像旋转 程序如下：

I=imread('pepper.jpg');

i=imrotate(I,-45,'bilinear');%用imrotate函数实现图片旋转的功能 subplot(121); imshow(I);

title('辣椒原图'); subplot(122); imshow(i);

title('旋转后的图');