Матэматычная мадэль сінхроннага генератара (Лабараторная работа № 10), страница 10

PSIf=Xf*If+Xad*Id                                                  183

PSIud=Xad*(Id+If)                                                  184

PSIuq=Xaq*Iq                                                       185

Mem=PSId*Iq-PSIq*Id                                                186

Mmex=Mem                                                           187

C        РАСЧЕТ НАЧАЛЬНЫХ УСЛОВИЙ Y(0) ДЛЯ СИСТЕМЫ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ

Y(1)=PSId                                                          200

Y(2)=PSIq                                                          201

Y(3)=PSIf                                                          202

Y(4)=PSIud                                                         203

Y(5)=PSIuq                                                         204

Y(6)=0.                                                            205

Y(7)=GAM                                                           206

DO 2 I=1,NUR                                                       210

2 Y0(I)=Y(I)                                                         211

C        РАСЧЕТ ПРОМЕЖУТОЧНЫХ ПЕРЕМЕННЫХ W1-W13, ИСПОЛЬЗУЕМЫХ В ПОДПРОГРАММЕ

C        "PRAV" ПРИ РАСЧЕТЕ ТОКОВ ПО ИЗВЕСТНЫМ ПОТОКОСЦЕПЛЕНИЯМ

W1=Xuq*Xq-Xaq*Xaq                                                  220

W2=Xad*Xad                                                         221

W3=Xq/W1                                                           222

W4=Xaq/W1                                                          223

W5=Xuq/W1                                                          224

W6=Xf*Xud*Xd+2.*Xad*Xad*Xad-W2*(Xud+Xf+Xd)                         225

W7=(Xf*Xd-W2)/W6                                                   226

W8=(Xf*Xad-W2)/W6                                                  227

W9=(Xad*Xd-W2)/W6                                                  228

W10=(Xd*Xud-W2)/W6                                                 229

W11=(Xud*Xad-W2)/W6                                                230

W12=(Xf*Xud-W2)/W6                                                 231

W13=120./180.*3.141593                                             232

C        ПОДГОТОВКА К РЕАЛИЗАЦИИ АЛГОРИТМА РУНГЕ-КУТТА

H1=OMEGA*H

TAU=OMEGA*T0                                                       240

T=T0                                                               241

NREZ=0                                                             242

C       ВЫДАЧА ЗАГЛАВИЯ ТАБЛИЦЫ РЕЗУЛЬТАОВ В ВЫХОДНОЙ ФАЙЛ С ИМЕНЕМ AAAA.REZ

WRITE(3,*)'           РЕЗУЛЬТАТЫ РАСЧЕТА'                          250

WRITE(3,*)'      T    PSId    PSIq    PSIf   PSIud   PSIuq    PSIa 251

*    PSIb    PSIc'                                                  252

WRITE(3,*)'             Id      Iq      If     Iud     Iuq      Ia 253

*      Ib      Ic'                                                  254

WRITE(3,*)'              S     GAM    Mmex     Mem'                255

WRITE(3,*)' НОРМАЛЬНЫЙ РЕЖИМ'                                      256

C НАЧАЛО ЦИКЛИЧЕСКОЙ ЧАСТИ РЕШЕНИЯ СИСТЕМЫ УРАВНЕНИЙ ГЕНЕРАТОРА

C        НАДО ЛИ ЗАНОСИТЬ ЗНАЧЕНИЯ ТЕКУЩЕГО ШАГА В ВЫХОДНЫЕ ФАЙЛЫ?

4 IF(T-H*REAL(NREZ*INTWR).LT.0.)GO TO 5                              260

C        РАСЧЕТ ПАРАМЕТРОВ, ВЫВОДИМЫХ В ВЫХОДНЫЕ ФАЙЛЫ, В ФАЗНЫХ

C        КООРДИНАТАХ А,В,С

Ia=Id*COS(GAM)+Iq*SIN(GAM)                                         270

Ib=Id*COS(GAM-W13)+Iq*SIN(GAM-W13)                                 271

Ic=Id*COS(GAM+W13)+Iq*SIN(GAM+W13)                                 272

PSIa=PSId*COS(GAM)+PSIq*SIN(GAM)                                   273

PSIb=PSId*COS(GAM-W13)+PSIq*SIN(GAM-W13)                           274

PSIc=PSId*COS(GAM+W13)+PSIq*SIN(GAM+W13)                           275

C        НАРАЩИВАНИЕ СЧЕТЧИКА ЗАПИСЕЙ, ВЫВОДИМЫХ В ВЫХОДНЫЕ ФАЙЛЫ

NREZ=NREZ+1                                                        280

C        ЗАНЕСЕНИЕ ИНФОРМАЦИИ ДЛЯ ВЫБРАННОГО УЗЛА В ВЫХОДНОЙ ФАЙЛ AAAA.REZ

WRITE(3,103)T,PSId,PSIq,PSIf,PSIud,PSIuq,PSIa,PSIb,PSIc            282

WRITE(3,104)Id,Iq,If,Iud,Iuq,Ia,Ib,Ic                              283

WRITE(3,105)S,GAM,Mmex,Mem                                         284

C        ЗАНЕСЕНИЕ ИНФОРМАЦИИ ДЛЯ ВЫБРАННОГО УЗЛА В ВЫХОДНОЙ ФАЙЛ AAAA.GRF