【50hz工频陷波器】

昨天做了一个50hz陷波器，记录一下~

单运算放大器双T陷波滤波器

matlab算参数

clc;
clear all;
close all;

fn=50;%stop
bw=2;%带宽
c=0.1*(10^-6);%uf

r=1/(4*pi*fn*c);%%实际取2r
fprintf('r = %.2f 欧姆\n', r);
fprintf('r = %.2f 欧姆\n', 2*r);

q=fn/bw;
k=1-1/(4*q);
fprintf('k = %.2f 欧姆\n', k);

syms x%x:r3 y:r4
y=10000;
eqn = k*x+k*y-y == 0;
sol = solve(eqn, x);
disp('方程的解为:');
disp(sol);

fndb=20*log10(1/q);%%预期缺口深度，db
fprintf('ln(10) = %.4f\n', fndb);