% ЛР №16. НЕПАРАМЕТРИЧЕСКИЙ СПЕКТРАЛЬНЫЙ АНАЛИЗ
%
%
% Введите ИСХОДНЫЕ ДАННЫЕ
Nb = 4
N = 128
Fs = 5000
A1 = 0.84
A2 = 1.26
f1 = 625
f2 = 1250
sigma = [0 0.84 1.68 3.36]
% Проверьте ПРАВИЛЬНОСТЬ ввода ИСХОДНЫХ ДАННЫХ
% При ПРАВИЛЬНЫХ ИСХОДНЫХ ДАННЫХ введите 1
% При НЕПРАВИЛЬНЫХ ИСХОДНЫХ ДАННЫХ введите 0 и ПОВТОРИТЕ ввод
--> 1
%
%
% Для продолжения нажмите <ENTER>
%
%
% п.1. ПРОВЕРКА ИНФОРМАТИВНОСТИ ПЕРИОДОГРАММЫ В ЗАВИСИМОСТИ ОТ УРОВНЯ ШУМА
%
%
% Для вывода ГРАФИКОВ ПОСЛЕДОВАТЕЛЬНОСТЕЙ и ПЕРИОДОГРАММ
% с различными СКО шума нажмите <ENTER>
%
%
% Для продолжения нажмите <ENTER>
%
%
% п.2. ПРОВЕРКА ИНФОРМАТИВНОСТИ ПЕРИОДОГРАММЫ В ЗАВИСИМОСТИ
% ОТ ПЕРИОДА ДИСКРЕТИЗАЦИИ ПО ЧАСТОТЕ
%
%
% Для вывода ГРАФИКОВ ПЕРИОДОГРАММ при различной размерности ДПФ нажмите <ENTER>
%
%
% Для продолжения нажмите <ENTER>
%
%
% п.3. ПРОВЕРКА ОЦЕНКИ СПМ НА АСИМПТОТИЧЕСКУЮ НЕСМЕЩЕННОСТЬ И СОСТОЯТЕЛЬНОСТЬ
%
%
% Для вывода ГРАФИКОВ зависимости СМЕЩЕНИЯ и СРЕДНЕГО КВАДРАТА
% ОТКЛОНЕНИЯ ИСТИННОЙ СПМ ОТ ЕЕ ОЦЕНКИ нажмите <ENTER>
%
%
% Для продолжения нажмите <ENTER>
%
%
% п.4. ФОРМИРОВАНИЕ СЛУЧАЙНОЙ ПОСЛЕДОВАТЕЛЬНОСТИ С ТРЕБУЕМОЙ АКФ
%
%
% Для вывода ГРАФИКА ТРЕБУЕМОЙ АКФ нажмите <ENTER>
%
%
% С помощью кнопки Zoom in определите и введите ЧЕТНЫЙ ПОРЯДОК КИХ-ФИЛЬТРА
%
    R = 1000
%
% Проверьте ПРАВИЛЬНОСТЬ ввода ИСХОДНЫХ ДАННЫХ
% При ПРАВИЛЬНЫХ ИСХОДНЫХ ДАННЫХ введите 1
% При НЕПРАВИЛЬНЫХ ИСХОДНЫХ ДАННЫХ введите 0 и ПОВТОРИТЕ ввод
-->
%
%
% Для вывода ГРАФИКА ТРЕБУЕМОЙ АКФ И ОЦЕНКИ АКФ РЕАКЦИИ нажмите <ENTER>
%
%
% Для вывода ГРАФИКОВ ИХ, НОРМАЛЬНОГО БЕЛОГО ШУМА И РЕАКЦИИ КИХ-фильтра нажмите <ENTER>
%
%
% Для продолжения нажмите <ENTER>
%
%
% п.5. ФИЛЬТРАЦИЯ СЛУЧАЙНОЙ ПОСЛЕДОВАТЕЛЬНОСТИ С ТРЕБУЕМОЙ АКФ
%
%
% Для вывода ГРАФИКОВ ИСХОДНОЙ и ПРЕОБРАЗОВАННОЙ последовательностей нажмите <ENTER>
%
%
% Для продолжения нажмите <ENTER>
%
%
% п.6. РАСЧЕТ ПЕРИОДОГРАММЫ
%
%
% Для вывода ГРАФИКА ПЕРИОДОГРАММЫ нажмите <ENTER>
%
%
% Для продолжения нажмите <ENTER>
%
%
% п.7. РАСЧЕТ ПЕРИОДОГРАММЫ ДАНЬЕЛЛА
%
%
% Для вывода ГРАФИКОВ ПЕРИОДОГРАММ ДАННЬЕЛЛА нажмите <ENTER>
%
%
% Для продолжения нажмите <ENTER>
%
%
% п.8. РАСЧЕТ ПЕРИОДОГРАММЫ БАРТЛЕТТА
%
%
% Для вывода ГРАФИКОВ ПЕРИОДОГРАММЫ БАРТЛЕТТА нажмите <ENTER>
%
%
% Для продолжения нажмите <ENTER>
%
%
% п.9. РАСЧЕТ ПЕРИОДОГРАММЫ УЭЛЧА
%
%
% Для вывода ГРАФИКОВ ПЕРИОДОГРАММЫ УЭЛЧА нажмите <ENTER>
%
%
% Для продолжения нажмите <ENTER>
%
%
% п.10. РАСЧЕТ ОЦЕНКИ СПМ ПО МЕТОДУ БЛЭКМАНА-ТЬЮКИ
%
%
% Для вывода ГРАФИКОВ ОЦЕНКИ СПМ, вычисленных по методу БЛЭКМАНА-ТЬЮКИ, нажмите <ENTER>
%
%
% Для продолжения нажмите <ENTER>
%
%
% п.11. ОПРЕДЕЛЕНИЕ ПОКАЗАТЕЛЕЙ КАЧЕСТВА ОЦЕНОК СПМ
%
%
% Для вывода ЗНАЧЕНИЙ СКО нажмите <ENTER>
%
   СКО периодограммы

STD_S =

     1.215579860642440e-06

   СКО периодограмм ДАНЬЕЛЛА при различном количестве усредняемых частот K

STD_SD =

   5.000000000000000   0.000000864292665
  10.000000000000000   0.000000769391079
  20.000000000000000   0.000000648468758

   СКО периодограмм БАРТЛЕТТА при различной длине фрагмента L

STD_SB =

   1.0e+02 *

   1.000000000000000   0.000000004102294
   0.500000000000000   0.000000005080934
   0.250000000000000   0.000000005789710

   СКО периодограммы УЭЛЧА при различной длине фрагмента L и величине перекрытия Q

STD_SW =

   1.0e+02 *

   1.000000000000000   0.130000000000000   0.000000006275760
   0.500000000000000   0.130000000000000   0.000000005655281
   0.250000000000000   0.130000000000000   0.000000005474376

   СКО оценок СПМ по методу БЛЭКМАНА-ТЬЮКИ при различных окнах

STD_SBT =

  3×2 cell array

    {'Rectangular Window'}    {[8.835478751280251e-07]}
    {'Hamming Window'    }    {[7.174645087221074e-07]}
    {'Chebyshev Window'  }    {[6.566124051042524e-07]}

%
%
% Для вывода ЗНАЧЕНИЙ ДОБРОТНОСТИ нажмите <ENTER>
%
   ДОБРОТНОСТЬ периодограммы

Q_S =

    0.4925

   ДОБРОТНОСТЬ периодограмм ДАНЬЕЛЛА при различном количестве усредняемых частот K

Q_SD =

    5.0000    0.9742
   10.0000    1.2293
   20.0000    1.7305

   ДОБРОТНОСТЬ периодограмм БАРТЛЕТТА при различной длине фрагмента L

Q_SB =

  100.0000    4.3242
   50.0000    2.8189
   25.0000    2.1709

   ДОБРОТНОСТЬ периодограмм УЭЛЧА при различной длине фрагмента L и величине перекрытия Q

Q_SW =

  100.0000   13.0000    1.5221
   50.0000   13.0000    1.8988
   25.0000   13.0000    2.5194

   ДОБРОТНОСТЬ оценок СПМ по методу БЛЭКМАНА-ТЬЮКИ при различных окнах

Q_SBT =

  3×2 cell array

    {'Rectangular Window'}    {[0.9322]}
    {'Hamming Window'    }    {[1.4137]}
    {'Chebyshev Window'  }    {[1.6879]}

% Для продолжения нажмите <ENTER>
%
%
% п.12. ПОСТРОЕНИЕ СПЕКТРОГРАММЫ
%
%
% Для вывода ГРАФИКА СПЕКТРОГРАММЫ дискретного гармонического сигнала нажмите <ENTER>
%
%
% РАБОТА ЗАВЕРШЕНА
>>