Я имею дело с решением волнового уравнения в двух различных случаях, представленных на рисунке как случай A и случай B. Эти два рисунка были получены путем освещения плиты падающей волной, идущей справа (случай A) и слева (случай A). Случай Б). Я проверил решение для случая A. Оно кажется правильным. Таким образом, я использовал ту же методологию для получения решения для случая B. Однако я подозреваю, что решение для случая B неверно, даже если я не могу понять, почему. Независимо от правильного/неправильного результата волнового уравнения в случае B, я хотел бы знать, есть ли способ вывести решение волнового уравнения в случае B из решения случая A, применяя своего рода симметрию. Обратите внимание, чтоη( х )
не обязательно имеет какую-либо симметрию относительноИкс
.
Это мое решение. Учитывая волновое уравнение:
г2у( х )гИкс2+к2( 1 + п( х ) ) у( х ) = 0
Ее решение можно найти, используя метод функций Грина. Детали метода приведены для случая B. Применительно к случаю A волновое уравнение имеет решение:
у( х ) =с1е− j k х+с2ек х _+j k2∫л0ек | _ х - т |η( т ) у( т ) дт
где
Дж2= - 1
и
с1
и
с2
– постоянная интегрирования. Для определения констант необходимо рассмотреть граничные условия.
Если обратиться к случаю А на рисунке, то можно различить падающую волну и отраженную волну с правой стороны и прошедшую волну с левой стороны. В этом случае выполняются следующие граничные условия:
у( л ) = 1 +рлу′( L ) знак равно - j k ( 2 - y( Л ) )у( 0 ) =Тлу′( 0 ) знак равно - j k y( 0 )
Применяя эти граничные условия, различают два случая
т - х > 0
(левая сторона рисунка) и
т - х < 0
(правая сторона рисунка) можно найти дифференциальное уравнение в
Тл
.
Теперь я хотел бы повторить те же вычисления на второй картинке (Случай B). На этом рисунке падающая волна идет с левой стороны, а прошедшая — с правой. В случае А волна, распространяющаяся в том же направлении оси х, имеет положительный экспоненциальный член. У меня есть сомнения относительно рисунка B. В этом случае волна, бегущая в том же направлении оси x, зависит от− х
. Это верно? Я полагаю, что это просто вопрос условностей. Граничные условия для случая на рисунке B, конечно, переключаются:
у( 0 ) = 1 +рлу′( 0 ) знак равно - j k ( 2 - y( 0 ) )у( л ) =Тлу′( L ) знак равно - j k y( л )
Что касается случая A, я сначала дифференцирую общее решение:
у′( Икс ) знак равно - j kс1е− j k х+ к _с2е+ к х _+j k2∫л0ггИкс(ек | _ х - т |) η( т ) у( т ) дт
Теперь я рассмотрел два разных случая: whel
т - х > 0
(слева) и
т - х < 0
(верно). Для
т - х > 0
Я могу написать:
у′( Икс ) знак равно - j kс1е− j k х+ к _с2е+ к х _− j k U( х )
где
U( х ) =∫л0ггИкс(е- j k ( Икс - т )) η( т ) у( т ) дт
Затем
у′( Икс ) знак равно - j k (с1е− j k х−с2е+ к х _+ У( х ) )
В то время когда
т - х < 0
У меня есть:
у′( Икс ) знак равно - j k (с1е− j k х−с2е+ к х _− К( х ) )
где
К( х ) =∫л0ггИкс(еj k ( Икс - т )) η( т ) у( т ) дт
Граничные условия для
х = 0
соответствуют случаю
т - х > 0
, а граничные условия в
х = L
соответствуют случаю
т - х > 0
. Для
х = 0
у нас есть
у( 0 ) =с1+с2+ У( 0 )у′( 0 ) знак равно - j k (с1−с2+ У( 0 ) )
Так что
с2= 0
. Повторяя те же вычисления в
х = L
можно определить, что
с1"="едж к л
. Тогда решение:
у( х ) =еj k ( L - Икс )+j k2∫л0ек | _ х - т |η( т ) у( т ) дт
Теперь дело на рисунке Б. Вот тут я не уверен. Судя по картинке
е− j k х
является распространяющейся вперед волной. Зеленые функции предназначены для
х > L
:
{у1( х ) =С1е− j k ху1( л ) =А′
В то время как для
х < 0
:
{у2( х ) =с2ек х _у2( 0 ) =Б′
Затем
{у1( л ) =с1е− j k L"="А′у2( 0 ) =с2"="Б′⇒у1( х ) =А′еj k ( L - Икс )⇒у2( х ) =Б′ек х _
Вронскиан это
Вт"="∣∣∣А′еj k ( L - Икс )− j kА′еj k ( L - Икс )Б′ек х _j kБ′ек х _∣∣∣= 2 j kА′Б′едж к л
Так что функция Грина
г ( Икс , т ) знак равно⎧⎩⎨⎪⎪А′еj k ( L - Икс )Б′еj k т2 й кА′Б′едж к лА′еj k ( L - т )Б′ек х _2 й кА′Б′едж к л
Что сводится к
г ( Икс , т ) знак равно12 й кек | _ т - х |
Так что решение дифференциального уравнения для случая на рисунке B равно
у( х ) =с1е− j k х+с2ек х _+j k2∫л0ек | _ т - х |η( т ) у( т ) дт
Упакс
Qмеханик
Упакс
Qмеханик