matlab 50hz陷波,新手求消除50HZ工频干扰陷波滤波器源程序

该楼层疑似违规已被系统折叠隐藏此楼查看此楼

%陷波器的设计

%陷波器的传输函数为

% 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为与陷波器深度相关的参数，a越大，深度越深。

%已知信号中50Hz工频干扰，信号为x=sin(2*pi*50*n*Ts)+sin(2*pi*125*n*Ts);

%要求通过陷波器滤除50Hz干扰信号

%参数设置：采样率Ts=0.001s，采样长度：512点

clf;clear;

%设置初值

f0=50;

Ts=0.001;

fs=1/Ts;

NLen=512;

n=0:NLen-1;

%陷波器的设计

apha=-2*cos(2*pi*f0*Ts);

beta=0.96;

b=[1 apha 1];

a=[1 apha*beta beta^2];

figure(1);

freqz(b,a,NLen,fs);%陷波器特性显示

x=sin(2*pi*50*n*Ts)+sin(2*pi*125*n*Ts);%原信号

y=dlsim(b,a,x);%陷波器滤波处理

%对信号进行频域变换。

xfft=fft(x,NLen);

xfft=xfft.*conj(xfft)/NLen;

y1=fft(y,NLen);

y2=y1.*conj(y1)/NLen;

figure(2);%滤除前后的信号对比。

subplot(2,2,1);plot(n,x);grid;

xlabel('Time (s)');ylabel('Amplitude');title('Input signal');

subplot(2,2,3);plot(n,y);grid;

xlabel('Time (s)');ylabel('Amplitude');title('Filter output');

subplot(2,2,2);plot(n*fs/NLen,xfft);axis([0 fs/2 min(xfft) max(xfft)]);grid;

xlabel('Frequency (Hz)');ylabel('Magnitude (dB)');title('Input signal');

subplot(2,2,4);plot(n*fs/NLen,y2);axis([0 fs/2 min(y2) max(y2)]);grid;

xlabel('Frequency (Hz)');ylabel('Magnitude (dB)');title('Filter output');

matlab 50hz陷波,新手求消除50HZ工频干扰陷波滤波器源程序相关推荐

50hz 60hz 级联陷波器,心电信号50Hz陷波器的FPGA实现
收稿日期: 2009 - 03 - 23 心电信号 50 Hz 陷波器的 FPGA 实现林霖 , 张志德 (南方医科大学生物医学工程学院 , 广东广州 510515) [中图分类号]TH772 ...
50Hz工频干扰消除
今天整理工频干扰消除算法. 我们知道,设计数字滤波器,和模拟滤波器的实质,其实就是求一组系数,逼近要求的频率响应. 模拟滤波器已经很成熟,因此,数字滤波器的设计,将S平面映射到Z平面就型.采用双线性变 ...
陷波器介绍_50Hz工频信号陷波器设计
文章目录学习目标: 基本概念: 基本原理: 参数的具体计算及选择: 具体计算 # 陷波器的意义学习目标: 了解陷波器的基本概念掌握50HZ工频陷波器的基本电路图基本概念: 提示:这里可以添加要 ...
怎么用matlab分析孔隙度,基于MATLAB软件的声波测井孔隙度求取
·99·2016年第 24 期基于MATLAB软件的声波测井孔隙度求取向旻 (新疆工程学院,新疆乌鲁木齐 830091) 摘要 :目前,声波测井是一种重要的孔隙度测井方法,被各大油田广泛的使用 ...
将z的第二列除以根号三MATLAB,matlab中如何用牛顿法求根号2的近似值
如何用matlab求出图中各条直线的斜率 1.瞬时斜率:求导.方法:d(k)=(u(k)-u(k-1))/T(k);d(0)=0;其中d就是瞬时斜率2.平均斜率:求拟合.方法:p=polyfit(x, ...
Matlab实现连通域标记算法求图像连通域
Matlab实现连通域标记算法求图像连通域连通域连通域标记算法连通域连通区域(Connected Component)一般是指图像中具有相同像素值且位置相邻的前景像素点组成的图像区域(Regi ...
ubuntu14.04扩展屏幕后，打开matlab就死机。求大神帮忙解决！！！万分感谢
ubuntu14.04扩展屏幕后,打开matlab就死机.求大神帮忙解决!!!万分感谢
matlab输出的特征向量,关于matlab中的eig函数(求特征值和特征向量)(最新整理)
<关于matlab中的eig函数(求特征值和特征向量)(最新整理)>由会员分享,可在线阅读,更多相关<关于matlab中的eig函数(求特征值和特征向量)(最新整理)(3页珍藏版)& ...
MATLAB输出四种常用波形(方波、正弦波、三角波、锯齿波)
1.我们在处理数学以及信号的时候经常会用到MATLAB,因为它可以使我们直观地观察到我们需要的信号的模型.因此掌握基本的MATLAB还是有必要的.这里我们来讲解一下如何用MATLAB产生四种常用波形( ...

matlab 50hz陷波,新手求消除50HZ工频干扰陷波滤波器源程序

matlab 50hz陷波,新手求消除50HZ工频干扰陷波滤波器源程序相关推荐

最新文章

热门文章