企业官网建设流程全解析

厦门上网站设计建设,图片网站的优化,企业文化的重要性和意义,门户网站建设发展趋势目录 1.butter函数,巴特沃斯滤波器 2. cheby1函数#xff0c;切比雪夫I型滤波器 3. cheby2函数#xff0c;切比雪夫II型滤波器 4.ellip函数#xff0c;椭圆滤波器 5.yulewalk函数#xff0c;Yule-Walker滤波器 6.各函数对比分析 IIR滤波器即无限长单位冲激响应滤波器…目录1.butter函数,巴特沃斯滤波器2. cheby1函数切比雪夫I型滤波器3. cheby2函数切比雪夫II型滤波器4.ellip函数椭圆滤波器5.yulewalk函数Yule-Walker滤波器6.各函数对比分析IIR滤波器即无限长单位冲激响应滤波器其输出不仅依赖当前输入还依赖历史输出核心特征是存在反馈回路冲激响应理论上无限长。IIR滤波器的差分方程为其中1.butter函数,巴特沃斯滤波器巴特沃斯滤波器的核心是最大平坦幅频特性通带内无纹波阻带单调衰减。在MATLAB中butter函数的程序如下Fs 1200; % 采样频率 Fc 100; % 截止频率 [n, Wn] buttord(Fc/(Fs/2), 120/(Fs/2), 1, 40); % 确定阶数和归一化截止频率 [b, a] butter(n, Wn); % 设计滤波器b分子系数a分母系数 [h, f] freqz(b, a, 1024, Fs); % 计算频率响应 % 绘图 figure; subplot(2,1,1); plot(f, 20*log10(abs(h))); title(巴特沃斯幅频响应); xlabel(频率(Hz)); ylabel(幅度(dB)); grid on; subplot(2,1,2); plot(f, unwrap(angle(h))); title(巴特沃斯相频响应); xlabel(频率(Hz)); ylabel(相位(rad)); grid on;其幅频响应如下图所示2. cheby1函数切比雪夫I型滤波器切比雪夫I型通带等纹波阻带单调衰减以通带纹波换取更陡的过渡带。在MATLAB中cheby1函数的程序如下% 设计切比雪夫I型低通滤波器通带纹波1dB阻带衰减40dB Fs 1200; Fc 100; Fs_stop 120; Rp 1; Rs 40; [n, Wn] cheb1ord(Fc/(Fs/2), Fs_stop/(Fs/2), Rp, Rs); [b, a] cheby1(n, Rp, Wn); % 核心函数阶数、通带纹波、归一化截止频率 [h, f] freqz(b, a, 1024, Fs); % 绘图 figure; subplot(2,1,1); plot(f, 20*log10(abs(h))); title(切比雪夫I型幅频响应); xlabel(频率(Hz)); ylabel(幅度(dB)); grid on; subplot(2,1,2); plot(f, unwrap(angle(h))); title(切比雪夫I型相频响应); xlabel(频率(Hz)); ylabel(相位(rad)); grid on;其幅频响应如下图所示3. cheby2函数切比雪夫II型滤波器切比雪夫II型通带单调阻带等纹波纹波出现在阻带通带保持平坦。其中ωs​为阻带截止频率其余参数同切比雪夫I型。在MATLAB中cheby2函数的程序如下% 设计切比雪夫II型低通滤波器阻带纹波40dB通带平坦 Fs 1200; Fc 100; Fs_stop 120; Rp 1; Rs 40; [n, Wn] cheb2ord(Fc/(Fs/2), Fs_stop/(Fs/2), Rp, Rs); [b, a] cheby2(n, Rs, Wn); % 核心函数阶数、阻带衰减、归一化截止频率 [h, f] freqz(b, a, 1024, Fs); % 绘图 figure; subplot(2,1,1); plot(f, 20*log10(abs(h))); title(切比雪夫II型幅频响应); xlabel(频率(Hz)); ylabel(幅度(dB)); grid on; subplot(2,1,2); plot(f, unwrap(angle(h))); title(切比雪夫II型相频响应); xlabel(频率(Hz)); ylabel(相位(rad)); grid on;其幅频响应如下图所示4.ellip函数椭圆滤波器椭圆滤波器考尔滤波器通带和阻带均等纹波过渡带最陡阶数最低但相位非线性最强。其中Rn(x)为n阶椭圆有理函数由雅可比椭圆函数推导核心特征是通带/阻带均有纹波实现最小阶数下的最陡衰减。在MATLAB中ellip函数的程序如下Fs 1200; Fc 100; Fs_stop 120; Rp 1; Rs 40; [n, Wn] ellipord(Fc/(Fs/2), Fs_stop/(Fs/2), Rp, Rs); [b, a] ellip(n, Rp, Rs, Wn); % 核心函数阶数、通带纹波、阻带衰减、归一化截止频率 [h, f] freqz(b, a, 1024, Fs); % 绘图 figure; subplot(2,1,1); plot(f, 20*log10(abs(h))); title(椭圆滤波器幅频响应); xlabel(频率(Hz)); ylabel(幅度(dB)); grid on; subplot(2,1,2); plot(f, unwrap(angle(h))); title(椭圆滤波器相频响应); xlabel(频率(Hz)); ylabel(相位(rad)); grid on;其幅频响应如下图所示5.yulewalk函数Yule-Walker滤波器Yule-Walker滤波器是基于频域采样的IIR滤波器设计方法通过最小二乘法拟合期望的幅频响应无显式闭合公式核心步骤给定频点和对应幅度构造期望频率响应计算自相关函数解Yule-Walker方程Rar(R为自相关矩阵a为分母系数r为自相关向量)得到滤波器系数。在MATLAB中yulewalk函数的程序如下% 设计Yule-Walker滤波器拟合自定义幅频响应 Fs 1200; f [0, 80, 100, 250,600]/(Fs/2); % 归一化频点0~1 m [1, 1, 0.5, 0, 0]; % 对应频点的幅度 n 8; % 滤波器阶数 [b, a] yulewalk(n, f, m); % 核心函数阶数、频点、幅度 [h, freq] freqz(b, a, 1024, Fs); % 绘图 figure; plot(f*(Fs/2), m, ro-, freq, abs(h), b-); title(Yule-Walker滤波器期望vs实际幅频响应); xlabel(频率(Hz)); ylabel(幅度); legend(期望,实际); grid on;其幅频响应如下图所示6.各函数对比分析函数幅频特性过渡带陡峭度相位非线性阶数相同指标适用场景butter通带无纹波阻带单调平缓较小最高要求通带平坦、相位线性如音频cheby1通带等纹波阻带单调较陡中等中等允许通带纹波需较陡过渡带cheby2通带单调阻带等纹波较陡中等中等要求通带平坦允许阻带纹波ellip通带 阻带均等纹波最陡最大最低阶数受限、过渡带要求极高如通信yulewalk拟合自定义幅频响应灵活不确定按需设定非标准滤波器自定义频响