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）求模型的零输入响应

=+

Y=[1.9691 6.4493],初始条件为x(0=程序如下

a=[-0.5572 -0.7814;0.7814 0];

c=[1.969 6.4493];

x0=[1 ; 0]%初始状态

x0 =

1

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('旋转后的图';

实验小结：

通过做这次大作业以后才发现自己对matlab软件还不是很熟悉，看到别的同学做的东西都是很复杂而且功能很强大，在自己做的大作业却是相当的简单，有时想一下真的好惭愧。

通过这次大作业，自己又重新熟悉了matlab的一些功能，刚开始的时候本来想用systemview来仿真一个通信系统的，但那软件不太会用，最后就选择了直接用matlab来处理数字信号处理中的一些问题，做完以后又发现做的内容太少了，于是又在原来的基础上又加了几个图像处理的简单例子，所以看起来内容就不会那么单调了。感觉到学习某一软件就要勤加练习，这样才会有效果。