数字滤波器及切比雪夫滤波器matlab程序参照.pdfVIP

数字滤波器及切比雪夫滤波器 matlab 程序 数字信号频谱理解 数字信号没有固定的频率， 应该首先规定一个抽样频率， 即信号的最高频率 fm 的二倍， 然后求信号的频率。 %% 数字信号的频谱理解程序示例 clear clc x = ones(1,100); for i = 1:2:99 x(i) = -1; end fs = 100; %规定抽样频率即规定信号的最高频率 fm 的二倍 X=fft(x,512); %信号的傅里叶变换返回 512 个点两个周期即最高频率 fs f=(0:255)/256*(fs/2); %画图只画到 fm 即只画一个周期的频谱图 plot(f,abs(X(1:256))); % 频谱图 xlabel( 频率 (Hz)); ylabel( 幅值 ); title( 信号频率 ); grid; Filter 和 impz 简介 分别用 filter （滤波函数 filter(b,a,x) ）及 impz 函数获取指定离散系统（ b=[0 0.5 0.3 0.2] ）,a=a=[1 0.2 0.4 -0.8] ）的单位抽样响应； 用 freqz 函数获取系统的频率响应； 分别用 root 及 zplane 函数获取系统的零极点图及增益。 程序代码： %E2_2_SignalProcess.m 文件源代码 %Matlab 常用信号分析处理函数演示实例：编写一个 M 文件，分别用 filter 及 impz 函 数获取指定离散系统 % （b=[0.8 0.5 0.6] ，a=[1 0.2 0.4 -0.8] ）的单位抽样响应；用 freqz 函数获取系统的频率 响应；分 %别用 root 及 zplane 函数获取系统的零极点图及增益。 L=128; %单位抽样序列的长度 Fs=1000; %采样频率为 1KHz 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); title(filter() 获取单位抽样响应 ,fontsize,8); subplot(222);stem(ImpzOut); title(impz() 获取单位抽样响应 ,fontsize,8); subplot(223);plot(f,mag); xlabel( 频率 (Hz),fontsize,8); ylabel(