本发明涉及体积介质阻挡放电(dielectricbarrierdischarge:dbd)反应器控制技术领域，特别是一种功率密度调制(powerdensitymodulation:pdm)工作模式的高压电源激励体积dbd空气净化装置的供电频率调控系统和方法。
背景技术：
近几年来，大气压等离子体被广泛应用于环境修复，污染控制，生物医学，流动控制和材料处理等方面。体积dbd是有绝缘介质插入放电空间，并形成一定放电间隙体积的非平衡态气体放电，具有放电稳定和放电面积大等特点，是在常压下产生低温等离子体的有效方法之一，适合于大规模工业生产，具有广阔的应用前景。在实际应用中，dbd的供电电压、供电电流、供电能量等电学参量，是低温等离子体研究和应用中极其重要的参量，在放电过程中它们相互作用，共同影响多重微放电的产生、放电的效果和等离子体的特性。
体积dbd反应器通常采用高压或脉冲高压电源驱动，采用pdm工作模式的高压电源驱动dbd反应器具有功率密度调控方便、放电反应器温度低(等离子体非平衡态高)、工作稳定等优点。在介质阻挡放电中，随着放电过程的进行，dbd反应器的等效电容会发生改变。由于容性的反应器与感性的激励电源之间存在匹配关系，当反应器的等效电容发生变化时，反应器与电源的匹配特性下降。在pdm电源的输入电压相同的情况下，pdm电源产生的高压激励电压峰值会减小，即原本注入到反应器上的供电电压会减小，而只有当电源的供电频率和电极的谐振频率越接近，放电效果才越好。此时就需要调整反应器的供电频率，从而提高dbd反应器的供电电压，提高反应器的放电强度。
技术实现要素：
本发明要解决的技术问题为：利用pdm工作模式的高压激励电源对体积dbd反应器进行激励放电，实现dbd反应器的空气净化作用，同时在dbd反应器的运行过程中，通过调整供电频率，达到调整dbd反应器供电能量的目的，保证反应器的放电强度。
本发明采取的技术方案为：一种pdm工作模式的高压激励电源驱动dbd空气净化装置的供电频率调控系统，包括控制单元、可编程电源、pdm高压激励电源、体积dbd反应器、电气参数检测单元和光谱检测单元；
dbd反应器的地线上串接有积分电容cm；电气参数检测单元包括供电电压检测电路、放电电流检测电路和积分电压检测电路；供电电压检测电路检测体积dbd反应器的激励电源输入端的供电电压信号，放电电流检测电路检测体积dbd反应器的放电电流信号，积分电压检测电路检测积分电容cm上的积分电压信号，光谱检测单元检测体积dbd反应器反应区域的发射光谱强度信号，分别传输至控制单元；
可编程电源控制pdm电源的供电电压峰值输出，pdm电源的输出端连接体积dbd反应器的激励电压输入端；
控制单元根据接收到的供电电压信号、放电电流信号、积分电压信号和发射光谱强度信号，计算体积dbd反应器的供电能量和相对光量子产率能效比，进而根据计算结果控制pdm电源对体积dbd反应器的供电频率，使得体积dbd反应器在pdm电源的激励下放电，并工作在设定的最佳相对光量子产率能效比波动范围内。
本发明所述设定的相对光量子产率能效比波动范围即最佳eer附近的范围，最佳eer波动范围的端点值利用牛顿爬山算法得到，波动范围的区间长度值为设定值。
进一步的，体积dbd反应器包括进气口和出气口，进气口上设有气体流速调节单元和气体流速传感器，出气口上设有气体质量检测单元；气体流速调节单元的控制输入端连接控制单元，气体流速传感器检测进气口的气体流速信号，气体质量检测单元检测出气口的气体质量信号，分别传输至控制单元；控制单元根据接收到的气体流速信号，通过气体流速调节单元对气体流速进行调节。
优选的，控制单元通过气体流速调节单元控制体积dbd反应器内的初始气体流速为2.5m/s；dbd的运行过程中，当检测到气体流速在2.0m/s～3.5m/s范围之外时，对气体流速进行调节。
优选的，供电电压检测电路包括依次连接的分压电路、电压跟随器、整流电路和真有效值转换电路，真有效值转换电路的输出端连接控制单元；分压电路包括串联在dbd激励电源高压端与接地端之间的多个分压电阻，其中一个分压电阻r2的阻值为所有分压电阻阻值之和的1/1000；该分压电阻r2的电压输出端通过电容c1连接电压跟随器的输入端。电容c1用于消除直流分量。供电电压检测电路最终检测到的电压值即为分压电阻r2上的电压，由于r2的阻值为事先设置，因此根据其上电压与供电电压的关系可求供电电压。
优选的，放电电流检测电路包括依次连接的电流互感器、电压跟随器、信号放大器、积分电路、精密整流电路和真有效值转换电路，真有效值转换电路的输出端连接控制单元。
优选的，放电电流检测电路中，电流互感器为均匀密绕在一环形非磁性骨架上的空心螺线管，dbd感应器的地线垂直通过电流互感器线圈所在平面；
电流互感器输出至电压跟随器的输出电压e(t)为：
其中，i是流经放电回路(地线)的电流，t是时间，h为环形非磁性骨架高度，n为线圈匝数，μ0为真空磁导率，rb、ra分别为环形非磁性骨架内径和外径；
本发明放电电流检测电路的工作原理为：先用电流互感器感应系统中的放电电流，将电流信号转化为电压信号，然后通过电压跟随器对干扰进行隔离，再对电压信号进行放大，将放大的信号进行积分，信号经过积分电路后转变为一个与感应到的电流信号呈正比的电压信号，经精密整流电路和真有效值转换电路输入到内置adc模块的mcu接口，实现对数据的采集和处理。本发明中的供电电压、放电电流、光谱检测等信号均属于频率较高的信号，相应的信号放大和调整电路的运算放大器均采用电流型运算放大器。
优选的，积分电压检测电路包括依次连接的分压电路、电压跟随器、整流电路和真有效值转换电路；分压电路包括并联在积分电容cm上的多个分压电阻，其中一个分压电阻r8的阻值为所有分压电阻阻值之和的1/10；该分压电阻r8的电压输出端通过电容c4连接电压跟随器的输入端。
积分电压检测电路的电压检测原理与供电电压检测电路的检测原理相同，由于积分电容两端电压在几十伏范围之内，因此电阻分压部分只需将积分电压衰减10倍，即使得分压电阻r8两端的电压是积分电容电压的1/10。
优选的，光谱检测单元包括依次连接的紫外传感器、电压跟随器、隔离电容、放大电路、电压跟随器、整流电路和真有效值转换电路；紫外传感器的信号接收端上设有对除·oh光谱之外的光谱进行滤除的紫外带通滤光片。
本发明还公开上述pdm工作模式下的高压激励电源驱动体积dbd空气净化装置的供电频率调控系统和方法的供电频率调控方法，包括步骤：
s1，设定pdm电源的初始供电频率，按照初始供电频率控制体积dbd反应器开始运行；
s2，在体积dbd反应器的运行过程中，检测体积dbd反应器的供电电压、放电电流、积分电容电压和发射光谱强度信号；
s3，根据检测到的信号计算体积dbd反应器的供电能量和相对光量子产率能效比eer；
s4，利用牛顿爬山算法确定ddb反应器的最佳eer、eer的最佳波动范围以及供电频率的调节范围，包括：
s401，将供电频率的变化作为扰动变量，多次控制pdm电源改变供电频率输出，并利用步骤s2至s3，得到每次供电频率改变后对应的相对光量子产率能效比；
s402，设定eer的最佳波动范围的区间长度值
s403，将最佳eer波动范围对应的供电频率范围作为供电频率的调节范围
s5，控制pdm电源的供电频率输出在供电频率的调节范围内。
优选的，步骤s4包括步骤：
s411，将当前eer及其对应的供电频率作为参考点a，控制pdm电源以一设定步长改变供电频率输出，作为参考点a的扰动量，使得eer增加；计算供电频率改变后的eer，将该eer及其对应的供电频率作为第一参考点；
s412，控制pdm电源以设定步长沿步骤s411中供电频率的改变方向，再次改变供电频率输出，作为第二参考点的扰动量；计算供电频率改变后的eer，若eer减少则转至步骤s413，若eer增加，则重复步骤s412，直至eer减少，转至步骤s413；
s413，定义使得eer减少的扰动量为加在参考点m上的扰动量，eer增加到参考点m之前的参考点为参考点c，eer减少后的参考点为参考点d；以设定步长沿步骤s422中供电频率的相反改变方向，改变供电频率的输出，使得eer回到参考点m；
s414，设定eer的最佳波动范围的区间长度值
s415，取区间(c，m)对应的eer范围的中点值作为参考点p1，区间(m，d)对应的eer范围的中点值作为参考点p2；比较参考点p1和参考点p2对应eer值的大小，若p1的eer值大于p2的eer值，eer的最大值在区间(p1，m)内，舍去区间(m，p2)；若p1的eer值小于p2的eer值，则eer的最大值在区间(m，p2)内，舍去区间(p1，m)；若p1与p2的eer值相等，则eer最大值在区间(p1，p2)内；转至步骤s416；
s416，将当前包含eer最大值的参考点区间长度值与
有益效果
本发明利用pdm电源作为体积dbd反应器的激励电源，通过可编程电源控制pdm电源的供电电压。在供电电压不变的情况下，本发明通过检测dbd运行过程中的放电参数，计算体积dbd反应器的供电能量和相对光量子产率能效比，并利用牛顿爬山算法确定最佳eer波动范围及其对应的供电频率调节范围，进而通过控制供电频率的变化，使得体积dbd反应器始终工作在最佳eer范围内，保障较高较稳定的供电能量和供电电压，从而使得dbd的反应效率稳定并保持较高水平。
附图说明
图1所示为本发明调控系统原理示意框图；
图2所示为体积dbd反应器的结构示意图；
图3所示为体积dbd反应器的剖面结构示意图；
图4所示为本发明一种实施例的数据采集与控制功能实现示意图；
图5所示为供电电压检测电路结构示意图；
图6所示为放电电流检测电路结构示意图；
图7所示为积分电压检测电路结构示意图；
图8所示为光谱强度检测电路结构示意图；
图9所示为功率密度调制电源供电波形示意图；
图10所示为牛顿爬山算法示意图；
图11所示为供电频率调节原理图。
具体实施方式
以下结合附图和具体实施例进一步描述。
实施例1
参考图1所示，本实施例pdm工作模式下的高压激励电源驱动体积dbd空气净化装置的供电频率调控系统和方法，包括控制单元(数据采集与控制)、可编程电源(可编程ac/dc电源)、pdm电源、体积dbd反应器、电气参数检测单元和光谱检测单元；
图2和图3所示为体积dbd反应器的结构示意图。反应器具体结构如图3所示，由两个电介质阻挡层，固定于电介质阻挡层两侧的电极，固定两块电介质阻挡层和电极的固定块，形成两块电介质阻挡层与固定块之间的放电通道。作为一种优选实施例，电介质阻挡层为石英电介质阻挡层，电极为铜板电极。图2中，l1为介质阻挡层的长度，约为180到200mm，l2为电极的长度，约为160到180mm，d1为介质阻挡层的宽度，约为40到60mm，d2为电极的宽度，约为60到80mm，h为反应器的高度，约为2到7mm。
结合图1，本实施例中，体积dbd反应器的地线上串接有积分电容cm；电气参数检测单元包括供电电压检测电路、放电电流检测电路和积分电压检测电路；供电电压检测电路检测体积dbd反应器的激励电源输入端的供电电压信号，放电电流检测电路检测体积dbd反应器的放电电流信号，积分电压检测电路检测积分电容cm上的积分电压信号，光谱检测单元检测体积dbd反应器反应区域的发射光谱强度信号，分别传输至控制单元；
可编程电源控制pdm电源的供电电压峰值输出，pdm电源的电压输出端连接体积dbd反应器的激励电压输入端；
控制单元根据接收到的供电电压信号、放电电流信号、积分电压信号和发射光谱强度信号，计算体积dbd反应器的供电能量和相对光量子产率能效比，进而根据计算结果控制pdm电源对体积dbd反应器的供电频率，使得体积dbd反应器在pdm电源的激励下放电，并工作在设定的最佳相对光量子产率能效比波动范围内。
本发明所述设定的相对光量子产率能效比波动范围即最佳eer附近的范围，最佳eer波动范围的端点值利用牛顿爬山算法得到，波动范围的区间长度值为设定值。
体积dbd反应器包括进气口和出气口，进气口上设有气体流速调节单元和气体流速传感器，出气口上设有气体质量检测单元；气体流速调节单元的控制输入端连接控制单元，气体流速传感器检测进气口的气体流速信号，气体质量检测单元检测出气口的气体质量信号，分别传输至控制单元；控制单元根据接收到的气体流速信号，通过气体流速调节单元对气体流速进行调节。
气体流速调节单元可采用电磁阀，控制器芯片通过调节电磁阀的开度控制气体流速，为现有技术。气体流速调节单元可配合气体流速传感器实现流速的反馈控制，气体流速传感器设置于气体流速调节单元之后的dbd进气口内。气体质量检测单元可采用现有用于检测空气质量的传感器。
本实施例中，控制单元通过气体流速调节单元控制体积dbd反应器内的初始气体流速为2.5m/s；dbd的运行过程中，当检测到气体流速在2.0m/s～3.5m/s范围之外时，对气体流速进行调节。
参考图4所示，控制单元可采用现有集数据采集与控制功能于一体的微控制单元芯片，如stc公司的stc12c、stm32、stc89系列单片机芯片。推荐使用封装为sop-20，单片机内置了8位adc和通用i/o口，速度均可达到100khz，8路adc模块可以用作电压检测、电流检测、光谱检测等。信号采集并传送到mcu的adc接口，之后mcu根据内部设定的模型判断放电装置是否工作在最佳放电状态，从而控制整个系统运行时供电频率输出。
为了实现使用过程的可视化，本实施例系统还设置显示单元，显示单元的输入端连接控制器，用于各种检测数据以及计算结果数据的显示。显示单元可采用现有显示屏，如lcd1602，lcd12864，lcd16864，lcd12232等型号，显示屏与mcu的i/o口连接。
参考图4所示的端口连接，本发明对供电频率调节的原理为：pdm电源内部是通过图11所示的电路的输出方波信号改变供电周期控制信号的频率。施密特触发器u10用于脉冲波形的产生，mcu产生直流信号vdac作为施密特触发器的直流偏置，可以改变电容ct的充放电时间，从而改变施密特触发器生成的方波频率。将施密特触发器生成的脉冲波形通过d触发器u11产生一个稳定的占空比为50％的方波。voc口输出的信号即为pdm电源内部高压逆变电路的控制信号，从而改变pdm电源输出的供电电压频率。所以，改变vdac，供电周期控制信号的频率发生改变，就可以调节pdm高压激励电源输出供电频率，达到改变电路的匹配状态，改变放电效果的目的。
供电电压检测电路
参考图5所示，本实施例供电电压检测电路包括依次连接的分压电路、电压跟随器、整流电路和真有效值转换电路，真有效值转换电路的输出端连接控制器；分压电路包括串联在dbd激励电源高压端与接地端之间的多个分压电阻，其中一个分压电阻r2的阻值远小于其它分压电阻的阻值之和；该分压电阻r2的电压输出端通过电容c1连接电压跟随器的输入端。电容c1用于消除直流分量。供电电压检测电路最终检测到的电压值即为分压电阻r2上的电压，由于r2的阻值为事先设置，因此根据其上电压与供电电压的关系可求供电电压。
供电电压检测检测电路中，分压电路还包括保护tvs管，保护tvs管并接在分压点与接地端之间。可导入大电压，吸收浪涌功率，防止出现瞬时高压尖峰通过分压电路，造成检测电路短路，损坏元器件情况，从而保护测试电路。
分压电阻r2的阻值为供电电压的1/1000。即r1/r2＝999/1，一方面使得供电电压信号衰减1000倍，另一方面方便后续的计算。
参考图5所示，由于pdm电源输出的供电电压峰-峰值在20到40kv，不方便直接将供电电压信号送入示波器中读取。因此，本发明采用的测量供电电压的方法是电阻分压器法，即将高电压波形转换成低电压波形，分压电路由高压臂和低压臂组成，r1、r2分别代表高电压臂和低电压臂的电阻。反应器供电端电压加在整个分压电路上，而输出电压则取自其中的低压电阻r2，在此，我们将供电电压信号衰减1000倍，即r2两端的电压是反应器供电端电压的1/1000。将r2两端的电压信号通过电容c1，起到消除直流分量的作用。另外，由于该电路检测的是高压信号，为了防止出现瞬时高电压造成检测电路元器件损坏，在c1后端并联一个tvs管d1，用于导入大电压，吸收浪涌功率，保护测试电路。然后将采集到的信号送入电压跟随器u1，电压跟随器的作用为：使r2两端的电压信号不会受后级测量系统的阻抗影响，所以后级测量系统应为高阻抗状态，由于电压跟随器具有输入高输入阻抗，低输出阻抗的特点，在采集电路与测量系统之间接入电压跟随器，一方面可以对采集电路的影响，另一方面也降低了对测量系统的影响。采集到的电压信号经过精密整流电路进行整流，再通过真有效值(rms)转换电路对交流信号进行真有效值(rms)转换。真有效值转换芯片可选择ad637、ad736、ad737等。之后将交流信号变为直流信号，然后输入到图4所示的内置adc模块的mcu的gp1/adc1，进行信号的处理。
图5中，假定被测电压为u，r2两端电压为ur2。根据电流连续性原理，被测电压为：
放电电流检测电路
参考图6所示，本实施例放电电流检测电路包括依次连接的电流互感器、电压跟随器、信号放大器、积分电路、精密整流电路和真有效值转换电路，真有效值转换电路的输出端连接控制器。
放电电流检测电路中的电压跟随器都采用电流型运算放大器。
放电电流检测电路中，电流互感器为均匀密绕在一环形非磁性骨架上的空心螺线管，dbd感应器的地线垂直通过电流互感器线圈所在平面；电流互感器的输出电压正比于被测电流随时间的变化率，因此电流互感器输出至电压跟随器的输出电压e(t)为：
其中，i是流经放电回路(地线)的电流，t是时间，h为环形非磁性骨架高度，n为线圈匝数，μ0为真空磁导率，rb、ra分别为环形非磁性骨架内径和外径；
由于电流互感器的输出电压与被测电流的微分成正比，所以需要对其输出的电压信号进行积分处理才能得到与回路电流成正比的输出电压。电流互感器输出的电压信号先经过一个电压跟随器u2再进行后续处理，电压跟随器的作用是降低后续电路对采集信号的干扰。另外，由于互感器感应到的电压信号很小，因此将其输出的电压信号先经过一个放大电路(包括u3)进行信号的放大，再经过积分电路(包括u4以及c2)。信号经过积分电路后已转变为一个与感应到的电流信号呈正比的电压信号，再经整流电路和真有效值转换电路输入到图4所示内置adc模块的mcu的gp3/adc3接口，实现对数据的采集和处理。
本发明放电电流检测电路的工作原理为：先用电流互感器感应系统中的放电电流，将电流信号转化为电压信号，然后通过电流型电压跟随器对干扰进行隔离，再对电压信号进行放大，将放大的信号进行积分，信号经过积分电路后转变为一个与感应到的电流信号呈正比的电压信号，再经精密整流电路和真有效值转换电路输入到图4所示内置adc模块的mcu的gp3/adc3接口，实现对数据的采集和处理。
图6中，u4与周围器件构成的电路起到对电流互感器输出电压进行积分的作用，此电路中，积分电路原理及公式如下：
u4o＝-uc(4)
u4o为运放u4的输入端电压，uc为积分电容c2上的积分电压；
将(3)代入(5)可得积分电路的输出电压为：
根据公式(6)，电流互感器输出的电压信号经过积分电路后，得到的电压信号与其感应到的电流信号呈一次函数关系。将此信号送到mcu中，可以设置感应电流与电压的关系，进而得到放电回路中的放电电流。
积分电压检测电路
参考图7所示，积分电压检测电路包括依次连接的分压电路、电压跟随器、整流电路和真有效值转换电路；分压电路包括并联在积分电容cm上的多个分压电阻，其中一个分压电阻r8的阻值远小于其它分压电阻的阻值之和；该分压电阻r8的电压输出端通过电容c4连接电压跟随器的输入端。再经整流电路和真有效值转换电路输入到图2所示内置adc模块的mcu的gp2/adc3接口，实现对数据的采集和处理。
积分电压检测电路中的电压跟随器都采用电流型运算放大器。
积分电压检测电路的电压检测原理与供电电压检测电路的检测原理相同，由于积分电容两端电压在几十伏范围之内，因此电阻分压部分只需将积分电压衰减10倍，即使得分压电阻r8两端的电压是积分电容电压的1/10。积分电压检测电路的电压检测原理与供电电压检测电路的检测原理相同，由于积分电容两端电压在几十伏范围之内，因此电阻分压部分只需将积分电压衰减10倍，即使得分压电阻r8两端的电压是积分电容电压的1/10。
光谱检测单元
参考图8所示，光谱检测单元包括依次连接的紫外传感器、电压跟随器、隔离电容、放大电路、电压跟随器、整流电路和真有效值转换电路；紫外传感器的信号接收端上设有对除·oh光谱之外的光谱进行滤除的紫外带通滤光片。
光谱检测单元中的电压跟随器采用电流型运算放大器。
由于本发明使用pdm高压电源驱动介质阻挡放电结构，在高压电场的作用下，物质吸收能量发生电离。气体中的原子或分子与电子之间的碰撞，原子与原子之间碰撞，光与原子或分子的相互作用等，会产生各种各样的放电生成物，包括高能电子、正负离子、自由基、臭氧(o3)和紫外线(uv)，以及少量的氮氧化物(nox)等，使废气中的有机物分子、重金属元素、病菌等有毒害物种，发生强氧化反应，最终达到降解、杀菌、消毒等目的。
气体在介质阻挡放电的作用下吸收电能，发生电子跃迁，产生的激发态的活性物种，活性物种在由激发态变回基态的过程中会发出特定波长的光，因此产生活性物种的浓度可以通过发射光谱的强度来表示。在本装置的放电中产生紫外辐射，紫外线的波长范围为180-400nm，具体分为uvc辐射(180-280nm)、uvb辐射(280-320nm)、uva辐射(320-400nm)。由于人眼看不见紫外辐射，放电产生的紫外线等原子和分析光谱成分和辐射强度需要通过专门的光学探测器件测量。使用半导体紫外光学探测器，可以将光学信号的强度通过电压信号显现出来，可选用以下型号传感器，如gs-ab-0603e(检测范围300-370nm)、gs-abc-2835(检测范围210-370nm)、gt-uvv-l(检测范围200-440nm)等，在dbd通道中观察到一系列光谱，其来自一些活性物种如羟基(·oh)、氧自由基(·o)和氮氧化物(nox)。
本发明重点关注放电产生的·oh光谱(309-318nm)，因此为了测量放电区域的活性物质种类和浓度，本发明利用紫外传感器探测反应器放电区域的活性物中的光谱强度，将光学信号转化为电信号，用紫外传感器输出的电压信号表示活性物种产生的相对浓度。由于介质阻挡放电过程中产生的活性物种会发出多种波长的光谱，而本发明中主要关注·oh的光谱，因此在紫外传感器前放置滤光片对光谱进行初步滤除。紫外带通型滤光片可以选择et313/25bp(中心波长313nm，带宽25nm)，et325/20bp(中心波长325nm，带宽20nm)，xbpa310，zbpa310(中心波长310nm，带宽10nm)。
将紫外传感器的电压信号通过电压跟随器u7进行隔离，由于紫外传感器输出的信号的频率较高，所以u7使用响应时间较快的电流型运算放大器用作电压跟随器。由于后级运算放大器的输入阻抗无穷大，造成电荷的积累，因此需要与大地之间接一个电阻r12作为能的释放通道，为了防止高压振荡，r12与r13进行匹配，阻值相等。紫外传感器输出的电压信号含有直流分量，如果将其直接通过后级的放大电路，可能会造成放大饱和失真，因此使用电容c6隔离传感器输出信号中的直流分量。下一步使用反相放大电路u8对取得的信号幅值进行放大。之后经过精密整流电路对将交流信号变为直流信号，将真有效值转换后的信号输入如图4所示的内置adc模块的mcu的gp0/adc0接口，实现对数据进行采集和显示。
放大电路u8对应的方程式为：
式中，v表示紫外传感器输出的电压信号经过电压跟随器的电压值，u8o表示放大电路输出的电压值。
实施例2
本发明还公开上述pdm工作模式下的高压激励电源驱动体积dbd空气净化装置的供电频率调控系统和方法的供电频率调控方法，包括步骤：
s1，设定pdm电源的初始供电频率，按照初始供电频率控制体积dbd反应器开始运行；
s2，在体积dbd反应器的运行过程中，检测体积dbd反应器的供电电压、放电电流、积分电容电压和发射光谱强度信号；
s3，根据检测到的信号计算体积dbd反应器的供电能量和相对光量子产率能效比eer；
s4，利用牛顿爬山算法确定ddb反应器的最佳eer、eer的最佳波动范围以及供电频率的调节范围，包括：
s401，将供电频率的变化作为扰动变量，多次控制pdm电源改变供电频率输出，并利用步骤s2至s3，得到每次供电频率改变后对应的相对光量子产率能效比；
s402，设定eer的最佳波动范围的区间长度值
s403，将最佳eer波动范围对应的供电频率范围作为供电频率的调节范围
s5，控制pdm电源的供电频率输出在供电频率的调节范围内。
优选的，步骤s4包括步骤：
s411，将当前eer及其对应的供电频率作为参考点a，控制pdm电源以一设定步长改变供电频率输出，作为参考点a的扰动量，使得eer增加；计算供电频率改变后的eer，将该eer及其对应的供电频率作为第一参考点；
s412，控制pdm电源以设定步长沿步骤s411中供电频率的改变方向，再次改变供电频率输出，作为第二参考点的扰动量；计算供电频率改变后的eer，若eer减少则转至步骤s413，若eer增加，则重复步骤s412，直至eer减少，转至步骤s413；
s413，定义使得eer减少的扰动量为加在参考点m上的扰动量，eer增加到参考点m之前的参考点为参考点c，eer减少后的参考点为参考点d；以设定步长沿步骤s422中供电频率的相反改变方向，改变供电频率的输出，使得eer回到参考点m；
s414，设定eer的最佳波动范围的区间长度值
s415，取区间(c，m)对应的eer范围的中点值作为参考点p1，区间(m，d)对应的eer范围的中点值作为参考点p2；比较参考点p1和参考点p2对应eer值的大小，若p1的eer值大于p2的eer值，eer的最大值在区间(p1，m)内，舍去区间(m，p2)；若p1的eer值小于p2的eer值，则eer的最大值在区间(m，p2)内，舍去区间(p1，m)；若p1与p2的eer值相等，则eer最大值在区间(p1，p2)内；转至步骤s416；
s416，将当前包含eer最大值的参考点区间长度值与
本发明的调控思想和调控标准为：整个系统中，随着放电条件的改变，放电特性会有很大不同，因此我们要通过调节的放电条件使整个系统工作在最佳的放电状态。随着放电条件的改变，活性物种产生的浓度、反应器供电能量、相对光量子产率能效比会跟着改变。这些参量随着放电条件的改变都具有变化规律，随着放电强度的增加，活性物种的浓度和反应器的供电能量都相应增加，但是并不意味着相对光量子产率能效比也是这个趋势。根据放电过程中的因变量(活性物种的浓度、反应器供电能量、相对光量子产率能效比)随自变量(pdm电源的供电频率)的变化情况，设计一种评估方法，判断在何种放电条件下得到最佳的放电效果。
在系统运行过程中，通过调节pdm电源的供电频率，可以改变激励电源与反应器的匹配效果，因此，反应器的供电电压峰值，放电强度等会发生改变，影响放电效果。根据反应器的供电能量和放电产生活性物种的光谱强度，获得此时反应器相对光量子产率能效比(eer)，根据牛顿爬山算法改变供电频率，使反应器的工作条件在最佳eer附近震荡。
结合图9所示的典型供电波形图可以看出，pdm电源是通过改变一个调功周期内的放电时间ton,j来调整功率密度，因此采用以下计算方法获得供电能量：
公式(8)是单个供电周期供电能量的计算公式，采用李萨如图形算法时，可表示为积分电容和李萨如图形面积的乘积即ed,i＝cmsd,i，cm为积分电容的电容值，sd,i为相应供电周期对应的李萨如图形面积。其中，ed,i是单个供电周期的供电能量，td,i表示单个供电周期的时间，u(t)是供电电压，i(t)是供电电流，此外，下标d,i表示供电周期的次数。
pdm高压激励电源的供电频率f表示为：
一个调功周期tm,j中供电时间的占空比是供电持续时间与调功周期的比值，通过下式表示：
下标m,j表示功率密度调整的周期次数，占空比调整范围为dj＝0～1。
当dbd系统工作在连续状态时，运行时间内总供电周期个数可以表示为：
tt是放电系统的运行总时间，其中ndm,j是第j个调功周期的供电周期个数，
系统运行时间内的总供电能量et表示如下：
pdm高压激励电源与反应器的匹配关系如下：
pdm电源和体积dbd反应器组成串联谐振回路。假设高压激励电源中的谐振电感值为l，体积dbd反应器的等效电容值为c，回路的阻抗z为：
其中r为回路的电阻值，ωl为电源中的感抗，
由于：
ω0＝2πf0(15)
f0表示谐振频率:
在介质阻挡放电中，随着放电过程的进行，体积dbd反应器的等效电容会发生改变。由于容性的反应器与感性的激励电源之间存在匹配关系，当反应器的等效电容发生变化时，反应器与电源的匹配特性下降。在pdm电源的输入电压相同的情况下，pdm电源产生的高压激励电压峰值会减小，即原本注入到反应器上的供电电压减小，而只有当电源的供电频率和电极的谐振频率越接近，放电效果才越好。此时就需要mcu调节pdm电源的供电频率，从而提高体积dbd反应器的供电电压，提高反应器的放电强度。
根据式(16)：当反应器等效电容c发生改变时，系统的固有谐振频率f0发生改变，此时需调节电源的供电频率，当pdm电源的供电频率与系统的谐振频率相同时，供电电压最大，供电能量最大。为了评估是否达到谐振匹配，定义pdm高压激励电源输出端能量与输入端能量的比值为pdm高压激励电源的能量效率p＝et/ein。其中，ein是可编程电源提供给pdm电源输入端的能量，et是pdm电源提供给体积dbd反应器的供电能量。当能量效率p最高时，系统达到谐振匹配。在本系统中，供电频率的改变可以直接通过td,i的变化观察到，因此，当c发生改变时，通过图4所示的gp4和gp5端口调控图11所示vdac接口的信号，即可改变pdm高压激励电源的供电频率，即改变td,i，就可以改变电源的匹配效果，从而改变放电状态。
相对光量子产率能效比计算
本发明利用光谱检测电路对体积dbd反应器的光谱进行分析，以紫外光电传感器探测到的·oh的相对光谱强度为标准。因此，来自放电区域的活性物质的相对光量子产率能效比(eer)可以通过下式求出：
其中i是放电区域内·oh发射光谱的相对强度，em,j为单个调功周期的供电能量：
等离子体区域中活性物质的相对浓度与其发射光谱的相对强度呈正相关。因此，可以通过相对光量子产率的eer来评价反应物质产率。
相对光量子产率随着供电能量的改变而改变，供电能量与供电电压有关。供电电压幅值是随着pdm电源与反应器的匹配特性改变的，pdm电源供电周期td,i决定了系统的匹配特性。本发明中，根据相对光量子产率能效比，确定供电频率的最佳调节范围。
本发明通过牛顿爬山算法得到最佳的eer波动范围和供电频率范围，通过不断调节放电反应系统的供电频率来比较调整前后dbd产生eer的变化情况，再根据eer的变化来调整td,i，即供电频率。
如图10所示为一种实施例的牛顿爬山算法的xx-eer图，横轴表示一个调功周期的供电能量em,j，纵轴表示相对光量子产率eer。牛顿爬山法又称扰动观察法，本文中是通过不断调节放电反应系统的供电频率来比较调整前后dbd产生eer的变化情况，再根据变化情况来调整供电频率，使放电反应器工作在最佳能效比波动范围内。
图10所示的牛顿爬山算法具体步骤为：
(1)在a点加一个扰动变量，改变供电频率，使反应器的eer达到b点；
(2)此时检测到供电能量的提高使得反应器的相对光量子产率能效比增加，则继续原来的方向增加扰动变量，使反应器工作在c点；
(3)继续原来的方向加扰动变量，使反应器工作在m点；
(4)继续原来的方向加扰动变量，使反应器工作在d点；
(5)此时检测到(4)加的扰动变量使得反应器的相对光量子产率能效比减小，则改变原来的方向继续加扰动变量，使反应器的eer重新达到m点；反应器可在c点，m点，d点三个工作点间波动。
上述扰动量的改变方向即增加或者减少。c点、d点对应的供电频率可作为供电频率调节范围的端点值，区间(c，d)可成为较佳调节范围。
在牛顿爬山算法中，根据扰动变量的步长所确定的m点不一定是eer的最高点，在确定了较佳参数调节范围c到d后，再用优选法对调节范围进行重新设定，根据下列步骤找出eer的最高点：
(1)在(c,m)区间内取中点p1(eer值的中间值对应的点)，(m,d)区间内取中点p2；
(2)当p1对应的eer值大于p2对应的eer值时，相对光量子产率能效比的极大值在(p1,m)区间内，(m,p2)的区间舍去；
(3)反之，极大值在(m,p2)区间内，舍去(p1,m)区间；
(4)当p1和p2对应的函数值相等时，极大值在(p1,p2)范围内，舍去(c,p1)和(p2,d)的区间；
(5)在剩余的区间内重新取中点，找到p3，p4，……，pn，pn+1。以步骤(1)到(4)的方式进行迭代计算；直到剩余区间范围小于设定值
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和变形，这些改进和变形也应视为本发明的保护范围。