Кривые контура равной громкости из ISO 226:2003. У кого есть первоисточник?

[Nikolay_Po]

Добрый день (или др. время суток)!

Сабж.

Во всём Интренет - лишь этот график: http://en.wikipedia.org/wiki/File:Lindos1.svg.
Он векторный и мог бы быть точным, но человек его просто нарисовал в OpenOffice Draw на основании данных стандарта.

Вопрос: В каком виде эти контуры определены в стандарте? Полагаю, должны быть определены математически. Кто видел документ ISO 226:2003 целиком?

Мне данные необходимы для экспериментов. Я уже сделал аппроксимацию кривой для громкости 80фон нормальными функциями (Гаусса). С неплохой точностью, при этом хватило 7 функций для среднеквадратичной ошибки в 0.13дБ и пиковой не более 0.25дБ относительно графика по пикселям.
Но сама форма кривых на графике из Википедии предполагает невысокую точность - они тоже были введены вручную, даже просто нарисованы...

[Astronavt]

точно настроенный параметрик

Offтопик:
Сейчас нам ИГВИН задаст параметрик...

[Nikolay_Po]

Я использую как минимум для трех вещей - АЧХ на слух, моды комнаты, и одинаковость АЧХ для двух АС (КИЗ посередине).

Вот про проверку чёткости и положения КИЗ я и не подумал. Спасибо!
Всё верно. Для этих целей и хочу более качественный сигнал создать. И не придумал ничего лучшего, чем третьоктавные полосы, тем более, что они на эквалайзере есть.

[gross]

IMHO, подобрать на слух небольшую тонкомпенсацию на НЧ и, если динамический диапазон сигнала достаточно большой, обработать _допороговым_ децилинейным компрессором, коэффициент сжатия около 1,3...1,6. В цепи управления компрессора должен стоять фильтр с АЧХ по типу МЭК-А или МЭК-В, чтобы он работал ближе к субъективной громкости. Еще максимально задавить комнатные моды, чтобы их гул не забивал все. В первую очередь - резонансные поглотители и оптимизация расстановки АС, во вторую (для наведения окончательного лоска) - точно настроенный параметрик. Работа реализуемая, но очень муторная.

Да уж, работы много - не то слово. С ТК все понятно, она у меня работает уже много лет, вопрос только в ее величине - под каждый диск она нужна своя, либо иметь одну, но мириться с неидеальностью. Компрессор нужен. Я слушаю в основном записи прежних лет, когда писали еще с приличным ДД, и на небольших громкостях тихие уровни слышны очень плохо. С компрессором вопрос - чем регулировать, чтобы не искажать при этом? Комнатные моды - это самая мутная тема, тут вариантов много. Единственный, неприемлемый для меня - превращать комнату в звуковую студию. Это правильно с т.з. акустики, но с т.з. жизни в комнате - сомнительно.

[Nikolay_Po]

Пока ограничился созданием теста АЧХ на основе КРГ.
Встала проблема расчёта опорного SPL для частоты 1кГц. Пытаюсь подготовить сигнал - не получается. Либо на НЧ (25Гц) перегрузка, либо не хватает SPL на 1кГц - мала громкость.
Дело в том, что SPL НЧ для достижения равной громкости крайних НЧ и 1кГц весьма и весьма велик, к тому же слабо зависит от SPL на 1 кГц. Считаю фильтр, фильтрую, нормализую - ничего не выходит. Понял, что для достижения равной громкости хотя бы в 50-60 фон во всём звуковом диапазоне, требуется мощный усилитель и чувствительная мощная акустика. В итоге придётся ограничиваться более узкой полосой тестирования.
Универсальный тест 25-16кГц сделать не получилось из-за проблем с мощностью на крайних НЧ и неопределённостью КРГ в области крайних ВЧ. Поэтому, полагаю, придётся готовить сигнал для частных случаев, в зависимости от частотного диапазона и максимального SPL связки усилитель-акустика.

[Nikolay_Po]

Текста стандарта так и не нашёл. Для некоммерческой задачи покупать стандарт не стал.
За всё время нашёл только один код на MathLab Central, автор Jeff Tackett. Этот код и переписал для "R", добавив сплайновую интерполяцию для произвольных значений частот:

## Equal-loudness contour function.
## May complain ISO 226:2003.
## Contains sole and nearly official numbers I*ve found in the Internet at the moment.
## Site: "http://www.mathworks.com/matlabcentral", author Jeff Tackett.
## The function valid if 20<=Freq<=12500 (Hz) and 1<=Loudness<=100 (phon).
## Adopted to R by Nikolay_Po.

EqSpl <- function(Freq,Loudness)
{
    f <- c(20,25,31.5,40,50,63,80,100,125,160,200,250,315,400,500,630,800,1000,1250,1600,2000,2500,3150,4000,5000,6300,8000,10000,12500)
    af <- c(0.532,0.506,0.480,0.455,0.432,0.409,0.387,0.367,0.349,0.330,0.315,0.301,0.288,0.276,0.267,0.259,0.253,0.250,0.246,0.244,0.243,0.243,0.243,0.242,0.242,0.245,0.254,0.271,0.301)
    Lu <- c(-31.6,-27.2,-23.0,-19.1,-15.9,-13.0,-10.3,-8.1,-6.2,-4.5,-3.1,-2.0,-1.1,-0.4,0.0,0.3,0.5,0.0,-2.7,-4.1,-1.0,1.7,2.5,1.2,-2.1,-7.1,-11.2,-10.7,-3.1)
    Tf <- c(78.5,68.7,59.5,51.1,44.0,37.5,31.5,26.5,22.1,17.9,14.4,11.4,8.6,6.2,4.4,0,2.2,2.4,3.5,1.7,-1.3,-4.2,-6.0,-5.4,-1.5,6.0,12.6,13.9,12.3)
    A <- 4.47E-3*(10^(0.025*Loudness)-1.15)+(0.4*10^(((Tf+Lu)/10)-9))^af
    Lp <- ((10/af)*log10(A))-Lu+94
    Level <- splinefun(f,Lp,method="monoH.FC")
    return (Level(Freq))
}

## Example
EqSpl(1000,80)

Результаты хорошо согласуются с данными из разных источников и очень похожи на цифры из стандарта. Автор кода-первоисточника явно ссылается на ISO 226, похоже, видел сам стандарт. Да, точно видел - ссылается на конкретные разделы.
Для себя считаю тему закрытой.