求设计一个带通滤波器到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')%处理后频谱

%写上标题

%设计低通滤波器：

[N,Wc]=buttord()

%估算得到Butterworth低通滤波器的最小阶数N和3dB截止频率Wc

[a,b]=butter(N,Wc) %设计Butterworth低通滤波器

[h,f]=freqz() %求数字低通滤波器的频率如派响应

figure(2) % 打开窗口2

subplot(221) %图形显示分割窗口

plot(f,abs(h)) %绘制Butterworth低通滤波器的幅频响应图

title(巴氏低通滤波器燃旅'')

grid %绘制带网格的图像

sf=filter(a,b,s) %叠加函数S经过低通滤波器以后的新函数

subplot(222)

plot(t,sf) %绘制叠加函数S经过低通滤波器以后的时域图形

xlabel('时间 (seconds)')

ylabel('时间按幅度')

SF=fft(sf,256) %对叠加函数S经过低通滤波器以后的新函数进行256点的基—2快速傅立叶变换

w= %新信号角频率

subplot(223)

plot())%绘制叠加函数S经过低通滤波器以后的频谱图

title('低通滤波后的频谱图')

%设计高通滤波器

[N,Wc]=buttord()

%估算得到Butterworth高通滤波器的最小阶数N和3dB截止频率Wc

[a,b]=butter(N,Wc,'high')%设计Butterworth高通滤波器

[h,f]=freqz()%求数字高通滤波器的频率响应

figure(3)

subplot(221)

plot())%绘制Butterworth高通滤波器的幅频响应图

title('巴氏高通滤波器')

grid%绘制带网格的图像

sf=filter()%叠加函数S经过高通滤波器以后的新函数

subplot(222)

plot(t,sf)%绘制叠加函数S经过高通滤波器以后的时域图形

xlabel('Time(seconds)')

ylabel('Time waveform')

w%新信号角频率

subplot(223)

plot())%绘制叠加函数S经过高通滤波器以后的频谱图

title('高通滤波后的频谱图')

%设计带通滤波器

[N,Wc]=buttord([)

%估算得到Butterworth带通滤波器的最小阶数N和3dB截止频率Wc

[a,b]=butter(N,Wc)%设计Butterworth带通滤波器

[h,f]=freqz()%求数字带通滤波器的频率响应

figure(4)

subplot(221)

plot(f,abs(h))%绘制Butterworth带通滤波器的幅频响应图

title('butter bandpass filter')

grid%绘制带网格的图像

sf=filter(a,b,s)%叠加函数S经过带通滤波器以后的新函数

subplot(222)

plot(t,sf)%绘制叠加函数S经过带通滤波器以后的时域图形

xlabel('Time(seconds)')

ylabel('Time waveform')

SF=fft()%对叠加函数S经过带通滤渣段贺波器以后的新函数进行256点的基—2快速傅立叶变换

w=( %新信号角频率

subplot(223)

plot('))%绘制叠加函数S经过带通滤波器以后的频谱图

title('带通滤波后的频谱图')