Міністерство освіти і науки, молоді та спорту України

МІНІСТЕРСТВО ОСВІТИ І НАУКИ, МОЛОДІ ТА СПОРТУ УКРАЇНИ

НАЦІОНАЛЬНІЙ ТЕХНІЧНИЙ УНІВЕРСИТЕТ УКРАЇНИ

«КИЇВСЬКИЙ ПОЛІТЕХНІЧНИЙ ІНСТИТУТ»

ВИДАВНИЧО–ПОЛІГРАФІЧНИЙ ІНСТИТУТ

Кафедра репрографії

Практична робота №2

З дисципліни: Проблемно-орієнтовані засоби управління

На тему Перевірка розподілу вимірювань (у відповідності з законом Гаусса та корекція дефектного значення

Виконала:

Студентка групи СРп–81м

Шалашинська В. Г.

Перевірив:

Д-р техн.наук, професор

Морфлюк В.Ф.

Київ 2012

Мета роботи: за допомогою програмної мови ПОЯ проаналізувати, чи нормально розподілені значення вимірювань по Гаусу, та виконати корекцію дефектного значення.

Текст програми:

МОДУЛЬ ВИЗНАЧЕННЯ СЕРЕДНЬОГО ТА СТАНДАРТНОГО ВIДХИЛЕННЯ

НАЧ " ВИЗНАЧЕННЯ XCP ТА SIG";

ПРИ M=0;

99 ПРИ L=1;

10 ПРИ I=1;

ПРИ S=0;

1 ПРИ S=S+X(I);

ПРИ I=I+1;

ЕСЛ (I<=N) ПЕР 1;

ПРИ XCP=S/N;

ПРИ I=1;

ПРИ SI=0;

2 ПРИ XI=X(I);

ПРИ SI=SI+(XI-XCP)^2;

ПРИ I=I+1;

ЕСЛ (I<=N) ПЕР 2;

ПРИ N1=N-1;

ПРИ SIG1=SI/N1;

ПРИ GN=SIG1;

ПРИ G=0;

91 ПРИ G=G+0.01;

ПРИ G1=G*G;

ЕСЛ (G1>GN) ПЕР 90;

ПЕР 91;

90 ПРИ SIG=G;

ПРИ G1=XCP-SIG;

ПРИ G2=XCP+SIG;

ПРИ Q(1)=0;

ПРИ Q(2)=0;

ПРИ Q(3)=0;

ПРИ Q(4)=0;

ПРИ P(1)=0.16;

ПРИ P(2)=0.34;

ПРИ P(3)=0.34;

ПРИ P(4)=0.16;

ПРИ I=1;

4 ЕСЛ (X(I)<=G1) ПЕР 16;

ЕСЛ (X(I)<=XCP) ПЕР 17;

ЕСЛ (X(I)<=G2) ПЕР 18;

ПЕР 19;

16 ПРИ Q(1)=Q(1)+1;

ПЕР 20;

17 ПРИ Q(2)=Q(2)+1;

ПЕР 20;

18 ПРИ Q(3)=Q(3)+1;

ПЕР 20;

19 ПРИ Q(4)=Q(4)+1;

20 ПРИ I=I+1;

ЕСЛ (I<=N) ПЕР 4;

ПРИ K=1;

ПРИ SIK=0;

15 ПРИ E(K)=N*P(K);

ПРИ SIK=SIK+(Q(K)-E(K))^2/E(K);

ПРИ K=K+1;

ЕСЛ (K<=4) ПЕР 15;

ЕСЛ (SIK<=4) ПЕР 9;

ВЫВ " ВИМIРЮВАННЯ РОЗПОДIЛЕНI НЕ НОРМАЛЬНО";

ПРИ I=1;

ПРИ XI=X(I);

ПРИ PA=XI-XCP;

ЕСЛ (PA>=0) ПЕР 50;

ПРИ PA=-PA;

50 ПРИ X1=PA;

ПРИ J=I;

ПРИ SUS=X(J);

ПРИ I=2;

6 ПРИ XI=X(I);

ПРИ RA=XI-XCP;

ЕСЛ (RA>=0) ПЕР 60;

ПРИ RA=-RA;

60 ЕСЛ (RA>X1) ПЕР 5;

ПЕР 14;

5 ПРИ J=I;

ПРИ X1=RA;

ПРИ SUS=X(J);

14 ПРИ I=I+1;

ЕСЛ (I<=N) ПЕР 6;

ПРИ SA=SUS-XCP;

ЕСЛ (SA>=0) ПЕР 70;

ПРИ SA=-SA;

70 ПРИ T=SA/SIG;

ПРИ P=0.031*T^3-0.3007*T^2+0.963*T-0.014;

ЕСЛ (P<1) ПЕР 80;

ПРИ P=0.9999;

80 ПРИ T1=(1-P)*N;

ЕСЛ (T1<0.5) ПЕР 8;

ПРИ I=1;

8 ПРИ X(J)=XCP; корректирует дефектное значение (присваивает ему среднее значение(прошлое)

ВЫВ J,X(J);

ПРИ L=L+1;

ЕСЛ (L<=10) ПЕР 10;

ПАУ "ПРИЙМАЙТЕ РIШЕННЯ - Є ПОМИЛКА?????";

ПЕР 24;

9 ПРИ I=1;

ВЫВ XCP,SIG;

200 ПРИ W=X(I);

ЗАП W>M.IZM;

ПРИ I=I+1;

ЕСЛ (I<=N) ПЕР 200;

ВЫВ " ВИМIРЮВАННЯ РОЗПОДIЛЕНI НОРМАЛЬНО";

ПАУ " ПАУЗА ";

24 КОН " ПРОЦЕС ОБРОБКИ ЗАКIНЧЕНО";

Перевірка роботи програми:

1. При ненормальному розподіленні невеликих значень (з корекцією дефектного значення):

Відкоректовані дані, що записалися в файл А1.IZM:

.13000E+01

.18300E+01

.15347E+01

.13600E+01

.12800E+01

.14700E+01

.15100E+01

.14900E+01

.14800E+01

.15200E+01

.13600E+01

.16400E+01

.17800E+01

.14400E+01

.11600E+01

.14800E+01

2. При ненормальному розподіленні великих значень (з корекцією дефектного значення).

Введення значень нашого масиву:

НАЧ " СТАТИСТИЧНЕ ВИМIРЮВАННЯ ПАРАМЕТРIВ ";

ПРИ N=17;

ПРИ X(1)=72.30;

ПРИ X(2)=71.83;

ПРИ X(3)=71.55;

ПРИ X(4)=71.36;

ПРИ X(5)=71.28;

ПРИ X(6)=71.47;

ПРИ X(7)=71.51;

ПРИ X(8)=71.49;

ПРИ X(9)=71.48;

ПРИ X(10)=71.52;

ПРИ X(11)=71.36;

ПРИ X(12)=71.64;

ПРИ X(13)=71.78;

ПРИ X(14)=71.44;

ПРИ X(15)=71.16;

ПРИ X(16)=71.48;

ПРИ X(17)=71.44;

ВЫВ " ВИМIРЮВАННЯ ЗАКIНЧЕНО";

КОН " ПЕРЕХIД ДО ОБРОБКИ ТА ВИЗНАЧЕННЯ";

Результат роботи програми:

Висновок: в ході виконання даної роботи з використанням програмної мови ПОЯ було проаналізовано розподіл виміряних значень та проведено корекцію дефектного значення

Дата добавления: 2015-08-28