用於估計蓄電池充電率的裝置的製作方法
2023-07-25 15:40:41 4
專利名稱：用於估計蓄電池充電率的裝置的製作方法
技術領域：
本實用新型涉及用於估計蓄電池充電率(縮寫為SOC)的裝置和方法。
背景技術：
2000年11月24日公開的首次公開號為2000-323183、2000年9月29日公開的首次公開號為2000-268886的日本專利申請和日本電子工程會社(T.IEEE Japan)1992年公開的卷號為112-C，No.4、名稱為「Estimation of Open Voltage and Residual Values for pbBattery by Adaptive Digital Filter」的日本文件預先舉例說明過所提出的用於蓄電池的SOC估計裝置。那就是說，由於蓄電池的充電率(或稱作充電狀態，即，SOC)與斷路電路電壓V0(當電池的電源供給切斷時的電池端電壓，也稱作電動勢或斷路電壓)相關，所以當獲得斷路電壓V0時能估計出充電率。然而，在切斷電源供給(結束充電和放電)後直到端電壓穩定需要相當長的時間。因此，從充電和放電結束起需要預定的時間間隔以確定精確的斷路電壓V0。因此，從充電/放電時間期間或充電和放電之後，立即確定出準確的斷路電壓是不可能的，並且使用上述方法也不可能獲得充電率。不過，為了確定出斷路電壓V0，使用上述首次公開號為2000-323183的日本專利申請中公開的方法可以估計出斷路電壓V0。
發明內容
然而，在上述公開號為2000-323183的日本專利申請中所公開的方法中，斷路電壓V0是從非循環(非衰退類型)電池模型(該模型的輸出值僅根據輸入值的當前值和過去值確定)中計算出的，該電池模型的特性完全不同於使用自適應數字濾波器(順序類型模型參數鑑別算法)的電池的物理特性。該充電率SOC從該值中被使用。因此，當該方法被應用到實際的電池特性(輸入電流，輸出電壓)時，根據該電池特性，估計計算完全收斂到實際值，或者根本不收斂。所以，準確的估計出充電率SOC是困難的。
因此，本實用新型的目的是提供用於準確的估計蓄電池的充電率(SOC)和準確的估計出其它與充電率(SOC)有關的參數的裝置和方法。
根據本實用新型的一個方面，提供一種用於蓄電池的充電率估計裝置，包括能夠測量出流經蓄電池的電流的電流檢測部分；能夠測量出蓄電池兩端電壓的端電壓檢測部分；參數估計部分，使用在公式(1)顯示的時間連續序列中的電池模型計算自適應數字濾波並同時估計所有的參數，這些參數對應於是公式(1)中的偏移量的斷路電壓和是瞬時量的係數A(s)，B(s)和C(s)；以及充電率估計部分，從前面得出的斷路電壓V0和使用斷路電壓V0的充電率SOC之間的關係中估計充電率，V=B(s)A(s)I+1C(s)V0---(1)]]>其中s表示拉普拉斯變換算子，A(s)，B(s)和C(s)表示s的多項式函數。
根據本實用新型的另一方面，提供一種用於蓄電池的充電率估計方法，包括測量流過蓄電池的電流；測量蓄電池的端電壓；利用在公式(1)表示的時間連續序列中存在的電池模型計算適應性數字濾波；同時估計所有的參數，這些參數對應於是公式(1)的偏移量的斷路電壓和是瞬時量的係數A(s)，B(s)和C(s)；並且從前面得出的斷路電壓V0和使用斷路電壓V0的充電率SOC之間的關係中估計充電率，V=B(s)A(s)I+1C(s)V0,]]>其中s表示拉普拉斯變換算子，A(s)，B(s)和C(s)表示s的多項式函數。
根據本實用新型的另一目的，提供一種用於蓄電池的充電率估計方法，包括測量流經蓄電池的電流I(k)；測量蓄電池的端電壓V(k)；當電流調零時將端電壓V(k)存儲作為端電壓的初始值ΔV(k)＝V(k)-V_ini；從公式(19)中確定出瞬時電流值I0(k)、I1(k)和I2(k)以及瞬時端電壓V1(k)、V2(k)和V3(k)，I0=1G1(s)I.]]>I1=sG1(s)I.]]>V1=sG1(s)V.]]>I2=s2G1(s)I.]]>V2=s2G1(s)V.]]>I3=s3G1(s)I.]]>V3=s3G1(s)V.]]>1G1(s)=1(p1s+1)3---(19)]]>其中p1表示確定G1(s)響應特性的常數；將瞬時電流值I0(k)、I1(k)和I2(k)以及瞬時端電壓V1(k)、V2(k)和V3(k)代入公式(18)，
(k)=3(k)1+3(k)T(k)P(k-1)(k)]]>θ(k)＝θ(k-1)-γ(k)·P(k-1)·ω(k)·[ωT(k)·θ(k-1)-y(k)]P(k)=11(k){P(k-1)-3(k)P(k-1)(k)T(k)P(k-1)1+3(k)T(k)P(k-1)(k)}=P(k)1(k)]]>1(k)={trace{P(k)}U:1trace{P(k)}U]]>{1:trace{P(k)}U1trace{P(k)}L]]>{trace{P(k)}L:trace{P(k)}L1]]>其中θ(k)表示在k(k＝0，1，2，3...)時間點的參數估計值，λ1，λ3(k)，γu和γL表示初始設定值，b＜λ1＜1，0＜λ3(k)＜∞，P(0)是一個足夠大的值，θ(0)提供一個非零但非常小的初始值，軌跡{P}表示矩陣P的軌跡，其中y(k)＝V1(k)ωT(k)＝[V3(k)V2(k)I3(k)I2(k)I1(k)I0(k)](k)=-a(k)-b(k)c(k)d(k)e(k)f(k)---(20);]]>將參數估計值θ(k)中的a、b、c、d、e和f代入公式(22)用來計算V0的可替代值V0`，該V0相應於某一時刻起斷路電壓估計值的變量ΔV0(k)，在該時刻執行估計的計算開始值；V0=(T1s+1)G2(s)V0=aV6+bV5+V4-cI6-dI5-eI4---(22);]]>
以及根據斷路電壓估計值的變量ΔV0(k)和端電壓初始值V_ini計算斷路電壓估計值V0(k)。
本實用新型的概述不必描述所有必要特徵以便本實用新型也可以是所述特徵的子集。
附圖簡述
圖1是根據本實用新型，在優選實施例中用於估計蓄電池充電率(SOC)的裝置功能方框圖。
圖2是根據本實用新型，在優選實施例中用於估計蓄電池充電率的裝置的具體電路方框圖。
圖3是表示蓄電池的等效電路模型的模型視圖。
圖4是表示斷路電壓和充電率(SOC)之間的關係的關係圖。
圖5是解釋在圖1所示第一優選實施例中的充電率估計裝置的電池控制器的微型計算機的操作流程圖。
圖6A，6B，6C，6D，6E，6F，6G，6H和6I是表示圖1所示的實施例中的充電率估計裝置中電流，電壓和各種參數的仿真結果的特徵圖。
發明詳述為了更好地理解本實用新型，以下將參考標記使用在附圖中。
圖1示出了根據本實用新型第一優選實施例中的充電率估計裝置的功能方框圖。在該圖1中，參考標記1表示基於將斷路電壓V0(k)作為偏移量的電池模型的參數估計部分。另外，參考標記2表示用於計算斷路電壓V0(k)的斷路電壓計算部分，以及參考標記3表示用於從斷路電壓計算充電率的充電率估計部分。另外，參考標記4表示用於檢測給電池充電和從電池放電的電流I(k)的電流I測量模塊，以及參考標記5表示電池的s端電壓以測量端電壓V(k)。
圖2示出了表示第一實施例中的充電率估計裝置的具體結構的方框圖。在該實施例中，負載如馬達由蓄電池驅動，而充電率估計裝置安裝在用馬達(負載)再生功率給蓄電池充電的系統中。在圖2中，參考標記10表示蓄電池(簡稱電池)，參考標記20表示負載如DC馬達，參考標記30表示電池控制器(電子控制單元)，用於估計出該電池的充電率(充電狀態)，該電池裝有包括ROM(只讀存儲器)，RAM(隨機存儲器)，CPU(中央處理器)，和輸入/輸出接口以及其它電子電路的微型計算機。參考標記40表示電流表，用於檢測給電池充電或從電池放電的電流，參考標記50表示電壓表，用於檢測電池的端電壓，參考標記60表示溫度計，用於檢測電池的溫度。這些儀表都被連接到電池控制器30。電池控制器30對應於參數估計部分1，斷路電壓V0(k)和充電率估計部分3的一部分。電流表40相應於電流I(k)測量部分，以及電壓表50相應於端電壓V(k)測量部分5。
首先，在第一實施例中使用的「電池模型」將在下面進行描寫。圖3是表示蓄電池的等效電路模型的等效電路。蓄電池的等效電路模型能用下面的公式(7)(＝公式(6))表示。
V=K(T2s+1)T1s+1I+1T3s+1V0---(7)]]>在公式(7)中，模型輸入是電流I(A)(正值表示充電而負值表示放電)，模型輸出是端電壓V[V]，斷路電壓是V0，K表示內部阻抗，T1-T3表示時間常數(T1≠T2≠T3，T1＜＜T3)而s表示拉普拉斯變換算子。
基於公式(7)的模型是衰減模型(一階)，其中正極和負極沒有專門分開。然而，相對容易的表示出實際電池的充電-放電特性是可能的。公式(7)，在公式(1)V＝B(s)/A(s)·I+1/C(s)·V0中，A(s)＝T1·s+1，B(s)＝k·(T2)·s+1，C(s)＝T3·s+1。
以下，背離基於公式(7)的電池模型的自適應數字濾波器將首先在下面進行描述。假如電流I的值乘以變量係數A是根據特定初始狀態積分而來，則斷路電壓V0能用公式(8)描述。
那就是說，V0=AsI---(8)]]>注意到公式(8)中的係數A能用相應於公式(2)中敘述的h替代，即，V0＝h/s·I。
假如公式(8)代入公式(7)，能夠得出公式(9)。
V0=K(T2s+1)T1s+1I+1T3s+1AsI---(9)]]>公式(9)相應於公式(3)，(V=(B(s)A(s)+1C(s)hs)I=sB(s)C(s)+hA(s)sA(s)C(s)I---(3))]]>對於公式(3)中的A(s)，B(s)，C(s)，下面的公式以與公式(7)同樣的方式代入公式(9)。
A(s)＝T1·s+1，B(s)＝k·(T2·s)+1C(s)＝T3·s+1。換句話說，公式(3)是導出的公式並且這應用到一階模型就是公式(9)。假如公式(9)成立，公式(10)就是給定的。
S·(T1·s+1)(T3·s+1)·V＝K·(T2·s+1)(T3·s+1)s·I+A·(T1·s+1)·I{T1·T3·s3+(T1+T3)·s2+s}·V＝{K·T2·T3·s3+K·(T2+T3)·s2+(K+A·T1)·s+A}·I(a·s3+b·s2+s)·V＝(c·s3+d·s2+e·s+f)·I ---...(10)
可以注意到，在公式(10)的最後一個等式中，參數被重寫如下a＝T1·T3，b＝T1+T3，c＝K·T2·T3，d＝K·(T1+T3)，e＝K+A·T1，and f＝A ---(11)。
假如把穩態低通濾波器G1(s)引入到公式(10)的兩邊並且整理，得到下面的公式(12)。
1G1(s)(as3+bs2+s)V=1G1(s)(cs3+ds2+es+f)I---(12)]]>詳細的講，公式(10)與公式(7)相反，假如將T1·S+1＝A(s)，K·(T2·S+1)＝B(s)，T3·s+1＝C(s)代入公式(10)，將會得出s·A(s)·C(s)·V＝B(s)·C(s)·s·I+A·A(s)·I。重新排列如下s·A(s)·C(s)·V＝[B(s)·C(s)·s·I+A·A(s)]·I---(12`)。假如把低通濾波器(LPF)G1(s)引入到公式(12`)的兩邊，將得出公式(4)。那就是說，sA(s)C(s)G1(s)V=sB(s)C(s)+hA(s)sA(s)C(s)I---(4)]]>應當注意到，s表示拉普拉斯變換算子，A(s)，B(s)，和C(s)表示s的多項式函數，h表示變量，以及1/G1(s)表示具有低通濾波特性的傳輸函數。那就是說，公式(4)是導出的函數，公式(12)是公式(4)在一階模型中的應用。
假如p1表示用來確定G1(s)的響應特性的常數並且能根據設計者的意願來確定，那麼實際測量的電流I和端電壓V可通過低通濾波器(LPF)和帶通濾波器(BPF)來處理，這些濾波器通過下面的公式(13)可以定義。
I0=1G1(s)I]]>I1=sG1(s)I,]]>V1=sG1(s)V,]]>I2=s2G1(s)I,]]>V2=s2G1(s)V,]]>I3=s3G1(s)I,]]>V3=s3G1(s)V,]]>1G1(s)=1(P1s+1)3---(13)]]>假如公式(12)用公式(13)所示的變量重寫，並且進行變形公式(14)能被表示為下面的公式(15)。
a·V3+b·V2+V1＝c·I3+d·I2+e·I1+f·I0V1＝-a·V3-b·V2+c·I3+d·I2+e·I1+f·I0---(14)。
V1=V3V2I3I2I1I0=-a-bcdef---(15).]]>公式(15)是測量值和未知參數的積-和等式。因此，自適應濾波器的標準(一般)類型(公式(16))與公式(15)是一致的。可以注意到ωT表示轉置向量，其中向量ω的行和列可以相互交換。
y＝ωT·θ---(16)。可注意到上述公式(16)中的y，ωT，和θ可用下面的公式(17)表示。
Y=V1,T=V3V2I3I2I1I0,=-a-bcdef---(17)]]>因此，假如對電流I和端電壓V處理的信號濾波器使用在數字濾波器處理計算中，能估計出未知參數向量θ。
在該實例中，使用「雙極限軌道增益法」改善了簡單「自適應數字濾波器通過最小二乘法」產生的邏輯缺陷。這樣的話，一旦估計值被收斂，即便改變參數，也不能再進行精確估計。用公式(16)估計未知參數向量θ的參數估計運算法則作為先決條件在公式(18)中示出。可注意到在k時間點的參數估計值是θ(k)。
(k)=3(k)1+3(k)T(k)P(k-1)(k)]]>θ(k)＝θ(k-1)-γ(k)·P(k-1)·ω(k)·[ωT(k)·θ(k-1)-y(k)]P(k)=11(k){P(k-1)-3(k)P(k-1)(k)T(k)P(k-1)1+3(k)T(k)P(k-1)(k)}=P(k)1(k)]]>1(k)={trace{P(k)}U:1trace{P(k)}U]]>{1:trace{P(k)}U1trace{P(k)}L]]>{trace{P(k)}L:trace{P(k)}L1---(18).]]>在公式(18)中，λ1，λ3(k)，γu和γL表示初始設定值，b＜λ1＜1，0＜λ3(k)＜∞。P(0)是一個足夠大的值，θ(0)提供了非零但非常小的初始值。另外，軌跡{P}表示矩陣P的軌跡。如上所述，自適應數字濾波器來自電池模型。
圖5顯示了執行電池控制器30的微型計算機的操作流程圖。圖5所示的例行程序是對於每個時間T0的恆定周期進行的。例如，I(k)是當前值，I(k-1)表示I(k)的前一個值。在步驟S10，電池控制器30測量電流I(k)並且以I(k-1)表示I(k)的前一個值。在步驟S20，電池控制器30執行蓄電池的斷開中繼的通-斷決定。那就是說，電池控制器30執行蓄電池的斷開中繼的轉換通斷控制。當中繼關閉時(電流I＝0)，例行程序轉到步驟S30。在中繼接合期間，例行程序轉到步驟S40。在步驟S30，當中繼被佔用時，例行程序轉到步驟S540。在步驟S530，電池控制器30存儲端電壓V(k)作為端電壓的初始值V_ini。在步驟S40，電池控制器30計算端電壓的差分值ΔV(k)。ΔV(k)＝V(k)-V_ini。這是因為在自適應數字濾波器中的估計參數初始值是0，使得在估計計算啟動期間估計參數不收斂。從而，所有的輸入都為零。在輸入都為零期間，在中繼斷開期間，步驟S30已經通過，並且由於I為零所以估計參數都保持初始狀態同時估計參數還保持活躍。
在步驟S50，在公式(13)的基礎上，低通濾波或帶通濾波實現電流I(k)和端電壓差值ΔV(k)。從公式(19)中計算出I0(k)至I3(k)和V1(k)至V3(k)。在這種情況下，為了提高公式(18)的參數估計算法的估計精度，低通濾波器的響應特性G1(s)設置得很緩慢以便減少觀察噪音。然而，假如該特性比蓄電池的響應特性快(時間常數T1的大概值是已知的)，就不能準確的估計出電池模型的每個參數。應當注意到在公式(19)中引用的P1表示由G1(s)的響應特性確定的常數。
I0=1G1(s)I,]]>I1=sG1(s)I,]]>V1=sG1(s)V.]]>I2=s2G1(s)I,]]>V2=s2G1(s)V,]]>I3=s3G1(s)I,]]>V3=s3G1(s)V,and]]>1G1(s)=1(p1s+1)3---(19).]]>在步驟S60，在步驟S(50)計算出的I0(k)至I3(k)和V(k)至V3(k)代入到公式(18)。那麼，在自適應數字濾波器中的參數估計算法，即，公式(18)被執行以計算參數估計值θ(k)。在公式(20)中示出y(k)，ωT(k)，和θ(k)。
y(k)＝V1(k)ωT(k)＝[V3(k)V2(k)I3(k)I2(k)I1(k)I0(k)}(k)=-a(k)-b(k)c(k)d(k)e(k)f(k)---(20).]]>在步驟S70，將步驟S60計算出的參數估計值θ(k)的a到e代入下面的公式(22)中，其中對上述的電池模型公式(7)變形來計算V0′，該V0′是斷路電壓V0的可選擇值。由於斷路電壓V0變量是平穩的，V0′可被替代的使用。應當注意，從估計計算開始時間中導出斷路電壓的變量ΔV0(k)。
應當注意到公式(21)中的等式[1/C1(s)]I用相應公式(22)的公式(24)代替。也應當注意到，在公式(22)的推導中，公式(21)中的k嚴格的不同於公式(21)中的e。然而，原則上，k＞＞A.T1，e約等於k(e≌k)。那麼，公式(22)中的每個係數a到e是在公式(23)中顯示的內容。
1T3s+1V0=V-K(T3+s+1)T1s+1I]]>(T1·s+1)·V0＝(T1·s+1)(T3·s+1)V-K·(T2·s+1)(T3·s+1)·I(T1·s+1)·V0＝{T1·T3·s2+(T1+T3)·s+1}·V-{K·T2·T3·s2+K·(T2+T3)·s+K}·I(T1s+1)G2V0=1G2(s)(as2+bs+K)I---(21).]]>V0=(T1s+1)G2(s)V0=aV6+bV5+V4-cI6-dI5-eI4----(22).]]>應當注意a＝T1·T3，b＝T1+T3，c＝K·(T2+T3}，d＝K·{T2+T3}，e＝K+λ·T1＝K---(23)。
I4=1G2(s)I,]]>V4=1G2(s)V,]]>I5=sG2(s)I,]]>V5=sG2(s)V,]]>1G2(s)=1p2s+11T1s+1,]]>I6=s2G2(s)I,]]>andV6=s2G2(s)V---(24)]]>公式(24)中的p2表示用來確定G2(s)的響應特性的常數。電池參數的T1為幾秒是公知的。因此，公式(24)的T′1被設置的約等於T1。因此，由於保持在公式(22)中的分子(T1·s+1)能被補償，斷路電壓V0的估計準確性能被提高。應當注意到公式(21)相應於公式(5)。那就是說，公式(21)能從(T1·S+1)＝(T1·S+1)(T3·S+1)·V-K·(T2·S+1)(T3·S+1)·(T3·S+1)·I得出。假如下面的三個公式被代入上面描寫的公式(21)的變形中，T1·S+1＝A(s)，K·(T2·S+1)＝B(s)，和T3·S+1＝C(s)。那就是說A(s)·V0＝A(s)·C(s)·V-B(s)·C(s)·I。假如重新安排，該結果在V0＝C(s)·V-B(s)·C(s)·I/A(s)，V0＝C(s)·[V-B(s)·I/A(s)]中。假如低通濾波器G2(s)被引入該公式的兩邊，結果會在公式(5)中。詳細的說，公式(5)是歸納的公式，並且公式(5)應用該一階模型是公式(2)。
在步驟S80，電池控制器30將斷路電壓初始值，即端電壓初始值V_ini添加到斷路電壓V0的變量ΔV0(k)，以從下面的公式(25)得到斷路電壓估計值V0(k)。
V0(k)＝ΔV0(k)+V_ini---(25)。
在步驟90，電池控制器30使用圖4所示的斷路電壓比充電率的關係圖，根據在步驟S80中計算出的斷路電壓V0(k)計算出充電率SOC(k)。可注意到，在圖4中，VL表示相應於SOC＝0％的斷路電壓和VH表示相應於SOC＝100％的斷路電壓。在步驟S100，電池控制器30存儲在後續計算中需要的必要數值，並且該程序結束。所上所述，已經描寫了用於估計蓄電池充電率的裝置的操作。
(1)如上描述，在蓄電池電流I和端電壓V中的關係，和斷路電壓V0在作為一般公式(1)的轉移函數中被構造，在優選實施例中，公式(7)等於公式(6)。因此，應用自適應數字濾波器，如最小平方(眾所周知的估計算法)是可能的。所以，估計以批處理形式存在的公式(即，是偏移量的斷路電壓V0和多詞學名公式A(s)，B(s)和C(s))中的參數是有可能的。這些參數被充電率，周圍溫度、惡化和瞬時變化大大的影響。順序的估計出具有好的準確度的自適應數字濾波器是可能的。那麼，假如圖4所示的斷路電壓V0和充電率被存儲，所估計的斷路電壓能轉換成充電率。因此，以相同的方式連續的估計出充電率作為上述參數是可能的。
(2)在公式(1)約等於公式(4)的情況下，其中公式(1)反映的是蓄電池的電流I和端電壓V的關係等式，在該公式中沒有包括偏移量(即，斷路電壓V0)，可以得到被濾波器處理的測量電流I和被濾波器處理的端電壓V之間的乘積-和-相加公式和未知參數(多詞學名公式A(s)，B(s)，和C(s)的係數參數和h)。正常應用的自適應數字濾波器(最小平均平方和眾所周知的參數估計算法)能直接應用在連續的時間序列中。
作為該結果，未知參數能以批處理方式被估計並且估計的參數h被代入公式(2)，斷路電壓V0的估計值能很容易的被計算出。所有的這些參數是瞬間變化的，自適應數字濾波器可以高準確地估計出任何時刻的充電率。由於在斷路電壓V0和充電率SOC之間建立的常數關係顯示在圖4中，假如該關係被提前存儲，充電率SOC能從斷路電壓V0的估計值中被估計出。
圖6A到6I整體上示出了信號時刻表，當每個參數被估計時，電流I和端電壓V輸入到自適應濾波器並且表示出仿真圖的結果。就公式(6)中的一階延遲的時間常數而言，T1＜T0。由於所有的參數a到f(指公式(11))被順利的估計，可以認為斷路電壓V0的估計值與真實值相當一致。
應當注意到，在預示出斷路電壓的圖6C中，公式(6)的右邊的秒項被描寫的原因是為了預示出斷路電壓估計值與真實值是幾乎沒有延遲的一致，儘管時間常數T3的最新量在輸入到自適應濾波器的端電壓上被測量。詳細的說，由於在公式(6)中，電池模型的參數估計格式化了自適應數字濾波器，所有的參數a到f能順利的被估計，並且斷路電壓V0的估計值能與真實值相當一致。
(3)另外，如(2)項所述，在該結構中，其中斷路電壓V0從公式(2)中被計算出，在估計值h收斂到真實值之前會發生積分，它的錯誤不能消除。然而，在公式(5)的結構中，其中沒包括積分，在參數估計值收斂到真實值之前的差錯在收斂之後不會受影響。
可以理解，在圖6I的部分①中，在估計值f收斂到真實值之前，僅即刻執行錯誤估計。在公式(2)中，該值也被積分使得錯誤沒被消除。然而，即便該值被積分錯誤也不會消除。但是在公式(5)的結構中，斷路電壓V0從不包括積分的公式中計算出。因此，在參數估計值收斂到實際時刻，該錯誤的估計部分能被消除。
(4)而且，在公式(6)代替公式(1)使用的情況下，計算時間和程序容量能被壓縮到最小而具有上述的優點。
日本專利申請號2002-340803(在日本於2002年10月25日申請)的全部內容在此被引作參考。本實用新型範圍參照下面的權利要求定義。
權利要求1.一種用於蓄電池的充電率估計裝置，包括測量流過蓄電池的電流(I)的電流檢測部分；測量蓄電池的端電壓(V)的端電壓檢測部分，其特徵在於所述充電率估計裝置還包括利用公式(1)中所示的連續時間序列的電池模型計算自適應數字濾波並同時估計所有的參數的參數估計部分，所述參數對應於作為公式(1)的偏移項的斷路電壓(V0)和作為瞬時項的係數A(s)、B(s)和C(s)；以及利用斷路電壓(V0)根據先前得出的斷路電壓(V0)和充電率(SOC)之間的關係來估計充電率的充電率估計部分，V=B(s)A(s)I+1C(s)V0---(1)]]>其中s表示拉普拉斯變換算子，A(s)、B(s)和C(s)表示s的多項式函數。
2.如權利要求1所述的蓄電池的充電率估計裝置，其中在公式(1)所示的連續時間序列中的電池模型的斷路電壓(V0)通過公式(2)被估計以提供公式(3)，並且數字濾波計算使用公式(3)和等同的公式(4)來執行，至少在公式(4)中估計h，該估計的h代入公式(2)得出斷路電壓(V0)，並且從先前得出的斷路電壓(V0)和充電率(SOC)之間的關係估計充電率，v0=hsI---(2)]]>V=(B(s)A(s)+1C(s)hs)I=sB(s)C(s)+hA(s)sA(s)C(s)I---(3)]]>sA(s)C(s)G1(s)V=sB(s)C(s)+hA(s)G1(s)I---(4).]]>其中s表示拉普拉斯變換算子，A(s)、B(s)和C(s)表示s的多項式函數，h表示變量，以及1/G1(s)表示具有低通濾波器特性的傳輸函數。
3.如權利要求1所述的蓄電池的充電率估計裝置，其中在公式(2)中估計在連續時間序列中的電池模型的斷路電壓(V0)以計算公式(3)，自適應濾波計算通過等同於公式(3)的公式(4)來執行，A(s)、B(s)和C(s)可從公式(4)中估計出，該估計出的A(s)、B(s)和C(s)代入公式(5)能夠確定出V0/G2(s)，並且根據先前得出的斷路電壓(V0)和充電率(SOC)之間的關係估計充電率，該充電率(SOC)使用獲得的V0/G2(s)代替斷路電壓V0，V0=hsI---(2)]]>V=(B(s)A(s)+1C(s)hs)I=sB(s)C(s)+hA(s)sA(s)C(s)I---(3)]]>sA(s)C(s)G1(s)V=sB(s)C(s)+hA(s)G1(s)I---(4),]]>V0G2(s)=C(s)G2(s)(V-B(s)A(s)I)---(5).]]>其中s表示拉普拉斯變換算子，A(s)、B(s)和C(s)表示s的多項式(公式)的函數，h表示變量，1/G1(s)和1/G2(s)表示具有低通濾波器特性的傳輸函數。
4.如權利要求1所述的蓄電池的充電率估計裝置，其中根據公式(6)中計算該電池模型，V=K(T2s+1)T1s+1I+1T3s+1V0,]]>其中k表示蓄電池的內部阻抗，T1、T2和T3表示時間常數，1/G1(s)表示具有三階或更多階的低通濾波器，以及1/G2(s)表示另一個具有二階或更多階的低通濾波器。
5.如權利要求4所述的蓄電池的充電率估計裝置，其中A(s)＝T1.s+1，B(s)＝k.(T2).s+1，C(s)＝T3.s+1。
6.如權利要求5所述的蓄電池的充電率估計裝置，其中V=K(T2s+1)T1s+1I+1T3s+1AsI---(9)]]>(a·s3+b·s2+s)·V＝(c·s3+d·s2+e·s+f)·I---(10).
7.如權利要求6所述的蓄電池的充電率估計裝置，其中a＝T1·T3.b＝T1+T3，c＝K·T2·T3，d＝K·(T2+T3)，e＝K+A·T1，f＝A---(11).
8.如權利要求7所述的蓄電池的充電率估計裝置，其中穩定的低通濾波器G1(s)引入公式(10)的兩邊，得出下面的公式1G1(s)(as3+bs2+s)V=1G1(s)(cs3+ds2+cs+f)I---(12).]]>
9.如權利要求8所述的蓄電池的充電率估計裝置，其中通過低通濾波器處理的實際測量的電流I和端電壓V表示如下I0=1G1(s)I.]]>I1=sG1(s)I,V1=sG1(s)V.]]>I2=s2G1(s)I,V2=s2G2(s)V.]]>I3=s3G1(s)I,V3=s2G1(s)V,and]]>1G1(s)=1(P1s+1)3---(13).]]>
10.如權利要求9所述的蓄電池的充電率估計裝置，其中利用公式(13)，公式(12)被重寫和重新安排如下V1=V3V2I3I2I1I0=-a-bcdef---(15)]]>並且公式(15)對應於一般公式，該一般公式與公式(16)y＝ωT·θ的普通自適應數字濾波器的標準形式一致，其中y＝V1、ωT＝[V3V2I3I2I1I0]並且=-a-bcdef---(17).]]>
11.如權利要求10所述的蓄電池的充電率估計裝置，其中將公式(16)作為先決條件的參數估計算法定義如下(k)=3(k)I+3(k)T(k)P(k-1)(k)]]>θ(k)＝θ(k-1)-γ(k)·P(k-1)·ω(k)·[ωT(k)·θ(k-1-y(k)]P(k)=11(k){P(k-1)-3(k)P(k-1)(k)T(k)P(k-1)1+3(k)T(k)P(k-1)(k)}P(k)3(k)]]>1(k)={trace{Pk}u:1trace{P(k)}u]]>{2:trace{Pk}u1trace{P(k)}L]]>{trace{Pk}L:trace{P(k)}L1(18),]]>其中θ(k)表示在時間點k(k＝0，1，2，3...)的參數估計值，λ1，λ3(k)，γu和γL表示初始設定值，b＜λ1＜1，0＜λ3(k)＜∞，P(0)是一個足夠大的值，θ(0)提供一個非零但非常小的初始值，以及軌跡{P}是矩陣P的軌跡。
專利摘要在蓄電池的充電率估計裝置和方法中，測量流過蓄電池的電流，測量蓄電池的端電壓，使用以公式(1)所示的連續時間序列中的電池模型來執行自適應數字濾波，同時估計所有的參數，相應於斷路電壓的參數是公式(1)中的偏移項和是瞬時項的係數A(s)，B(s)和C(s)，並且，從提前得出的斷路電壓V
文檔編號H02J7/00GK2724218SQ200320116788
公開日2005年9月7日 申請日期2003年11月25日 優先權日2002年11月25日
發明者湯本大次郎, 中村英夫 申請人:日產自動車株式會社