Реализовать модель математического маятника с учетом трения в точке подвеса

@Maksimyshka · Регистрация: 10.04.2012

Студворк — интернет-сервис помощи студентам

Реализовать модель математического маятника с учетом трения в точке подвеса.

Кликните здесь для просмотра всего текста

dQ/dt = x
dx/dt = -2*n*x-w*w*sinQ

$https://www.cyberforum.ru/cgi-bin/latex.cgi?\begin{cases} \frac{dQ}{dt} = x \\ \frac{dx}{dt} = -2nx-w^2sinQ \end{cases}$
x- угловая скорость, Q- отклонение маятника, n- коэффициент трения, w- частота собственных малых колебаний
Составил вот такую программу, но результат, мягко говоря не очень. Может подскажете в чем проблема?

Pascal
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
program ex1;
var T: Text;
x, h, a, b, x1, Q, k1_1, k2_1, k3_1, k4_1, k1_2, k2_2, k3_2, k4_2: real;
function F(x, N, C: real; p: byte):real;
begin
case p of
 1: f:=x1;
 2: f:=-2*0.9*x1-2*2*sin(Q);
 end;
end;
begin
assign(T, 'D:\1.txt');
rewrite(T);
a:=1; b:=10;
h:=0.001; x1:=1; Q:=45;
x:=a;
while (x<=b) do
  begin
writeln(T, x:15:5, x1:15:5, Q:15:5);
k1_1:=h*F(x,x1,Q,1);
k1_2:=h*F(x,x1,Q,2);
k2_1:=h*F(x+h/2, x1+k1_1/2, Q+k1_1/2, 1);
k2_2:=h*F(x+h/2, x1+k1_2/2, Q+k1_2/2, 2);
k3_1:=h*F(x+h/2, x1+k2_1/2, Q+k1_1/2, 1);
k3_2:=h*F(x+h/2, x1+k2_2/2, Q+k1_2/2, 2);
k4_1:=h*F(x+h, x1+k3_1, Q+k3_1, 1);
k4_2:=h*F(x+h, x1+k3_2, Q+k3_2, 2);
x1:=x1+1/6*(k1_1+2*k2_1+2*k3_1+k4_1);
Q :=Q+1/6*(k1_2+2*k2_2+2*k3_2+k4_2);
x:=x+h;
  end;
  close(T);
end.

Добавлено через 20 часов 39 минут
Ну выручите, пожалуйста.

@Legolas · 22.11.2014, 16:58

результат, мягко говоря не очень.

что значит не очень?

P.S. ещё было бы хорошо, если б вы немного прокоментировали код, описали какие переменные за что отвечают....

@Maksimyshka · 22.11.2014, 17:11 **[ТС]**

маятник крутится все время по кругу, x1 - угловая скорость, Q - угол отклонения h- шаг

@Legolas · 23.11.2014, 13:13

ладно, что такое:
h - пишите что шаг, в чем он измеряется? (в секундах?)
x - ? (просто икс)
a, b - ?
k1_k и ему подобные и зачем их вычислять так много раз?
p -? (почему оно равно 1 или 2)
что за 1/6?
почему при вызове функции вы некоторые параметры делите на 2
почему вы в функцию передаете параметр x, а не используете его?

Pascal
1
результат, мягко говоря не очень.

Самое главное - это значит не очень?
не компилируется/зацикливается/выдает не правильные результаты (как проверяете?)

Новые блоги и статьи Все статьи Все блоги /
Загрузка PNG-файла с альфа-каналом с помощью библиотеки SDL3_image на Android 8Observer8 27.01.2026 Содержание блога SDL3_image - это библиотека для загрузки и работы с изображениями. Эта пошаговая инструкция покажет, как загрузить и вывести на экран смартфона картинку с альфа-каналом, то есть с. . .	влияние грибов на сукцессию anaschu 26.01.2026 Бифуркационные изменения массы гриба происходят тогда, когда мы уменьшаем массу компоста в 10 раз, а скорость прироста биомассы уменьшаем в три раза. Скорость прироста биомассы может уменьшаться за. . .	Воспроизведение звукового файла с помощью SDL3_mixer при касании экрана Android 8Observer8 26.01.2026 Содержание блога SDL3_mixer - это библиотека я для воспроизведения аудио. В отличие от инструкции по добавлению текста код по проигрыванию звука уже содержится в шаблоне примера. Нужно только. . .	Установка Android SDK, NDK, JDK, CMake и т.д. 8Observer8 25.01.2026 Содержание блога Перейдите по ссылке: https:/ / developer. android. com/ studio и в самом низу страницы кликните по архиву "commandlinetools-win-xxxxxx_latest. zip" Извлеките архив и вы увидите. . .
Вывод текста со шрифтом TTF на Android с помощью библиотеки SDL3_ttf 8Observer8 25.01.2026 Содержание блога Если у вас не установлены Android SDK, NDK, JDK, и т. д. то сделайте это по следующей инструкции: Установка Android SDK, NDK, JDK, CMake и т. д. Сборка примера Скачайте. . .	Использование SDL3-callbacks вместо функции main() на Android, Desktop и WebAssembly 8Observer8 24.01.2026 Содержание блога Если вы откроете примеры для начинающих на официальном репозитории SDL3 в папке: examples, то вы увидите, что все примеры используют следующие четыре обязательные функции, а. . .	моя боль iceja 24.01.2026 Выложила интерполяцию кубическими сплайнами www. iceja. net REST сервисы временно не работают, только через Web. Написала за 56 рабочих часов этот сайт с нуля. При помощи perplexity. ai PRO , при. . .	Модель сукцессии микоризы anaschu 24.01.2026 Решили писать научную статью с неким РОманом

@Maksimyshka 2 / 2 / 2 Регистрация: 10.04.2012 Сообщений: 147
	22.11.2014, 17:11 [ТС]
	маятник крутится все время по кругу, x1 - угловая скорость, Q - угол отклонения h- шаг 0