Испытание на электромагнитную совместимость аппаратуры программно-технического комплекса (Лабораторная работа № 5), страница 3

           Рис. 5.2 Блок-схема обработки результатов измерений напряжения радиопомех   аппаратуры ПТК

Таблица 5.1.

Частота, МГц	Измерение напряжения радиопомех, дБ (мкВ)
	№ 1	№ 2	№ 3	№ 4
0.009	32	32	33	33
0.01	32	32	33	33
0.015	29	30	30	31
0.025	20	20	21	22
0.04	22	22	23	24
0.06	13	13	14	15
0.07	20	20	21	22
0.1	53	53	55	56
0.16	61	63	63	64
0.21	51	51	52	53
0.24	70	70	72	73
0.55	64	65	66	66
1	58	59	60	61
1.4	55	57	57	58
2	44	45	45	45
3.5	40	40	41	41
6.0	39	40	41	41
30	44	45	46	46
45	35	35	36	37
65	37	38	38	38

   3.2. Реализация алгоритма обработки впакете инженерного проектирования MatLab

В редакторе для m-файлов пакета MatLab создадим программу, соответствующую блок-схеме на рис. 5.2., и произвести запуск на исполнение с последующим анализом полученных результатов. Листинг программы приведен ниже.

%--------------------------------------------------------------------------

% Инициализация используемых в программе локальных переменных:

n = 4;              % Количество замеров;

k = 1.69;           % Коэффициент, связанный с числом испытанных устройств и  доверительной вероятностью;

L = 6;              % Коэффициент затухания кабеля;

X=[0.009    32  32  33  33

    0.01    32  32  33  33

    0.015   29  30  30  31

    0.025   20  20  21  22

    0.04    22  22  23  24

    0.06    13  13  14  15

    0.07    20  20  21  22

    0.1     53  53  55  56

    0.16    61  63  63  64

    0.21    51  51  52  53

    0.24    70  70  72  73

    0.55    64  65  66  66

    1       58  59  60  61

    1.4     55  57  57  58

    2       44  45  45  45

    3.5     40  40  41  41

    6.0     39  40  41  41

    30      44  45  46  46

    45      35  35  36  37

    65      37  38  38  38]%измерений напряжения радиопомех и частота

%--------------------------------------------------------------------------

% Учет коэффициента затухания кабеля

dimension=size(X);

 for i = 1:dimension(1)

            for j = 2:dimension(2)

            Xstroke(i,j-1)=X(i,j)+ L;

            end;        % Цикл расчётов для каждого значения ИРП на рассматриваемой частоте;

    end;        % Цикл расчётов для каждой из частоты;

    %---------------------------------------------

% Определение математического ожидания для каждой частоты

    for i = 1:dimension(1)

            Xline(i)=0;

            for j = 2:dimension(2)

            Xline(i)=Xline(i)+Xstroke(i,j-1);

            end;        % Цикл расчётов для каждого значения ИРП на рассматриваемой частоте;

            Xline(i)=Xline(i)/n;

    end;        % Цикл расчётов для каждой из частоты;    %---------------------------------------------

    %---------------------------------------------

    %

  % Определение с.к.о. для каждой частоты

    for i = 1:dimension(1)

            S(i)=0;

            for j = 2:dimension(2)

               S(i) =S(i) + (Xstroke(i,j-1) - Xline(i))^2;

            end;        % Цикл расчётов для каждого значения ИРП на рассматриваемой частоте;

            S(i) = sqrt(S(i)/(n-1));

    end;        % Цикл расчётов для каждой из частоты;    %---------------------------------------------

    %---------------------------------------------

 % Определение статистического значения радиопомех для каждой контрольной

 % частоты

    for i = 1:dimension(1)

            A(i)=Xline(i)+S(i)*k;

    end;        % Цикл расчётов для каждой из частоты;    %---------------------------------------------

        %---------------------------------------------

        % Проверка превышения или не превышения нормируемого значения;