JP6171897B2 - Approximation function creation program, approximation function creation method, approximation function creation device, and charging rate estimation program - Google Patents

Approximation function creation program, approximation function creation method, approximation function creation device, and charging rate estimation program Download PDF

Info

Publication number
JP6171897B2
JP6171897B2 JP2013247935A JP2013247935A JP6171897B2 JP 6171897 B2 JP6171897 B2 JP 6171897B2 JP 2013247935 A JP2013247935 A JP 2013247935A JP 2013247935 A JP2013247935 A JP 2013247935A JP 6171897 B2 JP6171897 B2 JP 6171897B2
Authority
JP
Japan
Prior art keywords
characteristic curve
approximate
range
charging
straight line
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.)
Active
Application number
JP2013247935A
Other languages
Japanese (ja)
Other versions
JP2015105875A (en
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.)
Fujitsu Ltd
Original Assignee
Fujitsu Ltd
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 Fujitsu Ltd filed Critical Fujitsu Ltd
Priority to JP2013247935A priority Critical patent/JP6171897B2/en
Publication of JP2015105875A publication Critical patent/JP2015105875A/en
Application granted granted Critical
Publication of JP6171897B2 publication Critical patent/JP6171897B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
    • Y02E60/10Energy storage using batteries

Description

本発明は、近似関数作成プログラム、近似関数作成方法、近似関数作成装置および充電率推定プログラムに関する。   The present invention relates to an approximation function creation program, an approximation function creation method, an approximation function creation device, and a charging rate estimation program.

従来から、充電を行うことにより、繰り返し使用可能な二次電池が知られている。この二次電池としては、例えば、リチウムイオン二次電池が挙げられる。二次電池は、蓄えられた残りの電力がどの程度かを把握するため、充電率(SoC:State of Charge)を精度よく求めることが、運用する上で重要である。そこで、例えば、二次電池の端子電圧からカルマンフィルタを用いて二次電池の充電率を推定する技術が知られている。   Conventionally, a secondary battery that can be repeatedly used by charging is known. Examples of the secondary battery include a lithium ion secondary battery. In order to know how much of the remaining electric power is stored in the secondary battery, it is important for operation to accurately calculate the state of charge (SoC). Thus, for example, a technique for estimating the charging rate of the secondary battery from the terminal voltage of the secondary battery using a Kalman filter is known.

特開2009−234557号公報JP 2009-234557 A 特開2012−139025号公報JP 2012-139025 A

ところで、二次電池は、負荷をかけていない状態の端子の電圧、所謂、開回路電圧(OCV:Open circuit voltage)と充電率との間に相関関係がある。このOCVと充電率との相関関係は、特性曲線により表すことができる。この特性曲線を正確に関数化することにより、二次電池の充電率の推定精度が向上する。   By the way, in the secondary battery, there is a correlation between the voltage of the terminal in a state where no load is applied, so-called open circuit voltage (OCV) and the charging rate. The correlation between the OCV and the charging rate can be expressed by a characteristic curve. By accurately functionalizing this characteristic curve, the estimation accuracy of the charging rate of the secondary battery is improved.

一方、カルマンフィルタは、非線形の関数を扱う場合、計算量が大きく増加する。このため、カルマンフィルタを用いて二次電池の充電率を推定する場合、特性曲線として非線形の関数を用いると、計算の負荷が大きくなる。このため、例えば、技術者が試行錯誤によって特性曲線を折れ線により近似した近似関数を求めている。   On the other hand, when the Kalman filter handles a nonlinear function, the amount of calculation greatly increases. For this reason, when estimating the charging rate of a secondary battery using a Kalman filter, if a non-linear function is used as a characteristic curve, the calculation load increases. For this reason, for example, an engineer obtains an approximate function in which a characteristic curve is approximated by a broken line by trial and error.

しかし、技術者の試行錯誤では、特性曲線を精度よく近似できない場合がある。また、二次電池は、種類によって特性曲線が変化する。このため、技術者がそれぞれの特性曲線の近似関数を試行錯誤で作成することは困難である。   However, the characteristic curve may not be approximated accurately by trial and error by an engineer. Further, the characteristic curve of the secondary battery changes depending on the type. For this reason, it is difficult for an engineer to create an approximate function of each characteristic curve by trial and error.

1つの側面では、本発明は、特性曲線の近似を精度よく行うことができる近似関数作成プログラム、近似関数作成方法、近似関数作成装置および充電率推定プログラムを提供することを目的とする。   In one aspect, an object of the present invention is to provide an approximation function creation program, an approximation function creation method, an approximation function creation device, and a charging rate estimation program capable of accurately approximating a characteristic curve.

一態様の近似関数作成プログラムは、コンピュータに、二次電池の開回路電圧と当該二次電池の充電率との相関関係を示した特性曲線のデータを取得する処理を実行させる。また、近似関数作成プログラムは、コンピュータに、取得されたデータにより示される特性曲線の充電率に対する開回路電圧の変化率が所定の閾値以下となる範囲内で、中間点を移動させつつ中間点の位置毎に、類似度を算出する処理を実行させる。近似関数作成プログラムは、コンピュータに、中間点を境として前記範囲を分割した各分割範囲の特性曲線のデータをそれぞれ近似した近似直線と、当該近似直線の近似に用いた特性曲線のデータとの類似度を算出させる。また、近似関数作成プログラムは、コンピュータに、算出される類似度が最も高い中間点を境とした各分割範囲の近似直線を用いて特性曲線を近似した近似関数を生成する処理を実行させる。   The approximate function creation program according to one aspect causes a computer to execute a process of acquiring characteristic curve data indicating a correlation between an open circuit voltage of a secondary battery and a charging rate of the secondary battery. Further, the approximate function creation program causes the computer to move the intermediate point while moving the intermediate point within a range where the change rate of the open circuit voltage with respect to the charging rate of the characteristic curve indicated by the acquired data is equal to or less than a predetermined threshold. For each position, a process for calculating the similarity is executed. The approximate function creation program uses a computer to approximate the approximate line obtained by approximating the characteristic curve data of each divided range obtained by dividing the range with the intermediate point as the boundary, and the characteristic curve data used for approximating the approximate line. Let the degree be calculated. Further, the approximate function creation program causes the computer to execute a process of generating an approximate function that approximates the characteristic curve using the approximate straight line of each divided range with the calculated intermediate point as the boundary.

特性曲線の近似を精度よく行うことができる。   The characteristic curve can be approximated with high accuracy.

図1は、リチウムイオン電池が使用される流れの一例を模式的に示した図である。FIG. 1 is a diagram schematically illustrating an example of a flow in which a lithium ion battery is used. 図2は、近似関数作成装置の全体構成を示す図である。FIG. 2 is a diagram showing the overall configuration of the approximate function creation device. 図3は、特性曲線の一例を示す図である。FIG. 3 is a diagram illustrating an example of a characteristic curve. 図4は、分割範囲毎の近似直線の一例を模式的に示した図である。FIG. 4 is a diagram schematically illustrating an example of an approximate line for each division range. 図5は、低充電率部分の近似直線の一例を模式的に示した図である。FIG. 5 is a diagram schematically illustrating an example of an approximate line of the low charge rate portion. 図6は、高充電率部分の近似直線の一例を模式的に示した図である。FIG. 6 is a diagram schematically illustrating an example of an approximate line of the high charge rate portion. 図7は、切替点の一例を模式的に示した図である。FIG. 7 is a diagram schematically illustrating an example of the switching point. 図8は、生成された近似関数の一例を示した図である。FIG. 8 is a diagram illustrating an example of the generated approximate function. 図9は、放電用の近似関数と充電用の近似関数の一例を示した図である。FIG. 9 is a diagram illustrating an example of an approximate function for discharging and an approximate function for charging. 図10は、近似関数作成処理の手順を示すフローチャートである。FIG. 10 is a flowchart showing the procedure of the approximate function creation process. 図11は、近似関数作成プログラムを実行するコンピュータの一例を示す図である。FIG. 11 is a diagram illustrating an example of a computer that executes an approximate function creation program.

以下に添付図面を参照して本願に係る近似関数作成プログラム、近似関数作成方法、近似関数作成装置および充電率推定プログラムについて説明する。なお、この実施例は開示の技術を限定するものではない。そして、各実施例は、処理内容を矛盾させない範囲で適宜組み合わせることが可能である。   Hereinafter, an approximate function creating program, an approximate function creating method, an approximate function creating device, and a charge rate estimating program according to the present application will be described with reference to the accompanying drawings. Note that this embodiment does not limit the disclosed technology. Each embodiment can be appropriately combined within a range in which processing contents are not contradictory.

[二次電池]
最初に、二次電池について簡単に説明する。この二次電池としては、例えば、リチウムイオン二次電池が挙げられる。以下では、二次電池としてリチウムイオン二次電池を例に説明する。なお、以下では、リチウムイオン二次電池をリチウムイオン電池とも表記する。
[Secondary battery]
First, the secondary battery will be briefly described. Examples of the secondary battery include a lithium ion secondary battery. Hereinafter, a lithium ion secondary battery will be described as an example of the secondary battery. Hereinafter, the lithium ion secondary battery is also referred to as a lithium ion battery.

リチウムイオン電池は、例えば、パソコンや、携帯電話、電気自動車、電動バイクなど、様々な機器で利用されている。例えば、電気自動車や電動バイク等では、リチウムイオン電池などの二次電池をモジュールごと交換することで、充電時間による使用時間のロスを無くす運用が知られている。このように、リチウムイオン電池は、各種の機器で用いられるため、充電率を精度よく求めることが運用する上で重要である。   Lithium ion batteries are used in various devices such as personal computers, mobile phones, electric cars, and electric motorcycles. For example, in an electric vehicle, an electric motorcycle, or the like, it is known to replace a secondary battery such as a lithium ion battery together with a module so as to eliminate a loss of usage time due to charging time. As described above, since lithium ion batteries are used in various devices, it is important to obtain a charging rate with high accuracy.

図1は、リチウムイオン電池が使用される流れの一例を模式的に示した図である。図1の例では、リチウムイオン電池10の概略的な構成が示されている。リチウムイオン電池10は、セル11と、制御チップ12とを有する。セル11は、化学的な変化により電力の蓄積および放電が可能な部品である。制御チップ12は、リチウムイオン電池10を制御する制御IC(Integrated Circuit)である。制御チップ12は、不揮発性メモリなどの記憶部12Aを有する。記憶部12Aには、リチウムイオン電池10の制御に用いる各種のプログラムや各種の情報が格納される。例えば、記憶部12Aには、充電率を推定する充電率推定プログラムが格納される。また、例えば、記憶部12Aには、充電率の推定に用いる各種の情報が格納される。一例として、記憶部12Aには、特性曲線を近似した近似関数の情報が格納される。この近似関数の情報は、例えば、セル11の電気的な特性を実験等によって測定した測定データに基づいて、近似関数作成装置20により作成される。リチウムイオン電池10は、型番などのタイプごとに特性が異なる。また、リチウムイオン電池10は、製造されたユニットや個体ごとに特性が異なる場合がある。近似関数の情報は、リチウムイオン電池10のタイプごとに測定データを求めて作成してもよく、ユニットや個体ごとに測定データを求めて作成してもよい。リチウムイオン電池10は、各種の機器21に搭載され、使用される。制御チップ12は、各種のプログラムが動作することにより、各種の制御を行う処理部として機能する。例えば、制御チップ12は、充電率推定プログラムが動作することにより、処理部として推定部12Bを有する。制御チップ12は、リチウムイオン電池10の端子電圧と記憶部12Aに記憶された近似関数の情報などの各種の情報に基づき、充電率を推定する。例えば、制御チップ12は、カルマンフィルタを用いて充電率を推定する。制御チップ12は、各種の制御を行う。例えば、制御チップ12は、セル11の充電率を推定した結果、セル11の過充電や過放電が発生する場合、電力の充放電に用いる端子からの電流を遮断する制御を行う。   FIG. 1 is a diagram schematically illustrating an example of a flow in which a lithium ion battery is used. In the example of FIG. 1, a schematic configuration of the lithium ion battery 10 is shown. The lithium ion battery 10 includes a cell 11 and a control chip 12. The cell 11 is a component capable of storing and discharging electric power by a chemical change. The control chip 12 is a control IC (Integrated Circuit) that controls the lithium ion battery 10. The control chip 12 includes a storage unit 12A such as a nonvolatile memory. The storage unit 12A stores various programs and various information used for controlling the lithium ion battery 10. For example, the storage unit 12A stores a charging rate estimation program for estimating the charging rate. Further, for example, various kinds of information used for estimating the charging rate is stored in the storage unit 12A. As an example, the storage unit 12A stores information on an approximate function that approximates a characteristic curve. The approximate function information is created by the approximate function creating device 20 based on, for example, measurement data obtained by measuring the electrical characteristics of the cell 11 through experiments or the like. The characteristics of the lithium ion battery 10 are different for each type such as a model number. Further, the lithium ion battery 10 may have different characteristics for each manufactured unit or individual. The approximate function information may be created by obtaining measurement data for each type of lithium ion battery 10 or may be created by obtaining measurement data for each unit or individual. The lithium ion battery 10 is mounted and used in various devices 21. The control chip 12 functions as a processing unit that performs various controls by operating various programs. For example, the control chip 12 has an estimation unit 12B as a processing unit by operating a charging rate estimation program. The control chip 12 estimates the charging rate based on various information such as the terminal voltage of the lithium ion battery 10 and information on the approximate function stored in the storage unit 12A. For example, the control chip 12 estimates the charging rate using a Kalman filter. The control chip 12 performs various controls. For example, when the charge rate of the cell 11 is estimated and the cell 11 is overcharged or overdischarged, the control chip 12 performs control to cut off the current from the terminal used for charging and discharging power.

なお、リチウムイオン電池10などの二次電池は、経時的に特性が変化する場合がある。このため、記憶部12Aに記憶された情報を更新可能としてもよい。例えば、制御チップ12の記憶部12Aに近似関数の情報を作成可能なプログラムを記憶させ、制御チップ12がセル11の特性曲線を測定し、測定された測定データに基づいて近似関数の情報を作成してもよい。この場合、制御チップ12またはリチウムイオン電池10が近似関数作成装置20として機能する。また、ネットワーク23を介して近似関数の情報を更新可能としてもよい。例えば、セル11の特性曲線を測定した測定データを機器21がネットワーク23を介して近似関数作成装置20へ送信する。近似関数作成装置20が、測定データに基づいて作成した近似関数の情報をネットワーク23を介して返信して、記憶部12Aに記憶された情報を更新するようにしてもよい。   Note that characteristics of the secondary battery such as the lithium ion battery 10 may change over time. For this reason, the information stored in the storage unit 12A may be updatable. For example, the storage unit 12A of the control chip 12 stores a program that can generate information on the approximate function, the control chip 12 measures the characteristic curve of the cell 11, and generates information on the approximate function based on the measured measurement data. May be. In this case, the control chip 12 or the lithium ion battery 10 functions as the approximate function creation device 20. Further, the approximate function information may be updated via the network 23. For example, the device 21 transmits measurement data obtained by measuring the characteristic curve of the cell 11 to the approximate function creation device 20 via the network 23. The approximate function creation device 20 may return information on the approximate function created based on the measurement data via the network 23 to update the information stored in the storage unit 12A.

次に、実施例1に係る近似関数作成装置20について説明する。図2は、近似関数作成装置の全体構成を示す図である。図2に示すように、近似関数作成装置20は、外部I/F(interface)部30と、表示部31と、入力部32と、記憶部33と、制御部34とを有する。なお、近似関数作成装置20は、図2に示した機能部以外にも既知のコンピュータが有する各種の機能部を有してもよい。   Next, the approximate function creation device 20 according to the first embodiment will be described. FIG. 2 is a diagram showing the overall configuration of the approximate function creation device. As shown in FIG. 2, the approximate function creating apparatus 20 includes an external I / F (interface) unit 30, a display unit 31, an input unit 32, a storage unit 33, and a control unit 34. Note that the approximate function creation device 20 may have various functional units included in a known computer in addition to the functional units illustrated in FIG.

外部I/F部30は、他の装置と情報を入出力するインタフェースである。例えば、外部I/F部30は、他の装置との間で通信制御を行うインタフェースである。かかる外部I/F部30の一態様としては、LANカードなどのネットワークインタフェースカードを採用できる。例えば、外部I/F部30は、実験等によって測定されたセル11の電気的な特性を示す測定データ40を他の装置から受信する。なお、外部I/F部30は、例えば、USB(Universal Serial Bus)ポートなどのインタフェースであってもよい。   The external I / F unit 30 is an interface that inputs and outputs information with other devices. For example, the external I / F unit 30 is an interface that performs communication control with other devices. As an aspect of the external I / F unit 30, a network interface card such as a LAN card can be employed. For example, the external I / F unit 30 receives measurement data 40 indicating the electrical characteristics of the cell 11 measured by an experiment or the like from another device. The external I / F unit 30 may be an interface such as a USB (Universal Serial Bus) port.

表示部31は、各種情報を表示する表示デバイスである。表示部31としては、LCD(Liquid Crystal Display)やCRT(Cathode Ray Tube)などの表示デバイスが挙げられる。表示部31は、各種情報を表示する。   The display unit 31 is a display device that displays various types of information. Examples of the display unit 31 include display devices such as an LCD (Liquid Crystal Display) and a CRT (Cathode Ray Tube). The display unit 31 displays various information.

入力部32は、各種の情報を入力する入力デバイスである。例えば、入力部32としては、マウスやキーボードなどの入力デバイスが挙げられる。入力部32は、管理者などからの操作入力を受付け、受付けた操作内容を示す操作情報を制御部34に入力する。   The input unit 32 is an input device that inputs various types of information. For example, examples of the input unit 32 include input devices such as a mouse and a keyboard. The input unit 32 receives an operation input from an administrator or the like, and inputs operation information indicating the received operation content to the control unit 34.

記憶部33は、各種のデータを記憶する記憶デバイスである。例えば、記憶部33は、ハードディスク、SSD(Solid State Drive)、光ディスクなどの記憶装置である。なお、記憶部33は、RAM(Random Access Memory)、フラッシュメモリ、NVSRAM(Non Volatile Static Random Access Memory)などのデータを書き換え可能な半導体メモリであってもよい。   The storage unit 33 is a storage device that stores various data. For example, the storage unit 33 is a storage device such as a hard disk, an SSD (Solid State Drive), or an optical disk. Note that the storage unit 33 may be a semiconductor memory that can rewrite data, such as a random access memory (RAM), a flash memory, and a non-volatile static random access memory (NVSRAM).

記憶部33は、制御部34で実行されるOS(Operating System)や各種プログラムを記憶する。例えば、記憶部33は、近似関数の作成に用いる各種のプログラムを記憶する。さらに、記憶部33は、制御部34で実行されるプログラムで用いられる各種データを記憶する。例えば、記憶部33は、測定データ40を記憶する。なお、記憶部33は、上記に例示した情報以外にも、他の電子データを併せて記憶することもできる。   The storage unit 33 stores an OS (Operating System) executed by the control unit 34 and various programs. For example, the storage unit 33 stores various programs used for creating the approximate function. Furthermore, the storage unit 33 stores various data used in the program executed by the control unit 34. For example, the storage unit 33 stores the measurement data 40. The storage unit 33 can also store other electronic data in addition to the information exemplified above.

測定データ40は、セル11の電気的な特性を実験等によって測定したデータである。本実施例では、測定データ40として、リチウムイオン電池10の電力の充放電に用いる端子を流れる電流を切り替えた際の端子の電圧変化を測定した電圧変化データ41を用意する。例えば、リチウムイオン電池10の端子を流れる電流を、符号を変えずに値の大きさを変えた際のリチウムイオン電池10の端子の電圧の変化を測定した電圧変化データ41を用意する。また、リチウムイオン電池10の端子を流れる電流の符号を変えた際のリチウムイオン電池10の端子の電圧の変化を測定した電圧変化データ41を用意する。例えば、リチウムイオン電池10の端子を流れる電流を5Aから10Aに切り替えた場合と、電流を−10Aから10Aに切り替えた場合の電圧変化を測定した電圧変化データ41を用意する。なお、本実施例では、電流の符号がプラスの場合、リチウムイオン電池10へ電流が流れる充電状態であり、電流の符号がマイナスの場合、リチウムイオン電池10から電流が流れる放電状態であるものとする。   The measurement data 40 is data obtained by measuring the electrical characteristics of the cell 11 through experiments or the like. In the present embodiment, as the measurement data 40, voltage change data 41 obtained by measuring the voltage change of the terminal when the current flowing through the terminal used for charging / discharging the power of the lithium ion battery 10 is prepared. For example, voltage change data 41 is prepared by measuring changes in the voltage of the terminal of the lithium ion battery 10 when the value of the current flowing through the terminal of the lithium ion battery 10 is changed without changing the sign. In addition, voltage change data 41 is prepared by measuring changes in the voltage of the terminal of the lithium ion battery 10 when the sign of the current flowing through the terminal of the lithium ion battery 10 is changed. For example, voltage change data 41 is prepared by measuring the voltage change when the current flowing through the terminal of the lithium ion battery 10 is switched from 5A to 10A and when the current is switched from -10A to 10A. In this embodiment, when the sign of current is positive, it is in a charged state in which current flows to the lithium ion battery 10, and when the sign of current is negative, it is in a discharged state in which current flows from the lithium ion battery 10. To do.

また、本実施例では、測定データ40として、リチウムイオン電池10の充電率を0%から一定電流で100%まで充電した際のOCVと充電率の相関関係を示した特性曲線の特性データ42を用意する。また、本実施例では、測定データ40として、リチウムイオン電池10の充電率を100%から一定電流で0%まで放電した際のOCVと充電率の相関関係を示した特性曲線の特性データ42を用意する。   Further, in this embodiment, as the measurement data 40, the characteristic data 42 of the characteristic curve showing the correlation between the OCV and the charging rate when the charging rate of the lithium ion battery 10 is charged from 0% to 100% with a constant current is used. prepare. In the present embodiment, as the measurement data 40, characteristic data 42 of a characteristic curve showing the correlation between the OCV and the charging rate when the charging rate of the lithium ion battery 10 is discharged from 100% to 0% at a constant current is used. prepare.

図3は、特性曲線の一例を示す図である。図3の例は、特性データ42に記憶された、リチウムイオン電池10のOCVと充電率の相関関係を示す特性曲線を示したものである。図3に示すように、リチウムイオン電池10のOCVと充電率には、相関関係がある。リチウムイオン電池10は、充電時と放電時でヒステリシスがあり、充電時と放電時でOCVと充電率との相関関係が異なる場合がある。図3の例では、充電時の特性曲線60と、放電時の特性曲線61が示されている。また、特性曲線60、61は、充電率が高い範囲および充電率が低い範囲で、充電率の変化に対するOCVの変化が大きくなっている。   FIG. 3 is a diagram illustrating an example of a characteristic curve. The example of FIG. 3 shows a characteristic curve indicating the correlation between the OCV of the lithium ion battery 10 and the charging rate stored in the characteristic data 42. As shown in FIG. 3, there is a correlation between the OCV and the charging rate of the lithium ion battery 10. The lithium ion battery 10 has hysteresis during charging and discharging, and the correlation between the OCV and the charging rate may be different between charging and discharging. In the example of FIG. 3, a characteristic curve 60 at the time of charging and a characteristic curve 61 at the time of discharging are shown. In the characteristic curves 60 and 61, the change in the OCV with respect to the change in the charge rate is large in the range where the charge rate is high and the range where the charge rate is low.

制御部34は、各種の処理手順を規定したプログラムや制御データを格納するための内部メモリを有し、これらによって種々の処理を実行する。制御部34は、図2に示すように、取得部50と、第1算出部51と、第2算出部52と、第3算出部53と、生成部54と、出力部55とを有する。   The control unit 34 has an internal memory for storing programs defining various processing procedures and control data, and executes various processes using these. As shown in FIG. 2, the control unit 34 includes an acquisition unit 50, a first calculation unit 51, a second calculation unit 52, a third calculation unit 53, a generation unit 54, and an output unit 55.

取得部50は、各種情報を取得する処理部である。例えば、取得部50は、特性曲線のデータを取得する。例えば、取得部50は、記憶部33に記憶された特性データ42を読み出すことにより、特性曲線のデータを取得する。なお、本実施例では、特性データ42を記憶部33に記憶しているため、取得部50は、記憶部33から特性データ42を取得するものとするが、これに限定されるものではない。特性データ42が他の装置に記憶されている場合、取得部50は、ネットワーク23を介して当該他の装置から特性データ42を取得してもよい。   The acquisition unit 50 is a processing unit that acquires various types of information. For example, the acquisition unit 50 acquires characteristic curve data. For example, the acquiring unit 50 acquires characteristic curve data by reading the characteristic data 42 stored in the storage unit 33. In the present embodiment, since the characteristic data 42 is stored in the storage unit 33, the acquisition unit 50 acquires the characteristic data 42 from the storage unit 33. However, the present invention is not limited to this. When the characteristic data 42 is stored in another device, the acquisition unit 50 may acquire the characteristic data 42 from the other device via the network 23.

第1算出部51は、各種の算出を行う処理部である。例えば、第1算出部51は、電圧変化データ41から、充電時の特性曲線と放電時の特性曲線からヒステリシスによる誤差を算出する。例えば、第1算出部51は、電圧変化データ41から、リチウムイオン電池10から流れる電流を切り替えた際の電流の変化量と電圧の変化量を求める。例えば、電流を5Aから10Aに切り替えた場合、電流の変化量は、5Aである。この電流を5Aから10Aに切り替えた際の電圧の変化量をV1とする。また、電流を−10Aから10Aに切り替えた場合、電流の変化量は、20Aである。この電流を−10Aから10Aに切り替えた際の電圧の変化量をV2とする。電流変化に関して発生する抵抗をR0とし、流れる電流の符号が切り替ることによるヒステリシスによる誤差をHとすると、以下の(1)、(2)式の関係が成り立つ。 The first calculation unit 51 is a processing unit that performs various calculations. For example, the first calculator 51 calculates an error due to hysteresis from the voltage change data 41 from the characteristic curve during charging and the characteristic curve during discharging. For example, the first calculation unit 51 obtains the amount of change in current and the amount of change in voltage when the current flowing from the lithium ion battery 10 is switched from the voltage change data 41. For example, when the current is switched from 5A to 10A, the amount of change in the current is 5A. The amount of change in voltage when this current is switched from 5 A to 10 A is defined as V 1 . Further, when the current is switched from −10 A to 10 A, the amount of change in the current is 20 A. The amount of change in voltage when this current is switched from −10 A to 10 A is V 2 . When the resistance generated with respect to the current change is R 0 and the error due to the hysteresis due to the switching of the sign of the flowing current is H, the following relations (1) and (2) are established.

5×R0=V1 (1)
20×R0+H=V2 (2)
5 × R 0 = V 1 (1)
20 × R 0 + H = V 2 (2)

この(1)式でR0に乗算される係数「5」は、電流を5Aから10Aに切り替えた場合の電流の変化量である。(2)式でR0に乗算される係数「20」は、電流を−10Aから10Aに切り替えた場合の電流の変化量である。 The coefficient “5” multiplied by R 0 in this equation (1) is the amount of change in current when the current is switched from 5A to 10A. The coefficient “20” multiplied by R 0 in equation (2) is the amount of change in current when the current is switched from −10 A to 10 A.

この(1)、(2)式から、以下の(3)式が求められる。すなわち、リチウムイオン電池10の端子を流れる電流を、符号を変えて変化させた場合の電流の変化量および電圧の変化量と、端子を流れる電流を、符号を変えずに変化させた場合、電流の変化量および電圧の変化量からヒステリシスによる誤差Hを算出できる。   From the equations (1) and (2), the following equation (3) is obtained. That is, when the current flowing through the terminal of the lithium ion battery 10 is changed by changing the sign, the current change amount and the voltage change amount, and when the current flowing through the terminal is changed without changing the sign, The error H due to hysteresis can be calculated from the amount of change and the amount of change in voltage.

H=V2−4×V1 (3) H = V 2 -4 × V 1 (3)

第1算出部51は、リチウムイオン電池10の端子を流れる電流を、符号を変えて変化させた場合の電流の変化量および電圧の変化量と、端子を流れる電流を、符号を変えずに変化させた場合の電流の変化量および電圧の変化量と、から誤差Hを算出する。なお、本実施例では、端子を流れる電流を変化させた際の電流の変化量および電圧の変化量からヒステリシスによる誤差Hを算出するが、特性データ42からヒステリシスによる誤差Hを求めてもよい。例えば、第1算出部51は、図3に示す充電時の特性曲線と放電時の特性曲線の電圧の差をヒステリシスによる誤差Hとしてもよい。   The first calculation unit 51 changes the amount of change in current and voltage when the current flowing through the terminal of the lithium ion battery 10 is changed by changing the sign, and the amount of change in current without changing the sign. The error H is calculated from the amount of change in current and the amount of change in voltage. In this embodiment, the error H due to hysteresis is calculated from the amount of change in current and the amount of change in voltage when the current flowing through the terminal is changed, but the error H due to hysteresis may be obtained from the characteristic data 42. For example, the first calculation unit 51 may set the difference in voltage between the characteristic curve during charging and the characteristic curve during discharging shown in FIG.

ここで、図3に示すように、特性曲線は、充電率が高い高充電率部分および充電率が低い低充電率部分で充電率の変化に対するOCVの変化が大きく、中央部分の充電率の変化に対するOCVの変化が小さくなっている。このような特性曲線を直線近似する場合、高充電率部分、中央部分、低充電率部分を分けて直線近似を行う。   Here, as shown in FIG. 3, the characteristic curve shows a large change in the OCV with respect to the change in the charge rate in the high charge rate portion where the charge rate is high and the low charge rate portion where the charge rate is low, and the change in the charge rate in the center portion The change in OCV with respect to is small. When linearly approximating such a characteristic curve, linear approximation is performed by dividing a high charge rate portion, a central portion, and a low charge rate portion.

第2算出部52は、特性曲線の中央部分の直線近似を行う処理部である。まず、第2算出部52は、高充電率部分と中央部分、および、中央部分と低充電率部分の境界を特定する。例えば、第2算出部52は、特性曲線の充電率に対する開回路電圧の変化率が所定の閾値以下となる範囲を中央部分と特定する。この閾値は、充電率に対する開回路電圧の変化が安定した範囲と見なせる値に管理者等により定められ、例えば、500とする。例えば、第2算出部52は、特性曲線の中央である充電率50%から高充電率側へ、順に充電率に対する開回路電圧の変化率を求め、最初に変化率が閾値を超えた充電率を高充電率側の境界とする。また、第2算出部52は、充電率50%から低充電率側へ、順に充電率に対する開回路電圧の変化率を求め、最初に変化率が閾値を超えた充電率を低充電率側の境界とする。この境界の判定に用いる特性曲線は、充電時の特性曲線であってもよく、放電時の特性曲線であってもよく、充電時の特性曲線と充電時の特性曲線を平均化した特性曲線であってもよい。また、第2算出部52は、充電時の特性曲線と放電時の特性曲線の何れか一方で変化率が閾値を超えた充電率を境界としてもよく、両方で変化率が閾値を超えた充電率を境界としてもよい。この高充電率側の境界と低充電率側の境界の間が中央部分となる。   The second calculation unit 52 is a processing unit that performs linear approximation of the central portion of the characteristic curve. First, the 2nd calculation part 52 specifies the boundary of a high charge rate part and a center part, and a center part and a low charge rate part. For example, the second calculation unit 52 specifies a range in which the change rate of the open circuit voltage with respect to the charging rate of the characteristic curve is equal to or less than a predetermined threshold as the central portion. This threshold is determined by an administrator or the like to a value that can be regarded as a stable range of change in the open circuit voltage with respect to the charging rate, and is set to 500, for example. For example, the second calculation unit 52 obtains the change rate of the open circuit voltage with respect to the charge rate in order from the charge rate of 50%, which is the center of the characteristic curve, to the high charge rate side, and the charge rate at which the change rate exceeded the threshold first. Is the boundary on the high charge rate side. Further, the second calculation unit 52 obtains the change rate of the open circuit voltage with respect to the charge rate in order from the charge rate of 50% to the low charge rate side, and first calculates the charge rate at which the change rate exceeded the threshold value on the low charge rate side. Boundary. The characteristic curve used for determining the boundary may be a characteristic curve at the time of charging or a characteristic curve at the time of discharging, and a characteristic curve obtained by averaging the characteristic curve at the time of charging and the characteristic curve at the time of charging. There may be. In addition, the second calculation unit 52 may use a charge rate at which the rate of change exceeds a threshold value in either one of the characteristic curve at the time of charging or the characteristic curve at the time of discharge as a boundary. The rate may be the boundary. A central portion is between the boundary on the high charge rate side and the boundary on the low charge rate side.

第2算出部52は、中央部分の充電時の特性曲線および放電時の特性曲線に対して、ヒステリシスによる誤差Hに応じた補正を行う。例えば、第2算出部52は、充電時の特性曲線の各充電率の電圧からH/2を減算し、放電時の特性曲線の各充電率の電圧にH/2を加算して、ヒステリシスによる誤差Hを半分ずつ補正する。なお、第2算出部52は、充電時の特性曲線と放電時の特性曲線の何れか一方のみを補正してもよい。例えば、第2算出部52は、充電時の特性曲線の各充電率の電圧からHを減算する補正を行ってもよい。   The second calculation unit 52 corrects the characteristic curve at the time of charging and the characteristic curve at the time of discharging according to the error H due to hysteresis. For example, the second calculation unit 52 subtracts H / 2 from the voltage of each charging rate of the characteristic curve at the time of charging, adds H / 2 to the voltage of each charging rate of the characteristic curve at the time of discharging, The error H is corrected by half. The second calculation unit 52 may correct only one of the characteristic curve during charging and the characteristic curve during discharging. For example, the 2nd calculation part 52 may perform correction | amendment which subtracts H from the voltage of each charge rate of the characteristic curve at the time of charge.

第2算出部52は、補正後の充電時の特性曲線と放電時の特性曲線から平均の特性曲線のデータを算出する。例えば、第2算出部52は、低充電率側の境界から所定の刻み幅毎に、充電時の特性曲線および放電時の特性曲線において、刻み幅の充電率よりも大きい直近の充電率と、当該直近の充電率での電圧とを特定する。そして、第2算出部52は、刻み幅毎に、特定した充電率および電圧をそれぞれ平均して平均の特性曲線のデータを算出する。例えば、低充電率側の境界の充電率が0.1であり、刻み幅を0.05とした場合、充電時の特性曲線で0.15よりも大きい直近の充電率がaであり、電圧がcであり、放電時の特性曲線で0.15よりも大きい直近の充電率がbであり、電圧がdであるものとする。この場合、第2算出部52は、充電率を(a+b)/2とし、電圧を(c+d)/2として平均の特性曲線のデータを算出する。   The second calculator 52 calculates data of an average characteristic curve from the corrected characteristic curve during charging and the characteristic curve during discharging. For example, the second calculation unit 52, for each predetermined step width from the boundary on the low charge rate side, in the characteristic curve at the time of charging and the characteristic curve at the time of discharging, the latest charging rate larger than the charging rate of the step size, The voltage at the latest charging rate is specified. Then, the second calculating unit 52 calculates the average characteristic curve data by averaging the specified charging rate and voltage for each step size. For example, when the charge rate at the low charge rate side boundary is 0.1 and the step size is 0.05, the most recent charge rate greater than 0.15 in the characteristic curve at the time of charge is a, and the voltage , C, the latest charge rate greater than 0.15 in the characteristic curve during discharge is b, and the voltage is d. In this case, the second calculation unit 52 calculates average characteristic curve data by setting the charging rate to (a + b) / 2 and the voltage to (c + d) / 2.

第2算出部52は、中央部分の範囲内で中間点を移動させる。そして、第2算出部52は、中間点の位置毎に、中間点を境として中央部分の範囲を分割し、分割した各分割範囲の平均の特性曲線のデータをそれぞれ近似した近似直線と、当該近似直線の近似に用いた平均の特性曲線のデータとの類似度を算出する。例えば、第2算出部52は、低充電率側の境界から所定の刻み幅で中間点を移動させる。第2算出部52は、中間点が移動する毎に、中間点を境として中央部分の範囲を分割する。そして、第2算出部52は、分割された分割範囲毎に、分割範囲内の平均の特性曲線のデータをそれぞれ近似した近似直線を求める。   The second calculation unit 52 moves the intermediate point within the range of the central portion. Then, the second calculation unit 52 divides the range of the central portion for each intermediate point position with the intermediate point as a boundary, and approximate straight lines approximating the average characteristic curve data of each divided range, The degree of similarity with the data of the average characteristic curve used for approximation of the approximate line is calculated. For example, the second calculation unit 52 moves the intermediate point with a predetermined step size from the boundary on the low charge rate side. Each time the intermediate point moves, the second calculator 52 divides the range of the central portion with the intermediate point as a boundary. Then, for each divided range, the second calculation unit 52 obtains an approximate line that approximates the data of the average characteristic curve in the divided range.

図4は、分割範囲毎の近似直線の一例を模式的に示した図である。図4の例では、低充電率側の境界C1と、高充電率側の境界C2とが示されている。また、図4の例では、境界C1と境界C2の間の中央部分の範囲が中間点により2つの分割範囲に分割され、分割範囲毎に、分割範囲の特性曲線のデータを近似した近似直線L1、L2が示されている。   FIG. 4 is a diagram schematically illustrating an example of an approximate line for each division range. In the example of FIG. 4, a boundary C1 on the low charging rate side and a boundary C2 on the high charging rate side are shown. In the example of FIG. 4, the range of the central portion between the boundary C1 and the boundary C2 is divided into two divided ranges by an intermediate point, and an approximate straight line L1 that approximates the characteristic curve data of the divided ranges for each divided range. , L2 is shown.

第2算出部52は、分割範囲毎に、近似直線と平均の特性曲線との類似度を求める。例えば、第2算出部52は、近似直線と平均の特性曲線との類似度として決定係数R2を求める。なお、類似度として用いる指標は、決定係数R2に限定されるものではなく近似直線と平均の類似度合いを示すものであれば何れかあってもよい。 The second calculator 52 obtains the similarity between the approximate straight line and the average characteristic curve for each divided range. For example, the second calculation unit 52 obtains the determination coefficient R 2 as the similarity between the approximate line and the average characteristic curve. The index used as the degree of similarity is not limited to the determination coefficient R 2 , and may be any as long as it shows the degree of similarity between the approximate line and the average.

第2算出部52は、類似度が最も高い中間点を境とした各分割範囲の近似直線を、中央部分の範囲の近似直線と特定する。例えば、第2算出部52は、決定係数R2が最も高い中間点を境とした各分割範囲の近似直線を、中央部分の範囲の近似直線と特定する。 The second calculation unit 52 identifies the approximate straight line of each divided range with the middle point having the highest similarity as the boundary as the approximate straight line of the central portion range. For example, the second calculation unit 52 identifies an approximate straight line in each divided range with an intermediate point having the highest determination coefficient R 2 as an approximate straight line in the range of the central portion.

なお、本実施例では、図4に示すように、中間点を1つとして中央部分の範囲が中間点により2つの分割範囲に分割する場合を例に説明するが、中間点を複数設けても中央部分の範囲を3以上の分割範囲に分割してもよい。例えば、第2算出部52は、中央部分の範囲内で、複数の中間点を移動させて複数の中間点をそれぞれ境として中央部分の範囲を分割する。そして、第2算出部52は、複数の中間点の位置毎に、各分割範囲の平均の特性曲線のデータをそれぞれ近似した近似直線と、当該近似直線の近似に用いた平均の特性曲線のデータとの類似度を算出してもよい。例えば、中間点を第1の中間点と第2の中間点の2つ設ける場合、第2算出部52は、第1の中間点を低充電率側の境界から高充電率側の境界へ移動させる。そして、第2算出部52は、第1の中間点を移動させる毎に、第2の中間点を第1の中間点の位置から高充電率側の境界へ移動させて、中央部分の範囲を第1の中間点および第2の中間点を境として3つの分割範囲に分割する。そして、第2算出部52は、第1の中間点および第2の中間点の位置毎に、各分割範囲の平均の特性曲線のデータをそれぞれ近似した近似直線と、当該近似直線の近似に用いた平均の特性曲線のデータとの類似度を算出してもよい。   In the present embodiment, as shown in FIG. 4, an example is described in which the middle portion is one and the range of the central portion is divided into two divided ranges by the middle point, but a plurality of middle points may be provided. The range of the central portion may be divided into three or more divided ranges. For example, the second calculation unit 52 moves the plurality of intermediate points within the range of the central portion, and divides the range of the central portion with the plurality of intermediate points as boundaries. Then, the second calculation unit 52 approximates the average characteristic curve data of each divided range for each of the plurality of intermediate points, and the average characteristic curve data used to approximate the approximate line. The degree of similarity may be calculated. For example, when two intermediate points are provided, the first intermediate point and the second intermediate point, the second calculation unit 52 moves the first intermediate point from the low charging rate side boundary to the high charging rate side boundary. Let Then, every time the second calculation unit 52 moves the first intermediate point, the second calculation unit 52 moves the second intermediate point from the position of the first intermediate point to the boundary on the high charge rate side, thereby reducing the range of the central portion. The first intermediate point and the second intermediate point are used as a boundary to divide into three divided ranges. Then, the second calculator 52 uses, for each of the positions of the first intermediate point and the second intermediate point, an approximate straight line approximating the data of the average characteristic curve of each divided range, and an approximation of the approximate straight line. The similarity with the average characteristic curve data may be calculated.

第3算出部53は、高充電率部分および低充電率部分の直線近似を行う処理部である。例えば、第3算出部53は、充電率が高充電率側の境界以上の範囲のデータを近似した近似直線と、充電率が低充電率側の境界以下の範囲のデータを近似した近似直線を求める。   The third calculation unit 53 is a processing unit that performs linear approximation of the high charge rate portion and the low charge rate portion. For example, the third calculation unit 53 calculates an approximate straight line approximating data in a range where the charging rate is equal to or higher than the boundary on the high charging rate side, and an approximate straight line approximating data in the range where the charging rate is equal to or lower than the boundary on the low charging rate side. Ask.

ここで、充電時の特性曲線には、低充電率部分に充電を開始した際の電流変化による過渡特性の影響が含まれる場合がある。また、放電時の特性曲線には、高充電率部分に放電を開始した際の電流変化による過渡特性の影響が含まれる場合がある。   Here, the characteristic curve at the time of charging may include the influence of the transient characteristic due to the current change when charging is started in the low charging rate portion. In addition, the characteristic curve at the time of discharge may include the influence of the transient characteristic due to the current change when the discharge is started in the high charge rate portion.

そこで、第3算出部53は、低充電率部分については放電時の特性曲線のデータを近似した近似直線を算出する。また、第3算出部53は、高充電率部分については充電時の特性曲線のデータを近似した近似直線を算出する。   Therefore, the third calculation unit 53 calculates an approximate straight line that approximates the data of the characteristic curve at the time of discharging for the low charge rate portion. In addition, the third calculation unit 53 calculates an approximate straight line that approximates the characteristic curve data at the time of charging for the high charge rate portion.

図5は、低充電率部分の近似直線の一例を模式的に示した図である。図5の例では、放電時の特性曲線のデータを近似した近似直線L3が示されている。図6は、高充電率部分の近似直線の一例を模式的に示した図である。図6の例では、充電時の特性曲線のデータを近似した近似直線L4が示されている。   FIG. 5 is a diagram schematically illustrating an example of an approximate line of the low charge rate portion. In the example of FIG. 5, an approximate straight line L3 that approximates the data of the characteristic curve at the time of discharge is shown. FIG. 6 is a diagram schematically illustrating an example of an approximate line of the high charge rate portion. In the example of FIG. 6, an approximate straight line L4 that approximates the characteristic curve data during charging is shown.

生成部54は、特性曲線を近似した近似関数を生成する処理部である。例えば、生成部54は、第2算出部52により特定された近似直線と、第3算出部53により算出された低充電率部分の近似直線および高充電率部分の近似直線を用いて特性曲線を近似した近似関数を生成する。この際、生成部54は、範囲が隣り合う近似直線の交点を近似直線の切替点として近似関数を生成する。   The generation unit 54 is a processing unit that generates an approximate function approximating the characteristic curve. For example, the generation unit 54 generates a characteristic curve using the approximate straight line specified by the second calculation unit 52 and the approximate straight line of the low charge rate portion and the approximate straight line of the high charge rate portion calculated by the third calculation unit 53. Generate an approximate function. At this time, the generation unit 54 generates an approximate function using an intersection of approximate lines with adjacent ranges as a switching point of the approximate lines.

図7は、切替点の一例を模式的に示した図である。図7の例では、近似直線L1と、近似直線L2とが交差する交点を近似直線の切替点とする。   FIG. 7 is a diagram schematically illustrating an example of the switching point. In the example of FIG. 7, an intersection point between the approximate line L1 and the approximate line L2 is set as a switch point of the approximate line.

図8は、生成された近似関数の一例を示した図である。図8の例では、放電時の特性曲線と充電時の特性曲線から生成された近似関数による線L5が示さている。   FIG. 8 is a diagram illustrating an example of the generated approximate function. In the example of FIG. 8, a line L5 by an approximate function generated from the characteristic curve at the time of discharging and the characteristic curve at the time of charging is shown.

生成部54は、生成した近似関数に、誤差Hに応じた変更を行って放電時の特性曲線を近似した放電用の近似関数と、充電用の特性曲線を近似した充電用の近似関数とを生成する。例えば、生成部54は、近似関数の各充電率の電圧からH/2を減算して、放電用の近似関数を生成する。また、生成部54は、近似関数の各充電率の電圧にH/2を加算して、充電用の近似関数を生成する。この誤差Hに応じた変更は、第2算出部52による誤差Hの補正と対応する割合で変更とすることが好ましい。例えば、第2算出部52が生成された特性曲線をH/2減算し、放電時の特性曲線をH/2を加算する補正をしている場合、生成部54は、近似関数をH/2減算する変更を行って放電用の近似関数を生成し、近似関数をH/2加算する変更を行って充電用の近似関数を生成する。   The generating unit 54 changes the generated approximate function according to the error H to approximate an approximate function for discharging and approximates a characteristic curve for discharging, and an approximate function for charging that approximates the characteristic curve for charging. Generate. For example, the generation unit 54 generates an approximate function for discharging by subtracting H / 2 from the voltage of each charging rate of the approximate function. Further, the generation unit 54 adds H / 2 to the voltage of each charging rate of the approximate function to generate an approximate function for charging. The change corresponding to the error H is preferably changed at a rate corresponding to the correction of the error H by the second calculation unit 52. For example, when the second calculating unit 52 performs correction by subtracting H / 2 from the generated characteristic curve and adding H / 2 to the characteristic curve during discharge, the generating unit 54 sets the approximate function to H / 2. The subtracting change is performed to generate an approximate function for discharging, and the subtracting function is added to H / 2 to generate an approximate function for charging.

図9は、放電用の近似関数と充電用の近似関数の一例を示した図である。図9の例では、近似関数による線L5をH/2減算した放電用の近似関数の線L6が示されている。また、図9の例では、近似関数による線L5をH/2加算した充電用の近似関数の線L7が示されている。   FIG. 9 is a diagram illustrating an example of an approximate function for discharging and an approximate function for charging. In the example of FIG. 9, a line L6 of the approximate function for discharge obtained by subtracting H / 2 from the line L5 by the approximate function is shown. In the example of FIG. 9, an approximate function line L7 for charging obtained by adding H / 2 to the line L5 by the approximate function is shown.

出力部55は、各種の出力を行う処理部である。例えば、出力部55は、生成部54により生成された近似関数の情報を出力する。出力部55による出力は、近似関数の情報を記憶部33に記憶させるものであってもよい。また、出力部55による出力は、表示部31に情報を表示させるものであってもよい。また、出力部55による出力は、外部I/F部30を介して外部の装置へ情報を出力するものであってもよい。   The output unit 55 is a processing unit that performs various outputs. For example, the output unit 55 outputs information on the approximate function generated by the generation unit 54. The output by the output unit 55 may store information on the approximate function in the storage unit 33. Further, the output by the output unit 55 may display information on the display unit 31. The output by the output unit 55 may output information to an external device via the external I / F unit 30.

この出力された近似関数の情報は、リチウムイオン電池10の制御チップ12の記憶部12Aに格納される。制御チップ12の推定部12Bは、記憶部12Aから近似関数の情報を読み出すことにより取得する。そして、推定部12Bは、取得した近似関数の情報を用いてリチウムイオン電池10の充電率を推定する。例えば、推定部12Bは、リチウムイオン電池10の電流および端子電圧を測定する。この端子電圧には、ノイズ等が含まれている。そこれ、また、推定部12Bは、端子電圧から近似関数の情報を用いてカルマンフィルタにより充電率を推定する。このため、制御チップ12は、近似関数の情報を用いて充電率を精度よく推定できる。   Information on the output approximate function is stored in the storage unit 12 </ b> A of the control chip 12 of the lithium ion battery 10. The estimation unit 12B of the control chip 12 acquires the information of the approximate function from the storage unit 12A. And the estimation part 12B estimates the charging rate of the lithium ion battery 10 using the information of the acquired approximate function. For example, the estimation unit 12B measures the current and terminal voltage of the lithium ion battery 10. This terminal voltage includes noise and the like. Accordingly, the estimation unit 12B estimates the charging rate by the Kalman filter using the information of the approximate function from the terminal voltage. For this reason, the control chip 12 can estimate the charging rate with high accuracy using the information of the approximate function.

なお、制御部34には、各種の集積回路や電子回路を採用できる。また、制御部34が有する機能部の一部を別の集積回路や電子回路とすることもできる。例えば、集積回路としては、ASIC(Application Specific Integrated Circuit)が挙げられる。また、電子回路としては、CPU(Central Processing Unit)やMPU(Micro Processing Unit)などが挙げられる。   Note that various integrated circuits and electronic circuits can be employed for the control unit 34. Further, a part of the functional unit included in the control unit 34 may be another integrated circuit or an electronic circuit. For example, an ASIC (Application Specific Integrated Circuit) is an example of the integrated circuit. Examples of the electronic circuit include a central processing unit (CPU) and a micro processing unit (MPU).

[処理の流れ]
続いて、本実施例に係る近似関数作成装置20が近似関数を作成する近似関数作成処理の流れについて説明する。図10は、近似関数作成処理の手順を示すフローチャートである。この近似関数作成処理は、所定のタイミング、例えば、入力部32あるいは外部I/F部30を介して他の装置から近似関数の作成が指示されたタイミングで処理が開始される。
[Process flow]
Subsequently, a flow of an approximate function creation process in which the approximate function creation device 20 according to the present embodiment creates an approximate function will be described. FIG. 10 is a flowchart showing the procedure of the approximate function creation process. The approximate function creation process is started at a predetermined timing, for example, at a timing when the creation of an approximate function is instructed from another device via the input unit 32 or the external I / F unit 30.

図10に示すように、取得部50は、記憶部33に記憶された特性データ42を読み出すことにより、特性曲線のデータを取得する(S10)。第1算出部51は、充電時の特性曲線と放電時の特性曲線とのヒステリシスによる誤差Hを算出する(S11)。第2算出部52は、高充電率部分と中央部分および中央部分と低充電率部分の境界を特定する(S12)。第2算出部52は、中央部分の充電時の特性曲線および放電時の特性曲線に対して、ヒステリシスによる誤差Hに応じた補正を行い、平均の特性曲線を算出する(S13)。   As shown in FIG. 10, the acquisition unit 50 acquires characteristic curve data by reading the characteristic data 42 stored in the storage unit 33 (S10). The first calculator 51 calculates an error H due to hysteresis between the characteristic curve at the time of charging and the characteristic curve at the time of discharging (S11). The 2nd calculation part 52 specifies the boundary of a high charge rate part, a center part, and a center part and a low charge rate part (S12). The second calculator 52 corrects the characteristic curve at the time of charging and the characteristic curve at the time of discharging according to the error H due to hysteresis, and calculates an average characteristic curve (S13).

第2算出部52は、中央部分の範囲内で、中間点を移動させつつ中間点の位置毎に、中間点を境として前記範囲を分割した各分割範囲の特性曲線をそれぞれ近似した近似直線と、当該近似直線の近似に用いたデータとの類似度を算出する(S14)。そして、第2算出部52は、類似度が最も高い中間点を境とした各分割範囲の近似直線を特定する(S15)。第3算出部53は、高充電率部分のデータを近似した近似直線と、低充電率部分のデータを近似した近似直線を求める(S16)。   The second calculation unit 52 approximates each of the characteristic curves of the divided ranges obtained by dividing the range with the intermediate point as a boundary for each position of the intermediate point while moving the intermediate point within the range of the central portion. Then, the similarity with the data used for the approximation of the approximate straight line is calculated (S14). And the 2nd calculation part 52 specifies the approximate straight line of each division | segmentation range on the boundary of the intermediate point with the highest similarity (S15). The 3rd calculation part 53 calculates | requires the approximate straight line which approximated the data of the low charge rate part, and the approximate straight line which approximated the data of the high charge rate part (S16).

生成部54は、中央部分の近似直線、高充電率部分の近似直線および低充電率部分の近似直線を用いて特性曲線を近似した近似関数を生成する(S17)。生成部54は、生成した近似関数に、誤差Hに応じた変更を行って放電時の特性曲線を近似した放電用の近似関数と充電時の特性曲線を近似した充電用の近似関数とを生成し(S18)、処理を終了する。   The generation unit 54 generates an approximate function that approximates the characteristic curve using the approximate straight line in the central portion, the approximate straight line in the high charge rate portion, and the approximate straight line in the low charge rate portion (S17). The generating unit 54 changes the generated approximate function according to the error H to generate an approximate function for discharging that approximates a characteristic curve at the time of discharge and an approximate function for charging that approximates the characteristic curve at the time of charging. (S18), and the process ends.

[実施例1の効果]
上述してきたように、近似関数作成装置20は、二次電池の開回路電圧と当該二次電池の充電率との相関関係を示した特性曲線のデータを取得する。近似関数作成装置20は、取得された特性曲線のデータにより示される特性曲線の充電率に対する開回路電圧の変化率が所定の閾値以下となる範囲内で、中間点を移動させる。そして、近似関数作成装置20は、間点の位置毎に、中間点を境として前記範囲を分割した各分割範囲の特性曲線のデータをそれぞれ近似した近似直線と、当該近似直線の近似に用いた特性曲線のデータとの類似度を算出する。近似関数作成装置20は、算出される類似度が最も高い中間点を境とした各分割範囲の近似直線を用いて特性曲線を近似した近似関数を生成する。これにより、近似関数作成装置20は、特性曲線の近似を精度よく行うことができる。
[Effect of Example 1]
As described above, the approximate function creation device 20 acquires characteristic curve data indicating the correlation between the open circuit voltage of the secondary battery and the charging rate of the secondary battery. The approximate function creation device 20 moves the intermediate point within a range where the change rate of the open circuit voltage with respect to the charging rate of the characteristic curve indicated by the acquired characteristic curve data is equal to or less than a predetermined threshold. Then, the approximate function creation device 20 uses the approximate line obtained by approximating the characteristic curve data of each divided range obtained by dividing the range at the intermediate point for each position of the intermediate point, and used to approximate the approximate line. The similarity with the characteristic curve data is calculated. The approximate function creation device 20 generates an approximate function that approximates the characteristic curve using an approximate line in each divided range with the calculated intermediate point as the boundary. Thereby, the approximate function creating apparatus 20 can accurately approximate the characteristic curve.

また、近似関数作成装置20は、分割範囲が隣り合う近似直線の交点を近似直線の切替点とした近似関数を生成する。これにより、本実施例に係る近似関数作成装置20は、近似直線をスムーズに切り替えることができる近似関数を生成できる。   In addition, the approximate function creation device 20 generates an approximate function using an intersection of approximate lines with adjacent division ranges as a switching point of the approximate line. Thereby, the approximate function creation apparatus 20 according to the present embodiment can generate an approximate function that can smoothly switch the approximate line.

また、近似関数作成装置20は、二次電池の充電時の特性曲線と放電時の特性曲線とのヒステリシスによる誤差を算出する。近似関数作成装置20は、充電時の特性曲線のデータおよび放電時の特性曲線のデータに対して、誤差に応じた補正を行った後に平均化したデータを用いて近似直線と、当該近似直線の近似に用いたデータとの類似度を算出する。近似関数作成装置20は、近似関数に誤差に応じた変更を行って放電時の特性曲線を近似した近似関数と充電時の特性曲線を近似した近似関数とを生成する。これにより、近似関数作成装置20は、放電時の特性曲線を近似した放電用の近似関数と充電時の特性曲線を近似した充電用の近似関数とを生成できる。   Further, the approximate function creation device 20 calculates an error due to hysteresis between the characteristic curve at the time of charging the secondary battery and the characteristic curve at the time of discharging. The approximate function creation device 20 uses the data obtained by correcting the characteristic curve data at the time of charging and the data of the characteristic curve at the time of discharging after performing correction according to the error, and the approximate straight line and the approximate straight line. The similarity with the data used for the approximation is calculated. The approximate function creation device 20 changes the approximate function according to the error to generate an approximate function that approximates the characteristic curve during discharging and an approximate function that approximates the characteristic curve during charging. Thereby, the approximate function creating apparatus 20 can generate an approximate function for discharging that approximates a characteristic curve at the time of discharging and an approximate function for charging that approximates the characteristic curve at the time of charging.

また、近似関数作成装置20は、特性曲線の範囲よりも充電率が高い高充電率範囲については充電時の特性曲線のデータを近似した近似直線を算出する。近似関数作成装置20は、特性曲線の範囲よりも充電率が低い低充電率範囲については放電時の特性曲線のデータを近似した近似直線を算出する。近似関数作成装置20は、算出された高充電率範囲および低充電率範囲の近似直線をさらに用いて、近似関数を生成する。近似関数作成装置20は、高充電率範囲の近似直線と、最も高充電率範囲側の分割範囲の近似直線との交点、および、低充電率範囲の近似直線と、最も低充電率範囲側の分割範囲の近似直線との交点をそれぞれ近似直線の切替点とした近似関数を生成する。これにより、本実施例に係る近似関数作成装置20は、高充電率範囲および低充電率範囲も含んだ近似関数を精度よく生成できる。   Further, the approximate function creation device 20 calculates an approximate straight line that approximates the characteristic curve data during charging for a high charge rate range in which the charge rate is higher than the range of the characteristic curve. The approximate function creation device 20 calculates an approximate straight line that approximates the data of the characteristic curve at the time of discharge for the low charge rate range where the charge rate is lower than the range of the characteristic curve. The approximate function creation device 20 further generates an approximate function by further using the approximate straight lines of the calculated high charge rate range and low charge rate range. The approximate function creation device 20 intersects the approximate straight line of the high charge rate range with the approximate straight line of the divided range on the highest charge rate range side, and the approximate straight line of the low charge rate range and the closest charge rate range side. Approximate functions are generated with the intersection of the divided range with the approximate line as the switching point of the approximate line. Thereby, the approximate function creating apparatus 20 according to the present embodiment can accurately generate an approximate function including the high charge rate range and the low charge rate range.

さて、これまで開示の装置に関する実施例について説明したが、本発明は上述した実施例以外にも、種々の異なる形態にて実施されてよいものである。そこで、以下では、本発明に含まれる他の実施例を説明する。   Although the embodiments related to the disclosed apparatus have been described above, the present invention may be implemented in various different forms other than the above-described embodiments. Therefore, another embodiment included in the present invention will be described below.

例えば、上記の実施例では、二次電池としてリチウムイオン電池を用いた場合について説明したが、開示の装置はこれに限定されない。二次電池は、リチウムイオンポリマー二次電池、カルシウムイオン二次電池、ナノワイヤバッテリ、ニッケル水素二次電池、鉛蓄電池であってもよい。また、リチウムイオン電池としては、コバルト酸リチウムイオン電池、リン酸鉄リチウムイオン電池等にも適用できる。これにより、多種の電池の充電率を測定できる。   For example, in the above-described embodiment, the case where a lithium ion battery is used as the secondary battery has been described, but the disclosed apparatus is not limited thereto. The secondary battery may be a lithium ion polymer secondary battery, a calcium ion secondary battery, a nanowire battery, a nickel metal hydride secondary battery, or a lead storage battery. Moreover, as a lithium ion battery, it is applicable also to a lithium cobalt ion battery, a lithium iron phosphate battery, etc. Thereby, the charge rate of various batteries can be measured.

[分散および統合]
また、図示した各装置の各構成要素は、必ずしも物理的に図示の如く構成されていることを要しない。すなわち、各装置の分散・統合の具体的形態は図示のものに限られず、その全部または一部を、各種の負荷や使用状況などに応じて、任意の単位で機能的または物理的に分散・統合して構成することができる。例えば、上記実施例の取得部50、第1算出部51、第2算出部52、第3算出部53、生成部54および出力部55を近似関数作成装置20の外部装置としてネットワーク経由で接続するようにしてもよい。また、上記実施例の取得部50、第1算出部51、第2算出部52、第3算出部53、生成部54および出力部55を別の装置がそれぞれ有し、ネットワーク接続されて協働することで、上記の近似関数作成装置20の機能を実現するようにしてもよい。
[Distribution and integration]
In addition, each component of each illustrated apparatus does not necessarily need to be physically configured as illustrated. In other words, the specific form of distribution / integration of each device is not limited to that shown in the figure, and all or a part thereof may be functionally or physically distributed or arbitrarily distributed in arbitrary units according to various loads or usage conditions. Can be integrated and configured. For example, the acquisition unit 50, the first calculation unit 51, the second calculation unit 52, the third calculation unit 53, the generation unit 54, and the output unit 55 of the above embodiment are connected as an external device of the approximate function creation device 20 via a network. You may do it. Further, the acquisition unit 50, the first calculation unit 51, the second calculation unit 52, the third calculation unit 53, the generation unit 54, and the output unit 55 of the above-described embodiment are respectively provided by different devices, and are connected via a network to cooperate. Thus, the function of the above approximate function creation device 20 may be realized.

[近似関数作成プログラム]
また、上記の実施例で説明した各種の処理は、予め用意されたプログラムをパーソナルコンピュータやワークステーションなどのコンピュータで実行することによって実現することができる。そこで、以下では、図11を用いて、上記の実施例と同様の機能を有する近似関数作成プログラムを実行するコンピュータの一例について説明する。図11は、近似関数作成プログラムを実行するコンピュータの一例を示す図である。
[Approximate function creation program]
The various processes described in the above embodiments can be realized by executing a prepared program on a computer such as a personal computer or a workstation. In the following, an example of a computer that executes an approximate function creation program having the same function as that of the above-described embodiment will be described with reference to FIG. FIG. 11 is a diagram illustrating an example of a computer that executes an approximate function creation program.

図11に示すように、コンピュータ300は、CPU(Central Processing Unit)310、ROM(Read Only Memory)320、HDD(Hard Disk Drive)330、RAM(Random Access Memory)340を有する。これら310〜340の各部は、バス400を介して接続される。   As shown in FIG. 11, the computer 300 includes a central processing unit (CPU) 310, a read only memory (ROM) 320, a hard disk drive (HDD) 330, and a random access memory (RAM) 340. These units 310 to 340 are connected via a bus 400.

ROM320には上記実施例の各処理部と同様の機能を発揮する近似関数作成プログラム320aが予め記憶される。例えば、上記実施例の取得部50、第1算出部51、第2算出部52、第3算出部53、生成部54および出力部55と同様の機能を発揮する近似関数作成プログラム320aを記憶させる。なお、近似関数作成プログラム320aについては、適宜分離しても良い。   The ROM 320 stores in advance an approximate function creation program 320a that exhibits the same function as each processing unit of the above embodiment. For example, the approximate function creation program 320a that performs the same functions as the acquisition unit 50, the first calculation unit 51, the second calculation unit 52, the third calculation unit 53, the generation unit 54, and the output unit 55 of the above embodiment is stored. . Note that the approximate function creation program 320a may be separated as appropriate.

そして、CPU310が、近似関数作成プログラム320aをROM320から読み出して実行することで、上記実施例と同様の動作を実行する。すなわち、近似関数作成プログラム320aは、上記実施例の取得部50、第1算出部51、第2算出部52、第3算出部53、生成部54および出力部55と同様の動作を実行する。   Then, the CPU 310 reads the approximate function creation program 320a from the ROM 320 and executes it, thereby executing the same operation as in the above embodiment. That is, the approximate function creation program 320a performs the same operations as the acquisition unit 50, the first calculation unit 51, the second calculation unit 52, the third calculation unit 53, the generation unit 54, and the output unit 55 of the above embodiment.

なお、上記した近似関数作成プログラム320aについては、必ずしも最初からROM320に記憶させることを要しない。近似関数作成プログラム320aはHDD330に記憶させてもよい。   The approximate function creation program 320a described above does not necessarily need to be stored in the ROM 320 from the beginning. The approximate function creation program 320 a may be stored in the HDD 330.

また、ROM320には、充電率を推定する充電率推定プログラムが記憶されてもよい。この場合、上記実施例の推定部12Bと同様の機能を発揮する充電率推定プログラムを記憶させる。なお、充電率推定プログラムについても、適宜分離しても良い。   The ROM 320 may store a charging rate estimation program for estimating the charging rate. In this case, a charging rate estimation program that exhibits the same function as the estimation unit 12B of the above embodiment is stored. Note that the charging rate estimation program may be appropriately separated.

この場合、CPU310が、充電率推定プログラムをROM320から読み出して実行することで、上記実施例と同様の動作を実行する。すなわち、充電率推定プログラムは、推定部12Bと同様の動作を実行する。   In this case, the CPU 310 reads the charge rate estimation program from the ROM 320 and executes it, thereby executing the same operation as in the above embodiment. That is, the charging rate estimation program performs the same operation as that of the estimation unit 12B.

なお、上記した充電率推定プログラムについても、必ずしも最初からROM320に記憶させることを要しない。充電率推定プログラムはHDD330に記憶させてもよい。   Note that the above-described charging rate estimation program is not necessarily stored in the ROM 320 from the beginning. The charging rate estimation program may be stored in the HDD 330.

また、例えば、コンピュータ300に挿入されるフレキシブルディスク(FD)、CD−ROM、DVDディスク、光磁気ディスク、ICカードなどの「可搬用の物理媒体」にプログラムを記憶させておく。そして、コンピュータ300がこれらからプログラムを読み出して実行するようにしてもよい。   Further, for example, the program is stored in a “portable physical medium” such as a flexible disk (FD), a CD-ROM, a DVD disk, a magneto-optical disk, an IC card or the like inserted into the computer 300. Then, the computer 300 may read and execute the program from these.

また、例えば、コンピュータ300に挿入されるフレキシブルディスク(FD)、CD−ROM、DVDディスク、光磁気ディスク、ICカードなどの「可搬用の物理媒体」にプログラムを記憶させておく。そして、コンピュータ300がこれらからプログラムを読み出して実行するようにしてもよい。   Further, for example, the program is stored in a “portable physical medium” such as a flexible disk (FD), a CD-ROM, a DVD disk, a magneto-optical disk, an IC card or the like inserted into the computer 300. Then, the computer 300 may read and execute the program from these.

さらには、公衆回線、インターネット、LAN、WANなどを介してコンピュータ300に接続される「他のコンピュータ(またはサーバ)」などにプログラムを記憶させておく。そして、コンピュータ300がこれらからプログラムを読み出して実行するようにしてもよい。   Furthermore, the program is stored in “another computer (or server)” connected to the computer 300 via a public line, the Internet, a LAN, a WAN, or the like. Then, the computer 300 may read and execute the program from these.

10 リチウムイオン電池
20 近似関数作成装置
33 記憶部
34 制御部
40 測定データ
41 電圧変化データ
42 特性データ
50 取得部
51 第1算出部
52 第2算出部
53 第3算出部
54 生成部
55 出力部
DESCRIPTION OF SYMBOLS 10 Lithium ion battery 20 Approximate function creation apparatus 33 Memory | storage part 34 Control part 40 Measurement data 41 Voltage change data 42 Characteristic data 50 Acquisition part 51 1st calculation part 52 2nd calculation part 53 3rd calculation part 54 Generation part 55 Output part

Claims (7)

コンピュータに
二次電池の開回路電圧と当該二次電池の充電率との相関関係を示した特性曲線のデータを取得し、
取得された特性曲線のデータにより示される特性曲線の充電率に対する開回路電圧の変化率が所定の閾値以下となる範囲内で、中間点を移動させつつ中間点の位置毎に、中間点を境として前記範囲を分割した各分割範囲の特性曲線のデータをそれぞれ近似した近似直線と、当該近似直線の近似に用いた特性曲線のデータとの類似度を算出し、
算出される類似度が最も高い中間点を境とした各分割範囲の近似直線を用いて特性曲線を近似した近似関数を生成する
処理を実行させることを特徴とする近似関数作成プログラム。
The computer obtains the characteristic curve data showing the correlation between the open circuit voltage of the secondary battery and the charging rate of the secondary battery,
Within the range where the change rate of the open circuit voltage with respect to the charging rate of the characteristic curve indicated by the acquired characteristic curve data is below a predetermined threshold, the intermediate point is delimited at each intermediate point position while moving the intermediate point. And calculating the similarity between the approximate line obtained by approximating the characteristic curve data of each divided range obtained by dividing the range and the characteristic curve data used for approximation of the approximate line,
An approximate function creation program that executes processing for generating an approximate function that approximates a characteristic curve using an approximate straight line of each divided range with an intermediate point having the highest similarity as a boundary.
前記生成する処理は、分割範囲が隣り合う近似直線の交点を近似直線の切替点とした近似関数を生成する
ことを特徴とする請求項1に記載の近似関数作成プログラム。
The approximation function creating program according to claim 1, wherein the generating process generates an approximate function using an intersection of approximate lines with adjacent division ranges as a switching point of the approximate line.
前記特性曲線のデータは、二次電池の充電時の特性曲線のデータと放電時の特性曲線のデータを含み、
コンピュータに、
二次電池の充電時の特性曲線と放電時の特性曲線とのヒステリシスによる誤差を算出する処理をさらに実行させ、
前記類似度を算出する処理は、充電時の特性曲線のデータおよび放電時の特性曲線のデータに対して、前記誤差に応じた補正を行った後に平均化したデータを用いて前記類似度を算出し、
前記生成する処理は、前記近似関数に前記誤差に応じた変更を行って放電時の特性曲線を近似した近似関数と充電時の特性曲線を近似した近似関数とを生成する
ことを特徴とする請求項1または2に記載の近似関数作成プログラム。
The data of the characteristic curve includes data of a characteristic curve at the time of charging the secondary battery and data of a characteristic curve at the time of discharging,
On the computer,
Further executing a process of calculating an error due to hysteresis between the characteristic curve at the time of charging the secondary battery and the characteristic curve at the time of discharging,
In the process of calculating the similarity, the similarity is calculated using data obtained by correcting the characteristic curve data during charging and the characteristic curve data during discharging after performing correction according to the error. And
The generating process includes generating an approximate function approximating a characteristic curve at the time of discharging and an approximate function approximating a characteristic curve at the time of charging by changing the approximate function according to the error. Item 3. An approximation function creation program according to item 1 or 2.
コンピュータに、
前記特性曲線の前記範囲よりも充電率が高い高充電率範囲については充電時の特性曲線のデータを近似した近似直線を算出し、前記範囲よりも充電率が低い低充電率範囲については放電時の特性曲線のデータを近似した近似直線を算出する処理をさらに実行させ、
前記生成する処理は、算出された高充電率範囲および低充電率範囲の近似直線をさらに用いて、高充電率範囲の近似直線と、最も高充電率範囲側の分割範囲の近似直線との交点、および、低充電率範囲の近似直線と、最も低充電率範囲側の分割範囲の近似直線との交点をそれぞれ近似直線の切替点とした近似関数を生成する
ことを特徴とする請求項3に記載の近似関数作成プログラム。
On the computer,
For the high charging rate range where the charging rate is higher than the range of the characteristic curve, calculate an approximate straight line approximating the characteristic curve data during charging, and for the low charging rate range where the charging rate is lower than the range during discharging To further execute a process of calculating an approximate straight line that approximates the characteristic curve data of
The processing to be generated further uses the approximate straight lines of the calculated high charge rate range and low charge rate range, and the intersection of the approximate straight line of the high charge rate range and the approximate straight line of the divided range on the highest charge rate range side And generating an approximate function with the intersection of the approximate straight line in the low charge rate range and the approximate straight line in the divided range on the lowest charge rate range side as the switching point of the approximate straight line, respectively. The approximate function creation program described.
コンピュータが
二次電池の開回路電圧と当該二次電池の充電率との相関関係を示した特性曲線のデータを取得し、
取得された特性曲線のデータにより示される特性曲線の充電率に対する開回路電圧の変化率が所定の閾値以下となる範囲内で、中間点を移動させつつ中間点の位置毎に、中間点を境として前記範囲を分割した各分割範囲の特性曲線のデータをそれぞれ近似した近似直線と、当該近似直線の近似に用いた特性曲線のデータとの類似度を算出し、
算出される類似度が最も高い中間点を境とした各分割範囲の近似直線を用いて特性曲線を近似した近似関数を生成する
処理を実行することを特徴とする近似関数作成方法。
The computer obtains characteristic curve data showing the correlation between the open circuit voltage of the secondary battery and the charging rate of the secondary battery,
Within the range where the change rate of the open circuit voltage with respect to the charging rate of the characteristic curve indicated by the acquired characteristic curve data is below a predetermined threshold, the intermediate point is delimited at each intermediate point position while moving the intermediate point. As an approximation straight line approximating the characteristic curve data of each divided range obtained by dividing the range, and the similarity between the characteristic curve data used for approximating the approximate straight line is calculated,
An approximation function creating method characterized by executing a process of generating an approximation function approximating a characteristic curve using an approximation straight line in each divided range with a calculated intermediate point as the boundary.
二次電池の開回路電圧と当該二次電池の充電率との相関関係を示した特性曲線のデータを取得する取得部と、
前記取得部により取得された特性曲線のデータにより示される特性曲線の充電率に対する開回路電圧の変化率が所定の閾値以下となる範囲内で、中間点を移動させつつ中間点の位置毎に、中間点を境として前記範囲を分割した各分割範囲の特性曲線のデータをそれぞれ近似した近似直線と、当該近似直線の近似に用いた特性曲線のデータとの類似度を算出する算出部と、
前記算出部により算出される類似度が最も高い中間点を境とした各分割範囲の近似直線を用いて特性曲線を近似した近似関数を生成する生成部と、
を有することを特徴とする近似関数作成装置。
An acquisition unit for acquiring data of a characteristic curve indicating a correlation between an open circuit voltage of the secondary battery and a charging rate of the secondary battery;
Within the range where the change rate of the open circuit voltage with respect to the charging rate of the characteristic curve indicated by the characteristic curve data acquired by the acquisition unit is equal to or less than a predetermined threshold, for each position of the intermediate point while moving the intermediate point, A calculation unit that calculates the similarity between the approximate line obtained by approximating the characteristic curve data of each divided range obtained by dividing the range with the intermediate point as a boundary, and the characteristic curve data used for approximation of the approximate line;
A generation unit that generates an approximate function that approximates a characteristic curve using an approximate straight line of each divided range with a middle point having the highest similarity calculated by the calculation unit;
An approximate function creation device characterized by comprising:
コンピュータに
二次電池の開回路電圧と当該二次電池の充電率との相関関係を示した特性曲線のデータにより示される特性曲線の充電率に対する開回路電圧の変化率が所定の閾値以下となる範囲内で、中間点を移動させつつ中間点の位置毎に、中間点を境として前記範囲を分割した各分割範囲の特性曲線のデータをそれぞれ近似した近似直線と、当該近似直線の近似に用いた特性曲線のデータとの類似度を算出し、算出される類似度が最も高い中間点を境とした各分割範囲の近似直線を用いて生成された特性曲線を近似した近似関数の情報を取得し、
取得した近似関数の情報を用いて前記二次電池の充電率を推定する
処理を実行させることを特徴とする充電率推定プログラム。
The rate of change of the open circuit voltage with respect to the charging rate of the characteristic curve indicated by the characteristic curve data indicating the correlation between the open circuit voltage of the secondary battery and the charging rate of the secondary battery is less than a predetermined threshold. For each position of the intermediate point while moving the intermediate point within the range, an approximate straight line approximating the characteristic curve data of each divided range obtained by dividing the range with the intermediate point as a boundary, and for approximating the approximate straight line The degree of similarity with the characteristic curve data was calculated, and information on the approximation function that approximated the characteristic curve generated using the approximate straight line of each divided range with the highest calculated similarity as the boundary is obtained. And
A charge rate estimation program for executing a process of estimating a charge rate of the secondary battery using information of the obtained approximate function.
JP2013247935A 2013-11-29 2013-11-29 Approximation function creation program, approximation function creation method, approximation function creation device, and charging rate estimation program Active JP6171897B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2013247935A JP6171897B2 (en) 2013-11-29 2013-11-29 Approximation function creation program, approximation function creation method, approximation function creation device, and charging rate estimation program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2013247935A JP6171897B2 (en) 2013-11-29 2013-11-29 Approximation function creation program, approximation function creation method, approximation function creation device, and charging rate estimation program

Publications (2)

Publication Number Publication Date
JP2015105875A JP2015105875A (en) 2015-06-08
JP6171897B2 true JP6171897B2 (en) 2017-08-02

Family

ID=53436060

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2013247935A Active JP6171897B2 (en) 2013-11-29 2013-11-29 Approximation function creation program, approximation function creation method, approximation function creation device, and charging rate estimation program

Country Status (1)

Country Link
JP (1) JP6171897B2 (en)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6561620B2 (en) * 2015-06-23 2019-08-21 日本電気株式会社 Measuring apparatus, measuring system, measuring result collecting method and program
JP6499939B2 (en) * 2015-07-09 2019-04-10 川崎重工業株式会社 Secondary battery charge state estimation method
WO2018025350A1 (en) * 2016-08-03 2018-02-08 富士通株式会社 Estimation device, estimation program, and charging control device
KR101972521B1 (en) * 2016-08-26 2019-04-25 주식회사 엘지화학 Apparatus and method for performance testing of a battery cell
JP6725380B2 (en) * 2016-09-20 2020-07-15 株式会社東芝 Optimization calculation device, optimization calculation method, and program
CN106772064B (en) * 2016-11-25 2019-05-28 广东电网有限责任公司电力科学研究院 A kind of health state of lithium ion battery prediction technique and device
JP7116886B2 (en) * 2017-02-20 2022-08-12 株式会社Gsユアサ state estimator
CN110991527B (en) * 2019-12-02 2023-05-05 国网湖北省电力有限公司电力科学研究院 Similarity threshold determining method considering average fluctuation rate of voltage curve

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3253207B2 (en) * 1993-01-25 2002-02-04 松下電器産業株式会社 Temperature compensated crystal oscillator
JPH08221388A (en) * 1995-02-09 1996-08-30 Nec Corp Fitting parameter decision method
JPH09243431A (en) * 1996-03-07 1997-09-19 Matsushita Electric Ind Co Ltd Flow rate sensor and flowmeter utilizing the same
JPH1132237A (en) * 1997-07-10 1999-02-02 Sony Corp Gamma correction device
KR101863036B1 (en) * 2011-11-30 2018-06-01 주식회사 실리콘웍스 Method for estimating the state of charge of battery and battery management system
US9201119B2 (en) * 2011-12-19 2015-12-01 Qualcomm Incorporated Battery fuel gauge
JP6065561B2 (en) * 2012-03-08 2017-01-25 日産自動車株式会社 Secondary battery control device and SOC detection method

Also Published As

Publication number Publication date
JP2015105875A (en) 2015-06-08

Similar Documents

Publication Publication Date Title
JP6171897B2 (en) Approximation function creation program, approximation function creation method, approximation function creation device, and charging rate estimation program
KR102452548B1 (en) Apparatus for determination battery degradation, system having the same and method thereof
US9952289B2 (en) Storage battery evaluating apparatus and method
JP6221884B2 (en) Estimation program, estimation method, and estimation apparatus
JP6299187B2 (en) Estimation program, estimation method, and estimation apparatus
US10564223B2 (en) Power storage system, control system of secondary battery, and control method of secondary battery
JP6330605B2 (en) Estimation program, estimation method, and estimation apparatus
JP2015038444A (en) Secondary battery remaining capacity estimation method and secondary battery remaining capacity estimation apparatus
JP2015184217A (en) Estimation program, estimation method, and estimation device
KR20210028476A (en) Method and apparatus charging battery
CN111707955B (en) Method, apparatus and medium for estimating remaining life of battery
JP2015215169A (en) Parameter calculation device, parameter calculation method, and program
WO2018025350A1 (en) Estimation device, estimation program, and charging control device
JP2017223536A (en) Device for estimating battery state and method for estimating battery state
WO2018211824A1 (en) Battery control device and vehicle system
JP6421411B2 (en) Estimation program for estimating battery charging rate, estimation method for estimating battery charging rate, and estimation device for estimating battery charging rate
CN113484763A (en) Method, device, equipment and storage medium for determining residual electric quantity of battery
JP2018194357A (en) Charge rate estimating device
JP2018009939A (en) Simulation method and simulation device
JP6287125B2 (en) Function creation program, function creation method, function creation device, and charging rate estimation program
CN106970328B (en) SOC estimation method and device
WO2018025306A1 (en) Estimation device, estimation program, and charging control device
CN114200315B (en) Method and device for predicting remaining charging time, electronic equipment and storage medium
CN115792665A (en) Method and device for determining residual capacity of battery, electronic equipment and storage medium
WO2018029849A1 (en) Estimation device, estimation program, and charging control device

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20160804

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20170526

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20170606

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20170619

R150 Certificate of patent or registration of utility model

Ref document number: 6171897

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150