注意：此函数现在可以从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]”）
