Skip to content

Commit 16e62a4

Browse files
authored
Update Lab7.md
1 parent 4056a78 commit 16e62a4

File tree

1 file changed

+40
-2
lines changed

1 file changed

+40
-2
lines changed

reports/Lab7.md

Lines changed: 40 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -8,7 +8,7 @@
88
- [Решение волнового уравнения при помощи явной разностной схемы крест](#решение-волнового-уравнения-при-помощи-явной-разностной-схемы-крест)
99
- [Решение волнового уравнения при помощи схемы с весами](#решение-волнового-уравнения-при-помощи-схемы-с-весами)
1010
- [Построение спектров отклонения струны, скорости и энергии](#построение-спектров-отклонения-струны-скорости-и-энергии)
11-
5. [Заключение и выводы 📝](#4-заключение-и-выводы-)
11+
4. [Заключение и выводы 📝](#4-заключение-и-выводы-)
1212

1313
## 1. Введение и цели работы 🎓
1414

@@ -264,10 +264,48 @@ Cпектральный признак дает следующие резуль
264264

265265
![Зависимость ошибки для схемы с весом 0.35](../images/Lab7/DependenceOfMaxErrorOnTimeStableWithSigma0.35UnstableWithSigma0.20.png)
266266

267-
**Рисунок 7 - Зависимость ошибки численного решения неявной схемой с весом $\sigma = 0.35$ от времени**
267+
**Рисунок 9 - Зависимость ошибки численного решения неявной схемой с весом $\sigma = 0.35$ от времени**
268268

269269
Данная схема является безусловно устойчивой, что и подтверждает полученный график ошибки.
270270

271271
---
272272

273+
### Построение спектров отклонения струны, скорости и энергии
273274

275+
Используя быстрое преобразование Фурье, построим спектры отклонения струны, скорости и энергии, а также исследуем, как изменится их внешний вид при изменении $\delta$:
276+
277+
При $\delta = 0.05$:
278+
279+
![Спектры при дельта 0.05](../images/Lab7/SpectrumOfDeflectionVelocityEnergyWithSigma0.25Delta0.05.png)
280+
281+
**Рисунок 10 - Спектры отклонения при $\delta = 0.05$**
282+
283+
При $\delta = 0.10$:
284+
285+
![Спектры при дельта 0.05](../images/Lab7/SpectrumOfDeflectionVelocityEnergyWithSigma0.25Delta0.10.png)
286+
287+
**Рисунок 10 - Спектры отклонения при $\delta = 0.10$**
288+
289+
Во всех случаях спектр отклонения u остаётся существенно ниже спектров скорости и энергии, что связано с тем, что в рассматриваемые моменты большая часть энергии всё ещё сосредоточена в кинетической форме, а амплитуда смещений мала по сравнению со скоростями.
290+
291+
Уменьшение $\delta$ (узкий импульс) приводит к заметному расширению спектров по модам k: максимумы спектра скорости и энергии становятся выше, а вклад высокочастотных мод (большие k) выражен сильнее. Узкий прямоугольный импульс по x ближе к дельта‑функции и содержит много высокочастотных гармоник, поэтому энергия распределяется по широкому диапазону мод.
292+
293+
При увеличении $\delta$ до `0.10` начальный профиль скорости становится шире и сглаженнее. В этом случае спектры скорости и энергии концентрируются преимущественно в низких модах: амплитуда у малых k растёт, а на высоких k заметно падает по сравнению со случаем $\delta = 0.05$.
294+
295+
Таким образом, рост $\delta$ приводит к перераспределению энергии от высокочастотных мод к длинноволновым нормальным колебаниям струны, что хорошо согласуется с физической интерпретацией более мягкого и протяжённого по пространству удара.
296+
297+
---
298+
299+
## 4. Заключение и выводы 📝
300+
301+
В работе была реализована и исследована полная цепочка постановки, численного решения и спектрального анализа одномерного волнового уравнения для натянутой струны с прямоугольным импульсом начальной скорости.
302+
303+
Главным результатом стало практическое подтверждение корректности разложений по собственным модам: реализованное точное решение в виде тригонометрического ряда хорошо совпадает с результатами численных схем при разумных шагах сетки. Явная крестовая схема показала ожидаемое поведение: при выполнении условия Куранта решение устойчиво, фронты импульса воспроизводятся достаточно точно, а максимальная ошибка остаётся на уровне $10^{-5}- 10^{-4}$ и медленно нарастает во времени за счёт накопления погрешностей и отражений волны от границ.
304+
305+
Схема с весами при $\sigma = 0.25$ подтвердила свою устойчивость без дополнительного ограничения на шаг по времени, но при тех же шагах дала несколько большую ошибку в максимальной норме и более сглаженные фронты, что связано с усиленным демпфированием высокочастотных мод.
306+
307+
Отдельно был исследован спектральный состав колебаний струны. С помощью БПФ построены спектры отклонения, скорости и локальной энергии по модам. Показано, что при заданном импульсе начальная энергия в основном сосредоточена в скорости, поэтому спектр смещения оказывается существенно ниже спектров скорости и энергии. Численные эксперименты с изменением ширины импульса $\delta$ продемонстрировали, что её уменьшение приводит к более широкому и высокочастотному спектру (возбуждается больше высоких мод), тогда как увеличение ширины импульса концентрирует энергию в низких модах и подавляет высокочастотные компоненты.
308+
309+
Наконец, варьирование параметра $\sigma$ в схеме с весами выявило ограниченность выбора: при больших $\sigma$ (например, 0.8) отрицательный вес среднего временного слоя приводит к сильным искажениям профиля и фактической потере физического смысла решения.
310+
311+
В целом выполненная работа показала, что для задач волнового типа важно не только соблюдение формальных условий устойчивости, но и разумный выбор численных параметров (шага, веса, ширины импульса) с учётом спектральных свойств решения. Реализованная модульная архитектура на C++ с точным решением, двумя конечно‑разностными схемами, решателем трёхдиагональных систем и блоком спектрального анализа позволяет легко проводить такие исследования и может служить основой для дальнейших экспериментов с более сложными начальными условиями и многомерными обобщениями.

0 commit comments

Comments
 (0)