Создайте функции взвешивания в непрерывном времени путем определения низкочастотного усиления, высокочастотного усиления и величины усиления на некоторой промежуточной частоте.

Например, создайте функцию взвешивания с усилением 40 дБ на низкой частоте, скатываясь до -20 дБ на высокой частоте. Укажите далее, что коэффициент усиления составляет около 10 дБ при 1 рад/с путем помещения этих значений в вектор [freq,mag]. Задайте коэффициент усиления в абсолютных модулях.

Wl = makeweight(100,[1,3.16],0.1);

Создайте весовую функцию с усилением -10 дБ на низкой частоте, поднявшись до 40 дБ на высокой частоте. Задайте частоту среза 0 дБ 10 рад/с. Чтобы задать частоту среза 0 дБ, можно использовать частоту среза в качестве второго входного параметра вместо вектора [freq,mag].

Wh = makeweight(0.316,10,100);

Постройте график величин функций взвешивания, чтобы подтвердить, что они соответствуют спецификациям отклика.

bodemag(Wl,Wh)
legend
grid on

Создайте профиль усиления, который откатывается на высокой частоте без уплощения. Задайте коэффициент усиления 40 дБ на низкой частоте и частоту среза 10 рад/с.

W = makeweight(100,[10 1],0);

Установка высокочастотного усиления 0 гарантирует, что частотная характеристика откатывается на высоких частотах без смещения. Постройте профиль усиления, чтобы подтвердить эту форму.

bodemag(W)
grid on

Создайте функции взвешивания в дискретном времени путем определения низкочастотного усиления, высокочастотного усиления, величины усиления на некоторой промежуточной частоте и шаге расчета.

Создайте функцию взвешивания со шаг расчета 0,1 с. Задайте коэффициент усиления 40 дБ на низкой частоте, скатываясь до -20 дБ на высокой частоте. Укажите далее, что коэффициент усиления составляет около 10 дБ при 0,01 рад/с. Обеспечьте все усиления в абсолютных модулях.

Wl = makeweight(100,[0.01,3.16],0.1,0.1);

Создайте весовую функцию с усилением -10 дБ на низкой частоте, поднявшись до 40 дБ на высокой частоте. Задайте частоту среза 0 дБ 2 рад/с и шаг расчета 0,1 с. Чтобы задать частоту среза 0 дБ, можно использовать частоту среза в качестве второго входного параметра вместо вектора [freq,mag].

Wh = makeweight(0.316,2,100,0.1);

Постройте график величин функций взвешивания, чтобы подтвердить, что они соответствуют спецификациям отклика.

bodemag(Wl,Wh)
grid on

Высокочастотное выравнивание Wh искажается из-за близости его частоты среза к частоте Найквиста.

По умолчанию makeweight создает весовые функции первого порядка. Если вам нужен более резкий переход между низкочастотным и высокочастотным усилениями, можно задать порядок с последним входным параметром. Например, предположим, что вы хотите создать функцию взвешивания со шаг расчета 0,1 с. Функция имеет коэффициент усиления -10 дБ на низкой частоте, поднимаясь до 40 дБ на высокой частоте. Кроме того, коэффициент усиления проходит через 6 дБ при 1 рад/с. Для сравнения создайте как функцию третьего порядка, так и функцию первого порядка с этими спецификациями.

W3 = makeweight(0.316,[1 2],100,0.1,3);
W1 = makeweight(0.316,[1 2],100,0.1);
bodemag(W3,W1)
legend('location','northwest')
grid on

Для функции первого порядка высокочастотное выравнивание искажается из-за близости его частоты среза к частоте Найквиста. Использование более резкого, более высокого порядка перехода гарантирует, что функция выровнялась до достижения частоты Найквиста.

Чтобы создать функции взвешивания в непрерывном времени более высокого порядка, установите Ts = 0. Например, создайте функции взвешивания в непрерывном времени с такими же спецификациями усиления, как W1 и W3.

W3c = makeweight(0.316,[1 2],100,0,3);
W1c = makeweight(0.316,[1 2],100);
bodemag(W3c,W1c)
legend('location','northwest')
grid on