PDE Mathematica - PullRequest
       42

PDE Mathematica

0 голосов
/ 30 марта 2020

Я новичок ie с Mathematica, поэтому, вероятно, мою проблему очень легко решить. Я хочу решить PDE:

Click to check the script

Первая проблема заключается в том, что программа не подставляет значения внутри w [x, y] и его вторых производных, чтобы использовать их как граничные условия. Он говорит мне, что Tag Equal in ((w1^(2,0))[x,y]==0)[x,y] is Protected

Другая проблема заключается в том, что он дает мне еще один тип ошибки:

"{NDsolve [{(w1 ^ (0,4)) [x, у] + (w1 ^ (2,2)) [х, у] + (w1 ^ (4,0)) [х, у] == 0, {w1 [0, у] == 0, w 1 [5 , у] == 0, ((w1 ^ (2,0)) [х, у] == 0) [0, у] == 0, ((w1 ^ (2,0)) [х, у] == 0) [5], у == 0}, {w1 [х, 0] == 3, w1 [х, 5] == 3, ((w1 ^ (0,2)) [х, у] == 0) [х, 0] == 0, ((w1 ^ (0,2)) [х, у] == 0) [х, 10] == 0}}, w 1 [х, у], {x, 0,10}, {y, 0,5}]} не является ни списком правил замены, ни действительной таблицей отправки и поэтому не может использоваться для замены. "

Заранее спасибо.


Я внес некоторые изменения.

Прежде всего я установил граничные условия, потому что у них была проблема. Более того, я также изменил уравнение.

Теперь он все еще дает мне ошибку:

"{NDsolve [{(w ^ (0,4)) [x, y] + (w ^ (2,2)) [х, у] + (ш ^ (4,0)) [х, у] == 2, {ш [0, у] == 0, ш [5, у] == 0 , Chix [0, у] == 0, Chix [5, у] == 0}, {ш [х, 0] == 0, ш [х, 10] == 0, чий [х, 0] = = 0, chiy [x, 10] == 0}}, w, {x, 0,5}, {y, 0,10}]} не является ни списком правил замены, ни действительной таблицей отправки, и поэтому не может использовать для замены. "

Click to check the script

Если нам удастся решить эту проблему, Mathematica должна дать нам ненулевое решение.

Большое спасибо!


Переписав его, он не дает никакого решения. И это тоже не дает ошибок. Похоже, он даже не заменяет граничные условия. Я решил изменить их, решив таким образом проблему другого рода. Это все еще не работает. Я думаю, что это очень странно, поскольку, если граничные условия были неправильными, это должно было обеспечить либо неопределенное, либо невозможное условие. Есть идеи?

Большое спасибо за помощь!

Click here to check the script

...