«Программы к моделированию физических задач (Pascal)» - Практическая работа

Содержание

1. Траектория движения при столкновении упругих шаров.

2.3. Моделирование артиллерийской задачи

2.4. Движение подводной лодки

2.5. Движение тела с учетом сопротивления среды

2.7. Столкновение упругих шаров

2.8. Движение небесного тела в гравитационном поле

2.9. Движение материального тела в поле тяготения планеты

2.10. Полет сверхзвукового самолета

2.11. Одноступенчатая ракета

2.12. Многоступенчатая ракета

2.13. Стыковка космического корабля

2.15. Физический маятник

2.14. Кривошипно-шатунный механизм

2.16. Параметрический маятник

2.17. Маятник Фуко

2.19. Связанные маятники

2.20. Колебания пружинного маятника

2.22. Успокоители механических колебаний

2.23. Связанные осцилляторы

2.24. Решение задачи Ферми-Паста-Улама

2.27. Распространение волн на воде. Солитон.

2.28. Форма капли жидкости

2.29. Замерзание капли

2.30. Решение уравнения теплопроводности

3.1. Силовые линии электрического поля

3.2. Движение заряженных частиц в кулоновском поле

3.3. Движение заряженной частицы в магнитном поле

3.4. Пространственный осциллятор

3.5. Разделение изотопов

3.6. Возмущение орбиты электрона в однородном магнитном поле движущимся протоном

3.8. Движение заряженной частицы в скрещенных полях

3.9. Фигуры Лиссажу

3.10. Выпрямление с фильтрацией

3.12. Генерирование колебаний

3.13. Параметрон

3.15. Телеграфное уравнение

3.16. Спектральный анализ

3.17. Радуга

4.1. Дифракция микрочастиц на отверстии

4.2. Квантово-механическая модель атомов. Часть 1. Расчет сферической части атома водорода

4.3. Квантово-механическая модель атомов. Часть 2. Расчет ра-диальной части атома водорода

4.7. Усиление оптического излучения

4.8. Генерация лазерного излучения

---

Введение

2.3. Моделирование артиллерийской задачи

Программа 1.

program Pushka;

{$R-}

{$Q-}

const

k: Integer = 0;

dt : Extended = 0.01;

da : Real = Pi / 25;

E: Extended = 0.1;

g: Real = 9.81;

I: Integer = 0;

Flag: Boolean = False;

Beg: Boolean = False;

Path: string = 'C:\temp';

var

X0 : Extended;

Y0 : Extended;

X1, X2: Extended;

Y1, Y2: Extended;

a, V0, L: Real;

PathTemp: string;

Num: string;

First: Boolean;

F: Text;

begin

Writeln('задайте скорость и угол');

Readln(V0, L);

a := Pi / 20;

I := 1;

while a < PI / 2 do

begin

X0 := 0;

Y0 := 0;

First := True;

Beg := False;

Str(I, Num);

PathTemp := Path + Num + '.dat';

Writeln(PathTemp);

Readln;

Assign(F, PathTemp);

Rewrite(F);

X1 := V0 * Cos(a) * dt - k * dt;

Y1 := X1 * Sin(a) / Cos(a);

while not Beg do

begin

{ Writeln(F, X0: 4: 4, ' ', Y0: 4: 4);}

Writeln(X0: 4: 4, ' ', Y0: 4: 4);

X2 := (X1 * (2 + k * dt) - X0) / (1 + k * dt);

Y2 := (Y1 * (2 - k * dt) - Y0 - g * Sqr(dt)) / (1 - k * dt);

X0 := X1; Y0 := Y1;

X1 := X2; Y1 := Y2;

if not First then

begin

Beg := Y0 <= 0;

Flag := ((X0 < (L + E)) or (X0 > (L - E))) and Beg;

end

else First := False;

end;

if Flag then

begin

Writeln('a =', (a * 180 / Pi):3:2);

Beg := False;

First := True;

X0 := 0;

Y0 := 0;

X1 := V0 * Cos(a) * dt - k * dt;

Y1 := X1 * Sin(a) / Cos(a);

while not Beg do

begin

Writeln(F, X0: 4: 4, ' ', Y0: 4: 4);

X2 := (X1 * (2 + k * dt) - X0) / (1 + k * dt);

Y2 := (Y1 * (2 - k * dt) - Y0 - g * Sqr(dt)) / (1 - k * dt);

X0 := X1;

Y0 := Y1;

X1 := X2;

Y1 := Y2;

if not First then

begin

Beg := Y0 <= 0;

Flag := ((X0 < (L + E)) or (X0 > (L - E))) and Beg;

end

else First := False;

end;

{ Readln;}

end;

a := a + da;

Inc(I);

Close(F);

end;

end.

---

Выдержка из текста работы

2.7. Столкновение упругих шаров

Программа 1.

PROGRAM Shary;

VAR

Vcx,Vx1,Vcy,V2x,m2,V1y,V2y,Vx10,Vi,Vc,L,Vxy,X1,Y1,

X2,Y2,Vxi,V10z,Vxz,Vx10z,Vy10,Vy10z,Vyi,Vx20,Vx20z,Vy20,

Vy20z,V20z,Vy1z,Vy2,Vx1z,V10x,Vx2,Vx2z,Vy2z,m1,w:real;

begin

m1:=1; m2:=2; Vx1:=3; V2x:=4; V1y:=5; V2y:=6; X1:=7;

X2:=8; Y1:=9; Y2:=10; w:= ;

Vx10z:=Vx10*cos(w);

Vy10z:=Vy10*sin(w);

Vx20z:=Vx20*cos(w);

Vy20z:=Vy20*sin(w);

Vcx:=(m1*Vx1+m2*V2x)/(m1+m2);

Vcy:=(m1*V1y+m2*V2y)/(m1+m2);

Vx10:=Vx1-Vcx;

IF L=1 THEN

Vi:={=(Vc*L)/L}(Vxy*(X2-X1)+Vcy*(Y2-Y1))/sqrt(sqr(X2-X1)+sqr(Y2-

Y1));{Vi-Ї а ««Ґ«м­ п};

Vxi:={=Vi*(X2-X1)/sqrt(X2-X1)}((Vxy*(X2-X1)+Vcy*(Y2-Y1))/(sqr(X2-

X1)+sqr(Y2-Y1)))*(X2-X1);

Vyi:={=Vi*(Y2-Y1)/sqrt(Y2-Y1)}((Vxy*(X2-X1)+Vcy*(Y2-Y1))/(sqr(X2-

X1)+sqr(Y2-Y1)))*(Y2-Y1);

V10z:=Vx10-2*Vxi;

Vx1z:=Vx10z+Vxy;

Vx10:=Vx1 - Vcx;

Vy10:=Vx1 - Vcy;

Vx20:=Vx2 - Vcx;

Vy20:=Vx2 - Vcy;

Vx10z:=Vx10 - 2*Vxi;

Vy10z:=Vy10 - 2*Vyi;

Vx20z:=Vx20 - 2*Vxi;

Vx20z:=Vx20 + 2*Vxi;

Vx1z:= Vx10z + Vcx;

Vy1z:= Vy10z + Vcy;

Vx2z:= Vx20z + Vcx;

Vy2z:= Vy20z + Vcy;

WRITELN (Vx1z:1:2,' ',Vx2z:1:2,' ',Vy1z:1:2, ' ',Vy2z:1:2);

readln

end.

---

Заключение

4.8. Генерация лазерного излучения

Program Laser;

CONST t=20;x=2;tp=0.5;t0=0.5;

VAR i,j,nt,nx,ntp,nt0:integer;

N0,h,g,L,gf:real;

N1,N,a,e1:array[0.200] of real;

e,e0:array[0.2000] of real;

ain,ain1:text;

BEGIN

assign(ain,'c:\laser.DAT'); REWRITE(ain);

assign(ain1,'c:\laser0.DAT'); REWRITE(ain1);

h:=0.1;

nt:=500; nx:=200; ntp:=30; nt0:=100;

g:=0.001; gf:=0.02; L:=0.1;

for j:=0 to nx do

begin

N[j]:=0;

a[j]:=0.02;

end;

for i:=0 to nt do

begin

{ if i<5*nt0 then L:=2*exp(-((i-nt0)*(i-nt0))/(ntp*ntp)/2) else L:=0;}

{ if i<6*nt0 then L:=1 else L:=0;}

for j:=0 to nx-1 do

begin

N1[j]:=N[j]-h*a[j]*N[j]-h*g*N[j];

e1[j+1]:=a[j]+h*a[j]*N[j]-h*gf*a[j]+h*L*a[j];

end;

for j:=0 to nx do

begin

N[j]:=N1[j];

a[j]:=e1[j];

end;

e[i]:=a[nx];

e0[i]:=N[nx-1];

end;

for i:=0 to nt do

begin writeln(ain,i*h,' ',e[i]);

writeln(ain1,i*h,' ',e0[i]);

end;

close (ain);

close (ain1);

end.

---

Примечания

К работе прилагается все исходники. Есть приложения.