注意:此函数现在可以从IoSR Matlab工具箱中作为IoSR.auditory.iso226使用。
-------------------------
ISO 226:2003标准等亮度等高线
[spl,f]=iso226(电话)
[spl,f]=iso226(电话,fq)
…=iso226(电话,fq,平方)
…=iso226(电话,[],平方)
[SPL,F]=ISO226(PHON)返回声压级
纯音频率F(Hz)的(SPL)(dB)
响度等级电话。计算值
根据ISO 226:2003,使用参考
标准中规定的频率。根据
标准,PHON仅在以下情况下在所有频率有效
20<=PHON<80(尽管函数将返回SPL
值超出此范围)。
PHON可以是任何大小的阵列;SPL和F将为
大小[1,29,M,N,P,…],其中M,N,P,。。。是尺寸
电话号码。
[SPL,F]=ISO226(PHON,FQ)返回纯
指定响度下FQ中的音调频率
级别。根据标准,FQ仅有效
低于12.5 kHz(尽管该功能将外推
SPL值超出此范围)。
FQ可以是任何大小的数组;SPL和F将为
尺寸[Q,R,S,…,M,N,P,…]其中Q,R、S,。。。是
FQ的尺寸。
…=iso226(phon,fq,sq)指定单例
尺寸将从输出中删除。使用
sq=false,将保留单个维度
(默认),否则将删除它们。
…=iso226(phon,[],sq)使用标准参考
SPL计算的频率。
例子
%绘制20到80个发声单元之间的等丰度轮廓
%计算SPL
声音=20:10:80;
[spl,f]=iso226(phons,[],true);
%绘图
数字;半对数(f,spl)
集合(gca,'xlim',[最小(f(:)))
图例(num2str(phons),‘位置’,‘西南’);
title(“不同响度水平的等响度轮廓(以音箱为单位)”)
xlabel(“频率[Hz]”)
ylabel(“声压级[dB]”)