matlab 倍频,1/3倍频程及Matlab程序实现

倍频程是声学里人的可听频率范围内，将声音的频谱进行一定规则的集中，变成有限的几个频点对应的强度，这样描述比较起来容易，是一种公约的描述形式。

使用1/3倍频程主要是因为人耳对声音的感觉，其频率分辨能力不是单一频率，而是频带，而1/3倍频程曾经被认为是比较符合人耳特性的频带划分方法，不过现在心理声学里提出了Critical

Band这么个频带划分方法，听说更符合人耳特性，但1/3倍频程仍在广泛使用。

分析频谱时，对于连续谱而言，分析某频率点上的声功率是没有意义的，因此有必要统计某一频带内的声功率。对于频带划分，倍频程和1/3倍频程是常用的划分方法之一，它们都是相对恒定带宽，例如1/3倍频程的带宽是中心频率的23%。

声学及振动测量仪器中的倍频程及1/3倍频程滤波主要是用于对噪声或振动进行频谱分析用的,它们是一种等百分比带宽滤波器,与人耳的频谱分析特性相似。在噪声测量中，使用1/3oct主要是将噪声的频率分布情况更直观的表示出来。便于今后的工作开展。

百分比=(2^(m/2)-2^(-m/2))*100%

其中m就是几倍频程，1/3倍频程m等于1/3。

先要知道1/3倍频程的划分方法，相关的书和国标都有公式和现成的数据表格，然后，将时间域的声信号fft变换到频率域，对定义的每个1/3倍频带的声压计算等效连续声压级。这就是1/3倍频程声压级。

function [g,f] =

oct3spec(B,A,Fs,Fc,s,n);

% OCT3SPEC Plots a

one-third-octave filter characteristics.

%

OCT3SPEC(B,A,Fs,Fc) plots the attenuation of the filter defined

by

%

B and A at sampling frequency Fs. Fc is the center frequency

of

%

the one-third-octave filter. The plot covers one decade on both

sides

%

of Fc.

%

%

OCT3SPEC(B,A,Fs,Fc,'ANSI',N) superposes the ANSI Order-N

analog

%

specification for comparison. Default is N = 3.

%

%

OCT3SPEC(B,A,Fs,Fc,'IEC',N) superposes the characteristics of

the

%

IEC 61260 class N specification for comparison. Default is N =

1.

%

%

[G,F] = OCT3SPEC(B,A,Fs,Fc) returns two 512-point vectors

with

%

the gain (in dB) in G and logarith