(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
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.实现两图像相减
程序如下：
用imread函数读入图片
Z=imsubtract(I,K;%用imsubtract函数实现两图片相减的功能
imshow(Z%显示相减后的图片
2.实现辆图像相加
程序如下：
用imread函数读入图片
K=imadd(I,J,'uint8';%用imadd函数实现两图片相加的功能
imshow(K
title('两图像相加'
3.图像切割
程序如下：
i=imcrop(I,[40,30,50,60];%用imcrop函数实现切割图片的功能
subplot(121;
imshow(I;
title('辣椒原图';
subplot(122;
imshow(i;
title('切割后的图';
4.图像旋转
程序如下：
i=imrotate(I,-45,'bilinear';%用imrotate函数实现图片旋转的功能
subplot(121;
imshow(I;
title('辣椒原图';
subplot(122;
imshow(i;
title('旋转后的图';
实验小结：
通过做这次大作业以后才发现自己对matlab软件还不是很熟悉，看到别的同学做的东西都是很复杂而且功能很强大，在自己做的大作业却是相当的简单，有时想一下真的好惭愧。
通过这次大作业，自己又重新熟悉了matlab的一些功能，刚开始的时候本来想用systemview来仿真一个通信系统的，但那软件不太会用，最后就选择了直接用matlab来处理数字信号处理中的一些问题，做完以后又发现做的内容太少了，于是又在原来的基础上又加了几个图像处理的简单例子，所以看起来内容就不会那么单调了。感觉到学习某一软件就要勤加练习，这样才会有效果。下载本文
