实验三 带阻滤波抑制心电50HZ工频干扰
【doc】一种抑制50Hz工频干扰的方法——基于FTF算法的自适应滤波
一种抑制50Hz工频干扰的方法——基于FTF算法的自适应滤波掌杂志一㈨,341~344f一种抑制50Hz工频干扰的方法-FTF算法的自适应滤波曾海坚王华卿曾祥森_———,'一(华中瞿工大学生袖工程未.武汉430074)林家瑞基于[L内容摘要在长时间动态心电记录(HOLTER)系统中.50Hz工频干扰会影响ECG信号的质量和可诊断性.从而对QRS复波的tl动识别及节律分析产生不良影响因此消除50Hz 工频干扰在ECG信号予处理中占有重要地位.50Hz工颇干扰的频率和幅度是随着电同参数的被动而变化的从本质上讲.50Hz工频干扰的波动过程是一种非平稳随机过程采用固定的窄带陷波器对50]-Iz工撷干扰进行去除,不可能完全抑制50Hz工频干扰本文提出丁一种基于F1下(快速横向tl适应滤波型)的tl适应方法对50Hz干扰进行跟踪,拟合.有效地抑制丁SOHz干扰j最后给出了tl适应滤波结果.并对此方法做了进一步讨非平稳适应咭渡m寸乏关键词非平稳白适应l咭渡丌r』一l(\r—々1导言在生物医学信号的采集过程中,有效地抑制50Hz工额干扰是人们长期研究的问题.人们提出了各种各样的解决方法.最为常用,简便的方法是设计一中心频率固定的带阻滤波器来消除50Hz工频干扰.但是,由于电网负荷的变化.使5OHz交流电在频率及幅度上缓慢波动.从本质上讲,50Hz交流电在频率及幅度上的缓慢波动是一非平稳随机过程.因此,采用固定的窄带陷波器不可能完全抑制50Hz工额干扰.自适应滤波能够自动跟踪输入信号的频率变化,同时能最低限度地减少有用信息(ECG信号)的损失,这对于全信息动态心电记录系统是非常重要的.目前的自适应滤波方法在实时性,收敛性方面可以得到保证(如FTF算法).因而是一种理想的抑制50Hz工额干扰的方法.自适应信息处理算法繁多,从其实质而言可分为遵循最小均方(LMS)准则及最小平方(IS)准则两大类….在70年代中期,Withrow等推出并完善了LMS算法,80年代.人们又从减小输入信号组合的相关性及将Kalman在最佳检测中的结论转化自适应优化等角度.修整补充了LMS算法.同时,人们建立了更完善利用输入信息的处理优化准则及相应算法即RLS型算法.在80年代后期90年代初提出的FTF(快速横向自适应滤波型)算法.使快速自适应信息处理达到实用阶段2基于FTF算法的自适应滤波由自适应原理可知:必须获取一个与干扰信号n(n)有关,但与信号s(n)(ECG信号)无关的参考信号N(n—N)(如图1所示).根据长时间动态心电记录盒(HOLTER系统)的特点,我们设计了一带宽2.5Hz,中心频率50Hz 的带通滤波器,用以从信遭中提取参考信号N (n)(如图1所示).为了尽量减少N(n—N)与s (n)的相关性.带通滤波器的带宽△B不能过大.根据有关资料表明:正常情况下50Hz交流电的频率变化△f<1Hz(一般小于0,2Hz).只要50Hz工额干扰频率波动在1Hz范围之内. 带通滤波器便可保证所需的50Hz干扰(即N (n)滤出,而有用信号s(n)基本上不会出现在参考信号N(n--N)中.图1为自适应滤波原理框图:.～V●f.n生物医学工程学杂志第13誊一田l自适应撼嬗原理框图F1B[oekDiagramofAdaptiveFilter其中:S(n){ECG信号;13.(n),n(13.一N):50Hz干扰信号;13.(n—N):拟合后的参考信号N(n—N):通过带通滤波器获得的参考信号;N(n—N 一1):通过移相器后的正交参考信号我们设计的带通滤波器为一整系数带通滤波器,且没有相位失真.带通滤波器的传递函数为:H(z{等)?(1一Zj(1)延时为N一40,通带增益为80.由前所述可知:我们只需对N(n—N)的幅度,频率两个变量进行调整,可采用两路加权迭加的方法对信道中50Hz干扰进行拟合.为了使自适应算法收敛速度快,我们在其中一路加入一个90'C的移相器,使图1中的两路参考信号(N(n—N),N(n—N一1))互相正交.我们知道当用几个互相正交的分量经加权组合来作自适应处理时,收敛速快.移相器传输函数为H(Z)一l—ZI1(2)由图1可知{n(n--N)=W.(n)*(1--Z.)*N(n—N)+W】(n)*N(n—N)(3)夸wo(n)一w0(n)+wl(n),W1(n)一一w1(n)则n(n--N)一Wo(n)*N(n—N)+w(n)*N(n—N一1)(4)从而完成对50Hz工额干扰的拟台.拟合信号13.(n—N)用于后面处理.(见下节)3FTF自适应算法自适应算法的目的在于寻求最佳加权矢量霄w(n)的迭代公式.最佳加权矢量能够实现对d(n)的最小二乘估计,我们知道实现对d(n) 的最小二乘估计的最佳加权矢量霄(n):WN(n)一k'?d(n)(5)所以必须寻找k的迭代公式.FTF算法通过利用向量空间的迭代(而不是矩阵的迭代)来组成整个算法步骤,如表1所示.算子k将迭代式中变量与有关进入d(n),(n),(n)联系起来,通过几何关系建立k.迭代式.然后转化为所需的迭代式.寰1FTF算法选代过程Table1FTFAlgorithmR~ersiveStet~第4期胡弦山等.一种抑.050Hz工鞭干扰的方法基于FTF算法的自适应滤谴值得注意的是:FTF算法引入了角参量(n).其更新迭代公式为:;…(n)一?(n一1)k'(6)H(n)=[1一一1(n)仃l(n)e(n!n1)]?(n)(7)其中为FPE模方量,e为BPE标量;m(n)为一中间参量.通过以上计算,迭代后得:霄(n)一霄N(n一1)+(n)e(nIn—1)(8)其中:e(nIn1)--X(n)一N(n)?W(n—1)(9)n(n)一N(n)?W(n一1)至此,便完成了霄w(n)的迭代过程.以LS为准则的RLS算法.也是一种横向结构算法,但由于其码元问的运算量大,基本上正比于N,而我们所采用的FTF(快速横向自适应滤波型)算法由于利用了向量空间的迭代(而不是阵的迭代)来组成整个算法步骤,通过向量空间法的分析转化,在一些相应部分,以向量或标量运算代替了矩阵运算.从而使计算量大为减少",计算量为7N—ION.应用FTF算法时,参量的初始化频骤很重要,初始化设置的合理与否,直接关系副收敛速度的快慢.通常情况下可令:(o)一(Oj一霄(0)一己(0)o圈2自适应滤谴结果叵根系般收敛曲拽Fig2Res.hsolAdaptireFEber8LConverge.ceCu~es0fWeighTs(A)涅有50Hz工期干扰的ECG信号(B)自适应浪'瞌结果{c)权系投w,w.收敛曲拽生物医学工程学杂志第18卷(O)一1.0.(o)一E(O)一8其中8为一正小值此外,FTF算法中参变量r可作为一个监测算法正常与否的标志.由于rC080,所以0<r<1.若在计算过程中r>1或r<O,则说明计算结果已经不正常,这通常是由于数值计算误差累积过大造成的.可采用相应措施加以补救,如当r>1时,令r=1,再继续进行迭代.4结果与讨论FTF型算法的主要优点是收敛速度快及对输入信号相关性不敏感.根据我们的实验结果,权系数收敛速度一般n<15,权系数收敛曲线如图2(c)所示而对于LMS算法一般约需要n=40,FTF型算法收敛速度快得多.滤波效果是令人满意的.如果考虑到FTF算法运算效率及对运算结果的监测有效与否.F1,F算法仍有可改进之处:1.FTF算法中间变量比较多,如果对其进行适当组合,随之将有关算式作相应变化,还可以略减少计算量.比如采用增益归一FTF算法,便是引入了g(n)/rw(n)(称之为cw)的变量,使运算量减少约1N.2,FTF算法中补救措施还不够有力主要问题在于r(T)的计算精度不是很高,如果引入式:1F10>来计算r(T).可使补救措施更为有效,及时参考文献1阵尚勤丰晓峰.快速自适应信息处理.民邮电出版社.1993l58WidrowB,GloverJR.口AdapliveNoiseCanceling: PrinciplendApplicatlon—PIEEE.1978i6316988榜福生.高上觊生物医学信号处理.高等教育出版社l989:84唐靖,赶干青简单整系数递归数字滤波器厦其在生物医学中的应用t三).中国医疗器械杂志.1989;18t8】5弦衄寰自适应滤波电子工业出版社.1989?4(收稿日期1998年10月8日)AMethodforRejecting50HzInterferenceAdaptive FilteringBasedonFTFAlgorithm HuZhenshanzengHaijianWangHuaqingZengXiangsenLinJiarui (DeptofB…g…rlgHUST.Wuhan-430074,Abstract50HzinterferenceharmfuItoautomaticdetectionofQRScomplexesandtheclassifi cationofar—rhythmiainHOLTERsystemrecordingTherefore—itisessentialtoreject50HzinterferenceinthepreprocessofECGsigna1.Frequencyandamplitudeof50Hzinterferencevarywiththevariationofthepara metersofpowernet—work.Fixedpropertyofanalogordigitalfilteringisdifficulttotrace50Hznoise.AnAdaptiveFilteringMethodispresentedtotrace50HzinterferenceinECGsignalSOastOeliminateit.Fa,stTransverseAda ptiveFiltering(FTF)algorithmisemployedtoimproveconvergenceandreducetimeofcalculation.Theres ultisgivenaswellaSfurtherdiscussionKeywordsNonstabilityAdaptivefilterFTFalgorithm。
工频50Hz的滤除电路
3.5 工频50Hz 的滤除电路工频干扰是经络信号的主要干扰，虽然前置放大电路对共模干扰具有较强的抑制作用，但有部分工频干扰是以差模信号方式进入电路的，且频率处于经络信号的频带之内，加上电极和输入回路不稳定等因素，前级电路输出的经络信号仍存在较强的工频干扰，所以必须专门滤除。
采用如下图所示是有源双T 带阻滤波器，该电路的Q 值随着反馈系数β（0<β<1）的增高而增高，Q 值与β的关系如下：()141Q β=-，调节R 16和R 17的比值可改变Q值。
图3-103.5.1参数计算：1314R R R ==，15/2R R =，C 7=C 8=C ，C 6=2C 先取0.15C F μ=，由公式012R f Cπ=计算得21.23R =KΩ，在软件上模拟后，调整为： 121422R R ==KΩ，8100.15C C F μ==，1312/211R R ==KΩ，9820.33C C F μ==。
3.5.2 Q 值讨论50Hz 陷波器的传递函数为：2202200()(/)s H s Q s ωωω+++p K （s ）= （1） 幅频特性为：2222220()()(/)A Q ωωωωωωω-=-+p K （2），p K ＝1，0ω＝100πrad 。
国家允许交流供电频率在49.5～50.5Hz 范围内，所以50Hz 陷波器的Q 值并不是越高越好，太高时，阻带过窄，若工频干扰频率发生波动，则根本达不到滤除工频干扰的目的。
而Q 值太小时，又可能会滤掉有用信号。
选择3dB 处截止频率为47.5Hz ，52.5Hz ，将1247.5rad ωπ⨯＝，2252.5rad ωπ⨯＝分别代入222222001()2()(/)A Q ωωωωωωω-==-+p K 中计算得，Q 1＝9.74，Q 2＝10.24，所以取 11214()R R R Q =+＝10，R 17＝22M ， R 16＝510K 。
50Hz工频抑制器设计与仿真
50Hz工频抑制器设计与仿真杨揆【摘要】介绍了一种基于双T网络和运算放大器的工频信号抑制器的设计,用以消除测试试验中50 Hz干扰信号,包括工频信号抑制器的设计原理与参数计算,选择.利用信号发生器与示波器对电路进行滤波试验,并且与DeweSofi数据采集仪的带阻功能进行比较.【期刊名称】《装备制造技术》【年(卷),期】2016(000)004【总页数】4页(P104-107)【关键词】工频干扰;抑制器;双T网络;同相比例运算放大器;带宽BW;带阻抑制器电路(BEF)【作者】杨揆【作者单位】海军装备部,陕西西安 710043【正文语种】中文【中图分类】TP274火炮电气系统低压信号测试中（如：后坐位移测试、依托霍尔器件电流测试、射速测试等）经常遇到50 Hz工频干扰信号通过地线、电磁场串扰进入信号电缆，使得终端测试数据叠加进50 Hz的干扰信号，测试曲线“毛刺”偏多，误差偏大。
图1为火炮射击时的后坐位移信号，由于有工频干扰信号的叠加，对位移信号峰值的读取易产生误差，特别在光标1和2（两条水平虚线）之间为误差易产生区域。
数据采集仪的滤波器多为低通滤波器，无法滤除工频这种特殊干扰，所以设计一种工频陷波器，对50 Hz单一频率进行衰减和滤除，在当前的测试工作中显得特别重要。
工频干扰陷波器对火炮测试过程中经常遇到的50 Hz干扰进行特定的滤除，测试曲线更加光滑、准确，便于测试数据后期分析，提高测试数据的精度。
1.1 抑制器原理图及系统组成如图2所示，工频抑制器由信号带阻电路、信号幅度调理电路、运算放大器调零电路、电源保护电路组成。
运算放大器使用前应先进行调零工作，将输入信号Vi接地，调节R6、R11变阻器阻值，直到运算放大器的输出端V6、V0为零（如图2所示）。
1.2 信号带阻电路1.2.1 信号带阻电路特征频率点的设计图2虚线左侧部分为信号带阻电路，在双T网络后加一级同相比例运算电路就构成了基本的二阶有源带阻滤波器（BEF）[1]。
如何去除心电仪器中的干扰信号的解决方法
如何去除心电仪器中的干扰信号的解决方法作者：张婷婷田桂林来源：《商情》2010年第17期[摘要]心电信号是一种微弱的信号,但是在心电信号中有很多干扰信号,在干扰信号的作用下往往使检查的结果发生微妙的变化,所以本文结合实际生活,介绍了滤除心电信号中干扰信号的一种方法。
[关键词]心电信号滤波 50Hz干扰信号一、前言人体的心电信号是一种淹没在许多较强的噪声和干扰之中,如肌电干扰,放大器噪声及其电极与皮肤之间的噪声,尤其是50Hz工频干扰等。
这些干扰对心电信号的检测影响很大,因此,在检测过程中,如何减少或消除这些干扰,突出所需要的心电信号,就是心电放大设计中要考虑的关键问题。
二、心电信号检测中的干扰分析由于心电信号的微弱,并且存在于强大的干扰信号中,因此对心电放大器的设计,主要的任务就是抗干扰。
所谓干扰,就是除有用信号以外的噪声或受环境影响而使信号变化部分的总称。
随着电子技术的发展,检测装置小型化,而大量的信号电路采用低电平,工作控制环境现场条件恶劣,所以装置的抗干扰技术是我们在设计检测系统时所必须要考虑的问题。
干扰对电子系统产生了很大的危害:在检测元件中,混入的干扰信号会使测量的数据产生误差,造成判断上的失误;在其它领域,如控制系统中,干扰信号可能错误的操作。
因此,为了使电子系统能够正常工作,必须在电子系统中采用抗干扰技术。
三、设计分析在心电信号经过前置级放大器后,然后需要将放大后的信号进行滤波,因为经过分析可以知道,所设计的电路需要完成的主要任务是将信号中的50Hz干扰信号进行滤除,所以所选滤波电路带阻滤波电路,如果从输入信号中减去带通滤波电路汇总处理过的信号,就可以得到带阻信号,而现在我将采用的是双T带阻滤波电路,对50Hz的工频干扰进行抑制,最后通过主放大器将心电信号进一步放大,便于后面波形的显示。
四、所选带阻滤波器的分析电路如图所示,由节点导纳方程不难导出电路的传递函数为:A(s)=V2(s)/V1(s)=Avf[1+(s/wo)2]/[1+2(2-Avf)s/wo+(s/wo)2]式中wo=1/RC,既是特征角频率,也是带阻滤波电路的中心角频率;Avf=Ao=1+Rf/R1为带阻滤波电路的通带电压增益;Q=1/2(2-Ao)。
滤除50Hz工频干扰的滤波电路设计_鲁连钢
收稿日期 ：２ ０ １ １—１ ２—２ ０ ， 男 ， 辽宁铁岭市人 ， 高级讲师 ， 主要从事物理学研究 ． 作者简介 ： 鲁连钢 （ １ ９ ５ ９－ ）
鲁连钢
滤除 ５ ０Ｈ ｚ工频干扰的滤波电路设计
９ １
ＲＷ ． 其中 ， 电位器 ＲＷ 具有调节电路的品质因数 Ｑ 的作用 ，Ｑ 值越大 ， 阻带宽度越窄 ， 选频特性越好 ． 滤除 ５ ０Ｈ ｚ工频干扰的滤波电路如图 ３ 所示 ， 针对 ５ ０Ｈ ｚ工频及其二次谐波 ， 设计了两个窄带阻带中心频率 ， 分 别为５ 工频干 ０Ｈ ｚ和１ ０ ０Ｈ ｚ的独立可调 Ｑ 值双 Ｔ 有源带阻滤波电路 ， 把二者先后串接即构成具有滤除 ５ ０Ｈ ｚ 扰作用的带阻滤波电路 ． 在 ５ ０Ｈ ｚ带阻滤波电路设计中 ， 令 Ｃ Ｃ ０ ． ０ １μ Ｆ， 则 Ｒ Ｒ ３ １ ８ ． ４ｋ Ω，Ｃ １＝ ２＝ １＝ ２≈ ３＝ ， ， ， ， 在 带阻滤波电路设计中 令 则 ０ ． ０ ２μ Ｆ Ｒ １ ５ ９ ． ２ｋ ０ ０Ｈ ｚ Ｃ Ｃ ０ ． ０ ２μ Ｆ Ｒ Ｒ ７ ９ ． ６ｋ Ω． １ Ω Ｃ ３≈ ４＝ ５＝ ５＝ ６≈ ６＝ ０ ． ０ ４μ Ｆ，Ｒ ３ ９ ． ８ｋ ０Ｈ ｚ及 １ ０ ０Ｈ ｚ处 ， 并且阻带带宽 Ω． 调整电位器 ＲＷ１ 、ＲＷ２ ， 使阻带中心频率分别在 ５ ７≈ 在 ２Ｈ ｚ以内 ， 以达到在输出电压ＶＯ 中完全滤除 ５ ０Ｈ ｚ工频信号及其二次谐波干扰的目的 ．
一种抑制心电信号50Hz工频干扰的新方法
・
（ ５ ）
率 ， 中： 其 ｄ＝ （ 一Ｉ ＿ ） １． Ｖ ｔ ／０ ｒ０ （） ７ ｅ＝ （ ＋一 ／０ ｌ １． ０ （） ８ 由于整个 Ｅ Ｇ波形在一个工频干扰周期 内都 Ｃ
＋＿ ｌ“ （ ０ ’６ 一＝ ＿） ） 只Ⅳ ｓ “４ Ｊ
式（ ） ６ 表明， 每计算一个采样点滤波后的值需 要利用该点前 ４ 点和后 ５点的原始值。该算法是把 这 些点 的连线 近似 看 成 一条 直 线 ， 处 理远 离 奇 异 在 点的波形时误差较小 ， 滤波效果也很好。但是在处 理靠近奇异点比如 Ｒ波波峰时 ， 由于这些点的连线
用工频干扰采样点 的幅值代数和为零 ， 在折线段滤 波时同样每次取一个周期的 １ 个 点。每个周期计 Ｏ
Ｋ ｙ ｗｏ ｄ ｅ ｒ ｓ：ＥＣ ｄｓｕｂ ｎ ｅｏ ｄ ｓｒ ￣ ｑ ｅ ｃ ； ｌ ｒ Ｇ； ｉｔｒａ ｃ ｆｉ ｕ ｔ ｎ ｙ ｅ ｕ ｎ ｙ ｆｔ ｉｅ
’
心电信号是心脏 电活动在人体体表的表现 ， 对
心 电信号的分析是 医生诊断心脏疾病的重要依 据。 信号一般 比较微 弱 ， 幅度为毫伏级 ， 率为 ０ ０ — 频 ．５
须 是工 频干 扰 ｆ 的整 数倍 。 （ ） 一 个 工 频 干扰 周 ２在
了一定的效果。但是方法一 中的 Ｍ值 因人 因时间
而异 ， 以确定 ； 法二需 要得 到采样 点 的工 频干扰 难 方
幅值 ， 要计算出幅值须知工频干扰正弦波的幅度、 初 相和频率 ， 最后经三角计算得出， 这无疑大大增加了 算法的复杂度 。
样有 Ｒ波峰值被削弱 ； 值小则滤波效果差。另一 种典型的方法【是首先定位 Ｒ波波峰， ３ 然后对于 Ｒ 值 ， ｅｋｖ 结果 即为滤 波后 的 值 。此方法 的不 足之处 在于噪声值 的获取 比较 困 难， 因为 ５ ｚ ０Ｈ 工频干扰是正弦波 ， 各个采样点的工 频干扰值是变化的 ， 准确地得 到需要进行大量浮点
消除ECG中基频漂移和工频干扰的算法
消除ECG 中基频漂移和工频干扰的算法在心电信号采集过程中,常常受到工频干扰、肌电干扰、基线漂移等各种干扰信号。
本系统中，我们采用了一种简单整系数陷波器，对工频干扰及基线漂移进行滤除。
我们设计的滤波器用到了减法技术，即用一个全通网络减去一个具有相同传输延迟和增益的窄带带通滤波器的输出，得到一个具有尖锐陷波特性的陷波器(NOTCH) ,其原理如下图所示。
该陷波滤波器实际上是一个梳状滤波器。
其传递函数为 211(()1R P R n n n Pz H z Q z -------）=z 式中Q=R/P 。
假设采样率Fs 为1000Hz ，工频干扰Fc 为50Hz ，所以P=Fs/Fc=20，n 和R 决定了滤波器的阻带带宽。
经过实验分析，本系统选取n=2，R=320，则有Q=16,由此得到梳状滤波器的传递函数为320300222011(()161z H z z------）=z 。
上述滤波器的差分方程为：y(n)=y1(n)-y2(n)/256;y1(n)=x(n-300);y2(n)=x(n)-2x(n-320)+x(n-640)+2y2(n-20)-y2(n-40);y1(n)为全通网络输出，y2(n)为窄带带通滤波器输出，y(n)为NOTCH 梳状滤波器输出。
NOTCH 滤波器的幅频特性、相频特性和群延迟特性分别如下图所示。
带阻滤波器原理图05010015020025030035040045050000.51Magnitude Response050100150200250300350400450500-1000-500Phase Response 00.51 1.52 2.53 3.5299.9995300300.0005Group Delay从图中可以看出，上述滤波器对50Hz 工频和基频干扰有较好的滤除效果。
此处采用的简单整系数带阻滤波器的系数都是2的整数幂，在滤波过程中不需要使用乘法，并可以用少量次数的递归计算实现非递归的线性相位FIR 滤波。
LMS算法滤除心电工频干扰
信号处理基础ECG心电图信号干扰消除实验报告姓名：xxxxxx班级：xx班学号：xxxxxxxxxx一、实验要求interference cancellationGiven a received signal r(k) which consist of a source signal s(k) and a sinusoidal interference with known frequency. The task is to extract s(k) from r(k). Notice that the amplitude and phase of the sinusoid is unknown. A well-known application is to remove 50/60 Hz power line interference in the recording of the electrocardiogram (ECG).二、Matlab仿真结果三、结果分析维纳滤波器参数是固定的，适合于平稳随机信号。
卡尔曼滤波器参数是时变的，适合于非平稳随机信号。
然而，只有在信号和噪声的统计特性先验已知的情况下，这两种滤波技术才能获得最优滤波。
在实际应用中，常常无法得到信号和噪声统计特性的先验知识。
在这种情况下，自适应滤波技术能够获得极佳的滤波性能，因而具有很好的应用价值。
本程序是基于LMS算法的自适应滤波，从上图可以看出，滤波效果也是很不错的。
四、源代码av=1;f0=0.2;[s,d]=generateDemoECG(av,f0);% 生成模拟ECG信号；N=length(s);f0=50;fs=500;t=[1:N];n1=100*sin(2*pi*f0*t/fs+20);d=n1+s; % 生成加50HZ干扰模拟ECG信号；%*************************LMS算法滤波*******************************x1=sin(2*pi*f0*t/fs);x2=cos(2*pi*f0*t/fs);t=(0:0.001:1);w0(1)=0;w1(1)=0;mu=0.007;y(1)=w0(1)*x1(1)+w1(1)*x2(1);e(1)=d(1)-y(1);w0(2)=w0(1)+2*mu*e(1)*x1(1);w1(2)=w1(1)+2*mu*e(1)*x2(1);for n=2:Ny(n)=w0(n)*x1(n)+w1(n)*x2(n);e(n)=d(n)-y(n);w0(n+1)=w0(n)+2*mu*e(n)*x1(n);w1(n+1)=w1(n)+2*mu*e(n)*x2(n);endr=d-w0(600)*x1-w1(600)*x2; %滤波后心电信号er=e.*e;%误差功率figure(1)subplot(2,1,1)plot(s((1:2000)))title('模拟ECG信号');xlabel('时间t/ms');ylabel('电压');subplot(2,1,2)plot(r(1:2000))title('滤波后模拟ECG信号');xlabel('时间t/ms');ylabel('电压');figure(2)subplot(2,1,1)plot(w0(1:1000))title('w0 逼近曲线');subplot(2,1,2)plot(w1(1:1000))title('w1 逼近曲线');figure(3)subplot(1,1,1)semilogy(er(1:2000));title('误差功率');generateDemoECG.mECG信号产生文件，将本文件代码放到MATLAB工作路径下function [x2,x3]=generateDemoECG(av,af)pv=0.1*208;qv=0.1*208;rv=1*208;sv=0.2*208;tv=0.2*208;st=0.1*208;x2(1:100)=0;x2(101:150)=pv-pv*cos([0:49] * 2*pi/50);x2(151:170)=0;x2(171:175) = -[1:5]*qv/5;x2(176:190)=-qv+[1:15]*(rv+qv)/15;x2(191:210)=rv- [1:20]*(rv+sv)/20;x2(210:219)=-sv + [1:10]*(sv+st)/10;x2(220:259)=st;x2(260:359)=tv-tv*cos([1:100]*2*pi/100) + st- [1:100]*st/100;x2(360:500)=0;x2(501:1000)=x2(1:500);x2(1001:2000)=x2(1:1000);x2(2001:4000)=x2(1:2000);x2(4001:10000)=[x 2(1:4000),x2(1:2000)];x3 = x2 + 208*av*sin([1:length(x2)]*af/500 *2*pi);%end of generateDemoECG。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
矿产资源开发利用方案编写内容要求及审查大纲
矿产资源开发利用方案编写内容要求及《矿产资源开发利用方案》审查大纲一、概述
㈠矿区位置、隶属关系和企业性质。
如为改扩建矿山, 应说明矿山现状、
特点及存在的主要问题。
㈡编制依据
(1简述项目前期工作进展情况及与有关方面对项目的意向性协议情况。
(2 列出开发利用方案编制所依据的主要基础性资料的名称。
如经储量管理部门认定的矿区地质勘探报告、选矿试验报告、加工利用试验报告、工程地质初评资料、矿区水文资料和供水资料等。
对改、扩建矿山应有生产实际资料, 如矿山总平面现状图、矿床开拓系统图、采场现状图和主要采选设备清单等。
二、矿产品需求现状和预测
㈠该矿产在国内需求情况和市场供应情况
1、矿产品现状及加工利用趋向。
2、国内近、远期的需求量及主要销向预测。
㈡产品价格分析
1、国内矿产品价格现状。
2、矿产品价格稳定性及变化趋势。
三、矿产资源概况
㈠矿区总体概况
1、矿区总体规划情况。
2、矿区矿产资源概况。
3、该设计与矿区总体开发的关系。
㈡该设计项目的资源概况
1、矿床地质及构造特征。
2、矿床开采技术条件及水文地质条件。