三分之一倍频程matlab

三分之一倍频程计算以及相关滤波器的matlab程序

function [BA] = adsgn（Fs）; 
% ADSGN  Design of a A-weighting filter.
%    [BA] = ADSGN（Fs） designs a digital A-weighting filter for 
%    sampling frequency Fs. Usage: Y = FILTER（BAX）. 
%    Warning: Fs should normally be higher than 20 kHz. For example 
%    Fs = 48000 yields a class 1-compliant filter.
%
%    Requires the Signal Processing Toolbox. 
%
%    See also ASPEC CDSGN CSPEC. 

% Author: Christophe Couvreur Faculte Polytechnique de Mons （Belgium）
%         couvreur@thor.fpms.ac.be
% Last modification: Aug. 20 1997 10:00am.

% References: 
%    [1] IEC/CD 1672: Electroacoustics-Sound Level Meters Nov. 1996. 


% Definition of analog A-weighting filter according to IEC/CD 1672.
f1 = 20.598997; 
f2 = 107.65265;
f3 = 737.86223;
f4 = 12194.217;
A1000 = 1.9997;
pi = 3.14159265358979;
NUMs = [ （2*pi*f4）^2*（10^（A1000/20）） 0 0 0 0 ];
DENs = conv（[1 +4*pi*f4 （2*pi*f4）^2][1 +4*pi*f1 （2*pi*f1）^2]）; 
DENs = conv（conv（DENs[1 2*pi*f3]）[1 2*pi*f2]）;

% Use the bilinear transformation to get the digital filter. 
[BA] = bilinear（NUMsDENsFs）; 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1340  2009-05-24 19:20  license.txt
     目录           0  1997-12-29 16:57  octave\
     文件        3474  1997-08-27 11:43  octave\README
     文件        1069  1997-08-26 18:01  octave\adsgn.m
     文件        1938  1997-08-26 18:01  octave\aspec.m
     文件        1579  1997-08-27 07:44  octave\bankdisp.m
     文件         984  1997-08-26 18:01  octave\cdsgn.m
     文件        1915  1997-08-26 18:01  octave\cspec.m
     文件        4420  1997-08-27 11:38  octave\filtbank.m
     文件         860  1997-08-26 18:01  octave\leq.m
     文件        2786  1997-08-26 18:01  octave\oct3bank.m
     文件        1466  1997-08-26 18:01  octave\oct3dsgn.m
     文件        3456  1997-09-04 08:47  octave\oct3spec.m
     文件        1458  1997-08-26 18:01  octave\octdsgn.m
     文件        3275  1997-09-04 08:47  octave\octspec.m