Лістинг basic-програми розрахунку центральних параметрів

CLS

1 REM Блок вводу вихідних даних

5 OPEN “KR.TXT” FOR OUTPUT AS #1

10 INPUT “ВВЕДІТЬ КІЛЬКІСТЬ ТОЧОК”; M

20 INPUT “ВВЕДІТЬ КРОК ДИСКРЕТИЗАЦІЇ СЕК (наприклад 10 мкс- 10_Е-60)”; Т

30 INPUT “ВВЕДІТЬ МЕРІЛО ПО ОСІ НАПРУГИ”; КО

40 DIM H (M):DIM SHO! (4):DIM SH! (4)

45 DIM 8! (4): W=0

50 INPUT “ВВЕДІТЬ КІЛЬКІСТЬ ТОЧОК АПРОКСИМАЦІЇ”; P

60 LM= (P+1) * (M-1) +1

65 DIM HPL (4, LM): DIM HP (P, M)

70 FOR I=1 TO M

80 INPUT “ВВЕДІТЬ ВИСОТУ ОСЦИЛОГРАММ У ЧЕРГОВІЙ ТОЧЦІ В ПОДІЛКАХ БЕЗ ЗНАКУ”; H (I)

85 PRINT “H (”;I;”)=”;H(I):IF W=1 GOTO 100

90 NEXT I

100 INPUT “ЧИ БУЛА ПОМИЛКА ПЫД ЧАС ЗВЕДЕННЯ ЗНАЧЕНЬ (ТАК=1, НІ=0) “; W

110 IF W=0 THEN GOTO 120

111 PRINT “Введить номер точки, яка треба виправити”

112 INPUT I: INPUT H (I)

113 GOTO 85

114 CLS

115 REM Розрахунок значень сигналу у точках апроксимації

120 L=1

130 FOR N=1 TO 4

140 HPL (N, 1)=H (1)^N

145 NEXT N

150 FOR I=2 TO M

160 FOR J=0 TO P

161 NEXT N

165 FOR N=1 TO 4

170 HP=H (I) – H (I)- H (I-1))* (P-J)/(P+1)

180 HPL (1, L)=HP

200 HPL (N, L) =HPL (1, L)^N

210 NEXT N

212 NEXT J

213 NEXT I

214 REM Розрахунок інтегральних величин

220 FOR N=1 TO 4

221 L=LM

224 SH (N) = (HPL (N, 1)+HPL (N, L))/2

226 NEXT N

230 FOR L=2 TO (LM-1)

231 FOR N=1 TO 4

240 SH (N) =SH (N) +HPL (N, L)

245 NEXT N

250 NEXT L

251 REM Виведення таблиці на екран

260 PRINT “I”, “H(I)^2”, “H (I)^3”, “H (I)^4

265 C=0

270 FOR L=1 TO LM STEP (P+1)

272 PRINT (L-C*P),

275 PRINT USING “######.##^^^^”; HPL (1, L), HPL (2, L), HPL (3, L), HPL (4, L)

282 INPUT Z

283 C=C+1

285 NEXT L

290 PRINT “ = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =”_________

300 FOR N=1 то 4

315 NEXT N

320 PRINT “S(N)”,

321 PRINT USING “######.##^^^^”; S(1), S(2), S(3), S(4)

325 REM Розрахунок та виведення на друкування узагальнених параметрів еквівалентного П-імпульсу

330 rem

340 UOB= S(2) / S(1) :TOB= ^2/ S(2)

350 F1= S(1) * S(3) / S(2) ^2

360 F2= S(1) ^2 *(S(4))/((S(2))^3

361 PRINT USING “UOB= ##.###^^^^”; UOB

362 PRINT USING “TOB= ##.###^^^^”; TOB

363 PRINT USING “F1= ##.###^^^^”; F1

364 PRINT USING “F2= ##.###^^^^”; F2

365 GOTO 398

366 REM Аналіз узагальнених параметрів та вибір моделі

367 IF F1 1 GOTO 396 ELSE GOTO 368

368 IF F1 =1.125 GOTO 396 ELSE GOTO 368

376 REM Розрахунок узагальнених параметрів трапецеїдальної моделі

377 PRINT «узагальнені параметри трапецеїдальної моделі»

378 UM=UOB* (3-SQR(9-8*F1))/2:PRINT USING “ UM=##.###^^^^”; UM

379 TI=TOB* (8/(3-SQR (9-8*F1;)-12/((3-SQR (9-8*F1))^2)):PRINT USING “ TI=##.###^^^^’;TI

380 TA=3*(1-SQR(9-8*F1))/(3-2*SQR(9-8*F1)):PRINT USING “ TA=##.####”; TA

394 GOTO 396

395 PRINT “Для уточнення значення амплітуди використовуйте графік (рис. 2 МВ) “

396 PRINT “Чи бажаєте повторити розрахунки? Так -1”: INPUT RP

397 IF RP=1 GOTO 100

398 PRINT “Чи бажаєте надрукувати розрахунки? Так -1, Ні – 0: INPUT PR

399 IF PR=1 GOTO 400 ELSE GOTO 424

400 REM Вивід таблиці у файл кр тхт для друку

401 PRINT #1, “I”, “H(I)”, “H(I)^2”, “H(I)^3”, “H(I)^4

402 C=0

403 FOR L=3 TO LM STEP (P+1)

404 PRINT #1, (L-C*P),

405 PRINT #1, USING “#####.##^^^^”; HPL (1, L), HPL (2, L), HPL (3, L), HPL (4, L)

406 C=C+1

407 NEXT L

408 PRINT #1, “====================================================”

409 FOR N=1 TO 4

410 S (N) =SH (N) * (KO^N) *T/(P+1)

411 NEXT N

412 PRINT #1, S (N)”,

413 PRINT #1, USING “#####.##^^^^”; S(1), S(2), S(3), S(4)

414 PRINT # 1, “Параметри еквівалентного П-імпульсу”

415 PRINT #1, USING “ UOБ=##.###^^^^”; UOB

416 PRINT #1, USING “ TOБ=##.###^^^^”; TOB

417 PRINT #1, USING “ F1=##.###^^^^”; F1

418 PRINT #1, USING “ F2=##.###^^^^”; F2

419 GOTO 424

420 PRINT #1, “ Узагальнені параметри трапецеїдальної моделі”