毛真真的回答：% 用切比雪夫最佳一致逼近设计线性相位FIR带通滤波器； %信号为0.5hz, 0.9hz, 1.1hz和1.5hz的正统信号叠加组成 %通带为[0.9,1.1] %频谱分辨率与信号实际长度N成正比 clear all; f1=0.5;f2=0.9;f3=1.1;f4=1.5;t=0:1203;N=length(t);fs=10;M=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)));grid;xlabel('hz');%处理前频谱 wc1=2*f2/fs;wc2=2*f3/fs;wc3=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));grid;title('带通'); 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)));grid;xlabel('hz');%处理后频谱 王晓晖的回答：% 用切比雪夫最佳一致逼近设计线性相位FIR带通滤波器； %信号为0.5hz, 0.9hz, 1.1hz和1.5hz的正统信号叠加组成 %通带为[0.9,1.1] %频谱分辨率与信号实际长度N成正比 clear all; f1=0.5;f2=0.9;f3=1.1;f4=1.5;t=0:1203;N=length(t);fs=10;M=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)));grid;xlabel('hz');%处理前频谱 wc1=2*f2/fs;wc2=2*f3/fs;wc3=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));grid;title('带通'); 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)));grid;xlabel('hz');%处理后频谱 绿水蓝天的回答：程序： wlp = 03*pi; whp = 04*pi; wc = [wlp/pi,whp/pi];n = 33; n=0:n-1; hn=fir1(n-1,wc, hanning (n)); %基于窗函数的 fir 滤波器[h2,w]=freqz(hn,1,512); %freqz是用离散傅里叶变换的标准公示计算的 fft使用快速傅里叶变换的公式计算的，求设计一个hamming窗设计的带通滤波器的matlab程序。 |
