发明名称 一种日照温差采集样本概率密度的测定方法
摘要 本发明公开了一种日照温差采集样本概率密度的测定方法,包括如下步骤:步骤10):确定日照温差采集样本;步骤20):确定正温差样本和负温差样本;步骤30):确定正温差样本和负温差样本中各温差值对应的累加概率值;步骤40):对正温差样本的累加分布进行拟合;步骤50):对负温差样本的累加分布进行拟合;步骤60):确定正温差样本的概率密度参数;步骤70):确定负温差样本的概率密度参数;步骤80):确定日照温差采集样本的概率密度。该测定方法可以准确测定日照温差采集样本的概率密度。
申请公布号 CN103048056B 申请公布日期 2014.09.10
申请号 CN201210551566.4 申请日期 2012.12.18
申请人 东南大学 发明人 王高新;丁幼亮;宋永生
分类号 G01K3/10(2006.01)I 主分类号 G01K3/10(2006.01)I
代理机构 南京苏高专利商标事务所(普通合伙) 32204 代理人 柏尚春
主权项 一种日照温差采集样本概率密度的测定方法,其特征在于,该测定方法包括如下步骤:步骤(1):确定日照温差采集样本:将多个温度传感器配接到同一个温度采集系统中,多个温度传感器对不同测点的日照温度同时进行采集,获取不同测点的日照温度采集样本,对其中两个测点的日照温度采集样本在同一时刻的温度值相减,得到日照温差采集样本,日照温差采集样本包含不同时刻对应的温差值;步骤(2):确定正温差样本和负温差样本:对步骤(1)得到日照温差采集样本,温差值大于或等于0℃的为正温差值,正温差值对应的样本作为正温差样本,温差值小于0℃的为负温差值,负温差值对应的样本作为负温差样本;步骤(3):确定正温差样本和负温差样本中各温差值对应的累加概率值:利用式(1)对正温差样本的累加分布特性进行分析,确定正温差样本中各温差值对应的累加概率值;<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>P</mi><mrow><mo>(</mo><msup><mi>T</mi><mo>+</mo></msup><mo>&le;</mo><msup><mi>t</mi><mo>+</mo></msup><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>q</mi><mn>1</mn></msub><mrow><mo>(</mo><msup><mi>T</mi><mo>+</mo></msup><mo>&le;</mo><msup><mi>t</mi><mo>+</mo></msup><mo>)</mo></mrow></mrow><msub><mi>l</mi><mn>1</mn></msub></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000508240820000011.GIF" wi="1084" he="142" /></maths>式中,T<sup>+</sup>表示正温差变量,t<sup>+</sup>为正温差样本中某一温差值,P(T<sup>+</sup>≤t<sup>+</sup>)表示t<sup>+</sup>对应的累加概率值,q<sub>1</sub>(T<sup>+</sup>≤t<sup>+</sup>)表示正温差样本中小于等于t<sup>+</sup>的温差值个数,l<sub>1</sub>为正温差样本中温差值的总数;利用式(2)对负温差样本的累加分布特性进行分析,确定负温差样本中各温差值对应的累加概率值:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>P</mi><mrow><mo>(</mo><msup><mi>T</mi><mo>-</mo></msup><mo>&le;</mo><msup><mi>t</mi><mo>-</mo></msup><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>q</mi><mn>2</mn></msub><mrow><mo>(</mo><msup><mi>T</mi><mo>-</mo></msup><mo>&le;</mo><msup><mi>t</mi><mo>-</mo></msup><mo>)</mo></mrow></mrow><msub><mi>l</mi><mn>2</mn></msub></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000508240820000012.GIF" wi="1062" he="131" /></maths>式中,T<sup>‑</sup>表示负温差变量,t<sup>‑</sup>为负温差样本中某一温差值,P(T‑≤t<sup>‑</sup>)表示t<sup>‑</sup>对应的累加概率值,q<sub>2</sub>(T‑≤t<sup>‑</sup>)表示负温差样本中小于等于t<sup>‑</sup>的温差值个数,l<sub>2</sub>为负温差样本中温差值的总数;步骤(4):对正温差样本的累加分布进行拟合:对步骤(2)得到的正温差样本,利用式(3)对正温差样本的累加分布特性进行拟合,式(3)表达如下:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>F</mi><mrow><mo>(</mo><msup><mi>T</mi><mo>+</mo></msup><mo>)</mo></mrow><mo>=</mo><msub><mi>a</mi><mn>0</mn></msub><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>n</mi><mn>1</mn></msub></munderover><mo>[</mo><msub><mi>a</mi><mi>i</mi></msub><mi>cos</mi><mrow><mo>(</mo><mi>iw</mi><msup><mi>T</mi><mo>+</mo></msup><mo>)</mo></mrow><mo>+</mo><msub><mi>b</mi><mi>i</mi></msub><mi>sin</mi><mrow><mo>(</mo><mi>iw</mi><msup><mi>T</mi><mo>+</mo></msup><mo>)</mo></mrow><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000508240820000021.GIF" wi="1384" he="158" /></maths>式中,T<sup>+</sup>表示正温差变量,F(T<sup>+</sup>)表示正温差变量的累加分布拟合函数,n<sub>1</sub>为≥4的整数,a<sub>0</sub>表示F(T<sup>+</sup>)的常数项,a<sub>0</sub>、w、a<sub>i</sub>和b<sub>i</sub>为待估参数,其中i为整数,且i=1、2、…、n<sub>1</sub>;基于最小二乘法,利用正温差值及式(1)得到的各正温差值对应的累加概率值,对函数F(T<sup>+</sup>)进行拟合,得到待估参数a<sub>0</sub>、w、a<sub>i</sub>和b<sub>i</sub>;步骤(5):对负温差样本的累加分布进行拟合:对步骤(2)得到的负温差样本,利用式(4)对负温差样本的累加分布特性进行拟合,式(4)表达如下:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>F</mi><mrow><mo>(</mo><msup><mi>T</mi><mo>-</mo></msup><mo>)</mo></mrow><mo>=</mo><msub><mi>c</mi><mn>0</mn></msub><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>n</mi><mn>2</mn></msub></munderover><mo>[</mo><msub><mi>c</mi><mi>j</mi></msub><mi>cos</mi><mrow><mo>(</mo><mi>j&lambda;</mi><msup><mi>T</mi><mo>-</mo></msup><mo>)</mo></mrow><mo>+</mo><msub><mi>d</mi><mi>j</mi></msub><mi>sin</mi><mrow><mo>(</mo><mi>j&lambda;</mi><msup><mi>T</mi><mo>-</mo></msup><mo>)</mo></mrow><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000508240820000022.GIF" wi="1392" he="136" /></maths>式中,T<sup>‑</sup>表示负温差变量,F(T<sup>‑</sup>)表示负温差变量的累加分布拟合函数,n<sub>2</sub>为≥4的整数,c<sub>0</sub>表示F(T<sup>‑</sup>)的常数项,c<sub>0</sub>、λ、c<sub>j</sub>和d<sub>j</sub>为待估参数,其中j为整数,且j=1、2、…、n<sub>2</sub>;基于最小二乘法,利用负温差值及式(2)得到的各负温差值对应的累加概率值,对函数F(T<sup>‑</sup>)进行拟合,得到待估参数c<sub>0</sub>、λ、c<sub>j</sub>和d<sub>j</sub>;步骤(6):确定正温差样本的概率密度参数:利用式(5)对F(T<sup>+</sup>)求导,得到f(T<sup>+</sup>),将正温差样本的温差值代入f(T<sup>+</sup>),得到正温差对应的概率密度值,利用式(6)对正温差样本的概率密度特性进行拟合:f(T<sup>+</sup>)=F′(T<sup>+</sup>) (5)<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>g</mi><mrow><mo>(</mo><msup><mi>T</mi><mo>+</mo></msup><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>m</mi><mn>1</mn></msub></munderover><msub><mi>v</mi><mi>k</mi></msub><mo>[</mo><msub><mi>&alpha;</mi><mi>k</mi></msub><msubsup><mi>&beta;</mi><mi>k</mi><mrow><mo>-</mo><msub><mi>&alpha;</mi><mi>k</mi></msub></mrow></msubsup><msup><mi>T</mi><mrow><mo>+</mo><msub><mi>&alpha;</mi><mi>k</mi></msub><mo>-</mo><mn>1</mn></mrow></msup><msup><mi>e</mi><msup><mrow><mo>(</mo><mo>-</mo><mfrac><msup><mi>T</mi><mo>+</mo></msup><msub><mi>&beta;</mi><mo>+</mo></msub></mfrac><mo>)</mo></mrow><msub><mi>&alpha;</mi><mi>k</mi></msub></msup></msup><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000508240820000023.GIF" wi="1215" he="179" /></maths>式中,g(T<sup>+</sup>)表示正温差变量的概率密度函数,m<sub>1</sub>为≥2的整数,v<sub>k</sub>表示威布尔分布函数的权重,且<img file="FDA0000508240820000024.GIF" wi="208" he="129" />α<sub>k</sub>表示威布尔分布函数的形状参数,β<sub>k</sub>表示威布尔分布函数的尺寸参数,其中,k为整数,且k=1、2、…、m<sub>1</sub>,v<sub>k</sub>、α<sub>k</sub>和β<sub>k</sub>为待估的正温差样本的概率密度参数;基于最小二乘法,利用正温差值及其概率密度值,对函数g(T<sup>+</sup>)进行拟合,确定正温差样本的概率密度参数v<sub>k</sub>、α<sub>k</sub>和β<sub>k</sub>;步骤(7):确定负温差样本的概率密度参数:利用式(7)对F(T<sup>‑</sup>)求导,得到f(T<sup>‑</sup>),将负温差样本的温差值代入f(T<sup>‑</sup>),得到负温差对应的概率密度值,然后将负温差取相反数,变为对应的正温差,用‑T<sup>‑</sup>表示负温差对应的正温差,最后利用式(8)对‑T‑的概率密度特性进行拟合:f(T<sup>‑</sup>)=F'(T<sup>‑</sup>) (7)<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>g</mi><mrow><mo>(</mo><mo>-</mo><msup><mi>T</mi><mo>-</mo></msup><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>p</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>m</mi><mn>2</mn></msub></munderover><msub><mi>&rho;</mi><mi>p</mi></msub><mo>[</mo><msub><mi>&gamma;</mi><mi>p</mi></msub><msubsup><mi>&eta;</mi><mi>p</mi><mrow><mo>-</mo><msub><mi>&gamma;</mi><mi>p</mi></msub></mrow></msubsup><msup><mrow><mo>(</mo><mo>-</mo><msup><mi>T</mi><mo>-</mo></msup><mo>)</mo></mrow><mrow><msub><mi>&gamma;</mi><mi>p</mi></msub><mo>-</mo><mn>1</mn></mrow></msup><msup><mi>e</mi><msup><mrow><mo>(</mo><mo>-</mo><mfrac><msup><mi>T</mi><mo>-</mo></msup><msub><mi>&eta;</mi><mi>p</mi></msub></mfrac><mo>)</mo></mrow><msub><mi>&gamma;</mi><mi>p</mi></msub></msup></msup><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000508240820000031.GIF" wi="1250" he="194" /></maths>式中,g(‑T<sup>‑</sup>)表示‑T<sup>‑</sup>的概率密度函数,m<sub>2</sub>为≥2的整数,ρ<sub>p</sub>表示威布尔分布函数的权重,且<img file="FDA0000508240820000032.GIF" wi="223" he="135" />γ<sub>p</sub>表示威布尔分布函数的形状参数,η<sub>p</sub>表示威布尔分布函数的尺寸参数,其中,p为整数,且p=1、2、…、m<sub>2</sub>,ρ<sub>p</sub>、γ<sub>p</sub>和η<sub>p</sub>为待估的负温差样本的概率密度参数;基于最小二乘法,利用‑T<sup>‑</sup>值及其概率密度值对函数g(‑T<sup>‑</sup>)进行拟合,确定负温差样本的概率密度参数ρ<sub>p</sub>、γ<sub>p</sub>和η<sub>p</sub>;步骤(8):确定日照温差采集样本的概率密度:根据步骤(6)得到的正温差样本的概率密度参数和步骤(7)得到负温差样本的概率密度参数,利用式(9)得到日照温差采集样本的概率密度y(T):<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mi>y</mi><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mfrac><msub><mi>l</mi><mn>1</mn></msub><mi>l</mi></mfrac><mi>g</mi><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow></mtd><mtd><mi>T</mi><mo>&GreaterEqual;</mo><mn>0</mn></mtd></mtr><mtr><mtd><mfrac><msub><mi>l</mi><mn>2</mn></msub><mi>l</mi></mfrac><mi>g</mi><mrow><mo>(</mo><mo>-</mo><mi>T</mi><mo>)</mo></mrow></mtd><mtd><mi>T</mi><mo>&lt;</mo><mn>0</mn></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000508240820000033.GIF" wi="1162" he="278" /></maths>式中,T为温差变量,当T<0时,T为T<sup>‑</sup>;当T≥0时,T为T<sup>+</sup>;l<sub>1</sub>为正温差样本中温差值的总数,l<sub>2</sub>为负温差样本中温差值的总数,l表示日照温差采集样本中温差值的总量,且l<sub>1</sub>+l<sub>2</sub>=l。
地址 211189 江苏省南京市江宁开发区东南大学路2号