Нестандартная конечно-разностная схема - Nonstandard finite difference scheme

Нестандартные конечно-разностные схемы это общий набор методов в числовой анализ что дает численные решения дифференциальные уравнения путем построения дискретной модели. Общие правила для таких схем точно не известны.[1][2]

Обзор

Конечно-разностная (FD) модель дифференциального уравнения (DE) может быть сформирована простой заменой производных приближениями FD. Но это наивный «перевод». Если мы буквально переводим с английского на японский, устанавливая взаимно однозначное соответствие между словами, исходное значение часто теряется. Точно так же наивная модель FD DE может сильно отличаться от исходной DE, потому что модель FD представляет собой разностное уравнение с решениями, которые могут сильно отличаться от решений DE. Для более технического определения см. Mickens 2000.[1]

Нестандартная (NS) конечно-разностная модель представляет собой бесплатный и более точный «перевод» дифференциального уравнения. Например, параметр (назовите его v) в DE может принимать другое значение ты в модели НС-ФД.

Пример

В качестве примера смоделируем волновое уравнение,

Наивная конечно-разностная модель, которую мы теперь называем стандартной (S) FD-моделью, находится путем аппроксимации производных с помощью FD-приближений. Центральная FD-аппроксимация первой производной второго порядка имеет вид

Применяя указанное выше приближение FD к , мы можем получить приближение ФД для ,

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

Если вставить раствор волнового уравнения (с ) в модель S-FD вы обнаружите, что

В целом потому что решение FD-приближения волнового уравнения не то же самое, что само волновое уравнение.

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

Таким образом, точная нестандартная конечно-разностная модель волнового уравнения имеет вид

Дополнительные подробности и расширения до двух и трех измерений, а также уравнений Максвелла можно найти в Cole 2002.[2]

Рекомендации

  1. ^ а б Миккенс, Р. (2000). Приложения нестандартных конечно-разностных схем. World Scientific.
  2. ^ а б Дж. Б. Коул, Высокоточный алгоритм Йи, основанный на нестандартных конечных различиях: новые разработки и проверки, IEEE Trans. по антеннам и распространению, т. 50, нет. 9, стр. 1185-1191 (2002).