%% FIR 陷波器设计
%% ======等波纹最佳逼近法 Equiripple Best Approximation
Method=========%%
% 目的：设计一个陷波器阻带在 50+-1.5Hz 以内,采样频率为 400Hz 的
滤波器，
% 并要求通带最大衰减为 0.1dB，阻带最小衰减为 80dB。
clc;
clear;close all;
rp=0.1; % 通带波纹最大衰减为 0.1dB
rs=80; % 阻带衰减为 80dB
fs=400; % 采样频率为 400Hz
f=[48.5,49.9,50.1,51.5]; % 根据阻带：49.9Hz~50.1Hz，通带：
0~48.5Hz 和 51.5Hz~pi，w=mo_w*T=mo_w/fs=2*pi*f/fs
m=[1 0 1]; % 期望幅度
dat1=(10^(rp/20)-1)/(10^(rp/20)+1); % 根据（7.4.4）和
（7.4.5）式计算通带和阻带波纹幅度 dat1，dat2 和 rip
dat2=10^(-rs/20);
rip=[dat1 dat2 dat1];
[M,fo,mo,wc]=remezord(f,m,rip,fs);
hn=remez(M,fo,mo,wc); % 得到滤波器单位脉冲响应
[hw,w]=freqz(hn);
%% ==========================产生信号====
======================== %%
Nlen=2038400;
n=0:Nlen-1;
Ts=1/fs;
xn=sin(2*pi.*48*n*Ts)+sin(2*pi.*48.5*n*Ts)+sin(2*pi.*49*n*Ts)
+sin(2*pi*50.*n*Ts)+sin(2*pi.*51*n*Ts)+sin(2*pi.*51.5*n*Ts)
祚一叶刀客
- 粉丝: 7
- 资源: 6
我的内容管理
收起
- 我的资源
快来上传第一个资源
我的收益 登录查看自己的收益
我的积分
登录查看自己的积分
我的C币
登录后查看C币余额
我的收藏 
我的下载 
下载帮助
评论0