К вопросу аппроксимации эмпирических зависимостей

Исследование результатов математического моделирования эмпирических зависимостей. Аппроксимация частотной характеристики акустической камеры. Количественные оценки параметров модели функции аппроксимации амплитудно-модулированного затухающего сигнала.

Рубрика Коммуникации, связь, цифровые приборы и радиоэлектроника
Вид статья
Язык русский
Дата добавления 30.04.2018
Размер файла 238,0 K

Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже

Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.

Размещено на http://www.allbest.ru/

К вопросу аппроксимации эмпирических зависимостей

Середа С.Н. кандидат технических наук, Доцент, Владимирского государственного университета имени Александра Григорьевича и Николая Григорьевича Столетовых

Аннотация

В статье представлены результаты математического моделирования эмпирических зависимостей на примере аппроксимации частотной характеристики акустической камеры. Предложена оригинальная модель функции аппроксимации амплитудно-модулированного затухающего сигнала. Приводятся количественные оценки параметров модели с использованием линеаризации нелинейных функций. Результаты исследования могут быть использованы для разработки шумозащитных экранов и решении прикладных задач математической статистики.

Ключевые слова: функция аппроксимации, линеаризация, математическая модель, частотная характеристика, шумозащита.

The article presents the results of mathematical modeling of empirical dependencies on the example of the approximation of acoustic chamber frequency characteristic. An original model of the approximation function of an amplitude-modulated damped signal is proposed. Quantitative estimates of the model parameters using the linearization of nonlinear functions are given. The results of the research can be used to develop noise shielding screens and solve applied problems of mathematical statistics.

Keywords: approximation function, linearization, mathematical model, frequency response, noise protection.

Во многих областях науки и техники возникает задача обработки результатов эксперимента и идентификации зависимостей между параметрами процесса или системы, которая может сводиться к поиску некоторой аналитической функции аппроксимации табличных эмпирических данных, что позволяет в дальнейшем проводить расчеты путем интерполяции и экстраполяции функции и прогнозировать динамику процесса. Хотя разработаны различные методы решения задачи аппроксимации, как-то кусочно-линейная аппроксимация, аппроксимация тригонометрическими функциями, полиномами Чебышева и др., в некоторых случаях не удается с помощью таких подходов получить приемлемо хорошее приближение, либо аналитическая функция не позволяет получить общую математическую модель исследуемого процесса как, например, при использовании сплайнов. При этом выдвижение гипотезы о виде функции аппроксимации должно учитывать дополнительную информацию о характере процесса, что скорее является искусством. Кроме того, определение значений параметров функции аппроксимации, обеспечивающих наилучшее приближение в смысле минимизации критерия среднеквадратического отклонения, может проводиться как аналитически, так и алгоритмически [1, C. 292], [2, C. 145], [8], [10].

Рассмотрим задачу поиска функции аппроксимации частотной характеристики акустической камеры, используемой в качестве лабораторного макета для оценки эффективности звукопоглощающих экранов [5, C. 18].

Цель исследования заключалась в определении собственной амплитудно-частотной характеристики (АЧХ) акустического экрана на основе резонатора Гельмгольца, установленного в лабораторную акустическую камеру [6]. Измерения уровня шума в камере проводились с помощью шумомера ВШВ-003. В качестве источника звука использовался генератор звуковой частоты TR-0157/002, усилитель звуковой частоты и акустическая система из двух колонок “ОДА”. Измерения проводились в треть октавных полосах частот, в диапазоне от 25 Гц до 8000 Гц. В данном диапазоне частот преобладает энергия большинства производственных, бытовых и транспортных шумов [7, C. 19].

Методика измерений амплитудно-частотной характеристики исследуемой системы включала следующие этапы [4, C. 28]:

Сначала замерялась интенсивность отраженного от глухой стенки звука (без резонатора). Производили три замера уровня звукового давления на каждой частоте перед глухой стенкой, затем находили среднее значение Lг.

Далее в камеру помещался экран, направленный щелями к источнику звука, и проводились измерения звукового давления перед экраном при разных настройках геометрических параметров резонатора Lр. Для каждого состояния было произведено несколько замеров и определялось средний статистический уровень звукового давления, чтобы минимизировать ошибки измерений.

Эмпирическая АЧХ звукопоглощающего экрана определялась как разность результатов измерений двух предыдущих этапов

Lэ=Lг - Lр (1)

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

Частотная характеристика звукопоглощающего экрана Lэ(x) (рис.1) получена при следующих значениях геометрических параметров резонатора: глубина короба h=75 мм, ширина щели b=0,5 мм, расстояние между щелями l=90,75 мм, глубина щели i=15 мм.

Рис. 1 - Эмпирическая частотная характеристика звукопоглощающего экрана

Для поиска аналитической функции, аппроксимирующей данные эксперимента, сделаем допущение о характере изменения анализируемой кривой. Во-первых, функция не является гладкой и содержит как максимумы, так и провалы на частотах, кратных резонансной частоте экрана. Во-вторых, как видно из рисунка, линия тренда TL(x), вычисленная путем усреднения, представляет собой монотонно убывающую кривую с ростом частоты. Проведем предварительную обработку путем нормирования эмпирической характеристики на ее максимум, и в дальнейшем будем проводить аппроксимацию нормированной характеристики L(x). Поскольку в эксперименте измерения проводились на наборе одно-тональных гармонических сигналов звукового диапазона частот, то можно принять гипотезу о частотном отклике звукопоглощающего экрана в форме амплитудно-модулированного затухающего сигнала соответствующей частоты [9]. Это позволяет определить математическую модель анализируемой характеристики в виде

S(x)=S1(x)+S2(x)*S3(x) (2)

где S1(x) - функция тренда; S2(x) - функция амплитуды; S3(x) - несущая функция; x - частотная координата.

Относительно вида функций, являющихся составными частями модели могут быть различные предположения. Предварительный анализ данных эксперимента, а также линии тренда позволяет поставить задачу оценки параметров функции аппроксимации вида

S(x)=exp (a1•xb1)+a2•exp (b2•x) • sin (щ•x+ц) (3)

где a1, a2, b1, b2 - параметры функции; щ - частота резонатора; ц - фазовый сдвиг.

Для определения значений параметров искомой функции воспользуемся методом приведения зависимостей к линейному виду для её составляющих. Так можно получить формулы линеаризации для функции тренда S1(x) (4) и функции амплитуды (5) соответственно.

Y1= ln(S1(x))= a1•xb1; U1=A+b1•Z1; A=lg(a1); Z1=lg(x); U1=lg(-Y1) (4)

Y2= S2(x)= a2•exp (b2•x); U2=ln(a2)+b2•x; U2=ln(Y2) (5)

В качестве критерия приближения принимаем минимум среднеквадратического отклонения

(6)

Достижению минимума целевой функции (6) в случае линейной зависимости с двумя параметрами ai, bi соответствует равенство нулю первых частных производных [1, C. 345], [3, C. 153]

(7)

Раскрытие выражения (7) в отношении функций S1(x) и S2(x) приводит к системам нормальных уравнений (7) и (8) соответственно

(8)

(9)

Решение указанных систем уравнений в программе математического моделирования MathCad с помощью встроенной функции lsolve(ai,bi) дает значения параметров a1=-0.015, b1=0.505, a2=0.288673, b2=-0.004712. При этом, значение ошибки составляет Q1=2.7*10-7 и Q2=3.8*10-7.

Оценка параметров функции S3(x) проводилась алгоритмически методом итераций, где критерием остановки итерационного процесса было совпадение глобального максимума и минимума на частоте провала функции аппроксимации S(x) с данными эксперимента. В результате поиска были получены значения параметров щ=420, ц=6.41, что позволяет установить их прямую взаимосвязь с резонансной частотой акустического экрана fp, определяемой по его геометрическим характеристикам [6]

(10)

где fm - частота глобального максимума эмпирической функции L(x);

fр - резонансная частота акустического экрана, определяемая по [6]

где b - ширина щели; h - глубина короба; i - глубина щели; l - расстояние между щелями; с=331 м/с - скорость звука в воздухе при температуре воздуха 20°С.

Результаты математического моделирования функции аппроксимации с установленными значениями параметров показаны на рис. 2. В качестве меры адекватности модели примем коэффициент корреляции, показывающий степень совпадения кривых, который составил 97,4%. Как видно из рисунка предложенная аналитическая функция достаточно точно локализует частоты, в которых наблюдаются глобальный максимум и основной провал при приемлемом приближении в остальном диапазоне частот. Полученная модель позволяет в дальнейшем рассчитать частотную характеристику шумозащитного экрана при изменении его конструктивных параметров и оценивать степень звукопоглощения.

Рис. 2 - Аппроксимация нормированной частотной характеристики L(x) функцией S(x)

Список литературы

1. Амосов А.А. Вычислительные методы для инженеров: Учеб. пособие / А.А. Амосов, Ю.А. Дубинский, Н.В. Копченова. - М. : Высш. шк., 1994. - 544 с.

2. Большаков А.А. Методы обработки многомерных данных и временных рядов: Учеб. пособие для вузов / А.А. Большаков, Р.Н. Каримов. - М. : Горячая линия - Телеком, 2007. - 522 с.

3. Гумеров А.М. Математическое моделирование химико-технологических процессов: Учеб. пособие / А.М. Гумеров. - 2-е изд. - СПб. : Издательство «Лань», 2014. - 176 с.

4. Гуськов П.М. Выбор уровня звукового давления в лабораторной акустической камере при исследовании характеристик шумозащитных экранов / П.М. Гуськов // Методы и устройства передачи и обработки информации. - 2017. - №1. - С. 28-31.

5. Калиниченко М.В. Некоторые аспекты применения резонансных поглотителей на урбанизированных территориях / М.В. Калиниченко // Машиностроение и безопасность жизнедеятельности. - 2013. - № 4. - С. 18-24.

6. Патент ПМ №139578 Российская Федерация, E01F8/00. Шумопоглотитель / Булкин В. В. Калиниченко М. В.; заявитель и патентообладатель Владимирский государственный университет имени Александра Григорьевича и Николая Григорьевича Столетовых. - №2013121443; заявл. 07.05.2013; опубл. 20.04.2014.

7. Bulkin V. V. Preliminary results of development and test of the lab anechoic chamber / V. V. Bulkin, M. V. Kalinichenko // Noise Theory and Practice, vol. 2. - 2016. - №2. - P. 19-26.

8. Powell M. J. D. Approximation Theory and Methods / M.J.D. Powell. - Cambridge University Press, 1981. - 339 p.

9. Sereda S.N. Mathematical Model of Acoustic Baffle based on Helmholtz Resonator / S.N. Sereda, V.V Bulkin, M.V. Kalinichenko, P.M. Guskov P.M. // IEEE Conference 2017 Dynamics of Systems, Mechanisms and Machines (Omsk, 2017).

10. Steffens Karl-Georg. The History of Approximation Theory. From Euler to Bernstein / Karl-Georg Steffens. - Boston: Birkhauser, 2006. - 219 p.

Размещено на Allbest.ru

...

Подобные документы

  • Определение импульсной характеристики фильтра. Расчет амплитудно- и фазово-частотной характеристик и методами разложения в ряд Фурье, наименьших квадратов и частотной выборки. Построение графиков и оценка точности аппроксимации (абсолютной погрешности).

    курсовая работа [677,0 K], добавлен 21.12.2012

  • Этапы процесса синтеза электрической схемы. Требования к частотной характеристике фильтра. Аппроксимация заданной амплитудно-частотной характеристики. Порядок расчета и соображения по методике настройки активных фильтров. Расчет величин элементов схемы.

    курсовая работа [490,3 K], добавлен 27.01.2010

  • Исследование влияния параметров проводной линии на характеристики ее выходного сигнала. Методика измерения параметров выходного импульса. Искажение сигнала вследствие частотной дисперсии. Описание интерфейса взаимодействия с пользователем модели.

    лабораторная работа [398,0 K], добавлен 06.07.2009

  • Проведение анализа системы, содержащей идеальный операционный усилитель. Определение вида выходного сигнала при известном напряжении на входе во временной области. Построение графика амплитудно-частотной (АЧХ) и фазо-частотной (ФЧХ) характеристики.

    курсовая работа [552,2 K], добавлен 14.02.2013

  • Расчет передаточных функций разомкнутой и замкнутой системы автоматического регулирования при отрицательной единичной обратной связи. Исследование характеристик САР: амплитудно-фазовой частотной, АЧХ, ФЧХ, логарифмической амплитудно-частотной и ЛФЧХ.

    контрольная работа [709,2 K], добавлен 06.12.2010

  • Аппроксимация кривой разгона объекта управления уравнением звена второго порядка с запаздыванием. Величина достоверности аппроксимации, передаточные функции датчика, преобразователя и исполнительного механизма. Проверка полученных систем на устойчивость.

    курсовая работа [779,2 K], добавлен 18.03.2014

  • Аппроксимация полиномом седьмой степени экспериментальной зависимости коэффициента усиления заданного усилительного каскада на полевом транзисторе типа 2П905А(119J). Определение параметров нелинейности третьего порядка и выбор режима работы каскада.

    курсовая работа [467,6 K], добавлен 01.04.2013

  • Графическое и аналитическое решение трансцендентного уравнения. Выполнение аппроксимации вольтамперной характеристики диодов различных видов методом полинома третьего порядка. Определение реакции цепи на входное воздействие при помощи интеграла Дюамеля.

    контрольная работа [3,3 M], добавлен 15.08.2012

  • Аппроксимация амплитудно-частотной характеристики фильтра. Определение передаточной функции фильтра нижних частот в области комплексной частоты. Схемотехническое проектирование устройства и его конструкторская реализация в виде узла с печатным монтажом.

    курсовая работа [330,8 K], добавлен 09.06.2015

  • Проектирование устройств фильтрации по рабочим параметрам. Виды аппроксимации частотных характеристик. Моделирование разрабатываемого фильтра на функциональном уровне в MathCAD, в частотной и временной областях, в нормированном и денормированном виде.

    курсовая работа [2,8 M], добавлен 28.06.2011

  • Исследование принципа действия и устройства коаксиального фильтра СВЧ диапазона. Построение амплитудно-частотной характеристики в заданном диапазоне частот. Проведение снятия зависимости амплитуды напряжения от частоты сигнала при отключенном фильтре.

    лабораторная работа [16,8 K], добавлен 28.10.2013

  • Назначение, типы и аппроксимация характеристик цифровых и аналоговых фильтров. Разработка на языке MATLAB программы моделирования ФВЧ методом Баттерворта, построение графиков амплитудно- и фазо-частотной характеристик; построение Simulink – модели.

    курсовая работа [883,8 K], добавлен 17.06.2011

  • Расчёт энергетических характеристик сигналов и информационных характеристик канала. Определение кодовой последовательности. Характеристики модулированного сигнала. Расчет вероятности ошибки оптимального демодулятора. Граничные частоты спектров сигналов.

    курсовая работа [520,4 K], добавлен 07.02.2013

  • Определение ширины полосы пропускания и допустимого коэффициента шума приёмника. Выбор числа поддиапазонов и их границ. Вычисление входной цепи, шумоподавления, регулировки громкости, стабилизатора и неравномерности амплитудно-частотной характеристики.

    курсовая работа [336,0 K], добавлен 30.10.2013

  • Определение параметров линейной схемы на резонансной частоте. Нахождение передаточной функции цепи по напряжению. Процесс построения управляющего сигнала. Отклик схемы на спектр амплитудно-модулированного колебания. Импульсная характеристика схемы.

    курсовая работа [1,1 M], добавлен 08.10.2012

  • Проектирование амплитудно–модулированного СВЧ–передатчика с частотной модуляцией. Расчет задающего генератора на диоде Ганна и выходного усилителя на лавинно–пролетном диоде. Выбор конструкции и эквивалентной схемы, определение электронного режима.

    курсовая работа [160,4 K], добавлен 20.09.2011

  • Рассмотрение методов измерения параметров радиосигналов при времени измерения менее и некратном периоду сигнала. Разработка алгоритмов оценки параметров сигнала и исследование их погрешностей в аппаратуре потребителя спутниковых навигационных систем.

    дипломная работа [3,6 M], добавлен 23.10.2011

  • Три схемы модуляции: амплитудная, угловая и импульсная. Особенности и подходы к реализации данных схем модуляции, предъявляемые к ним требования. Схемы перемножителей и направления исследования их элементов. Спектр амплитудно-модулированного сигнала.

    контрольная работа [735,4 K], добавлен 13.06.2012

  • Способы аппроксимации кривой разгона апериодическим звеном первого порядка с запаздыванием. Оптимальные настройки регулятора (метод Копеловича). Нахождение передаточной функции замкнутой системы. Моделирование АСР с использованием программы 20-sim.

    контрольная работа [418,7 K], добавлен 11.05.2012

  • Моделирование процесса дискретизации аналогового сигнала, а также модулированного по амплитуде, и восстановления аналогового сигнала из дискретного. Определение системной функции, комплексного коэффициента передачи, параметров цифрового фильтра.

    курсовая работа [1,5 M], добавлен 07.01.2014

Работы в архивах красиво оформлены согласно требованиям ВУЗов и содержат рисунки, диаграммы, формулы и т.д.
PPT, PPTX и PDF-файлы представлены только в архивах.
Рекомендуем скачать работу.