Доброе Время Суток, дорогие форумчане.
Очень  хотелось бы услышать совет(ы) от гуру форума сего. Ну или предложить фриланс ;) ( это уж как получится).
Итак,дано:
Непрерывная нестационарная математическая модель, описывающая распространение шума с учетом некоторых из основных факторов диссипативной среды. Предложенная модель построена на основе систем уравнений Максвелла и телеграфных  уравнений c учетом того факта, что сила тяжести и стратификация геофизических сред сильно модифицируют распространяющиеся в них звуковые волны.

ClicToZoom



Условия согласования ( рассчитаны ( согласование начальных и граничных):

ClicToZoom




ClicToZoom





Описывю я это все как : 

ClicToZoom






В Результате упрощtний и перехода к дискретной форме ( конечно-разностная аппроксимация) ,получилось :

ClicToZoom




Соб-сно,вопрос. Наиболее оптимальное (эффективное и не сильно трудоемкое) решение этой задачи численными методами (к конечно-разностной форме то уже пришли) с построением  графиков распределения.
Матлаб ? Фемлаб? C++ ? Delphi ? . 
Опыта работы с пакетами прикладными не имею.
Готовые примеры (по тематике) , ссылки, информация ( буду очень признателен).
 В Теории понятно, а как это переложить на Матлаб\Фемлаб , да запрограммировать ,  да чтоб с моделью распределения.
В случае аутсорсинга ( что очень предпочтительно и если кто то готов взяться) , -  материальное вознаграждение.

Фемлаб решает методм конечных элементов, что в случае прямоугольной сетки и простой геометрии модели порождает дополнительные сложности по сравнению с методом конечных разностей во временной области (FDTD). Ну и поскольку у вас две компоненты будет учавстсовать, то желательно использовать staggered grid (смещённую на полуширину ячейки сеть). Смысл подхода легко гуглится.

Поэтому я бы FDTD и посоветовал. Тем более ваша задача видимо будет решаться простым методом прогонки, т.е. внешний цикл по t и три внутренних по координатам. Начальные условия есть, граничные условия заданы. Всё просто до безобразия. 

С++, Delphi -- конечно хорошо, но в этом случае придётся вам писать всё с нуля. Лишние временные и трудовые затраты, взамен, возможно, получите небольшой прирост скорости относительно Matlab. 

Поэтому Matlab и FDTD. Реализаций в открытом доступе предостаточно. Гуглится по словам: "Matlab FDTD"

В аттаче прилагаю небольшой пример для Matlab решения уравнений Максвелла для квазистационарных полей в 2D.

Присоединённый файл ( Кол-во скачиваний: 29 )
 TwoD_FDTD.m 8,22 Kb

2W4FhLF

Спасибо за наводку и пример, скачал мануал по Матлабу.
Буду учиться, пробовать , о результатах отпишусь в этой теме. ПО крайней мере  Дискретная модель у меня уже есть, ее в матлаб "сунуть" можно.


2 Модератор : Может ли в Теме дальше вестись дискуссия,обсуждение вопроса ( Ответ на вопрос вроде уже дали, но пока еще не реализованно, могут возникнуть сложности). 


Оффтоп, но это уже на будущее  пока наверное не осилю.
Возможно ли в Матлабе ( как в Фемлабе) импортировать из CAD фигуру-чертеж? ( или описать его, это же собственно граничные условия, - ведь параллелепипед воокруг здания терминала аэропорта( уже есть в модели, смотрим пост выше) + в идеале само здание терминала внутри этого параллелепипеда ).
Ну это наверное  совсем сложно.

Это сообщение отредактировал(а) MrChe - 20.12.2010, 02:34

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

<< Непрерывная нестационарная математическая модель, описывающая распространение шума с учетом некоторых из основных факторов диссипативной среды. Предложенная модель построена на основе систем уравнений Максвелла и телеграфных  уравнений c учетом того факта, что сила тяжести и стратификация геофизических сред сильно модифицируют распространяющиеся в них звуковые волны. >>


Давненько, конечно, я все это проходил и изрядно запамятовал; и тем не менее осмеливаюсь смиренно спросить: каким образом уравнения Максвелла, описывающие распространение электромагнитных волн, могут помочь в описании распространения звука, подчиняющегося уравнениям гидродинамики?