陷波滤波器的以下为matlab的实现代码：
% 用切比雪夫最佳一致逼近设计线性相位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滤波器，但是原理需要根据你的论文了，可能是论文提出的什么方法，这是我的理解：
首先是生成一个信号，采样率1000hz，采样点数1000点；
wp和ws为滤波器参数，具体是低通等类型需要根据下面的分析；
Aw为窗口宽度，这里取奇数；
后面是设计hanning窗，窗幅度根据他提供的公式了，根据hanning窗设计滤波器，应该是低通了；
然后绘制出来设计的滤波器幅频和相频相应；
然后滤波；
会后绘制出来滤波前后波形，这里将250hz滤掉了。
你那个程序标点你复制过来的时候出问题了应该，下面的是可以运行的：
N=1000;
Fs=1000;
n=0:N-1;
t=n/Fs;
x=sin(2*pi*50*t)+sin(2*pi*250*t);
wp=0.5*pi;
ws=0.66*pi;
wdelta=ws-wp;
A=ceil(8*pi/wdelta);
if rem(A,2)==0
A=A+1;
end
Aw=A;
wc=(wp+ws)/2;
a=0:A-1;
alpha=(A-1)/2;
m=a-alpha+0.00001;
hd=sin(wc*m)/(pi*m);
win=hanning(Aw);
h=hd*win;
b=h;
freqz(b,1,512)
yfft=fftfilt(b,x,256);
n1=81:241;
t1=t(n1);
x1=x(n1);
y1=yfft(n1);
subplot(221);
plot(t1,x1);
title('Original Signal')
subplot(222);
plot(t1,y1);
title('Signal after the filter')
另外你问的，end的问题，end这里和if配对的，这是语法内容。
你可以看看这个，这是我以前的实验资料，可以参考，经典的做法：
http://bbs.21eic.com/IndivGroup_Index.asp?GroupID=5&GroupBoardid=3
% Matlab code for plotting the frequency response of digital notch filter
% Implemented as all pass filter section clear; close all fs = 1e6; fn = 200e3; fb = 50e3; omega0T =fn/(fs/2)*pi; deltaT = fb/(fs/2)*pi; a2 = (1-tan(deltaT/2))./(1+tan(deltaT/2)); a1 = (1+a2).*cos(omega0T); B = [1 -a1 a2]; A = [a2 -a1 1]; [H1 W1] = freqz(B,A,1024,'whole'); [H2 W2] = freqz(1,1,1024,'whole'); H3 = (H1+H2)/2; h = figure(1); subplot(2,1,1); plot([-512:511]/1024*fs/1e6,20*log10(fftshift(abs(H3))),'b-','LineWidth',4); grid on; ylabel('amplitude, dB'); title('notch filter, fs=1MHz, fn=200kHz, fb=50kHz'); axis([-0.5 0.5 -50 10]); subplot(2,1,2); plot([-512:511]/1024*fs/1e6,(fftshift(angle(H3)*180/pi)),'m-','LineWidth',4); grid on; xlabel('freq, MHz'); ylabel('angle, deg'); title('phase response'); axis([-0.5 0.5 -180 180]);
【滤波】的意思是什么?【滤波】是什么意思?
11.软件措施包括数字滤波、置看门狗定时时钟及程序陷阱等。 12.适应性保边滤波器 13.适用于三相三线制供电设备。可根据用户需求定制军用级滤波器。 14.梳齿陷波滤波器 15.数字式多通带滤波器 16.衰减的斜率是由滤波器的阶数所决定的。 17.虽然它们远不是理想滤波器，但仍然能把两个族分离得相当好。
陷波滤波器的工作原理以及电路如何实现
只滤除或衰减特定的频率时，可使用陷波滤波器，例如用它滤除电源频率引起的交流声、滤除基波后测量波形失真率等。采用双T电路时，如果采用大的Q值，无用的频率附近的信号也会跟着衰减，因此陷波器的Q值要求可变。双T电路由3个电阻、3个电容组成，基本上是双对称型的。单个无源滤波元件其衰减特性Q=0...
滤波器怎么接线
5、贵的可以加前面，便宜的从可靠性出发建议加后面。6、两根是输入线，两根是输出线，还有一根是地线。滤波器怎么选？1、倘若一滤波器的构成部分，较K常数型具有较尖锐的截止频率（即对频率范围选择性强），而同时对此截止频率以外的其他频率只有较小的衰减率者，称为m常数滤波器。所谓截止频率，亦即与...
什么是手机滤波器???
手机接收和发射时需要一定的频率，滤波器是把频率限定在一定的范围，滤波器有好多种，带通滤波器、高、中、低通滤波器等，市面上没有卖的，因为手机的滤波器大部分都集成在芯片里了。手机的滤波器大部分都是无源滤波器，一般在信号接收部分有高通滤波，一般双频手机有一个900M和一个1800M的高通滤波器...
陷波滤波器原理
最简单的（二阶）陷波滤波器是RLC串联电路，如下图 这个电路滤除的频率可以通过这个公式计算:陷波滤波器的数字实现：下图是数字式陷波滤波器原理图：以下为matlab的实现代码：Matlab code for plotting the frequency response of digital notch filter% Implemented as all pass filter sectionclear; ...
有源滤波与无源滤波有什么区别?
两者的具体区别有以下几点：1、有源滤波器是电子装置，而无源滤波器是机械的。2、有源滤波器是在检测到某一设定好的谐波次数后消除它，无源滤波器是通过电抗器与电容器的相互配合形成某次谐波通道来吸收谐波。3、采用无源滤波器是由于电容器的原因，所以同时可提高功率因素。采用有源滤波器只是消除谐波而...
FIR滤波器
H(z)=∑h(n)*z^-n 就是说，它有（N—1）阶极点在z = 0处，有（N—1）个零点位于有限z平面的任何位置。FIR滤波器有以下几种基本结构：2.1、 横截型（卷积型、直接型）（7.10）式的系统的差分方程表达式为 y(n)=∑h(m)x(n-m) ( 7.11）很明显，这就是线性移不变系统的卷积和...
SAWF是什么器件?有什么作用?使用时要注意什么问题?
声表面波滤波器的英文缩写为SAWF,声表面波滤波器具有体积小,重量轻、性能可靠、不需要复杂调整。在有线电视系统中实现邻频传输的关键器件。声表面波滤波器的特点是:(1)频率响应平坦,不平坦度仅为±0.3-±0.5dB,群时延±30-±50ns。(2)SAWF矩形系数好,带外抑制可达40dB以上。(3)插入损耗虽高达25-30dB,但...
滤波器的一般电路组成?
滤波器是指在复合的信号中滤除某些不需要的频率。比如在整流电路中要的是稳定的直流分量，不需要交流分量，所以要滤掉交流分量。在电视中要从复合的信号中去掉亮度信号取出色度信号，就用梳状滤波器滤掉亮度信号。要取出倦意信号就用6.5MHz的LS滤波器滤除其他信号取出6.5M的倦意信号。根据不同的要求，...
减少噪声的匹配滤波算法
海上时移地震油藏监测技术 上两式中：RY2Y2（m－n）为时间延迟为m－n的时期Y2地震记录在设计窗口中的自相关；RY1Y2（n）为时间延迟为n的时期Y1与时期Y2地震记录在设计窗口中的互相关，于是方程（4.8）可以进一步写成：海上时移地震油藏监测技术 求解方程组（4.11）得到匹配滤波器算子｛P（m），...
端木荷乳块： <p>陷波滤波器(英文叫:notch filter)指的是一种可以在某一个频率点迅速衰减输入信号,以达到阻碍此频率信号通过的滤波效果.从通过信号的频率范围的角度讲,陷波滤波器属于带阻滤波器的一种,只是他的阻带非常狭窄.既然陷波滤波器...
宁波市13959974628： 如何提高陷波电路的衰减倍数 - ？
端木荷乳块： 1,你需要提供你陷波器的输入信号频率2,你陷波器的BW(带宽是多少?)3,感觉你的线路似乎成了一个高通滤波器了 要么就是全部信号都衰减了 参考一下我的空间陷波器的设计 http://hi.baidu.com/solank/blog/item/963906834c4e31ae0df4d271.html
宁波市13959974628： 数字陷波器MATLAB函数 - ？
端木荷乳块： 好像没有直接的函数,你直接傅里叶变换,看着傅里叶变换的频谱图,把不需要的频率挖掉(即使那附近的频率值等于0即可).
宁波市13959974628： 如何用matlab设计一个50Hz陷波器 - ？
端木荷乳块： 亲没有给出采样频率之类的 随便编了一个 好久没用MATLAB了,应该是对的,亲把分给我呗 wp1=45;wp2=55;ws1=49.9; ws2=50.1; rp=0.3;rs=40;Fs=200; wp=[wp1,wp2]/(Fs/2);ws=[ws1,ws2]/(Fs/2); [N,wc]=ellipord(wp,ws,rp,rs,'z'); [num,den]=ellip(N,rp,rs,wc,'stop'); [H,W]=freqz(num,den); plot(W*Fs/(2*pi),abs(H));grid on; xlabel('频率/Hz');ylabel('幅值');
宁波市13959974628： MATLAB/SIMULINK里面的陷波器是什么? - ？
端木荷乳块： 陷波器是一种谐振电路,或者说是一种自动开关的感应器,在天线工程上应用它可以根据信号的频率,自动延长或缩短天线的长度.无线电接收机中专门用于消除某些无用信号以减小对有用信号的干扰的滤波器. 用在电路上滤除不需要的频率...
宁波市13959974628： 反陷波器怎么设计 用matlab？
端木荷乳块： %陷波器的设计 %陷波器的传输函数为 % B(1/z) (z-exp(j*2*pi*f0))*(z-exp(-j*2*pi*f0)) %H(z) = -------- = -------------------------------------------- % A(1/z) (z-a*exp(j*2*pi*f0))*(z-a*exp(-j*2*pi*f0)) %其中f0为陷波器要滤除信号的频率,a为与陷波器深度相关的参...
宁波市13959974628： MATLAB/SIMULINK里面的陷波器是什么?？
端木荷乳块： 陷波器是利用压电效应制成的带阻滤波器,它的作用是阻止或滤掉信号中有害分量对电路的影响
宁波市13959974628： 求Butterworth陷波滤波器MATLAB代码 - ？
端木荷乳块： 据我所知,matlab好像没有专门的陷波滤波器函数.matalb里的滤波器函数都是正规的模拟和FIR、IIR滤波器.而模拟电路里的陷波器往往是双T型的,靠相位抵消而实现的.估计Matlab的设计者不会对这些“怪招”感兴趣吧.正规的做法是用带阻滤波器函数来实现吧,它的优缺点我就不说了.
宁波市13959974628： 联系matlab用双线性变换法设计Butterworth低通滤波器m - ？
端木荷乳块： %下面程序可能对你有用,需根据你的指标修改一下前5行相关e79fa5e98193e58685e5aeb931333264656133参数就行,至于最后显示部分,相应变化一下即可.%1.用Butterworth低通模拟滤波器原型%%%%%%%%%%%%%%%%%%%%%...
宁波市13959974628： matlab butterworth滤波怎么用 - ？
端木荷乳块： 在MATLAB环境下IIR数字滤波器的典型设计和完全设计等方法. 典型设计是先按一定规则将给出的数字滤波器的技术指标转换成模拟低通滤波器的技术指标,据此产生模拟滤波器原型,然后把友橡哗模拟低通滤波器原型转换成模拟低通、高通...
