matlab与数字信号处理

matlab功能强大，能够处理图像和数字信号，matlab中提供了很多处理数字信号的函数，比如filter(),impz()等

L=128;					        %单位抽样序列的长度
Fs=1000;					    %采样频率为1 kHz
b=[0.8 0.5 0.6];				%系统函数的分子系数向量
a=[1 0.2 0.4 -0.8];			    %系统函数的分母系数向量
delta=[1 zeros(1,L-1)];			%生成长度为L的单位抽样序列

FilterOut=filter(b,a,delta);	%filter函数获取单位抽样响应
ImpzOut=impz(b,a,L);			%impz函数获取单位抽样响应 
[h,f]=freqz(b,a,L,Fs);			%freqz函数求频率响应
mag=20*log(abs(h))/log(10);  	%幅度转换成dB单位    
ph=angle(h)*180/pi;			    %相位值单位转换
zr=roots(b)					    %求系统的零点，并显示在命令窗口
pk=roots(a)				        %求系统的极点，并显示在命令窗口
g=b(1)/a(1)					    %求系统的增益，并显示在命令窗口

%绘图
figure(1);
subplot(221);stem(FilterOut);
subplot(222);stem(ImpzOut);
subplot(223);plot(f,mag);
subplot(224);plot(f,ph);
figure(2);
freqz(b,a);					%用feqz函数绘制系统频率响应
figure(3);  
zplane(b,a);				%用zplane函数绘制系统零极点图