х елементів),
В· оператор циклу while - перевіряє значення виразу (вираз логічного типу), якщо воно дорівнює True (умова виконується), то виконуються інструкції, що знаходяться між begin і end (інструкції циклу). Потім знову перевіряється значення виразу умова, і так продовжується до тих пір, поки значення виразу умова не стане рівним False. Таким чином, після слова while записується умова виконання інструкцій циклу.
В· оператор виведення writeln
В· оператор введення readln.
На представленій нижче блок-схемі (Рис.1.2) описана послідовність дій, виконаних при написанні програми.
Як видно з блок-схеми на початку програма значення необхідні для розрахунку добротності контуру:
Tmax: = 2 * Pi * sqrt (L * c2); Tmin: = 2 * Pi * sqrt (L * c1);
h: = 0.1 * Tmin;: = 0;: = sqrt (l/C);
Далі, використовуючи отримані дані, програма вважає добротність і вибирає залежно від неї необхідний час протікання перехідного процесу.
Q: = p/R; Q <0.5 then time: = Tmax else time: = 6 * Tmax;
Блок-схема
В
Рис. 1.2 Блок-схема алгоритму рішення системи диференціальних рівнянь у Pascal
Далі програма переходить до основного циклу:
t <= time Do
i: = i + di;: = uc + duc;
ur: = i * r;: = u-uc-ur;: = (ui * r-uc) * h/l;: = i * h/c;: = (u- (i + ki1/2) * r-(uc + ku1/2)) * h/l;: = (i + ki1/2) * h/c;: = (u-(i + ki2/2) * r-(uc + ku2/2)) * h/l;: = (i + ki2/2) * h/c;: = (u-(i + ki3) * r-(uc + ku3)) * h /l;: = (i + ki3) * h/c;: = (ku1 +2 * ku2 +2 * ku3 + ku4)/6;: = (ki1 +2 * ki2 +2 * ki3 + ki4)/6 ; (f, t: 7: 3, '', ki1: 7: 3, '', ki2: 7: 3, '',: 7: 3, '', ki4: 7: 3, '', ku1: 7: 3, '', ku2: 7: 3, '',: 7: 3, '', ku4: 7: 3, '', i: 7: 3, '', uc: 7: 3, '' , ur: 7: 3, '', ul: 7: 3);: = t + h;
В кінці роботи програма закриває файл з результатом.
Після завершення роботи програми користувач може ознайомитися з результатами переглянувши файл "tabpas. PAS . Повним текстом програми можна ознайомитися в додатку.
Розрахунок перехідного процесу в табличному редакторі Excel методом Рунге-Кутта 4-го порядку
За формулою (14) обчислимо період коливань і крок зміни часу, рівний однієї десятої періоду коливань:
В
В
= 2 * ПІ () * КОРІНЬ (B3 * C3) = 0,1 * I3
де C3 - адреса комірки зі значенням індуктивності, B3 - адреса комірки зі...
