求MATLAB freqz_m函数的源码

dongchao1207

百度Google都搜不到，好不容易找到个还要付费，拜托各位有的发个源码:lol 谢谢:loveliness:

ChaChing

搜寻一下matlab路径, 不就有freqz.m !?

beyondhxf

1. 
   
2. which freqz.m
   

复制代码

friendchj

可能楼主需要的 freqz_m.m，这个其实是能搜到的。

1. 
   
2. function [db,mag,pha,grd,w] = freqz_m(b,a);
   
3. % Modified version of freqz subroutine
   
4. % ------------------------------------
   
5. % [db,mag,pha,grd,w] = freqz_m(b,a);
   
6. %  db = Relative magnitude in dB computed over 0 to pi radians
   
7. % mag = absolute magnitude computed over 0 to pi radians
   
8. % pha = Phase response in radians over 0 to pi radians
   
9. % grd = Group delay over 0 to pi radians
   
10. %   w = 501 frequency samples between 0 to pi radians
    
11. %   b = numerator polynomial of H(z)   (for FIR: b=h)
    
12. %   a = denominator polynomial of H(z) (for FIR: a=[1])
    
13. %
    
14. [H,w] = freqz(b,a,1000,'whole');
    
15.     H = (H(1:1:501))'; w = (w(1:1:501))';
    
16.   mag = abs(H);
    
17.    db = 20*log10((mag+eps)/max(mag));
    
18.   pha = angle(H);
    
19. %  pha = unwrap(angle(H));
    
20.   grd = grpdelay(b,a,w);
    
21. %  grd = diff(pha);
    
22. %  grd = [grd(1) grd];
    
23. %  grd = [0 grd(1:1:500); grd; grd(2:1:501) 0];
    
24. %  grd = median(grd)*500/pi;
    

复制代码

dongchao1207

就是这个，谢谢:loveliness: 我搜到的都要money才让下:@L