ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 26.01.2021
Просмотров: 497
Скачиваний: 2
end;
NE=NE-1
%
оценка
показателя
(
реакции
)
по
выборке
наблюдений
L=l/NE
Yl(j)=L;
end;
%
определение
коэффициентов
регрессии
Cl=X*X';
b_l=inv(Cl)*X*Yl'
%
формирование
зависимости
реакции
системы
на
множестве
%
значений
факторов
Al=minf(1):0.0001:maxf(1);
Bl=minf(2):0.01:maxf(2);
[k N1]=size(Al);
[k N2]=size(Bl);
for i=1:N1,
for j=1:N2,
anl(i)=2*(Al(i)-minf(1))/(maxf(1)-minf(1))-1;
bnl(j)=2*(Bl(j)-minf(2))/(maxf(2)-minf(2))-1;
%
экспериментальная
поверхность
реакции
Yo(j,i)=b_l(1)+anl(i)*b_l(2)+bnl(j)*b_l(3)+anl(i)*bnl(j)*b_l(4);
end;
end;
Результаты
работы
этого
фрагмента
УП
,
выводимые
в
командном
окне
системы
,
отображаются
следующим
образом
:
Am = 0
R = 1
fraceks =
0.0010 1.0000
0.0010 2.0000
0.0100 1.0000
0.0100 2.0000
NP = 0.0010
mr = 1
NE = 334
L = 0.0078
NP = 0.0010
mr = 2
NE = 22
L = 0
NP = 0.0100
mr = 1
NE = 336
L = 0.0071
NP = 0.0100
mr = 2
NE = 22
L = 0
b_l =
0.0037
-0.0002
-0.0037
0.0002
Следует
отметить
,
что
в
этом
случае
при
тактическом
планировании
объем
испытаний
определяется
внутри
цикла
с
учетом
результатов
оценки
выборочной
дисперсии
.
Поэтому
число
испытаний
для
каждого
набора
исходных
данных
стратегического
плана
будет
различным
.
В
завершающей
части
m-
файла
УП
реализуется
графическое
отображение
результатов
моделирования
%
отображение
зависимостей
в
трехмерной
графике
[x,y]=meshgrid(A,B);
[xl,yl]=meshgrid(Al,Bl);
figure;
subplot(1,2,1),plot3(x,y,Yc),
xlabel('fact a'),
ylabel('fact b'),
zlabel('Yc'),
title('PO'),
grid on,
subplot(1,2,2),plot3(xl,yl,Yo),
xlabel('fact a'),
ylabel('fact b'),
zlabel('Yo'),
title('L'),
grid on;
Это
позволяет
получить
графики
зависимостей
реакции
системы
,
представленные
на
рис
. 6.18
а
)
б
)
Рис
. 6.18.
Графики
зависимостей
реакции
системы
Анализ
полученных
зависимостей
показывает
,
например
,
что
величина
L
при
выбранном
способе
задания
порога
не
зависит
от
NP
(
фактор
а
),
а
зависит
только
от
величины
mr,
характеризующей
относительное
превышение
порога
над
уровнем
входного
шума
.
Подобный
результат
выглядит
вполне
естественно
.
В
заключение
следует
отметить
,
что
наибольшую
сложность
при
реализации
аналогичного
типа
S-
моделей
для
исследования
систем
передачи
информации
импульсных
и
непрерывных
сигналов
с
различными
видами
модуляции
представляет
согласованная
настройка
блоков
модуляторов
и
демодуляторов
.
Для
этого
требуется
внимательно
изучить
соответствующие
разделы
помощи
MATLAB
и
Simulink.
Целесообразно
также
исследовать
стандартные
m-
функции
расчета
характеристик
фильтровых
устройств
,
если
их
параметры
фигурируют
в
окнах
настройки
блоков
.
Достаточно
наглядно
особенности
реализации
алгоритмов
модуляции
и
демодуляции
представлены
в
описании
пакета
расширения
Signal Processing Toolbox
системы
MATLAB.
Здесь
модуляция
осуществляется
на
основе
вызова
функций
y modulate(x,fc,fs,'method ',[.opt])
=
;
Эта
функция
генерирует
вектор
вещественных
отсчетов
модулированного
сигнала
с
несущей
частотой
fc
и
частотой
дискретизации
по
времени
fs.
Модулирующий
сигнал
задается
отсчетами
вектора
x.
Параметр
′
method
′
определяет
вид
модуляции
(amdsb – sc, amdsb – tc, amssb, fm, pm, pwm,
ptm
и
др
.).
Параметр
opt
определяет
относительное
соотношение
диапазона
изменения
модулирующего
сигнала
и
значения
соответствующего
параметра
несущего
радиосигнала
.
Например
,
при
задании
′
method
′
=
′
amdsb – tc
′
осуществляется
амплитудная
двухполосная
модуляция
сигнала
с
изменяемой
глубиной
модуляции
)
t
fc
pi
2
cos(
).
opt
x
(
y
∗
∗
∗
∗
+
=
;
Для
модуляции
используется
функция
x demod(y,fc,fs,'method ',[.opt])
=
;
Здесь
все
параметры
должны
соответствовать
используемым
при
модуляции
.
При
задании
′
method
′
=
′
amdsb – tc
′
при
демодуляции
реализуется
следующий
набор
операций
:
)
t
fc
pi
2
cos(
y
x
∗
∗
∗
∗
=
;
)
fs
2
*
fc
,
n
(
butter
]
denum
,
num
[
=
;
)
x
,
denum
,
num
(
filtfilt
x
=
;
В
этой
последовательности
функция
butter (…)
реализует
расчет
коэффициентов
передаточной
функции
низкочастотного
фильтра
Баттерворта
(Lowpass filter numerator - denominator)
порядка
n (n = 2
для
блока
демодуляции
,
представленного
настройками
на
рис
. 6.12)
с
нормированной
частотой
среза
1
fs
fc
2
Wn
≤
×
=
.
Функция
filtfilt (…)
реализует
фильтрацию
с
применением
параметров
рассчитанного
фильтра
с
исключением
фазовых
искажений
.
Аналогичным
образом
с
использованием
помощи
системы
могут
быть
определены
характеристики
фильтровых
устройств
для
детектирования
сигналов
с
другими
видами
модуляции
.