有窄带滤波方面的高手么！（有偿）

快闪

本人一份报告需要对测试数据进行窄带滤波处理，需要这方面高手帮忙！愿意付一定的酬劳，有意者请给zgdky@hotmail.com信箱留言！

dsp2008

MATLAB有现成的函数：fir1，设置好相应的参数即可。

vibmaster

仅供参考，如有不对，还请指正！！！
设计了一个带噪声的正弦信号

clc;
clear all;
close all;
N = 20;        % 每个周期的采样点数
a = 20;        % 周期数
m1 = 1;        % 输入信号的幅值

for n = 1:a*N
    signal(n) = m1*sin(2*pi*(n-1)/N ) ;           % 输入信号
end
snr = 30;         % 设定加入白噪声后的信噪比
signal1 = awgn(signal,snr,'measured');       % 加入噪声后的信号

现在要给这个正弦信号设计一个窄带滤波器，对其进行滤波，其中Fs = 10KHZ

M = 20;                       % 滤波器的阶数
Fs = 10000;                   % 采样频率
fs1 = 350; fp1 = 450; fp2 = 550; fs2 = 650;      % 通带和阻带边界频率
wp1 = 2*pi*fp1/Fs; ws1 = 2*pi*fs1/Fs; wp2 = 2*pi*fp2/Fs; ws2 = 2*pi*fs2/Fs;    % 将频率转换成相位   
wc1 = (ws1+wp1)/2/pi; wc2 = (wp2+ws2)/2/pi;     % 截止频率取通带阻带边缘频率的平均值
wc = [wc1,wc2];            
hn = fir1(M,wc,blackman(M+1));      % 设计FIR滤波器
y = filter(hn,1,signal1);         % 对带噪声的信号进行滤波
plot（y）

ZH----过客

你直接设计频域滤波程序不就完事了吗？？？频域滤波程序你可以随意设计需要的上下限频率，岂不是更简单！！！

Eminem

网上找的，希望给你带来帮助！！！

空间滤波和频域线性滤波的基础都是卷积定理（针对线性系统）：

在使用MATLAB中的fft2函数实现傅立叶滤波时，分两种情况：①无填充时②有填充。

①无填充。假设图像数据为f，直接调用fft2(f)进行FT滤波。

>> f=zeros(256,256);

>> f(1:128,:)=1;

>> imshow(f)%生成简单图像

>> [M,N]=size(f)%获得图像尺寸

>> F=fft2(f);%进行傅立叶滤波，无参数

>> sigma=10;

>> H=lpfilter('gaussian',M,N,sigma);%生成高斯低通滤波器，

>> G=H.*F;%在频域将图像的FT和生成的滤波器点乘（由公式（1），再进行IFT，就可得到滤波后的图像）

>> g=real(ifft2(G));%IFT，再取实部以还原图像

>> figure,imshow(g,[])

滤波后，图像边缘产生了模糊——这是预期的，低通滤波器就是能让图像模糊。但垂直边缘不模糊。其原因在于DFT计算中暗含的周期性。DFT将图像看成向四个方向（x,y轴正负向）延伸的无限周期序列。序列中，原图数据分别向四个方向延伸，可以想像其四周邻接的是什么图像。当滤波器下（同时存在）一块亮区域和暗区域时，滤波结果便得到模糊且变暗的图像。当滤波器位于序列图像的亮部时，它会遇到一个周期分量相同的区域，常数区域的平均值仍为常数，故不会产生模糊现象。

②有填充的情形

无填充时，会产生相邻周期之间的干扰（称为：折叠误差的干扰），可以通过使用0来填充函数的方法来避免。——冈萨雷斯《数字图像处理》

>> PQ=paddedsize(size(f))%求取填充尺寸，如f(x,y)和h(x,y)的大小分别为AxB,CxD，则P≥A+C-1，Q≥B+D-1

PQ =

       512   512

>> fp=fft2(f,PQ(1),PQ(2));%有填充的FT

>> hp=lpfilter('gaussian',PQ(1),PQ(2),2*sigma);%低通滤波器

>> gp=hp.*fp;%频域相乘

>> gp2=real(ifft2(gp));%IFT求取原图像

>> figure

>> imshow(gp2,[])

>> gpc=gp2(1:size(f,1),1:size(f,2)); %将左上部的矩形修剪为原始大小

>> figure

>> imshow(gpc,[])

红轮

直接用工具箱设计滤波器...

缱绻

红轮 发表于 2016-6-28 23:02
直接用工具箱设计滤波器...

怎么操作  具体说说呗  像楼上一样把程序贴出来