zoukankan      html  css  js  c++  java
  • 《DSP using MATLAB》示例Example7.18

    代码:

    M = 33; alpha = (M-1)/2; l = 0:M-1; wl = (2*pi/M)*l; 
    T1 = 0.1095; T2 = 0.598;
    Hrs = [zeros(1,11),T1,T2,ones(1,8),T2,T1,zeros(1,10)];    % Ideal Amp Res sampled
    Hdr = [0, 0, 1, 1]; wdl = [0, 0.6, 0.8, 1];               % Ideal Amp Res for plotting
    k1 = 0:floor((M-1)/2); k2 = floor((M-1)/2)+1:M-1;
    
    angH = [-alpha*(2*pi)/M*k1, alpha*(2*pi)/M*(M-k2)]; 
    H = Hrs.*exp(j*angH); h = real(ifft(H, M));
    
    [db, mag, pha, grd, w] = freqz_m(h, 1);
    [Hr, ww, a, L] = Hr_Type1(h);
    
    %Plot 
    
    figure('NumberTitle', 'off', 'Name', 'Exameple 7.18')
    set(gcf,'Color','white'); 
    
    subplot(2,2,1); plot(wl(1:17)/pi, Hrs(1:17), 'o', wdl, Hdr); axis([0, 1, -0.1, 1.1]); grid on;
    xlabel('frequency in pi nuits'); ylabel('Hr(k)'); title('Hignpass: M=33, T1=0.1095, T2=0.598');
    
    subplot(2,2,2); stem(l, h); axis([-1, M, -0.4, 0.4]); grid on;
    xlabel('n'); ylabel('h(n)'); title('Impulse Response');
    
    subplot(2,2,3); plot(ww/pi, Hr, wl(1:17)/pi, Hrs(1:17), 'o'); axis([0, 1, -0.1, 1.1]); grid on;
    xlabel('frequency in pi units'); ylabel('Hr(w)'); title('Amplitude Response');
    
    subplot(2,2,4); plot(w/pi, db); axis([0, 1, -100, 10]); grid on;
    xlabel('frequency in pi units'); ylabel('Decibels'); title('Magnitude Response');
    

      运行结果:

          幅度响应图看出,满足阻带衰减50dB的设计要求。

    牢记: 1、如果你决定做某事,那就动手去做;不要受任何人、任何事的干扰。2、这个世界并不完美,但依然值得我们去为之奋斗。
  • 相关阅读:
    小程序 视频
    b161: NOIP2007 4.Hanoi双塔问题
    命名规则、.gitignore、freopen()
    c++学习记录(九)
    c++学习笔记(八)
    2020面向对象程序设计寒假作业2
    c++学习记录(七)
    c++学习记录(六)
    c+学习记录(五)
    c++学习记录(四)
  • 原文地址:https://www.cnblogs.com/ky027wh-sx/p/6646999.html
Copyright © 2011-2022 走看看