Пятница, 29.11.2024
Мой сайт
Меню сайта
Статистика

Онлайн всего: 70
Гостей: 70
Пользователей: 0
Главная » 2017 » Февраль » 17 » Начальные условия волнового уравнения на примере косинуса
22:12
Начальные условия волнового уравнения на примере косинуса

Начальные условия волнового уравнения на примере косинуса

Тут:
http://akostina76.ucoz.ru/blog/2017-02-11-3842

… быстро придумалось решение волнового уравнения как, например:
u(x,t) = cos(6*t-2*x)
Константа «a» для этого уравнения  = 6/2=3.
Взяв начальные значения из этого решения можно увидеть какой вклад каждое из них вносит в итоговую волну.
Отсюда:
http://akostina76.ucoz.ru/blog/2017-02-12-3846


и:

Осталось всё аккуратно посчитать для конкретного косинуса:

Подстановка «a», вычисление интеграла:

Всё вместе:

Вывод до начального значения, т.е косинуса тут для проверки.
Интересно то какой вклад в решение дают обе части u(x,t), т.е слагаемые с φ и интеграл с ψ.
Обсчёт всего этого в Maple:

 

u_1(x,t) – функция от φ

u_1(x,t) – функция от ψ

Но если их просуммировать то будет ранее полученный вариант:

Лично у меня пространственного воображения не хватает чтобы понять почему вклады начальной формы и начальной скорости такие и как получается что при суммировании выходит такая форма. Вот они вместе, одно приподнято над другим:

Сразу видно насколько сильно поменяется картина если убрать одну из составляющих.
По такому графику, т.е по u(x) при t=0 и 0.05 видно чуть лучше:

За такой маленький промежуток времени вклад u_1() почти не поменялся, а вот вклад скорости – интеграла подрос довольно существенно. Он-то и сдвинул итоговую волну вправо на видимое расстояние. До «впадины» вклад «горба» положительный потому левая часть увеличивается, сдвигая волну вправо. А после «впадины» вклад отрицательный, но это очень удачно сдвигает волну на этом участке тоже вправо.
Три графика:

Все тексты в Maple:
omega:=6:
k:=2:
a:=omega/k:
u_cos:=(x,t)->cos(omega*t-k*x);
fi_t_0:=unapply(u_cos(x,0),x);
psi_t:=unapply(diff(u_cos(x,t),t),x,t);
psi_t_0:=(x)->psi_t(x,0);
u_1:=(x,t)->(fi_t_0(x-a*t)+fi_t_0(x+a*t))/2;
u_2:=(x,t)->int(psi_t_0(s),s=x-a*t..x+a*t)/(2*a);
u:=(x,t)->u_1(x,t)+u_2(x,t);
plot3d(u(x,t),x=0..6,t=0..1,axes=boxed);
plot3d(u_1(x,t),x=0..6,t=0..1,axes=boxed);
plot3d(u_2(x,t),x=0..6,t=0..1,axes=boxed);
plot3d({u_1(x,t),2+u_2(x,t),6+u_1(x,t)+u_2(x,t) },x=0..6,t=0..1,axes=boxed);
plot({u_1(x,0),2+u_2(x,0),4+u_1(x,0)+u_2(x,0) ,u_1(x,0.05),2+u_2(x,0.05),4+u_1(x,0.05)+u_2(x,0.05)},x=0..6,axes=boxed);
plot({u_1(x,0),2+u_2(x,0),4+u_1(x,0)+u_2(x,0) ,u_1(x,0.1),2+u_2(x,0.1),4+u_1(x,0.1)+u_2(x,0.1),u_1(x,0.2),2+u_2(x,0.2),4+u_1(x,0.2)+u_2(x,0.2)},x=0..6,axes=boxed);
p/s

Здесь у меня была ошибка:
http://akostina76.ucoz.ru/blog/2017-02-16-3862

Начальной формой в виде синуса функция задаётся с негласным выполнением условий в начале и в конце струны:

(скобка в формуле не там стояла).
Но это скорее везение и тоже своеобразное угадывание – подгонка. В общем случае эти концы запросто могут болтаться, т.к их пока, действительно не закрепили. 

Просмотров: 272 | Добавил: akostina76 | Рейтинг: 0.0/0
Всего комментариев: 0
Имя *:
Email *:
Код *:
Форма входа
Поиск
Календарь
Архив записей
Друзья сайта
  • Официальный блог
  • Сообщество uCoz
  • FAQ по системе
  • Инструкции для uCoz
  • Copyright MyCorp © 2024
    Бесплатный конструктор сайтов - uCoz