CN105404350B - 一种基于二分梯度法的热电发电mppt方法 - Google Patents

一种基于二分梯度法的热电发电mppt方法 Download PDF

Info

Publication number
CN105404350B
CN105404350B CN201510658970.5A CN201510658970A CN105404350B CN 105404350 B CN105404350 B CN 105404350B CN 201510658970 A CN201510658970 A CN 201510658970A CN 105404350 B CN105404350 B CN 105404350B
Authority
CN
China
Prior art keywords
power
point
power generation
thermoelectric
gradient method
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201510658970.5A
Other languages
English (en)
Other versions
CN105404350A (zh
Inventor
全书海
房伟
谢长君
黄亮
翟鹏程
唐新峰
曾春年
陈启宏
石英
张立炎
邓坚
孙晓明
肖朋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuhan University of Technology WUT
Original Assignee
Wuhan University of Technology WUT
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Wuhan University of Technology WUT filed Critical Wuhan University of Technology WUT
Priority to CN201510658970.5A priority Critical patent/CN105404350B/zh
Publication of CN105404350A publication Critical patent/CN105404350A/zh
Application granted granted Critical
Publication of CN105404350B publication Critical patent/CN105404350B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Control Of Eletrric Generators (AREA)

Abstract

本发明公开了一种基于二分梯度法的热电发电MPPT(Maximum Power Point Tracking，简称MPPT)方法，旨在利用汽车尾气废热进行热电转换并回收利用，对提高能源利用效率具有重要意义。该方法首先利用电子负载，测得热电发电模块输出电压与输出功率的实验值；然后采用最小二乘法拟合出热电模块的输出功率与输出电压的函数曲线；选定初始隔离区间，使用二分法缩小拟合出的功率电压曲线的最大功率点隔离区间，再利用梯度法快速寻找缩小隔离区间后的函数最大值，并以该点作为热电发电模块的最大功率点。相比传统方法，本发明提出的方法在跟踪精度和速度上有较大提升。

Description

一种基于二分梯度法的热电发电MPPT方法
技术领域
本发明属于热电发电领域，更具体地，属于一种基于二分梯度法的热电发电MPPT(Maximum Power Point Tracking:最大功率跟踪)方法。
背景技术
当前最大功率点跟踪技术在光伏发电领域应用广泛，而利用塞贝克效应实现汽车尾气废热发电，同样需要在不同条件下实时跟踪温差发电的最大功率。目前最为常用的MPPT方法有开路电压比例系数法、干扰观测法、电导增量法等。
由热电转换模块的输出特性可知，当热电模块的开路电压在外界条件变化时，热电模块的最高功率点电压也近似成比例变化。因此认为热电模块的最高功率和热电模块的开路电压之间存在近似的线性关系。开路电压比例系数法虽然结构简单，但是这种线性关系只是一种近似，工作点并不是最大功率，而且比例系数的选取对于系统影响很大，不易控制。
干扰观测法是在原有输出上增加一个扰动，通过实时采样计算得到瞬时功率，与上一次的输出功率进行比较，根据功率增减变化动态调节输出功率，使其逐步逼近最大功率点。但是初始值及跟踪步长的给定对跟踪精度和跟踪速度的影响较大。
其它常用的MPPT方法，如双步长扰动观察法，迭代比较法，电流固定参数法，变步长电导增量法等，均在以上两种主要跟踪方法上有所改进。但是在汽车实时动态工况快速变化的条件下，这些传统常用的方法便出现了它们本身的不足之处，无法满足热电发电模块最大功率跟踪速度快，精度高的实际需求。
发明内容
本发明的目的是针对现有技术的不足，而提供一种基于二分梯度法的热电发电MPPT方法。
为了实现上述目的，本发明所采用的技术方案是：
一种基于二分梯度法的热电发电MPPT方法，包括以下步骤：
(1)选取热电模块实测电压值作为自变量，功率值作为因变量，在发动机常见运行条件下测取功率电压值；
(2)使用最小二乘法对步骤(1)中测取的热电模块的功率电压值进行曲线拟合以获取其内在函数关系；
(3)使用二分法缩小由步骤(2)所拟合出的功率电压曲线的最大功率点的隔离区间；
(4)利用梯度法快速寻找经过步骤(3)缩小了隔离区间后的函数的最大值，并以该点作为热电模块的最大功率点的近似值输出给系统。
所述步骤(1)具体包括以下步骤：
(1-1)在汽车尾气发电系统中将热电模块以全串联连接方式分两层放置在气箱的上下表面；
(1-2)运行发电系统，当发动机工作常见的一种工况下，当尾气气箱温度稳定后，测量并选取功率电压值；
(1-3)改变发动机工况另一种常见的工况，重复(1-2)。
所述步骤(2)具体包括以下子步骤：
(2-1)假设拟合的曲线方程式为：
①
其中：y-输出功率；x-输出电压；ai：-第i次项系数,i＝0,1,…m；下同；
(2-2)由假设的拟合曲线方程式①得到误差平方和为下式：
②
(2-3)结合根据多元函数求极值原理，有如下等式：
③
(2-4)将写成矩阵形式求出最小二乘解：
矩阵形式为：Aα＝Y；
其中α＝[a0 a1 … am]T；
Y＝[y1 y2 … ym]T；
得到最小二乘解α＝(ATA)-1ATY ④
(2-5)拟合出其功率-电压曲线。
所述步骤(3)具体包括以下步骤：
(3-1)：取初始隔离区间[X1,X2]，隔离区间精度值e，并计算其中值点Xmid对应的功率f(Xmid)；
(3-2)：取步长点间距为Step>0，步长点功率值主要用来判断二分点与最大功率点的相对位置，则步长点Xgrad＝Xmid+Step；并计算此时步长点所对应的功率值f(Xgrad)；
(3-3)：比较中点功率f(Xmid)和步长点功率值f(Xgrad)的大小：
若f(Xmid)＞f(Xgrad)，说明Xmid在最大功率点右端，令F＝f(Xmid)，则有：
X1＝X1，X2＝Xgrad，Xmid＝(X1+X2)/2，Fmax＝f(Xmid)；
若f(Xmid)＜f(Xgrad)，说明Xmid在最大功率点左端。令F＝f(Xmid)，则有：
X1＝Xmid，X2＝X2，Xmid＝(X1+X2)/2；Fmax＝f(Xmid)；
(3-4)：重复过程(3-1)，(3-2)，(3-3)，直至fabs(F-Fmax)＜e，f(Xmid)＞f(Xgrad)，则新的隔离区间为[X1,Xmid]；f(Xmid)＜f(Xgrad)，则新的隔离区间为[Xmid,X2]。
步骤(4)具体包括以下子步骤：
(4-1)：初始点为xk＝Xmid，终止误差error，令k＝0；
(4-2)：计算若则停止，输出x0；否则转(4-3)；
(4-3)：取
(4-4)：进行一维搜索，求tk，使得：
令xk+1＝xk+tkpk，k:＝k+1，转(4-2)。
通过本发明所提出的上述MPPT方法，与现有方法相比，本发明具有以下的有益效果：
1、最大输出功率点跟踪精度高：由于采用了步骤(1)、(2)，测取的数值来自不同的工况，从而更具有代表性，此外，在使用最小二乘法对热电系统的功率电压值进行曲线拟合以获取其内在函数关系过程中，利用了多元函数求极值原理，将误差平方和降到了最小，从而使得最大输出功率点跟踪精度相对于传统方法大大提高。
2、最大输出功率点跟踪速度快：由于采用了步骤(3)、(4)，有效的将传统的二分法和梯度法结合，使其各自的优点融合，首先二分法的使用使最大功率点的隔离区间大大缩小，然后利用梯度法快速寻找函数最大值，此时即使较小的步长也可以收敛很快，从而使得在汽车实时动态工况快速变化的条件下依旧能够快速的跟踪最大功率点。
附图说明
图1为本发明的流程图。
图2为汽车尾气发电系统试验台架示意图。
图3为本发明中二分梯度法示意图。
图4为各类方法误差率对比图。
图5为本发明3100rpm@61NM工况下输出功率与电压关系曲线图。
图6为本发明3300rpm@69NM工况下输出功率与电压关系曲线图。
图7为本发明3100rpm@61NM工况输出功率与电压拟合曲线。
图8为本发明3300rpm@69NM工况输出功率与电压拟合曲线。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
一种基于二分梯度法的热电发电MPPT方法包括以下步骤：
(1)选取热电模块实测电压值作为自变量，功率值作为因变量，在发动机工况为3100rpm@61NM和3300rpm@69NM测取功率电压值。
具体包括以下子步骤：
(1-1)在汽车尾气发电系统中将热电模块以全串联连接方式分两层放置在气箱的上下表面。为更详细的解释本发明，进行试验如下：搭建汽车尾气发电系统试验台架，其示意图如图2所示，将热电模块以全串联连接方式分两层放置在气箱的上下表面。
(1-2)运行发电系统，当发动机工作在3100rpm@61NM工况下，当气箱温度稳定后，测量并选取8组功率电压值，描绘输出功率与电压关系曲线图。
(1-3)改变发动机工况为3300rpm@69NM，重复(1-2)。如图5所示是3100rpm@61NM工况下输出功率与电压关系曲线图，图6是3300rpm@69NM工况下输出功率与电压关系曲线图。
(2)使用最小二乘法对步骤(1)中测取的热电系统的功率电压值进行曲线拟合以获取其内在函数关系。
具体包括以下子步骤：
(2-1)假设拟合的曲线方程式为：
(y:输出功率；x:输出电压；ai:i次项系数,i＝0,1,…m；下同。) ①
(2-2)由假设的拟合曲线方程式①可得到误差平方和为下式：
②
(2-3)结合根据多元函数求极值原理，有如下等式：
③
(2-4)将将写成求出最小二乘解：
矩阵形式为：Aα＝Y；
α＝[a0 a1 … am]T；
Y＝[y1 y2 … ym]T；
④
则最小二乘解为：α＝(ATA)-1ATY； ⑤
(2-5)拟合出其功率-电压曲线：
如步骤(1-1)中所述汽车尾气发电系统试验台架，试验中，工况3100rpm@61NM时，如图5所示测得的电压与功率关系近似为多项式曲线，使用三次多项式进行曲线拟合：，拟合曲线如图7所示，功率电压曲线基本上为三次多项式形式。同样，当转速3300rpm@69NM时，这里也采用三次多项式作为其主曲线拟合形式，拟合曲线如图8所示。比较图7及图8，发现不同工况下，拟合曲线均为单峰曲线。因此可以通过求取其功率曲线的最高点作为其最大功率点。
(3)使用二分法缩小由步骤(2)所拟合出的功率电压曲线的最大功率点的隔离区间：
具体包括以下子步骤：
(3-1)：取初始隔离区间[X1,X2]，隔离区间精度值e，并计算其中值点Xmid对应的功率f(Xmid)；
(3-2)：取步长点间距为Step>0，步长点功率值主要用来判断二分点与最大功率点的相对位置。则步长点Xgrad＝Xmid+Step；并计算此时步长点所对应的功率值f(Xgrad)；
(3-3)：比较中点功率f(Xmid)和步长点功率值f(Xgrad)的大小：
若f(Xmid)＞f(Xgrad)，说明Xmid在最大功率点的右端。令F＝f(Xmid)，则有：
X1＝X1，X2＝Xgrad，Xmid＝(X1+X2)/2，
Fmax＝f(Xmid)；
若f(Xmid)＜f(Xgrad)，说明Xmid在最大功率点左端。令F＝f(Xmid)，则有：
X1＝Xmid，X2＝X2，Xmid＝(X1+X2)/2；
Fmax＝f(Xmid)；
(3-4)：重复过程(2-1)，(2-2)，(2-3)，直至fabs(F-Fmax)＜e。f(Xmid)＞f(Xgrad)，则新的隔离区间为[X1,Xmid]；f(Xmid)＜f(Xgrad)，则新的隔离区间为[Xmid,X2]。
(4)利用梯度法快速寻找经过步骤(3)缩小了隔离区间后的函数的最大值，并以该点作为系统最大功率点的近似值输出给系统。
具体包括以下子步骤：
(4-1)：初始点为xk＝Xmid，终止误差error，令k:＝0；
(4-2)：计算若则停止，输出x0，否则转(4-3)；
(4-3)：取
(4-4)：进行一维搜索，求tk，使得：
令xk+1＝xk+tkpk，k:＝k+1，转(4-2)。
如步骤(1-1)中所述热电实验台架下分别在2700r/min，2900r/min，3100r/min，3200r/min，3300r/min的工况条件下测试了传统二分法、梯度法、扰动观察法，以及本发明提出二分梯度法的最大功率点跟踪效果。将以上各种方法与实际值进行对比，可以得到误差情况，各类方法的误差率对比图如图4所示，可以明显看出二分梯度法的误差率比现有方法低。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (5)

1.一种基于二分梯度法的热电发电MPPT方法，其特征在于，包括以下步骤：
(1)选取热电系统实测电压值作为自变量，功率值作为因变量，在发动机常见运行条件下测取功率电压值；
(2)使用最小二乘法对步骤(1)中测取的热电系统的功率电压值进行曲线拟合，以获取其内在函数关系；
(3)使用二分法缩小由步骤(2)所拟合出的功率电压曲线的最大功率点的隔离区间；
(4)利用梯度法快速寻找经过步骤(3)缩小了隔离区间后的函数的最大值，并以该点作为系统最大功率点的近似值输出给系统。
2.根据权利要求1所述的基于二分梯度法的热电发电MPPT方法，其特征在于，步骤(1)具体包括以下子步骤：
(1-1)在汽车尾气发电系统中将热电模块以全串联连接方式分两层放置在气箱的上下表面；
(1-2)运行发电系统，当发动机工作常见的一种工况下，当尾气气箱温度稳定后，测量并选取多组功率电压值；
(1-3)改变发动机工况另一种常见的工况，重复步骤(1-2)。
3.根据权利要求1所述的基于二分梯度法的热电发电MPPT方法，其特征在于，步骤(2)具体包括以下子步骤：
(2-1)首先设定拟合的曲线方程式；
(2-2)由设定的拟合曲线方程式可得到误差平方和；
(2-3)结合根据多元函数求极值原理，得到最小二乘解；
(24)拟合出其功率-电压曲线得出其方程式。
4.根据权利要求1所述的基于二分梯度法的热电发电MPPT方法，其特征在于，步骤(3)具体包括以下子步骤：
(3-1)：取初始隔离区间，隔离区间精度值e，并计算其中值点对应的功率；
(3-2)：取步长点间距大于零，步长点功率值主要用来判断二分点与最大功率点的相对位置；则步长点为中点值加步长点间距，计算此时步长点所对应的功率值；
(3-3)：比较中点功率和步长点功率值的大小，确定此时最大功率点的位置；
(3-4)：重复步骤(3-1)，(3-2)，(3-3)，直至功率值满足精度要求。
5.根据权利要求1所述的基于二分梯度法的热电发电MPPT方法，其特征在于，步骤(4)具体包括以下子步骤：
(4-1)：初始点，给定终止误差，令k＝0；
(4-2)：计算误差，若误差小于终止误差，则停止，输出此时的功率点，否则转(4-3)；
(4-3)：进行一维搜索，令k＝k+1，转步骤(4-2)。
CN201510658970.5A 2015-10-12 2015-10-12 一种基于二分梯度法的热电发电mppt方法 Expired - Fee Related CN105404350B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510658970.5A CN105404350B (zh) 2015-10-12 2015-10-12 一种基于二分梯度法的热电发电mppt方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510658970.5A CN105404350B (zh) 2015-10-12 2015-10-12 一种基于二分梯度法的热电发电mppt方法

Publications (2)

Publication Number Publication Date
CN105404350A CN105404350A (zh) 2016-03-16
CN105404350B true CN105404350B (zh) 2017-03-01

Family

ID=55469887

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510658970.5A Expired - Fee Related CN105404350B (zh) 2015-10-12 2015-10-12 一种基于二分梯度法的热电发电mppt方法

Country Status (1)

Country Link
CN (1) CN105404350B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106437915B (zh) * 2016-09-06 2018-03-30 山东省科学院自动化研究所 Mppt控制器及方法、控制装置及方法、低温余热发电系统
CN109343648B (zh) * 2018-09-28 2020-06-09 江苏大学 一种基于新型开路电压法的mppt在线控制系统
CN109814652B (zh) * 2019-03-19 2020-04-21 武汉理工大学 一种基于非等分多分法的热电发电mppt方法
CN113110689B (zh) * 2021-04-25 2022-07-22 江苏大学 一种考虑反向温差的最大功率点追踪控制器及控制方法
CN114489228B (zh) * 2022-01-26 2022-11-01 四川大学 一种基于改进pso算法的mppt装置及方法
CN114943105B (zh) * 2022-07-26 2022-11-29 武汉理工大学 含螺旋纽带的环形热电发电机性能优化方法
CN115344078B (zh) * 2022-08-15 2023-09-22 深圳市冲锋兔新能源科技有限公司 Mppt功率快速追踪算法
CN116627204B (zh) * 2023-04-12 2023-10-20 团车互联网信息服务(北京)有限公司 一种太阳能汽车多峰mppt追踪控制方法及系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN202406064U (zh) * 2012-01-09 2012-08-29 浙江尖山光电股份有限公司 一种基于dc-dc变换的自整定光伏电源优化器
CN103399491A (zh) * 2013-08-06 2013-11-20 清华大学 光伏发电系统光伏组件机理模型参数辨识方法
CN104182407A (zh) * 2013-05-23 2014-12-03 中国科学院软件研究所 一种缩小数据查找范围的数据处理系统

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012181610A (ja) * 2011-02-28 2012-09-20 Seiko Instruments Inc エネルギー変換装置および最大電力変換回路および最大電力変換方法
KR101371226B1 (ko) * 2012-01-05 2014-03-10 신옥테크(주) 이동형 태양광 발전모듈을 위한 최대 전력 추종 제어방법

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN202406064U (zh) * 2012-01-09 2012-08-29 浙江尖山光电股份有限公司 一种基于dc-dc变换的自整定光伏电源优化器
CN104182407A (zh) * 2013-05-23 2014-12-03 中国科学院软件研究所 一种缩小数据查找范围的数据处理系统
CN103399491A (zh) * 2013-08-06 2013-11-20 清华大学 光伏发电系统光伏组件机理模型参数辨识方法

Also Published As

Publication number Publication date
CN105404350A (zh) 2016-03-16

Similar Documents

Publication Publication Date Title
CN105404350B (zh) 一种基于二分梯度法的热电发电mppt方法
CN100498636C (zh) 具有自适应搜索算法的最大功率跟踪捕获光伏控制方法
CN104778352B (zh) 基于stft算子的七参数光伏电池输出特性建模方法
CN108365822B (zh) 光伏电池在不同光照温度条件下工作特性的预测方法及系统
CN109974570B (zh) 一种差分电感式位移传感器及测量方法
CN105573400A (zh) 一种光伏发电系统最大功率跟踪控制方法
CN106885944A (zh) 一种用于下垂控制的直流微网线路阻抗检测方法及装置
CN101833602B (zh) 一种硅太阳能电池工程用数学模型的建模方法
CN102902298B (zh) 基于分段模型的光伏阵列 mppt 控制器及控制方法
CN102263527A (zh) 一种光伏发电系统最大功率点跟踪方法
CN104298296A (zh) 一种燃料电池最大功率跟踪控制方法
CN103092249B (zh) 一种太阳能电池最大功率点跟踪方法
CN105259971B (zh) 一种优化的mppt算法
CN105048962A (zh) 一种光伏电池参数测试系统
CN107168450A (zh) 一种太阳能光伏发电最大功率点跟踪控制系统与方法
CN106055017A (zh) 一种基于最大功率点追踪的太阳能功率优化方法和装置
CN108763645B (zh) 光伏组件5参数模型的模型参数求解算法
CN109767353B (zh) 一种基于概率分布函数的光伏发电功率预测方法
CN102637056B (zh) 一种维持光伏发电系统最大功率点的方法
CN107423257B (zh) 一种太阳能电池最大功率点参数求解方法
CN106468935A (zh) 一种折线拟合求太阳能电池最大功率点方法
CN107544610A (zh) 一种基于mpp电压规律与梯度寻优的光伏mppt控制方法
CN115343618B (zh) 一种新能源电机的测试设备及测试方法
CN110309956A (zh) 基于内部等效参数的光伏组件剩余寿命预测方法
CN113708722A (zh) 一种基于llc拓扑光伏发电系统的mppt控制方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170301

Termination date: 20211012