Фафурин В.А., Терюшов И.Н. Автоматизация технологических процессов и производств

COLOR 7

PRINT

PRINT "Решение определяем на интервале ( В, С )"

PRINT "Используем метод половинного деления"

PRINT "Точность метода равна N"

PRINT

PRINT "6. В - левая граница отрезка: ";: COLOR 11: INPUT "",

В

COLOR 7

PRINT "2. С - правая граница отрезка: ";: COLOR 11: INPUT "", С

COLOR 7

PRINT "3. N - точность расчета:

“

;: COLOR 11: INPUT"", N

PRINT

30 GOTO 80

40 LET P(1) = EXP(-A*T*L)

42 LET P(2) = EXP(-A * L)

44 LET P(3) = EXP(A * L)

46 LET P(4) = EXP(-A * M * T * L)

48 LET P(5) = EXP(-A * M * L)

50 LET P(6) = EXP(A * M * L)

52 LET Q(1) = (1 - P(2)) * (P(5) - 1) * (P(6) - 1)

54 LET Q(2) = 2 * (P(2) - 1) * (P(3) - 1) * (P(5) - 1) * (P(6) - 1)

56 LET Q(3) = К * (P(2) – 1) * (P(3) - 1) * (1 - P(2))

58 LET P(7) = 2 * T

^

2 * A * L * P(1) + 4 * T * P(1)

60 LET P(8) = T

^

2 * P(2) + T

^

2 * P{3) - 4 * T

62 LET P(9) = 2 * T

^

2 * A * M * L * P(4) - 4 * T

64 LET Q(4) = T

^

2 * P(5) + T

^

2 * P(6) + 4 * T * P(4)

68 LET F(1) = Q(1)* (P(7) - P(8))

70 LET F(2) = Q(2) * (T * A * L * P(1) - 1 * P(1)) * T

^

2

72 LET F(3) = Q(3) * (P(9) - Q(4))

74 DEF FNA (T) = F(1) - F(2) + F(3)

76 RETURN

80 LET T = В

90 GOSUB 40

100 LET F1 = FNA(T)

110 LET T = С

120 GOSUB 40

130 LET F2 = FNA(T)

140 LET D = (C + B)/2

145 LET T = D

150 GOSUB 40

160 LET F3 = FNA(T)

170 IF ABS(F3) <= N THEN 230

180 IF (SGN(F1) - SGN(F2)) = 0 THEN 210

190 LET C = D

200 GOTO 110

210 LET B = D

220 GOTO 80

230 GOTO 250

250 LET T = D

260 LET S(1) = 2 + T * P(2) - T * P(3) - 2 * P(1)

270 LET S(2) = T

^

2 * (P(2) - 1) * (P(3) - 1)

280 LET S(3) = 2 * (1 -P(1))

290 LET S(4) = T * (1 - P(2))

300 LET S(5) = К * (2 + T * P(5) - T * P(6) - 2 * P(4))

310 LET S(6) = T

^

2 * (P(5) - 1) * (P(6) - 1)

320 LET SO = S * (1 + S(1) / S(2) - S(3) / S(4) + S(5) / S(6))

330 PRINT "Интервал усреднения N(ext)=", T

335 PRINT "Дисперсия фильтра S(ext)=", S0

340 END

'Center:

' Center text on the given row.

SUB center (row%, text$)

LOCATE row %, 41 - LEN(text$) / 2

PRINT text$;

END SUB

"АСУТП KORF" - расчет статистических характеристик

стационарных случайных процессов: математического ожидания,

дисперсии и корреляционной функции

program korf;

uses crt;

const

n=223;

n1=20;

var

f: text;

st: string;

K,R,S,M,D:

REAL;

STROC, I, J, G: INTEGER;

X, T: ARRAY [1..400] OF REAL;

begin

CLRSCR;

writeln (‘СТАТИСТИЧЕСКИЕ ХАРАКТЕРИСТИКИ

СЛУЧАЙНОГО ПРОЦЕССА');

write ('Введите имя файла с данными: ');

readln(st);

assigN(f, st);

reset(f);

S: =0;

FOR I: =l TO N DO

begin

readln (f, X[I]);

S: =S+X[I];

end;

M: =S/n;

writeln (' МАТЕМАТИЧЕСКОЕ ОЖИДАНИЕ M=' , m:5:8);

D: =0;

FOR I: =l TO N DO

begin

X[I]: = X[I] - M;

D: = D + X[I] * X[I] / (N - 1) ;

end;

writeln(' ДИСПЕРСИЯ D=' , d:5:8);

readln;

writeln(' ЗНАЧЕНИЯ КОРРЕЛЯЦИОННОЙ ФУНКЦИИ ');

writeln;

writeln (' ЧИСЛО ШАГОВ ДИСКРЕ- ' ,'КОРРЕЛЯЦИОННАЯ ', ' НОРМИРОВАННАЯ');

writeln ('ТИЗАЦИИ МЕЖДУ СЕЧЕ- ', ' ФУНКЦИЯ СЛ. ', ' КОРРЕЛЯЦИОННАЯ');

writeln ('НИЯМИ СЛУЧАЙНОГО ', ' ПРОЦЕССА ', ' ФУНКЦИЯ');

writeln (' ПРОЦЕССА ');

writeln;

writeln (' J ', ' К ' ,

!

R ' ) ;

writeln (' _______________________________________________’);

writeln (' О ', D:5:8, ' ',

' 1 ');

STROC: =15;

FOR J: =l TO 20 DO

begin

G: =N-J;

K: =0;

FOR I: =l TO G DO

K: =K+X[I]*X[I+J];

T[J]: =K/G;

R: =T[J]/D;

writeln (' ', J, ‘ ‘,T[J]:5:8, ' ',

' R:5:8 ');

end;

IF STROC =24 THEN

REPEAT

UNTIL KEYPRESSED;

STROC: =STROC +1;

CLOSE(F);

READLN;

END.

"АСУТП KORF" - расчет статистических характеристик

стационарных случайных процессов: математического ожидания,

дисперсии и корреляционной функции

25 PRINT "СТАТИСТИЧЕСКИЕ ХАРАКТЕРИСТИКИ

СЛУЧАЙНОГО ПРОЦЕССА "

30 PRINT

40 READ n, n1

50 DIM x(310), t (50)

60 s = 0

70 FOR i = 1 TO n

80 READ x(i)

90 s = s + x(i)

100 NEXT i

110 M = s / n

120 PRINT " МАТЕМАТИЧЕСКОЕ ОЖИДАНИЕ М="; М

140 d = 0

150 FOR i = 1 TO n

160 x(i) = x(i) - M

170 d = d + x(i) * x(i) / (n - 1)

180 NEXT i

190 PRINT " ДИСПЕРСИЯ D="; d

210 PRINT

220 PRINT " ЗНАЧЕНИЯ КОРРЕЛЯЦИОННОЙ ФУНКЦИИ "

222 PRINT

224 PRINT"ЧИСЛО ШАГОВ ДИСКРЕ-","КОРРЕЛЯЦИОННАЯ","

НОРМИРОВАННАЯ"

225 PRINT "ТИЗАЦИИ МЕЖДУ СЕЧЕ-", " ФУНКЦИЯ СЛ. ", " КОРРЕЛЯЦИОННАЯ"

226 PRINT "НИЯМИ СЛУЧАЙНОГО ", " ПРОЦЕССА ", " ФУНКЦИЯ"

227 PRINT " ПРОЦЕССА "

228 PRINT

230 PRINT " J ", " К ", " R "

232 PRINT "______________________________________________"

235 PRINT " 0 ", ; " "; d, , " "; " 1 "

240 FOR J = 0 TO 20

250 G = n - J

260 k = 0

270 FOR i = 1 TO G

280 k = k + x(i) * x(i + J)

290 NEXT i

300 t(J) = k / G

310 R = t(J) / d

320 PRINT " "; USING "##"; J + 1;

PRINT " ";

PRINT USING "################# . ######"; t(J); R

330 NEXT J

340 GOTO 400

350 DATA 170, 20

352 DATA 160, 160.5, 161.6, 161.1, 160.9, 160.6, 159.5, 159.8, 160.5,

160.2, 159.8

354 DATA 160.6, 160.7, 160.8, 161.8, 163.1, 163.9, 163.5, 161.2, 161.4,

161, 160.8

356 DATA 160.5, 160.6, 160.2, 161, 161.9, 160.7, 159.5, 159.3, 158.3,

156.3, 156

358 DATA 157.1, 157.3, 158.8, 160, 161.1, 160.7, 157.7, 157.2, 158.9, 160,

161

360 DATA 158.7, 158.6, 159.8, 160.9, 159.6, 158.4, 157.8, 159.4, 160.9,

161.7, 160

362 DATA 157.7, 157, 157.9, 158.3, 158.7, 159, 159.4, 160.8, 161, 162,

161.5, 159.8

364 DATA 159.3, 159.4, 160, 158.9, 158.5, 158, 159, 159.6, 160.6,

160.9, 161, 161.7

366 DATA 160.8, 160.4, 159.9, 159.6, 159, 158, 157.9, 158, 158.6, 159,

159.6, 160

368 DATA 161, 160.5, 160.35, 159.2, 158.8, 158, 159.2, 160, 160, 160.4,

162.3, 163

370 DATA 161.7, 161.1, 160.5, 159.6, 159, 158.8, 158, 158.2, 158.5,

159, 159.6, 160

372 DATA 159.7, 159, 158.5, 159.2, 160, 161, 162, 162.5, 162, 161, 160,

159.3, 158.8

374 DATA 159, 159.7, 161, 161.8, 160.5, 160, 159, 158, 159.1, 159.9,

160, 160.9, 161

37 6 DATA 160.2, 159.5, 158.2, 159, 159.6, 161.5, 162, 161, 160, 159,

158.5, 158

378 DATA 160.7, 161.6, 161.9, 161.2, 160.8, 160.2, 159.5, 159, 158.4,

159.6, 160

399 PRINT

400 REM построение

PRINT

500 REM подготовка к построению

REM поиск мин фун Т( )

min = d: max = d

FOR i = 0 TO n1

IF min > t (i) THEN min = t(i): positMin = i

NEXT i

REM подготовка коэфициэнтов

maxx = positMin: IF maxx < nl THEN maxx = maxx + 1

minx = 0

maxy = max

miny = min

REM подготовка массива

maxx = maxx + 1

FOR i = maxx TO 1 STEP -1

t(i) = t(i - 1)

NEXT i

t(0) = d

REM расчет массштаба

mashX = 630 / (maxx - minx)

mashy = 460 / (maxy - miny)

REM поиск касательной

dy = t (0) - t (1)

i = 0

DO

dyold = dy

i = i + 1

dy = t(i) - t(i + 1)

LOOP WHILE dyold < dy

dy = dyold

REM расчет касательной

b = i * dy + t (i)

k = -dy

REM вывод уравнения касательной

PRINT "Уравнение касательной: Y="; k; "*Х+"; Ь; "т. ( "; i; ","; t(i); ")"

PRINT

PRINT "ГРАФИК КОРРЕЛЯЦИОННОЙ ФУНКЦИИ...";

G$ = ""

WHILE G$ =""

G$ = INKEY$

WEND

REM подготовка экрана

CLS

SCREEN 12

REM координатные оси

tmpy = (maxy - 0) * mashy

LINE (0, tmpy) - (640, tmpy)

LINE (0

f

0) - (0, 480)

REM построение

pointsx = (0 - minx) * mashX

pointsy = (maxy - t(0)) * mashy

FOR i = 1 TO maxx

pointx = (i - minx) * mashX

pointy = (maxy - t(i)) * mashy

LINE (pointsx, pointsy) - (pointx, pointy)

PSET (pointsx, pointsy), 4

pointsx = pointx: pointsy = pointy

NEXT i

REM построение касательной

i = 0

pointsx = (i - minx) * mashX

pointsy = (maxy - b) * mashy

FOR i = 1 TO maxx STEP .25

pointx = (i - minx) * mashX

pointy = (maxy - k * i - b) * mashy

IF ((pointy > 1) AND (pointy < 480)) AND ((pointsy > 0) AND (pointsy<

<479)) THEN LINE (pointsx, pointsy) - (pointx, pointy)

pointsx = pointx: pointsy = pointy

NEXT i

END

"DIGNOSTICA " - Контроль достоверности исходной

информации и диагностика отказов информационно-

измерительных каналов

DIM x(4), dx(4), sigma(4), p(4)

COLOR 7, 0

CLS

PRINT " Контроль достоверности исходной информации и"

PRINT "диагностика отказов информационно-измерительных

каналов"

PRINT

PRINT "Ввод исходной информации для расчета"

PRINT "___________________________________________"

PRINT "Численные значения кол-ва информации, поступающей с

каждого ИИК:"

FOR i = 1 ТО 4

PRINT "X("; i;") =";

INPUT "", x(i)

NEXT i

INPUT "Допустимая погрешность выполнения уравнения связи:", I

PRINT : PRINT "Допустимые погрешности работы ИИК:"

FOR i = 1 ТО 4

PRINT "dX("; i;") =";

INPUT "", dx(i)

NEXT i

PRINT: PRINT "Среднеквадратичные погрешности измерения:"

FOR i = 1 TO 4

PRINT "sigma("; i;") =";

INPUT "", sigma(i)

NEXTi

PRINT " **********************Расчет*********************** "

SX = x(1) + x(2) + x(3)-x(4)

PRINT USING "Погрешность выполнения уравнения связи =

### . #####"; SX

IF ABS(SX)<= I THEN PRINT "Все результаты измерения

достоверные": END

PRINT "Среди результатов измерения имеются недостоверные"

k = 0

FOR i = 1 TO 4

k = k + 1 / (sigma(i)

^

2)

NEXTi

k= 1 /k

PRINT USING "K= ### . #####"; k

PRINT "Весовые коэффициенты:"

FOR i = 1 TO 4

p(i) = k / (sigma(i)

^

2)

PRINT USING "P(#)= ### . #####"; i; p(i)

NEXT i

Фафурин В.А., Терюшов И.Н. Автоматизация технологических процессов и производств

Подождите немного. Документ загружается.