精华内容
下载资源
问答
  • 这是我课程设计的一部分,自己编写的程序,呵呵.
  • 切比雪夫II 型带通IIR 数字滤波器设计1.设计思路(1) 数字—模拟指标转换。利用双线性变换的频率预畸变公式2tan 2ωT =Ω,把所要求的数字滤波器)(z H 数字频率指标转换为相应的模拟滤波器)(s H 的模拟频率指标。(2) ...

    切比雪夫II 型带通IIR 数字滤波器设计

    1.设计思路

    (1) 数字—模拟指标转换。利用双线性变换的频率预畸变公式

    2

    tan 2ωT =Ω,把所要求的数字滤波器)(z H 数字频率指标转换为相应

    的模拟滤波器)(s H 的模拟频率指标。

    (2) 低通模拟指标转换。将模拟滤波器)(s H 的频率指标归一化原型低通

    滤波器)(p H LP 的频率指标。

    (3) 模拟滤波器的设计,得到低通滤波器的归一化传输函数)(p H LP 。 (4) 模拟频率变换,将模拟低通滤波器归一化传输函数)(p H LP 转换成所

    需要的模拟滤波器传输函数)(z H 。

    (5) 模拟—数字滤波器变换。利用双线性变换得到所要求的数字滤波器

    传输函数1

    1112)()(--+-==z z T s s H Z H 。

    2.设计要求及方案

    设计一带通切比雪夫II 型IIR 滤波器,要求如下:

    通带上下边沿频率分别为300 Hz 和400 Hz ,

    通带最大衰减dB p

    3=α,

    阻带上下边沿频率分别为200Hz 和500 Hz ,

    阻带最小衰减 dB s 18=α,采样频率Hz f s 2000=。

    2.1带通滤波器的边缘频率为

    Hz f p 3001=,Hz f p 4002= Hz f s 2001=,Hz f s 5002=

    给定的系统取样频率为Hz f s 2000=,相应的数字频率为

    ππω3.0211=?

    =s

    p p f f ππω4.0222=?

    =s

    p p f f

    展开全文
  • [Matlab]切比雪夫Ⅰ型滤波器设计:低通、高通、带通和带阻-------(2) ​ 切比雪夫滤波器,又名“车比雪夫滤波器”,是在通带或阻带上频率响应幅度等波纹波动的滤波器。切比雪夫滤波器来自切比雪夫分布,以...

    [Matlab]切比雪夫Ⅰ型滤波器设计:低通、高通、带通和带阻-------(2)

    ​ 切比雪夫滤波器,又名“车比雪夫滤波器”,是在通带或阻带上频率响应幅度等波纹波动的滤波器。切比雪夫滤波器来自切比雪夫分布,以“切比雪夫”命名,是用以纪念俄罗斯数学家巴夫尼提·列波维其·切比雪夫。

    切比雪夫Ⅰ型滤波器特点:

    ​ 切比雪夫滤波器在过渡带比巴特沃斯滤波器的衰减快,但频率响应的幅频特性不如后者平坦。切比雪夫滤波器和理想滤波器的频率响应曲线之间的误差最小,但是在通频带内存在幅度波动。

    1、幅度特性是在一个频带内(通带或阻带)范围内具有等波纹特性;
    2、Ⅰ型在通带范围内是等波纹的,在阻带范围内是单调的。

    根据频率响应曲线波动位置的不同,切比雪夫滤波器可以分为以下两种:

    I型切比雪夫滤波器

    通带(或称“通频带”)上频率响应幅度等波纹波动的滤波器称为“I型切比雪夫滤波器”;

    原始信号设定:

    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    %  Cheby1Filter.m
    %  切比雪夫Ⅰ型滤波器的设计
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    clear;
    close all;
    clc;
    
    fs = 1000; %Hz 采样频率
    Ts = 1/fs;
    N  = 1000; %序列长度
    t = (0:N-1)*Ts;
    delta_f = 1*fs/N;
    f1 = 50;
    f2 = 100;
    f3 = 200;
    f4 = 400;
    x1 = 2*0.5*sin(2*pi*f1*t);
    x2 = 2*0.5*sin(2*pi*f2*t);
    x3 = 2*0.5*sin(2*pi*f3*t);
    x4 = 2*0.5*sin(2*pi*f4*t);
    x = x1 + x2 + x3 + x4; %待处理信号由四个分量组成
     
    X = fftshift(abs(fft(x)))/N;
    X_angle = fftshift(angle(fft(x)));
    f = (-N/2:N/2-1)*delta_f;
     
    figure(1);
    subplot(3,1,1);
    plot(t,x);
    title('原信号');
    subplot(3,1,2);
    plot(f,X);
    grid on;
    title('原信号频谱幅度特性');
    subplot(3,1,3);
    plot(f,X_angle);
    title('原信号频谱相位特性');
    grid on;
    

    [外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-ZDDVsYax-1573631016234)(G:\研究生\项目小组任务\笔记\第四周笔记\qubify.bmp)]

    低通滤波器设计:

    %设计一个切比雪夫低通滤波器,要求把50Hz的频率分量保留,其他分量滤掉
    wp = 55/(fs/2);  %通带截止频率,取50~100中间的值,并对其归一化
    ws = 90/(fs/2);  %阻带截止频率,取50~100中间的值,并对其归一化
    alpha_p = 3; %通带允许最大衰减为 db
    alpha_s = 40;%阻带允许最小衰减为 db
    %获取阶数和截止频率
    [ N1 wc1 ] = cheb1ord( wp , ws , alpha_p , alpha_s);
    %获得转移函数系数
    [ b a ] = cheby1(N1,alpha_p,wc1,'low');
    %滤波
    filter_lp_s = filter(b,a,x);
    X_lp_s = fftshift(abs(fft(filter_lp_s)))/N;
    X_lp_s_angle = fftshift(angle(fft(filter_lp_s)));
    figure(2);
    freqz(b,a); %滤波器频谱特性
    figure(3);
    subplot(3,1,1);
    plot(t,filter_lp_s);
    grid on;
    title('低通滤波后时域图形');
    subplot(3,1,2);
    plot(f,X_lp_s);
    title('低通滤波后频域幅度特性');
    subplot(3,1,3);
    plot(f,X_lp_s_angle);
    title('低通滤波后频域相位特性');
    

    [外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-a6R66EbD-1573631016238)(G:\研究生\项目小组任务\笔记\第四周笔记\lowp.bmp)]

    高通滤波器设计:

    %设计一个高通滤波器,要求把400Hz的频率分量保留,其他分量滤掉
    wp = 350/(fs/2);  %通带截止频率,取200~400中间的值,并对其归一化
    ws = 380/(fs/2);  %阻带截止频率,取200~400中间的值,并对其归一化
    alpha_p = 3; %通带允许最大衰减为  db
    alpha_s = 20;%阻带允许最小衰减为  db
    %获取阶数和截止频率
    [ N2 wc2 ] = cheb1ord( wp , ws , alpha_p , alpha_s);
    %获得转移函数系数
    [ b a ] = cheby1(N2,alpha_p,wc2,'high');
    %滤波
    filter_hp_s = filter(b,a,x);
    X_hp_s = fftshift(abs(fft(filter_hp_s)))/N;
    X_hp_s_angle = fftshift(angle(fft(filter_hp_s)));
    figure(4);
    freqz(b,a); %滤波器频谱特性
    figure(5);
    subplot(3,1,1);
    plot(t,filter_hp_s);
    grid on;
    title('高通滤波后时域图形');
    subplot(3,1,2);
    plot(f,X_hp_s);
    title('高通滤波后频域幅度特性');
    subplot(3,1,3);
    plot(f,X_hp_s_angle);
    title('高通滤波后频域相位特性');
    

    [外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-LPt70b8v-1573631016257)(G:\研究生\项目小组任务\笔记\第四周笔记\hightp.bmp)]

    带通滤波器设计:

    %设计一个带通滤波器,要求把50Hz和400Hz的频率分量滤掉,其他分量保留
    wp = [65 385 ] / (fs/2);  %通带截止频率,50~100、200~400中间各取一个值,并对其归一化
    ws = [75 375 ] / (fs/2);  %阻带截止频率,50~100、200~400中间各取一个值,并对其归一化
    alpha_p = 3; %通带允许最大衰减为  db
    alpha_s = 20;%阻带允许最小衰减为  db
    %获取阶数和截止频率
    [ N3 wn ] = cheb1ord( wp , ws , alpha_p , alpha_s);
    %获得转移函数系数
    [ b a ] = cheby1(N3,alpha_p,wn,'bandpass');
    %滤波
    filter_bp_s = filter(b,a,x);
    X_bp_s = fftshift(abs(fft(filter_bp_s)))/N;
    X_bp_s_angle = fftshift(angle(fft(filter_bp_s)));
    figure(6);
    freqz(b,a); %滤波器频谱特性
    figure(7);
    subplot(3,1,1);
    plot(t,filter_bp_s);
    grid on;
    title('带通滤波后时域图形');
    subplot(3,1,2);
    plot(f,X_bp_s);
    title('带通滤波后频域幅度特性');
    subplot(3,1,3);
    plot(f,X_bp_s_angle);
    title('带通滤波后频域相位特性');
    

    [外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-3XBKXRhq-1573631016258)(G:\研究生\项目小组任务\笔记\第四周笔记\bandp.bmp)]

    带阻滤波器设计:

    %设计一个带阻滤波器,要求把50Hz和400Hz的频率分量保留,其他分量滤掉
    wp = [65 385 ] / (fs/2);  %通带截止频率?,50~100、200~400中间各取一个值,并对其归一化
    ws = [75 375 ] / (fs/2);  %阻带截止频率?,50~100、200~400中间各取一个值,并对其归一化
    alpha_p = 3; %通带允许最大衰减为  db
    alpha_s = 20;%阻带允许最小衰减为  db
    %获取阶数和截止频率
    [ N4 wn ] = cheb1ord( wp , ws , alpha_p , alpha_s);
    %获得转移函数系数
    [ b a ] = cheby1(N4,alpha_p,wn,'stop');
    %滤波
    filter_bs_s = filter(b,a,x);
    X_bs_s = fftshift(abs(fft(filter_bs_s)))/N;
    X_bs_s_angle = fftshift(angle(fft(filter_bs_s)));
    figure(8);
    freqz(b,a); %滤波器频谱特性
    figure(9);
    subplot(3,1,1);
    plot(t,filter_bs_s);
    grid on;
    title('带阻滤波后时域图形');
    subplot(3,1,2);
    plot(f,X_bs_s);
    title('带阻滤波后频域幅度特性');
    subplot(3,1,3);
    plot(f,X_bs_s_angle);
    title('带阻滤波后频域相位特性');
    

    [外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-ccbiuMxk-1573631016258)(G:\研究生\项目小组任务\笔记\第四周笔记\bands.bmp)]

    展开全文
  • [Matlab]切比雪夫Ⅱ型滤波器设计:低通、高通、带通和带阻-------(3) II型切比雪夫滤波器 在阻带(或称“阻频带”)上频率响应幅度等波纹波动的滤波器称为“II型切比雪夫滤波器”。 也称倒数切比雪夫滤波器,较不...

    [Matlab]切比雪夫Ⅱ型滤波器设计:低通、高通、带通和带阻-------(3)

    II型切比雪夫滤波器

    阻带(或称“阻频带”)上频率响应幅度等波纹波动的滤波器称为“II型切比雪夫滤波器”。

    也称倒数切比雪夫滤波器,较不常用,因为频率截止速度不如I型快,也需要用更多的电子元件。II型切比雪夫滤波器在通频带内没有幅度波动,只在阻频带内有幅度波动。

    切比雪夫Ⅱ型滤波器特点:

    1、幅度特性是在一个频带内(通带或阻带)范围内具有等波纹特性;

    2、Ⅱ型在通带范围内是单调的,在阻带范围内是等波纹的。

    原始信号设定:

    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    %  Cheby2Filter.m
    %  切比雪夫Ⅱ型滤波器的设计
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     
    clear;
    close all;
    clc;
     
    fs = 1000; %Hz 采样频率
    Ts = 1/fs;
    N  = 1000; %序列长度
    t = (0:N-1)*Ts;
    delta_f = 1*fs/N;
    f1 = 50;
    f2 = 100;
    f3 = 200;
    f4 = 400;
    x1 = 2*0.5*sin(2*pi*f1*t);
    x2 = 2*0.5*sin(2*pi*f2*t);
    x3 = 2*0.5*sin(2*pi*f3*t);
    x4 = 2*0.5*sin(2*pi*f4*t);
    x = x1 + x2 + x3 + x4; %待处理信号由四个分量组成
     
    X = fftshift(abs(fft(x)))/N;
    X_angle = fftshift(angle(fft(x)));
    f = (-N/2:N/2-1)*delta_f;
     
    figure(1);
    subplot(3,1,1);
    plot(t,x);
    title('原信号');
    subplot(3,1,2);
    plot(f,X);
    grid on;
    title('原信号频谱幅度特性');
    subplot(3,1,3);
    plot(f,X_angle);
    title('原信号频谱相位特性');
    grid on;
    

    [外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-lwY7FDzh-1573631797431)(G:\研究生\项目小组任务\笔记\第四周笔记\qubify1.bmp)]

    低通滤波器:

    %设计一个切比雪夫低通滤波器,要求把50Hz的频率分量保留,其他分量滤掉
    wp = 55/(fs/2);  %通带截止频率,取50~100中间的值,并对其归一化
    ws = 60/(fs/2);  %阻带截止频率,取50~100中间的值,并对其归一化
    alpha_p = 3; %通带允许最大衰减为 db
    alpha_s = 40;%阻带允许最小衰减为 db
    %获取阶数和截止频率
    [ N1 wc1 ] = cheb2ord( wp , ws , alpha_p , alpha_s);
    %获得转移函数系数
    [ b a ] = cheby2(N1,alpha_s,wc1,'low');
    %滤波
    filter_lp_s = filter(b,a,x);
    X_lp_s = fftshift(abs(fft(filter_lp_s)))/N;
    X_lp_s_angle = fftshift(angle(fft(filter_lp_s)));
    figure(2);
    freqz(b,a); %滤波器频谱特性
    figure(3);
    subplot(3,1,1);
    plot(t,filter_lp_s);
    grid on;
    title('低通滤波后时域图形');
    subplot(3,1,2);
    plot(f,X_lp_s);
    title('低通滤波后频域幅度特性');
    subplot(3,1,3);
    plot(f,X_lp_s_angle);
    title('低通滤波后频域相位特性');
    

    [外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-Aiv7LRqY-1573631797432)(G:\研究生\项目小组任务\笔记\第四周笔记\lowp1.bmp)]

    高通滤波器:

    %设计一个高通滤波器,要求把400Hz的频率分量保留,其他分量滤掉
    wp = 350/(fs/2);  %通带截止频率,取200~400中间的值,并对其归一化
    ws = 380/(fs/2);  %阻带截止频率,取200~400中间的值,并对其归一化
    alpha_p = 3; %通带允许最大衰减为  db
    alpha_s = 40;%阻带允许最小衰减为  db
    %获取阶数和截止频率
    [ N2 wc2 ] = cheb2ord( wp , ws , alpha_p , alpha_s);
    %获得转移函数系数
    [ b a ] = cheby2(N2,alpha_s,wc2,'high');
    %滤波
    filter_hp_s = filter(b,a,x);
    X_hp_s = fftshift(abs(fft(filter_hp_s)))/N;
    X_hp_s_angle = fftshift(angle(fft(filter_hp_s)));
    figure(4);
    freqz(b,a); %滤波器频谱特性
    figure(5);
    subplot(3,1,1);
    plot(t,filter_hp_s);
    grid on;
    title('高通滤波后时域图形');
    subplot(3,1,2);
    plot(f,X_hp_s);
    title('高通滤波后频域幅度特性');
    subplot(3,1,3);
    plot(f,X_hp_s_angle);
    title('高通滤波后频域相位特性');
    

    [外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-PzagABKE-1573631797433)(G:\研究生\项目小组任务\笔记\第四周笔记\hightp1.bmp)]

    带通滤波器:

    %设计一个带通滤波器,要求把50Hz和400Hz的频率分量滤掉,其他分量保留
    wp = [65 385 ] / (fs/2);  %通带截止频率,50~100、200~400中间各取一个值,并对其归一化
    ws = [75 375 ] / (fs/2);  %阻带截止频率,50~100、200~400中间各取一个值,并对其归一化
    alpha_p = 3; %通带允许最大衰减为  db
    alpha_s = 40;%阻带允许最小衰减为  db
    %获取阶数和截止频率
    [ N3 wn ] = cheb2ord( wp , ws , alpha_p , alpha_s);
    %获得转移函数系数
    [ b a ] = cheby2(N3,alpha_s,wn,'bandpass');
    %滤波
    filter_bp_s = filter(b,a,x);
    X_bp_s = fftshift(abs(fft(filter_bp_s)))/N;
    X_bp_s_angle = fftshift(angle(fft(filter_bp_s)));
    figure(6);
    freqz(b,a); %滤波器频谱特性
    figure(7);
    subplot(3,1,1);
    plot(t,filter_bp_s);
    grid on;
    title('带通滤波后时域图形');
    subplot(3,1,2);
    plot(f,X_bp_s);
    title('带通滤波后频域幅度特性');
    subplot(3,1,3);
    plot(f,X_bp_s_angle);
    title('带通滤波后频域相位特性');
    

    [外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-oNpbpAdi-1573631797434)(G:\研究生\项目小组任务\笔记\第四周笔记\bandp1.bmp)]

    带阻滤波器:

    %设计一个带阻滤波器,要求把50Hz和400Hz的频率分量保留,其他分量滤掉
    wp = [65 385 ] / (fs/2);  %通带截止频率?,50~100、200~400中间各取一个值,并对其归一化
    ws = [75 375 ] / (fs/2);  %阻带截止频率?,50~100、200~400中间各取一个值,并对其归一化
    alpha_p = 3; %通带允许最大衰减为  db
    alpha_s = 40;%阻带允许最小衰减为  db
    %获取阶数和截止频率
    [ N4 wn ] = cheb2ord( wp , ws , alpha_p , alpha_s);
    %获得转移函数系数
    [ b a ] = cheby2(N4,alpha_s,wn,'stop');
    %滤波
    filter_bs_s = filter(b,a,x);
    X_bs_s = fftshift(abs(fft(filter_bs_s)))/N;
    X_bs_s_angle = fftshift(angle(fft(filter_bs_s)));
    figure(8);
    freqz(b,a); %滤波器频谱特性
    figure(9);
    subplot(3,1,1);
    plot(t,filter_bs_s);
    grid on;
    title('带阻滤波后时域图形');
    subplot(3,1,2);
    plot(f,X_bs_s);
    title('带阻滤波后频域幅度特性');
    subplot(3,1,3);
    plot(f,X_bs_s_angle);
    title('带阻滤波后频域相位特性');
    

    [外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-hPjWzOis-1573631797436)(G:\研究生\项目小组任务\笔记\第四周笔记\bands1.bmp)]

    展开全文
  • 关于切比雪夫型IIR滤波器设计matlab代码
  • PAGE PAGE # / 7 基于的带通带阻滤波器设计实例 以下两个滤波器都是切比雪夫 I型数字滤波器不是巴特沃尔滤波器请使用者注意! 1.带通滤波器 (13) 滞通滤波 %吏用注意事项通带或阻带的截止频率与采样率的选取范围是不...
  • 切比雪夫Ⅰ型滤波器特点: 1、幅度特性是在一个频带内(通带或阻带)范围内具有等波纹特性;...% 切比雪夫Ⅰ型滤波器设计 % clear; close all; clc; fs = 1000; %Hz 采样频率 Ts = 1/fs; N = 10...

    切比雪夫Ⅰ型滤波器特点:

       1、幅度特性是在一个频带内(通带或阻带)范围内具有等波纹特性;

       2、Ⅰ型在通带范围内是等波纹的,在阻带范围内是单调的。

    测试代码:

    %  Cheby1Filter.m
    %  切比雪夫Ⅰ型滤波器的设计
    %
    
    clear;
    close all;
    clc;
    
    fs = 1000; %Hz 采样频率
    Ts = 1/fs;
    N  = 1000; %序列长度
    t = (0:N-1)*Ts;
    delta_f = 1*fs/N;
    f1 = 50;
    f2 = 100;
    f3 = 200;
    f4 = 400;
    x1 = 2*0.5*sin(2*pi*f1*t);
    x2 = 2*0.5*sin(2*pi*f2*t);
    x3 = 2*0.5*sin(2*pi*f3*t);
    x4 = 2*0.5*sin(2*pi*f4*t);
    x = x1 + x2 + x3 + x4; %待处理信号由四个分量组成
    
    X = fftshift(abs(fft(x)))/N;
    X_angle = fftshift(angle(fft(x)));
    f = (-N/2:N/2-1)*delta_f;
    
    figure(1);
    subplot(3,1,1);
    plot(t,x);
    title('原信号');
    subplot(3,1,2);
    plot(f,X);
    grid on;
    title('原信号频谱幅度特性');
    subplot(3,1,3);
    plot(f,X_angle);
    title('原信号频谱相位特性');
    grid on;
    
    %设计一个切比雪夫低通滤波器,要求把50Hz的频率分量保留,其他分量滤掉
    wp = 55/(fs/2);  %通带截止频率,取50~100中间的值,并对其归一化
    ws = 90/(fs/2);  %阻带截止频率,取50~100中间的值,并对其归一化
    alpha_p = 3; %通带允许最大衰减为 db
    alpha_s = 40;%阻带允许最小衰减为 db
    %获取阶数和截止频率
    [ N1 wc1 ] = cheb1ord( wp , ws , alpha_p , alpha_s);
    %获得转移函数系数
    [ b a ] = cheby1(N1,alpha_p,wc1,'low'); 
    %滤波
    filter_lp_s = filter(b,a,x);
    X_lp_s = fftshift(abs(fft(filter_lp_s)))/N;
    X_lp_s_angle = fftshift(angle(fft(filter_lp_s)));
    figure(2);
    freqz(b,a); %滤波器频谱特性
    figure(3);
    subplot(3,1,1);
    plot(t,filter_lp_s);
    grid on;
    title('低通滤波后时域图形');
    subplot(3,1,2);
    plot(f,X_lp_s);
    title('低通滤波后频域幅度特性');
    subplot(3,1,3);
    plot(f,X_lp_s_angle);
    title('低通滤波后频域相位特性');
    
    
    
    %设计一个高通滤波器,要求把400Hz的频率分量保留,其他分量滤掉
    wp = 350/(fs/2);  %通带截止频率,取200~400中间的值,并对其归一化
    ws = 380/(fs/2);  %阻带截止频率,取200~400中间的值,并对其归一化
    alpha_p = 3; %通带允许最大衰减为  db
    alpha_s = 20;%阻带允许最小衰减为  db
    %获取阶数和截止频率
    [ N2 wc2 ] = cheb1ord( wp , ws , alpha_p , alpha_s);
    %获得转移函数系数
    [ b a ] = cheby1(N2,alpha_p,wc2,'high');
    %滤波
    filter_hp_s = filter(b,a,x);
    X_hp_s = fftshift(abs(fft(filter_hp_s)))/N;
    X_hp_s_angle = fftshift(angle(fft(filter_hp_s)));
    figure(4);
    freqz(b,a); %滤波器频谱特性
    figure(5);
    subplot(3,1,1);
    plot(t,filter_hp_s);
    grid on;
    title('高通滤波后时域图形');
    subplot(3,1,2);
    plot(f,X_hp_s);
    title('高通滤波后频域幅度特性');
    subplot(3,1,3);
    plot(f,X_hp_s_angle);
    title('高通滤波后频域相位特性');
    
    
    %设计一个带通滤波器,要求把50Hz和400Hz的频率分量滤掉,其他分量保留
    wp = [65 385 ] / (fs/2);  %通带截止频率,50~100、200~400中间各取一个值,并对其归一化
    ws = [75 375 ] / (fs/2);  %阻带截止频率,50~100、200~400中间各取一个值,并对其归一化
    alpha_p = 3; %通带允许最大衰减为  db
    alpha_s = 20;%阻带允许最小衰减为  db
    %获取阶数和截止频率
    [ N3 wn ] = cheb1ord( wp , ws , alpha_p , alpha_s);
    %获得转移函数系数
    [ b a ] = cheby1(N3,alpha_p,wn,'bandpass'); 
    %滤波
    filter_bp_s = filter(b,a,x);
    X_bp_s = fftshift(abs(fft(filter_bp_s)))/N;
    X_bp_s_angle = fftshift(angle(fft(filter_bp_s)));
    figure(6);
    freqz(b,a); %滤波器频谱特性
    figure(7);
    subplot(3,1,1);
    plot(t,filter_bp_s);
    grid on;
    title('带通滤波后时域图形');
    subplot(3,1,2);
    plot(f,X_bp_s);
    title('带通滤波后频域幅度特性');
    subplot(3,1,3);
    plot(f,X_bp_s_angle);
    title('带通滤波后频域相位特性');
    
    
    %设计一个带阻滤波器,要求把50Hz和400Hz的频率分量保留,其他分量滤掉
    wp = [65 385 ] / (fs/2);  %通带截止频率?,50~100、200~400中间各取一个值,并对其归一化
    ws = [75 375 ] / (fs/2);  %阻带截止频率?,50~100、200~400中间各取一个值,并对其归一化
    alpha_p = 3; %通带允许最大衰减为  db
    alpha_s = 20;%阻带允许最小衰减为  db
    %获取阶数和截止频率
    [ N4 wn ] = cheb1ord( wp , ws , alpha_p , alpha_s);
    %获得转移函数系数
    [ b a ] = cheby1(N4,alpha_p,wn,'stop'); 
    %滤波
    filter_bs_s = filter(b,a,x);
    X_bs_s = fftshift(abs(fft(filter_bs_s)))/N;
    X_bs_s_angle = fftshift(angle(fft(filter_bs_s)));
    figure(8);
    freqz(b,a); %滤波器频谱特性
    figure(9);
    subplot(3,1,1);
    plot(t,filter_bs_s);
    grid on;
    title('带阻滤波后时域图形');
    subplot(3,1,2);
    plot(f,X_bs_s);
    title('带阻滤波后频域幅度特性');
    subplot(3,1,3);
    plot(f,X_bs_s_angle);
    title('带阻滤波后频域相位特性');
    

      

    效果:

    原始信号:

    生成的低通滤波器和滤波后的效果:

    生成的高通滤波器和滤波后的结果:

    生成的带通滤波器和滤波后的结果:

    生成的带阻滤波器和滤波后的结果:

    转载于:https://www.cnblogs.com/alimy/p/9140695.html

    展开全文
  • 切比雪夫Ⅱ型滤波器特点: 1、幅度特性是在一个频带内(通带或阻带)范围内具有等波纹特性;...% 切比雪夫Ⅱ型滤波器设计 % clear; close all; clc; fs = 1000; %Hz 采样频率 Ts = 1/fs; N = 10...
  • 滤波器的种类很多,从功能上可将滤波器分为低、带、高、带阻类型。从实现方法上可分为FIR、IIR类型 。从设计方法上可分为Chebyshev(切比雪夫),Butterworth(巴特沃兹)。而本次课程设计上要用到的切比雪夫滤波器...
  • 要求 带阻滤波器,中心频率为10kHz,在±250Hz(9.75kHz,1O.25kHz)衰减为3dB,在±100Hz(9.9kHz、1O.1kHz)最小衰减为30dB,信号源和负载阻抗为600Ω。  解 ①将所给条件变换为几何对称的设计指标。因为带宽...
  • 基于matlab的带通、带阻滤波器设计实例以下两个滤波器都是切比雪夫I型数字滤波器,不是巴特沃尔滤波器,请使用者注意!1.带通滤波器function y=bandp(x,f1,f3,fsl,fsh,rp,rs,Fs)%带通滤波%使用注意事项:通带或阻带...
  • matlab设计模拟带阻滤波器

    万次阅读 多人点赞 2018-10-24 22:13:11
    简单记录下在matlab上如何设计出模拟的带阻滤波器,包括:巴特沃斯滤波器、切比雪夫I型滤波器、切比雪夫II型滤波器、椭圆型滤波器。 %设计带阻滤波器 %巴特沃斯、切比雪夫I型、切比雪夫II型、椭圆型滤波器 clear ...
  • 下两个滤波器都是切比雪夫I型数字滤波器,不是巴特沃尔滤波器,请使用者注意!1.低通滤波器使用说明:将下列代码幅值然后以m文件保存,文件名要与函数名相同,这里函数名:lowp。function y=lowp(x,f1,f3,rp,rs,Fs) ...
  • Matlab——切比雪夫I型滤波器 cheb1ord 、 chey1 and freqz 切比雪夫Ⅰ型滤波器在通带或阻带上频率响应幅度等波纹波动的滤波器切比雪夫滤波器和理想滤波器的频率响应曲线之间的误差最小,但是在通频带内存在幅度...
  • 本资源包含一个matlab代码文件和两个音乐信号,通过切比雪夫一型滤波器滤除高频噪声,同时用均值滤波滤除白噪声
  • 有频率为5Hz、15Hz、30Hz的叠加余弦信号,设计低通滤波器保留5Hz的频率分量,设计高通滤波器保留30Hz的频率分量,设计带通滤波器保留20Hz的频率分量,设计带阻滤波器滤除5Hz和30Hz的频率分量。 一、pandas是什么? ...
  •  巴特沃思数字带阻滤波器设计: 抽样频率为10kHZ,,在-2 dB衰减处的边带频率是1.5 kHZ,4 kHZ, 在-13dB衰减处频率是2kHZ和3kHZ  分别绘制这三种数字滤波器的幅度响应曲线和相位响应曲线;  采用切比雪夫Ⅰ型...
  • 巴特沃思数字带阻滤波器设计: 抽样频率为10kHZ,,在-2 dB衰减处的边带频率是1.5 kHZ,4 kHZ, 在-13dB衰减处频率是2kHZ和3kHZ 分别绘制这三种数字滤波器的幅度响应曲线和相位响应曲线; 采用切比雪夫Ⅰ型滤波器为...
  • 利用模拟低通滤波器设计带通、带阻、高通模拟滤波器 参考:《数字信号处理 第四版》高西全、丁玉美 最近学到第六章了,利用模拟低通滤波器设计不同类型的模拟滤波器,我们重点在设计巴特沃斯滤波器。 巴特沃斯滤波器...
  • 测试代码: ...%% 使用双线性变换法设计带阻滤波器 %% 2018年6月13日 16:48:04 %% author:Alimy close all; clear; clc; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %代码正文 %%%%%%%%%%%%%%%%%%...

空空如也

空空如也

1 2 3 4 5 ... 20
收藏数 404
精华内容 161
关键字:

切比雪夫带阻滤波器设计