Цитата:
Сообщение от kot64rus
Блин, похоже я с математикой совсем раздружился. Я совсем не догоняю, как решать это =( Можно как нибудь на примере. Куда что подставить надо. Весь день сижу понять не могу(
|
Да нет вопросов! Пожалуйста.
Одно замечание.
Чтобы осциллятор начал реально осциллировать, исходную скорость следует задавать
отрицательной. В противном случае получится просто однонаправленное движение с выходом скорости на некую константу.
Итак:
Код:
Const
Dt=0.1;
V=-0.5;
Q=0.0;
Var
f:Text;
V_old,V_new,Q_old,Q_new,t:Real;
i:Word;
Begin
Assign(f,'D:\Result1.dat');
Rewrite(f);
t:=0;
Q_old:=Q;
V_old:=V;
Writeln(f,t:5:1,Q:15:3,V:15:3);
for i:=1 to 350 do
begin
t:=Dt*i;
Q_new:=Q_old+V_old*Dt;
V_new:=V_old+Exp(-Q_old)*(Exp(-Q_old)-1)*Dt;
Writeln(f,t:5:1,Q_new:15:3,V_new:15:3);
Q_old:=Q_new;
V_old:=V_new;
end;
Close(f);
End.
Естественно, имя файла и путь к нему поставьте свои.
Результат:
А вот как там "исследовать точки равновесия на устойчивость" - это я не в курсе.