|
|
  |
Представление не монотонной не периодической финитной функции, подобно ряду Фурье |
|
|
|
Apr 16 2009, 16:36
|
    
Группа: Свой
Сообщений: 1 928
Регистрация: 11-07-06
Пользователь №: 18 731

|
Цитата(EUrry @ Apr 16 2009, 19:17)  Это всё понятно! Вопрос в том, что воможна ли модель с конечным числом параметров для абсолютно неизвестной априори зависимости во всем диапазоне изменения аргумента? Конечно возможна! Вопрос только в адекватности модели - один период синуса можно смоделировать и прямой, и параболой.  Как Вы справедливо заметили, основная головная боль - выбор базиса (или базисов). После этого f1, f2, ...,fn становятся известными. Фраза "для абсолютно неизвестной априори зависимости" мне непонятна - Вы же имеете "экспериментальные данные MEAS", по которым можно оценить удобство использования различных базисов (хотя бы на уровне периодические/непериодические) или просто какого-либо безбазисного набора функций (который можно предположить, исходя из физики процесса). А дальше уточнять коэффиценты всякими корреляциями и ковариациями до достижения требуемой точности модели. Ваша задача в частной постановке - веселая многоходовочка (типа, шаг вперед, два назад), на которой люди (не я) в 80-х легко диссеры писали (благо, хватало желающих с экспериментальными данными и народнохозяйственной важностью), а в общей - тема для серьезных мужей, в математике премногосведущих (тоже не я  ).
|
|
|
|
|
Apr 16 2009, 16:46
|

Гуру
     
Группа: Свой
Сообщений: 3 218
Регистрация: 14-11-06
Из: Н. Новгород
Пользователь №: 22 312

|
Цитата(xemul @ Apr 16 2009, 20:36)  Конечно возможна! Вопрос только в адекватности модели - один период синуса можно смоделировать и прямой, и параболой.  Фраза "для абсолютно неизвестной априори зависимости" мне непонятна - Вы же имеете "экспериментальные данные MEAS", по которым можно оценить удобство использования различных базисов В том то и дело, что не MEAS надо аппроксимировать, а входящие туда зависимости, которые могут себя вести произвольно: могут быть и более менее периодичны, могут иметь резонанстный характер (хотя, это пока лучше вообще не трогать!!!  ) и т. д. Цитата Ваша задача в частной постановке - веселая многоходовочка (типа, шаг вперед, два назад), на которой люди (не я) в 80-х легко диссеры писали (благо, хватало желающих с экспериментальными данными и народнохозяйственной важностью), а в общей - тема для серьезных мужей, в математике премногосведущих (тоже не я  ). Да вот и я тоже с такими разделами математики не знаком особо, поэтому и решил у сообщества спросить.
--------------------
Все не могут только сеять разумное, доброе, вечное: кому-то надо и пахать! Природа не терпит пустоты: там, где люди не знают правды, они заполняют пробелы домыслом. © Бернард Шоу
|
|
|
|
|
Apr 16 2009, 17:50
|
    
Группа: Свой
Сообщений: 1 928
Регистрация: 11-07-06
Пользователь №: 18 731

|
Цитата(EUrry @ Apr 16 2009, 20:46)  В том то и дело, что не MEAS надо аппроксимировать, а входящие туда зависимости, которые могут себя вести произвольно: могут быть и более менее периодичны, могут иметь резонанстный характер (хотя, это пока лучше вообще не трогать!!!  ) и т. д. Бр-р-р... Что же в таком случае MEAS? Я решил, что это набор(ы) измерений нескольких величин, для одной части которых нужно придумать модели, зависящие от другой части. Сколькомерный у Вас MEAS, и сколько моделей Вы по нему хотите придумать? Если двумерный (как на первой картинке), то получается, что Вам нужно придумать всего-то y = U(F i(x, a, b, c,...)) (через U я обозначил набор функций F i с коэффициентами a, b, c,..., аргументом которых является единственная переменная x). Курс численных методов, для которого это стандартная задача, проходят, по-моему, на 3-ем курсе. Решение задачи сводится к определению набора функций и вычислению их коэффициентов. Или, н-р, иначе - к исключению незначащих функций (которые не улучшают корреляцию модели с экспериментальными данными) из заданного набора. А с периодическими и резонансными как раз проще всего.
|
|
|
|
|
Apr 16 2009, 18:24
|

Гуру
     
Группа: Свой
Сообщений: 3 218
Регистрация: 14-11-06
Из: Н. Новгород
Пользователь №: 22 312

|
Цитата(xemul @ Apr 16 2009, 21:50)  Бр-р-р... Что же в таком случае MEAS? Я решил, что это набор(ы) измерений нескольких величин, для одной части которых нужно придумать модели, зависящие от другой части. Сколькомерный у Вас MEAS, и сколько моделей Вы по нему хотите придумать? Если двумерный (как на первой картинке)... Да, зря я эту картинку прилепил! Она наводит на мысль, что ее и нужно аппроксимировать, в чём, в принципе, нет проблем. Давайте еще так попробуем: Имеется некая система A, описываемая матрицей [A(x)] с коэффициентами, зависящими от аргумента x. Система окружена некими "мешающими", а, может, и вспомогательными (тут как метод поставить) подсистемами, описываемыми матрицами [B1(x)], [B2(x)], ..., [Bn(x)]. Результат измерений MEAS(х) содержит в себе параметры как системы А, так и параметры "мешающих" подсистем. Измерение представляется моделью MOLEL(х) = F{[A(x)], [B1(x)], [B2(x)], ..., [Bn(x)]}, где F - известное соотношение между элементами указанных матриц. Собственно, неизвестными как априори, так и апостериори являются зависимости от х элементов матриц [A(x)], [B1(x)], [B2(x)], ..., [Bn(x)], которые и нужно представить параметрическими моделями. Далее в пространстве этих параметров минимизируется оптимизационная функция OPT(x)=MEAS(x) - MOLEL(х). И, если всё проходит удачно, то находятся, как параметры самой системы A, так и параметры окружающих "мешающих"/вспомогательных подсистем. P. S. Кстати, в сообщении №9 ошибка была (надо запретить злую комбинацию ctrl+с -> ctrl+v!!!  ). Было так (зачеркнутое не нужно!) OPT=MEAS - MODEL=F(f1, f2, ...,fn). Удалось исправить!
Сообщение отредактировал EUrry - Apr 16 2009, 18:31
--------------------
Все не могут только сеять разумное, доброе, вечное: кому-то надо и пахать! Природа не терпит пустоты: там, где люди не знают правды, они заполняют пробелы домыслом. © Бернард Шоу
|
|
|
|
Guest_TSerg_*
|
Apr 17 2009, 08:58
|
Guests

|
Цитата(EUrry @ Apr 15 2009, 14:25)  Известно, что периодическую функцию можно представить в виде ряда Фурье. Даже некоторые не периодические представляют, исскуственно создав им период, который устремляется к бесконечности. А не придумали ли математики подобное разложение для произвольной не периодической не монотонной финитной функции? Например, представление экспериментальной кривой на рисунке. Хоть какое-то приближение. Если существует что-то подобное, то где можно посмотреть? Понимаю, что замахнулся не на того, но что поделать!  Вам сюда Ивахненко, МГУА http://www.gmdh.net/gmdh.htmhttp://iissvit.narod.ru/rass/vip17.htm
|
|
|
|
|
Apr 17 2009, 10:36
|
Участник

Группа: Участник
Сообщений: 70
Регистрация: 15-10-07
Из: Ιθάκη
Пользователь №: 31 371

|
Цитата(EUrry @ Apr 16 2009, 16:03)  С целью "извлечения" из результатов измерений входящих в него параметров, например, S-параметров СВЧ устройства. Да, в принципе, чего угодно, в зависимости от области исследований. Я думаю, в таком случае, базис для разложения должен быть не произвольным, а позволяющим дать физическую интерпретацию результатов аппроксимации. Например, если у Вас частотная зависимость S11, может быть использован базис в виде полюсных функций, где значение полюса wp определяет резонанс в анализируемом объекте: R = a1/(w-w1) + ... + aM/(w-wM). Программы типа HFSS, CST, MWO и подобные используют эту модель (иногда с успехом, иногда с провалом) для интерполяции/экстраполяции расчетных точек для ускорения частотного свипирования (я, например, не рискую использовать этот подход для анализируемой структуры общего типа). В некоторых ситуациях, когда анализируемый объект является распределенной структурой и отклик от него есть сумма сигналов отражения с разными запаздываниями (и переотражениями можно пренебречь), модель может быть представлена суммой экспонент R = a1*exp(jwt1) + ... + aM*exp(jwtM) с очевидной физической интерпретацией: ai - амплитуда отражения от i-ой неоднородности, ti = 2*zi/c - временное местоположение этого отражателя. Базисные функции в обоих случаях зависят от параметров {wi} и {ti} - нелинейно, поэтому необходимо использовать методы нелинейной оптмизации, отсюда возможен ряд проблем. Но при использовании начальной информации об анализируемом объекте задача, в принципе, решаемая. Оптимизационная задача может быть решена методами локальной оптимизации при адекватном выборе начального приближения.
|
|
|
|
|
Apr 18 2009, 19:44
|

Гуру
     
Группа: Свой
Сообщений: 3 218
Регистрация: 14-11-06
Из: Н. Новгород
Пользователь №: 22 312

|
Цитата(Ulysses @ Apr 17 2009, 14:36)  Например, если у Вас частотная зависимость S11, может быть использован базис в виде полюсных функций, где значение полюса wp определяет резонанс в анализируемом объекте: R = a1/(w-w1) + ... + aM/(w-wM). Программы типа HFSS, CST, MWO и подобные используют эту модель (иногда с успехом, иногда с провалом) для интерполяции/экстраполяции расчетных точек для ускорения частотного свипирования (я, например, не рискую использовать этот подход для анализируемой структуры общего типа). Тоже пробовал в HFSS быстрый свип - результат не был адекватен. Поэтому, больше не использую - лучше подождать, но быть более уверенным. Цитата В некоторых ситуациях, когда анализируемый объект является распределенной структурой и отклик от него есть сумма сигналов отражения с разными запаздываниями (и переотражениями можно пренебречь), модель может быть представлена суммой экспонент R = a1*exp(jwt1) + ... + aM*exp(jwtM) с очевидной физической интерпретацией: ai - амплитуда отражения от i-ой неоднородности, ti = 2*zi/c - временное местоположение этого отражателя. В моей задаче переотражения есть и, причем, используются с делом. Возможно, такая модель более подходит для импульсной рефлектометрии с целью определения местоположения неоднородностей. Цитата(TSerg) Вам сюда Ивахненко, МГУА Спасибо за ссылки, в ближайшее время посмотрю.
--------------------
Все не могут только сеять разумное, доброе, вечное: кому-то надо и пахать! Природа не терпит пустоты: там, где люди не знают правды, они заполняют пробелы домыслом. © Бернард Шоу
|
|
|
|
|
Apr 20 2009, 07:28
|

Гуру
     
Группа: Свой
Сообщений: 3 218
Регистрация: 14-11-06
Из: Н. Новгород
Пользователь №: 22 312

|
Цитата(scifi @ Apr 20 2009, 10:57)  И моё мнение таково, что задача нерешаемая. Но не из-за недоразвитости человечества, а из-за неправильной постановки задачи. Давно известно, что правильно заданный вопрос - это половина ответа. По-моему, постановка задачи ясна. Если бы было всё так просто и была создана абсолютно универсальная модель, то не "извращались" бы с теми же "окнами". Повторюсь, целью вопроса было "прощупывание современной почвы" в данном направлении. За спрос, как говорится, не бьют. В данном случае считаю, что спрашивать надо обширнее, а уж что достанется в ответ, тому и будешь рад. Информации для анализа больше будет.
--------------------
Все не могут только сеять разумное, доброе, вечное: кому-то надо и пахать! Природа не терпит пустоты: там, где люди не знают правды, они заполняют пробелы домыслом. © Бернард Шоу
|
|
|
|
|
Apr 20 2009, 12:24
|
Группа: Новичок
Сообщений: 3
Регистрация: 19-04-09
Из: Украина
Пользователь №: 47 953

|
Я много времени провел за изучением функционального анализа. Представление функции в виде разложения по какому либо ряду - хорошо изучено. И вроде бы серьезной альтернативы еще не встречал.
|
|
|
|
|
  |
1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)
Пользователей: 0
|
|
|