求设计一个带通滤波器到matlab程序

用切比雪夫最佳一致逼近设计线性相位FIR带通滤波器；

%信号为0.5hz,

0.9hz,

1.1hz和1.5hz的正统信号叠加组成

%通带为[0.9,1.1]

%频谱分辨率与信号实际长度N成正比

clear

all

f1=0.5f2=0.9f3=1.1f4=1.5t=0:1203N=length(t)fs=10M=512

x1=sin(2*pi*(f1/fs)*t)+sin(2*pi*(f2/fs)*t)+sin(2*pi*(f3/fs)*t)+sin(2*pi*(f4/fs)*t)

figure(1)

subplot(211)plot(t,x1)title('原信号')

y=fft(x1)

f=(0:1/N:1/2-1/N)*fs

subplot(212)plot(f,abs(y(1:N/2)))gridxlabel('hz')%处理前频谱

wc1=2*f2/fswc2=2*f3/fswc3=2*f4/fs%归一化角频率，用于下面的f1

f1=[0

wc1-0.05

wc1

wc2

wc2+0.05

1]

A=[0

0

1

1

0

0]%设置带通或带阻，1为带通，0为带阻

weigh=[1

1

1

]%设置通带和阻带的权重

b=remez(60,f1,A,weigh)%传函分子

h1=freqz(b,1,M)%幅频特性

figure(2)

f=(0:1/M:1-1/M)*fs/2

subplot(211)plot(f,abs(h1))gridtitle('带通')

x2=filter(b,1,x1)

S1=fft(x2)

f=(0:1/N:1/2-1/N)*fs

subplot(212)plot(f,abs(S1(1:N/2)))gridxlabel('hz')%处理后频谱

实际上就是高通和

低通滤波器

的合体嘛，理解上不难。

150MHz信号带宽是6MHz，那你就设计一个147MHz~153MHz的带通滤波器呗；

145MHz信号带宽是30KHz，那你同理可以计算出上下限并设计出带通滤波器的导通和截至

频点

。

分离开这两个信号，再在不同的端口用

示波器

捕捉看就好啊。