一种基于二分梯度法的热电发电mppt方法
【技术领域】
[0001] 本发明属于热电发电领域，更具体地，属于一种基于二分梯度法的热电发电 MPPT(MaximumPowerPointTracking:最大功率跟踪）方法。
【背景技术】
[0002] 当前最大功率点跟踪技术在光伏发电领域应用广泛，而利用塞贝克效应实现汽车 尾气废热发电，同样需要在不同条件下实时跟踪温差发电的最大功率。目前最为常用的 MPPT方法有开路电压比例系数法、干扰观测法、电导增量法等。
[0003] 由热电转换模块的输出特性可知，当热电模块的开路电压在外界条件变化时，热 电模块的最高功率点电压也近似成比例变化。因此认为热电模块的最高功率和热电模块的 开路电压之间存在近似的线性关系。开路电压比例系数法虽然结构简单，但是这种线性关 系只是一种近似，工作点并不是最大功率，而且比例系数的选取对于系统影响很大，不易控 制。
[0004] 干扰观测法是在原有输出上增加一个扰动，通过实时采样计算得到瞬时功率，与 上一次的输出功率进行比较，根据功率增减变化动态调节输出功率，使其逐步逼近最大功 率点。但是初始值及跟踪步长的给定对跟踪精度和跟踪速度的影响较大。
[0005] 其它常用的MPPT方法，如双步长扰动观察法，迭代比较法，电流固定参数法，变步 长电导增量法等，均在以上两种主要跟踪方法上有所改进。但是在汽车实时动态工况快速 变化的条件下，这些传统常用的方法便出现了它们本身的不足之处，无法满足热电发电模 块最大功率跟踪速度快，精度高的实际需求。
【发明内容】
[0006] 本发明的目的是针对现有技术的不足，而提供一种基于二分梯度法的热电发电 MPPT方法。
[0007] 为了实现上述目的，本发明所采用的技术方案是：
[0008] -种基于二分梯度法的热电发电MPPT方法，包括以下步骤：
[0009] (1)选取热电模块实测电压值作为自变量，功率值作为因变量，在发动机常见运行 条件下测取功率电压值；
[0010] (2)使用最小二乘法对步骤（1)中测取的热电模块的功率电压值进行曲线拟合以 获取其内在函数关系；
[0011] (3)使用二分法缩小由步骤（2)所拟合出的功率电压曲线的最大功率点的隔离区 间；
[0012] (4)利用梯度法快速寻找经过步骤（3)缩小了隔离区间后的函数的最大值，并以 该点作为热电模块的最大功率点的近似值输出给系统。
[0013] 所述步骤（1)具体包括以下步骤：
[0014] (1-1)在汽车尾气发电系统中将热电模块以全串联连接方式分两层放置在气箱的 上下表面；
[0015] (1-2)运行发电系统，当发动机工作常见的一种工况下，当尾气气箱温度稳定后， 测量并选取功率电压值；
[0016] (1-3)改变发动机工况另一种常见的工况，重复（1-2)。
[0017] 所述步骤（2)具体包括以下子步骤：
[0018] (2-1)假设拟合的曲线方程式为：
[0020] 其中：y_输出功率；X-输出电压；ai:-第i次项系数，i= 0, 1，…m;下同；
[0021] (2-2)由假设的拟合曲线方程式①得到误差平方和为下式：
[0023] (2-3)结合根据多元函数求极值原理，有如下等式：
[0025] (2-4)将写成矩阵形式求出最小二乘解：
[0026] 矩阵形式为：Αα=Y;
[0027]其中α= [a。a!…am]T;
[0028] Y=hy2 …yJT;
[0030] 得到最小二乘解a= (ATA) 1TY④
[0031] (2-5)拟合出其功率-电压曲线。
[0032] 所述步骤（3)具体包括以下步骤：
[0033] (3-1):取初始隔离区间[Xi，X2]，隔离区间精度值e，并计算其中值点乂^对应的功 率D;
[0034] (3-2):取步长点间距为St印>0,步长点功率值主要用来判断二分点与最大功率 点的相对位置，则步长点Xgrad=X^+Step;并计算此时步长点所对应的功率值f(Xgrad);
[0035] (3-3):比较中点功率f(XJ和步长点功率值f(Xgrad)的大小：
[0036] 若fU>f(XgMd)，说明在最大功率点右端，令F=f(XJ，则有：
[0037] (Χ1+Χ2)/2，Ρ_=?·(Χ-);
[0038] 若fU<f(XgMd)，说明在最大功率点左端。令F=f(XJ，则有：
[0039]X1=Xnid,X2=X2,Xnid= (X^/2 ;Fnax=f(Xnid);
[0040] (3-4):重复过程（3-1)，（3-2)，（3-3)，直至fabs(F-F_) <e，f(X-) >f(Xgrad)， 则新的隔离区间为[ΙΧ^] ;f(XMd) <f(XgMd)，则新的隔离区间为[Χ"^，Χ2]。
[0041] 步骤（4)具体包括以下子步骤：
[0042] (4-1):初始点为xk=Xmid，终止误差error，令k=0 ;
[0043](4-2):计算ν/(χ,)，若 <?T〇r，则停止，输出x。;否则转（4-3);
[0044] (4-3):取凡=7/(4) ;
[0045] (4-4):进行一维搜索，求tk，使得：
[0048] 通过本发明所提出的上述MPPT方法，与现有方法相比，本发明具有以下的有益效 果：
[0049] 1、最大输出功率点跟踪精度高：由于采用了步骤（1)、（2)，测取的数值来自不同 的工况，从而更具有代表性，此外，在使用最小二乘法对热电系统的功率电压值进行曲线拟 合以获取其内在函数关系过程中，利用了多元函数求极值原理，将误差平方和降到了最小， 从而使得最大输出功率点跟踪精度相对于传统方法大大提高。
[0050] 2、最大输出功率点跟踪速度快：由于采用了步骤（3)、（4)，有效的将传统的二分 法和梯度法结合，使其各自的优点融合，首先二分法的使用使最大功率点的隔离区间大大 缩小，然后利用梯度法快速寻找函数最大值，此时即使较小的步长也可以收敛很快，从而使 得在汽车实时动态工况快速变化的条件下依旧能够快速的跟踪最大功率点。
【附图说明】
[0051] 图1为本发明的流程图。
[0052] 图2为汽车尾气发电系统试验台架示意图。
[0053] 图3为本发明中二分梯度法示意图。
[0054] 图4为各类方法误差率对比图。
[0055] 图5为本发明3100rpm@61匪工况下输出功率与电压关系曲线图。
[0056] 图6为本发明3300rpm@69匪工况下输出功率与电压关系曲线图。
[0057] 图7为本发明3100rpm@61匪工况输出功率与电压拟合曲线。
[0058] 图8为本发明3300rpm@69匪工况输出功率与电压拟合曲线。
【具体实施方式】
[0059] 为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实例，对本 发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不 用于限定本发明。
[0060] -种基于二分梯度法的热电发电MPPT方法包括以下步骤：
[0061] (1)选取热电模块实测电压值作为自变量，功率值作为因变量，在发动机工况为 3100rpm@61NM和 3300rpm@69NM测取功率电压值。
[0062] 具体包括以下子步骤：
[0063] (1-1)在汽车尾气发电系统中将热电模块以全串联连接方式分两层放置在气箱的 上下表面。为更详细的解释本发明，进行试验如下：搭建汽车尾气发电系统试验台架，其示 意图如图2所示，将热电模块以全串联连接方式分两层放置在气箱的上下表面。
[0064] (1-2)运行发电系统，当发动机工作在3100rpm@61NM工况下，当气箱