带通滤波基础问题 附简单的源程序

sdf4940

%目的是想进行10-30HZ的带通滤波，fft后应该只剩下sin(100*t)这一个波形（因为其频率为100/(2*pi)在10-30HZ之
%间），%但滤波结果不是保留10-30HZ，而是保留10-30角频率。
%请问我该如何改正。简单的源程序如下：



clc;
clear all;

fs=2*pi*512;% 采样频率(HZ)
wp1=2*pi*10;%允许最低频率(角频率)
wp2=2*pi*30;%允许最高频率(角频率)
ws1=2*pi*5; %最低截止频率(角频率)
ws2=2*pi*50;%最高截止频率(角频率)
wp=[wp1 wp2];
ws=[ws1 ws2];
rp=3 ;%波纹系数
rs=60 ;%波纹系数
[n,wn]=ellipord(wp/(fs/2),ws/(fs/2),rp,rs) ;%获取固有频率和阶
[b,a]=ellip(n,rp,rs,wn) ;%返回传递函数参数
%freqz(b,a,512,1000);%1000是频率点数
[h w]=freqz(b,a);
figure(1);
plot(w*fs/(2*pi),abs(h));grid;%频率响应

t=[0:0.01:100];
y1=sin(t)+sin(10*t)+sin(20*t)+sin(30*t)+sin(100*t)+sin(200*t);
y=filter(b,a,y1) ;%用椭圆滤波器滤波，y为滤波后的结果

z=fft(y);
N=length(t);
Fs=100;
%f0=1/(2*pi);
F=([0:N-1])*Fs/N;
out=abs(z);
figure(2)
plot(F(1:N/4),out(1:N/4));
title('幅频特性');
xlabel('频率');
ylabel('幅值');

sdf4940

问题已经知道，是滤波和信号的采样频率不一致造成的