Дискуссионный математический форумМатематический форум
Математический форум Math Help Planet

Обсуждение и решение задач по математике, физике, химии, экономике

Теоретический раздел
Часовой пояс: UTC + 3 часа [ Летнее время ]
новый онлайн-сервис
число, сумма и дата прописью

Часовой пояс: UTC + 3 часа [ Летнее время ]




Начать новую тему Ответить на тему  [ Сообщений: 2 ] 
Автор Сообщение
 Заголовок сообщения: Дисктертное синус-преобразование Фурье от второй производной
СообщениеДобавлено: 02 май 2017, 14:44 
Не в сети
Начинающий
Зарегистрирован:
21 фев 2017, 16:14
Сообщений: 5
Cпасибо сказано: 0
Спасибо получено:
0 раз в 0 сообщении
Очков репутации: 1

Добавить очки репутацииУменьшить очки репутации
Здравствуйте! Никак не могу уяснить для себя вывод формулы дискретного синус-преобразования от второй производной решетчатой функции. Вторая производная представлена аппроксимацией по конечно - разностной схеме. Начну по порядку. Имеется дифференциальное уравнение для функции [math]\upsilon \left( {x,z} \right)[/math]: [math]\frac{{{\partial ^2}\upsilon }}{{\partial {z^2}}} + 2 \cdot i \cdot k\frac{{\partial \upsilon }}{{\partial x}} = 0[/math] с граничными условиями [math]\upsilon \left( {x,0} \right) = \upsilon \left( {x,{z_{\max }}} \right) = 0[/math]. Это уравнение дискретизируется по z. Причем вторая производная по [math]z[/math] аппроксимируется по конечно - разностной схеме следующим образом:[math]\frac{{{\partial ^2}\upsilon }}{{\partial {z^2}}} \sim \frac{{\upsilon \left( {x,z + \Delta z} \right) + \upsilon \left( {x,z - \Delta z} \right) - 2 \cdot \upsilon \left( {x,z} \right)}}{{\Delta {z^2}}}[/math]. Если обозначить дискретный оператор взятия второй производной по [math]z[/math] по конечно - разностной схеме как [math]D^2[/math], то, применяя его к последовательности значений исходной функции [math]\upsilon[/math], дискретизированной с шагом [math]\Delta z[/math], получим последовательность значений:
[math]{w_1}\left( x \right) = \frac{{{\upsilon _2}\left( x \right) - 2 \cdot {\upsilon _1}\left( x \right)}}{{\Delta {z^2}}}[/math]
[math]{w_j}\left( x \right) = \frac{{{\upsilon _{j + 1}}\left( x \right) + {\upsilon _{j - 1}}\left( x \right) - 2 \cdot {\upsilon _j}\left( x \right)}}{{\Delta {z^2}}}[/math]
[math]{w_{L - 1}}\left( x \right) = \frac{{{\upsilon _{L - 2}}\left( x \right) - 2 \cdot {\upsilon _{L - 1}}\left( x \right)}}{{\Delta {z^2}}}[/math]
Таким образом дискретное представление дифференциального уравнения будет выглядеть так:
[math]{D^2}\widetilde \upsilon \left( x \right) + 2 \cdot i \cdot k \cdot \frac{{\partial \widetilde \upsilon }}{{\partial x}} = 0[/math]
Далее необходимо применить к уравнению дискретное синус-преобразование [math]\widetilde W[/math]следующего вида:
[math]{U_l} = \sqrt {\frac{2}{L}} \cdot \sum_{j = 1}^{L - 1} {{\upsilon _j} \cdot \sin \left( {\frac{{\pi \cdot j \cdot l}}{L}} \right)}[/math]
Как утверждается в книге, если применить данное преобразование к дискретизированной версии второй производной, то есть [math]\widetilde W{D^2}\widetilde \upsilon[/math], то результирующая последовательность значений будет удовлетворять следующему соотношению:
[math]{W_l} = - 4 \cdot {\sin ^2}\left( {\frac{{\pi \cdot l}}{{2 \cdot L}}} \right) \cdot {U_l}[/math]
где [math]{U_l}[/math] результат дискретного синус-преобразования от исходных значений дискретизированной функции [math]\widetilde \upsilon[/math].
Я никак не смог вывести данной соотношение. В книге перед данной формулой сказано следующее: "Rearranging terms, we write the discrete sine transform [math]\widetilde W[/math] of [math]{D^2}\widetilde \upsilon[/math] as". И далее конечная формула. Насколько я понимаю данная фраза говорит о том что мы меняем местами оператор дифференцирования и оператор дискретного синус-преобразования. При этом получается примерно следующее:
[math]{W_l} = \sum_{j = 1}^{L - 1} { - 4 \cdot {\upsilon _j} \cdot \sin (} \frac{{\pi \cdot j \cdot l}}{L}) \cdot {\sin ^2}\left( {\frac{{\pi \cdot j}}{{2 \cdot L}}} \right)[/math]
Похоже я что - то недопонимаю. Буду рад любой подсказке и помощи. Заранее спасибо.

Вернуться к началу
 Профиль  
Cпасибо сказано 
 Заголовок сообщения: Re: Дисктертное синус-преобразование Фурье от второй производной
СообщениеДобавлено: 03 май 2017, 15:01 
Не в сети
Начинающий
Зарегистрирован:
21 фев 2017, 16:14
Сообщений: 5
Cпасибо сказано: 0
Спасибо получено:
0 раз в 0 сообщении
Очков репутации: 1

Добавить очки репутацииУменьшить очки репутации
Провел численный эксперимент в Mathcad и вяснилось что если дискретизировать функцию [math]f\left( x \right) = \sin \left( x \right) \cdot \cos \left( x \right)[/math] в интервале [math][0..\pi ][/math] то ни приведенная ни полученная мною формулы не подтверждаются. Окончательно не понимаю что происходит..

Вернуться к началу
 Профиль  
Cпасибо сказано 
Показать сообщения за:  Поле сортировки  
Начать новую тему Ответить на тему      Страница 1 из 1 [ Сообщений: 2 ]

 Похожие темы   Автор   Ответы   Просмотры   Последнее сообщение 
Косинус и Синус преобразование Фурье,вопрос по интегралам

в форуме Ряды Фурье и Интегральные преобразования

hatefiles

0

364

03 май 2016, 20:53

Шаблон для второй обыкновенной производной

в форуме Численные методы

HJey

0

246

11 июн 2018, 11:18

Найти значение второй производной в точке

в форуме Дифференциальное исчисление

RETU

17

707

08 июл 2018, 00:06

Геометрический смысл второй частной производной

в форуме Дифференциальное исчисление

HitGirl

6

347

10 мар 2020, 11:22

Простое выражение для второй производной по индексу Бесселя

в форуме Дифференциальное исчисление

lan1967

0

204

11 апр 2015, 13:03

Исследование функции с помощью первой и второй производной

в форуме Пределы числовых последовательностей и функций, Исследования функций

marinyshka

5

516

24 окт 2015, 12:39

Правая разностная схема с использованием второй производной

в форуме Дифференциальные и Интегральные уравнения

stut

3

511

05 сен 2014, 08:50

Равенство нулю второй производной ограниченной функции

в форуме Задачи со школьных и студенческих олимпиад

anpe0681

3

968

26 янв 2017, 01:49

Найти синус, косинус преобразования Фурье

в форуме Ряды Фурье и Интегральные преобразования

kalina_vladi

1

1026

03 май 2014, 14:44

Преобразование Фурье

в форуме Ряды Фурье и Интегральные преобразования

timdeygun

0

496

14 дек 2016, 19:23


Часовой пояс: UTC + 3 часа [ Летнее время ]



Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и гости: 7


Вы не можете начинать темы
Вы не можете отвечать на сообщения
Вы не можете редактировать свои сообщения
Вы не можете удалять свои сообщения
Вы не можете добавлять вложения

Найти:
Перейти:  

Яндекс.Метрика

Copyright © 2010-2023 MathHelpPlanet.com. All rights reserved