Исследование характеристик аналоговой цепи

Дана аналоговая пассивная цепь:

Рис. 1. Схема аналоговой цепи

Для цепи заданы параметры:

R₁=1000 Ом

R₂=4000 Ом

С₁=1,5 мкФ

C₂=0,5 мкФ

Воспользуемся символическим методом расчёта электрических цепей.

Рис. 2. Схема для расчёта аналоговой цепи

Применим метод «двух узлов». Узел a – заземлим, его потенциал будет равен нулю. Найдём потенциал узла c.

Ток, протекающий через резистор R₂:

Далее найдём потенциал в точке b. Так как , а ,

Получим частотную передаточную функцию цепи:

Код в Matlab:

clc;clear;

R1=1000;

R2=4000;

C1=1.5*10^(-6);

C2=0.5*10^(-6);

Wa0=(R2*C2)^2;

Wa1=R2*C2;

Wb0=C1*R1*(R2*C2)^2;

Wb1=C2*R2*(C2*R2+2*C1*R1+C2*R1);

Wb2=2*C2*R2+C1*R1+C2*R1;

Wanalog=tf([Wa0 Wa1 0], [Wb0 Wb1 Wb2 1])

Подставив значения параметров R и C, получим:

Так как преобразования Фурье и Лапласа имеют схожий вид, можно выполнить замену p=jω. Тогда передаточная функция аналоговой цепи:

Найдём для этой передаточной функции нули и полюса.

Код в Matlab:

disp('Полюса:');

disp(pole(Wanalog));

disp('Нули:');

disp(zero(Wanalog));

Нули:

p=0; p=-500

Полюса:

p=-1000; p=-500; p=-1000/3

Исходя из найденных значений, можем записать передаточную функцию в виде:

Этой передаточной функции соответствует дифференциальное уравнение:

Получим логарифмические амплитудно-частотные и фазо-частотные характеристики для этого фильтра:

Амплитудно-частотная характеристика:

Фазо-частотная характеристика:

Построим графики этих характеристик:

Код в Matlab:

figure(1);

bode(Wanalog,'g');

grid on;

Рис. 3. Логарифмические частотные характеристики аналогового фильтра

Для проверки предыдущих вычислений можно сравнить полученную выше ЛАЧХ с ЛАЧХ, полученной при непосредственном моделировании электрической цепи в пакете Electronics Workbench:

Рис. 4. Моделирование электрической цепи в пакете Electronics Workbench

Получим импульсную переходную функцию этого фильтра.

Для этого разложим передаточную функцию на простые дроби:

Найдём коэффициенты:

C₀=1,5

C₁=-0,5

С помощью обратного преобразования Лапласа, найдём ИПФ:

Построим в системе Matlab график ИПФ:

Код в Matlab:

figure(2);

impulse(Wanalog,'g');

grid on;

Рис. 5. График импульсной переходной функции аналогового фильтра.