Майер роберт Валерьевич основы компьютерного моделирования: Учебное пособие. Глазов: ггпи, 2005. 25 с



Download 1,6 Mb.
Pdf ko'rish
bet9/23
Sana22.04.2022
Hajmi1,6 Mb.
#573218
TuriУчебное пособие
1   ...   5   6   7   8   9   10   11   12   ...   23
Bog'liq
А 3 МАЙЕР Роберт Валерьевич

3. Алгоритм.
Пусть 
в момент времени материальная точка имеет 
координаты 
и 
проекции скорости 
Запишем 
второй закон 
Ньютона:
Отсюда 
следует, 
что 
проекции 
ускорения 
точки 
в 
момент 
времени 
равны:
Определив координаты и проекции скорости точки в момент 
времени 
, можно повторить процедуру вычисления требуемое 
количество раз и построить траекторию движения точки. 
Построим алгоритм модели. 
1. Задают массу материальной точки 
, коэффициент вязкости , начальные 
координаты 
и 
проекции 
скорости 
силовое 
поле 
, а также шаг по времени 
2. Начало 
цикла 
по Дают 
приращение 
по 
времени: 
переменной присваивают значение 


3. Определяют ускорение, скорость и координату тела в следующий момент 
времени: 
4. Результаты вычислений 
выводят на экран в числовом 
виде либо строят соответствующие точки на координатной плоскости. 
5. Возвращение к операции 2. Если цикл по закончился, -- выход из цикла. 
4. Компьютерная программа.
Предлагаемая компьютерная программа 
позволяет изучить движение материальной точки в различных силовых полях 
с учетом действующей на точку силы трения. 
program PROGRAMMA3; 
uses crt, graph; 
var v, B, q, F, Fx, Fy : real; 
r, x, y, vx,vy,ax,ay : real; Gd, Gm, i: integer; 
const M=500; mm=100; dt=0.005; rr=0.1; k=2; 
Begin 
Gd:= Detect; InitGraph(Gd, Gm, 'c:\bp\bgi'); 
if GraphResult <> grOk then Halt(1); 
line(320,240,640,240); line(320,240,320,0); circle(320,240,5); 
x:=100; y:=120; vx:=1; vy:=-2; 
Repeat 
begin 
{--Заданние силового поля--} 
(* Fy:=3; Fx:=0; *) 
(* Fx:=-k*x; Fy:=-k*y; *) 
(* r:=sqrt(x*x+y*y); F:=M*mm/(r*r); 
Fx:=-F*x/r; Fy:=-F*y/r; *) 
B:=2; q:=1; F:=B*v*q; v:=sqrt(vx*vx+vy*vy); 
Fx:=F*vy/v; Fy:=-F*vx/v; 
(* B:=2; q:=1; F:=B*v*q; v:=sqrt(vx*vx+vy*vy); 
Fx:=F*vy/v; Fy:=-0.5-F*vx/v; *) 
{--Расчет скоростей и ускорений--} 
ax:=(Fx-rr*vx)/mm; ay:=(Fy-rr*vy)/mm; 
vx:=vx+ax*dt; vy:=vy+ay*dt; x:=x+vx*dt; y:=y+vy*dt; 
circle(round(x)+320,240-round(y),2); setcolor(12); 


circle(round(x)+320,240-round(y),1); setcolor(15); 
end; 
until KeyPressed; 
CloseGraph; 
END. 

Download 1,6 Mb.

Do'stlaringiz bilan baham:
1   ...   5   6   7   8   9   10   11   12   ...   23




Ma'lumotlar bazasi mualliflik huquqi bilan himoyalangan ©hozir.org 2024
ma'muriyatiga murojaat qiling

kiriting | ro'yxatdan o'tish
    Bosh sahifa
юртда тантана
Боғда битган
Бугун юртда
Эшитганлар жилманглар
Эшитмадим деманглар
битган бодомлар
Yangiariq tumani
qitish marakazi
Raqamli texnologiyalar
ilishida muhokamadan
tasdiqqa tavsiya
tavsiya etilgan
iqtisodiyot kafedrasi
steiermarkischen landesregierung
asarlaringizni yuboring
o'zingizning asarlaringizni
Iltimos faqat
faqat o'zingizning
steierm rkischen
landesregierung fachabteilung
rkischen landesregierung
hamshira loyihasi
loyihasi mavsum
faolyatining oqibatlari
asosiy adabiyotlar
fakulteti ahborot
ahborot havfsizligi
havfsizligi kafedrasi
fanidan bo’yicha
fakulteti iqtisodiyot
boshqaruv fakulteti
chiqarishda boshqaruv
ishlab chiqarishda
iqtisodiyot fakultet
multiservis tarmoqlari
fanidan asosiy
Uzbek fanidan
mavzulari potok
asosidagi multiservis
'aliyyil a'ziym
billahil 'aliyyil
illaa billahil
quvvata illaa
falah' deganida
Kompyuter savodxonligi
bo’yicha mustaqil
'alal falah'
Hayya 'alal
'alas soloh
Hayya 'alas
mavsum boyicha


yuklab olish