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

Например, создать взвешивающую функцию с коэффициентом усиления 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