Сборник описаний задач специального лабораторного практикума "Компьютерный эксперимент" темы 3, 4, 5. Будко В.Н. - 8 стр.

UptoLike

Составители: 

Рубрика: 

Алгоритм модели :
Ввод t0,r0 печать точки Z1,Z2 на экране
Ввод x0,y0 X:=Z1
X:=x0 y:=Z2
y:=y0 Vx:=Z3
печать точки x,y на графике экрана Vy:=Z4
Ввод «угол» идти «цикл»
Vx=V0*cos(угол) y t0=1000
Vy=V0*sin(угол) r0=10
X0=60
<цикл> C:=X0/X y0=170
qq:=0.5*q*q
Z1:=X+1.5*Vx
Z2:=y+1.5*Vy x
Z3:=Vx+(q+qq)*Vy
Z4:=Vy+(q-qq)*Vx
10 INPUT t0=; t0 : INPUT r0= ;r0
20 LET V0=2*пi*r0/t0
30 INPUT x0= ; x0 :INPUT y0= ;y0:LET x=x0 : LET y=y0 : PLOT x,y
40 INPUT угол= ;угол
50 LET Vx=V0*cos(угол) : LET Vy=V0*cos(угол)
60 LET C:=X0/X : LET q=2*пi*c/t0 : LET qq=0.5*q*q
70 LET Z1=x+1.5*Vx :LET Z2=y+1.5*Vy
80 LET Z3=Vx+(q+qq)*Vy
90 LET Z4=Vy+(q-qq)*Vx
100 PLOT Z1,Z2
110 LET x=Z1 : LET y=Z2 : LET Vx=Z3 : LET Vy=Z4
130 GOTO 60
ТЕМА4.3Нелинейная автоколебательная система на туннельном диоде.
I
R L
i I(U) U0/R
+ U A P
- C Im
I2 B
U
U1 U0 U2 U3
Если выбрать рабочую точку «р» на падающем участке характеристики (участок AB вольт-
амперной характеристики с отрицательным сопротивлением ), то в схеме возникнут автоколебания
и её работа будет описываться диф .уравнениями :
Вольт-амперную характеристику туннельного диода можно апроксимировать либо
экспонециальной зависимостью [2] :
0/0**2:0 trV
π
=
0/**2: tCq
π
=
2
π
−=угол
U 0
)4(
)(
,
0
C
UIi
dt
dU
L
URiU
dt
di
=
⋅−
=
)
5
(
)
1
)
)
exp(
)
(
+
=
u
D
u
u
A
u
I
β
α
2
,
26
,
1
,
,
/
1
,
I
:
m/
=
=
=
=
m
B
I
D
u
u
e
A
где
µ
ϕ
β
α
А лго р и тм мо д е ли :
             Вво д t0,r0                                         пе чать то чки Z1,Z2 на экр ане
              V 0 := 2 * π * r 0 / t 0
              Вво д x0,y0                                               X:=Z1
              X:=x0                                                     y:=Z2
              y:=y0                                                    Vx:=Z3
              пе чать то чки x,y на гр аф и ке экр ана                      Vy:=Z4
              Вво д «уго л»                                             и д ти «ци кл»
              Vx=V0*cos(у   го л)                                  y                                 t0=1000
              Vy=V0*sin(у  го л)                                                                     r0=10
                                                                                                    X0=60
<ци кл>        C:=X0/X                                                                                y0=170
               q := 2 * π * C / t 0                                                                           угол = − π
                                                                                                                            2
              qq:=0.5*q*q
              Z1:=X+1.5*Vx
              Z2:=y+1.5*Vy                                                                   x
              Z3:=Vx+(q+qq)*Vy
              Z4:=Vy+(q-qq)*Vx

 10 INPUT “t0=”; t0 : INPUT “r0=” ;r0
 20 LET V0=2*пi*r0/t0
 30 INPUT “x0=” ; x0 :INPUT “y0=” ;y0:LET x=x0 : LET y=y0 : PLOT x,y
 40 INPUT “у го л=” ;уго л
 50 LET Vx=V0*cos(у    го л) : LET Vy=V0*cos(у
                                             го л)
 60 LET C:=X0/X : LET q=2*пi*c/t0 : LET qq=0.5*q*q
 70 LET Z1=x+1.5*Vx :LET Z2=y+1.5*Vy
 80 LET Z3=Vx+(q+qq)*Vy
 90 LET Z4=Vy+(q-qq)*Vx
 100 PLOT Z1,Z2
 110 LET x=Z1 : LET y=Z2 : LET Vx=Z3 : LET Vy=Z4
 130 GOTO 60



                        ТЕ М А 4.3Н е ли не йная авто ко ле б ате льная си сте ма на ту
                                                                                      нне льно м д и о д е .
                                                                               I
                   R                    L
                               i                  I(U)                  U0/R
          +     U0                                                        U                 A       P

          -                                              C                         Im
                                                                                   I2                               B

                                                                                                                        U
                                                                                           U1           U0 U2 U3

    Е сли выб р ать р аб о чу  ю то чку«р » на пад ающ е м у    частке х ар акте р и сти ки (у часто к AB во льт-
ампе р но й х ар акте р и сти ки с о тр и цате льным со пр о ти вле ни е м ), то в сх е ме во зни кнут авто ко ле б ани я
и е ёр аб о та б уд е т о пи сываться д и ф .у р авне ни ями :

     di   U 0−i⋅R −U                     dU    i − I (U )
        =                           ,        =                                      (4)
     dt        L                          dt       C

  Во льт-ампе р ну           ю х ар акте р и сти куту  нне льно го д и о д а мо ж но апр о кси ми р о вать ли б о
экспо не ци ально й зави си мо стью [2] :
                                                                  1
  I (еu ): A= =A e⋅ uI m/⋅ exp(
  гд                        u ,α−α=⋅ 1u/) u+ ,DD⋅ (exp(
                                                   = I ,ββ⋅ u=) − 1) , ϕ ≈ 26    ( 5 )µ B , m ≈ 2