Filtros usando matlab

%Rechaza banda filter
w=logspace(0,9,1000);
denom=(j*w+17846.1+74988.6*j).*(j*w+17846.1-74988.6*j);
num=(j*w-77082.9*j).*(j*w+77082.9*j);
H=num./denom;
Hdb=20*log10(abs(H));
semilogx(w,Hdb)
colormap(gray);grid;
xlabel('frecuencia(rad/s)');
ylabel('H(jw)(db)');

%pasa altas filter
w=logspace(0,9,1000);
denom=(j*w+17846.1+74988.6*j).*(j*w+17846.1-74988.6*j);
num=(j*w*12.973/1000000).*(j*w*12.973/1000000);
H=num./denom;
Hdb=20*log10(abs(H));
semilogx(w,Hdb)
colormap(gray);grid;
xlabel('frecuencia(rad/s)');
ylabel('H(jw)(db)');

%filtro pasa banda
w=logspace(0,9,1000);
denom=(j*w+9803.9-76456.9*j).*(j*w+9803.9+76456.9*j);
num=(j*w*3.3/1000000);
H= num./denom;
Hdb=20*log10(abs(H));
semilogx(w,Hdb)
colormap(gray);grid;
xlabel('frecuencia(rad/s)');
ylabel('H(jw)(db)');