gaotong.m

数字信号处理 各种滤波器的设计

%高通(haxes1,haxes2)
function []=gaotong(haxes1,haxes2)
haxes1=axes('position',[0.20,0.10,0.65,0.40]);
haxes2=axes('position',[0.20,0.55,0.65,0.40]);
fp=1,fs=2,ap=1;
fc=fp*(10^(0.1*ap)-1)^(-0.1);
[B1,A1]=butter(5,fc,'s');
[hf1,f1]=freqs(B1,A1,1024);
i=1;
while f1(i)<1
    h1(i)=hf1(i);
    fg1(i)=f1(i);
    i=i+1;
end
j=i-1;
while (f1(i)>=1)&&(f1(i)<2)
    h2(i-j)=hf1(i);
    fg2(i-j)=f1(i);
    i=i+1;
end
k=i-1;
while (f1(i)>=2)&&(f1(i)<3)
    h3(i-k)=hf1(i);
    fg3(i-k)=f1(i);
    i=i+1;
end

[B2,A2]=butter(5,fc,'high','s');
[hf2,f2]=freqs(B2,A2,1024);
i=1;
while f2(i)<1
    m1(i)=hf2(i);
    mn1(i)=f2(i);
    i=i+1;
end
j=i-1;
while (f2(i)>=1)&&(f2(i)<2)
    m2(i-j)=hf2(i);
    mn2(i-j)=f2(i);
    i=i+1;
end
k=i-1;
while (f2(i)>=2)&&(f2(i)<3)
    m3(i-k)=hf2(i);
    mn3(i-k)=f2(i);
    i=i+1;
end

axes(haxes2);
plot(fg1,abs(h1),'r');
axis([0.1,3,0,1.2]);
hold on;
%title('低通滤波特性');

axes(haxes1);
plot(mn1,abs(m1),'r');
axis([0.1,3,0,1.2]);
hold on;
%title('高通滤波特性');
pause(1)
axes(haxes2);
plot(fg2,abs(h2),'k');
axis([0.1,3,0,1.2]);
hold on;
%title('低通滤波特性');
axes(haxes1);
plot(mn2,abs(m2),'k');
axis([0.1,3,0,1.2]);
hold on;
%title('高通滤波特性');
pause(1)
axes(haxes2);
plot(fg3,abs(h3));
axis([0.1,3,0,1.2]);
hold on;
%title('低通滤波特性');
axes(haxes1);
plot(mn3,abs(m3));
axis([0.1,3,0,1.2]);
hold on;
%title('高通滤波特性');