"充放電を行う蓄電素子と、前記蓄電素子の外部における温度を取得するための温度センサと、前記蓄電素子の外部における温度と、熱の移動を表す式とを用いて、前記蓄電素子の内部抵抗に対応した温度を示す基準点の温度を算出するコントローラと、を有し、前記コントローラは、前記蓄電素子の充放電を再開するとき、前記蓄電素子の内部における温度バラツキが解消しているか否かを判別し、前記温度バラツキが解消しているときには、前記蓄電素子の表面における表面温度を前記基準点の温度として用い、前記温度バラツキが解消していないときには、前記基準点の温度を算出することを特徴とする蓄電システム。前記コントローラは、前記温度バラツキが解消していないとき、前記蓄電素子の外面における温度および熱伝導方程式を用いて、前記基準点の温度を算出することを特徴とする請求項1に記載の蓄電システム。前記コントローラは、前記蓄電素子の充放電を行っていない間の経過時間を計測し、前記経過時間が所定時間よりも長いときに、前記温度バラツキが解消していると判別することを特徴とする請求項1又は2に記載の蓄電システム。前記コントローラは、前記蓄電素子の充放電を行っていない間において、所定周期で前記表面温度を取得し、前記所定周期における前記表面温度の変化量が閾値よりも小さいときに、前記温度バラツキが解消していると判別することを特徴とする請求項1又は2に記載の蓄電システム。前記コントローラは、前記蓄電素子の充放電を停止したときから、前記蓄電素子の充放電を再開させるときまでの間における温度変化を用いて、前記基準点の温度を算出することを特徴とする請求項1から4のいずれか1つに記載の蓄電システム。前記コントローラは、前記基準点の温度を用いて、前記蓄電素子の内部状態を推定することを特徴とする請求項1から5のいずれか1つに記載の蓄電システム。前記内部状態は、SOHであることを特徴とする請求項6に記載の蓄電システム。前記蓄電素子は、車両の走行に用いられるエネルギを出力することを特徴とする請求項1から7のいずれか1つに記載の蓄電システム。充放電を行う蓄電素子の外部における温度を取得するステップと、前記蓄電素子の充放電を再開するとき、前記蓄電素子の内部における温度バラツキが解消しているか否かを判別するステップと、前記温度バラツキが解消しているとき、前記蓄電素子の内部抵抗に対応した温度を示す基準点の温度として、前記蓄電素子の表面における表面温度を用い、前記温度バラツキが解消していないとき、前記蓄電素子の外部における温度と、熱の移動を表す式とを用いて、前記基準点の温度を算出するステップと、ことを特徴とする蓄電素子の温度推定方法。前記温度バラツキが解消していないとき、前記蓄電素子の外面における温度および熱伝導方程式を用いて、前記基準点の温度を算出することを特徴とする請求項9に記載の蓄電素子の温度推定方法。前記蓄電素子の充放電を行っていない間の経過時間を計測し、前記経過時間が所定時間よりも長いときに、前記温度バラツキが解消していると判別することを特徴とする請求項9又は10に記載の蓄電素子の温度推定方法。前記蓄電素子の充放電を行っていない間において、所定周期で前記表面温度を取得し、前記所定周期における前記表面温度の変化量が閾値よりも小さいときに、前記温度バラツキが解消していると判別することを特徴とする請求項9又は10に記載の蓄電素子の温度推定方法。前記蓄電素子の充放電を停止したときから、前記蓄電素子の充放電を再開させるときまでの間における温度変化を用いて、前記基準点の温度を算出することを特徴とする請求項9から12のいずれか1つに記載の蓄電素子の温度推定方法。前記基準点の温度は、前記蓄電素子の内部状態を推定するために用いられることを特徴とする請求項9から13のいずれか1つに記載の蓄電素子の温度推定方法。前記内部状態は、SOHであることを特徴とする請求項14に記載の蓄電素子の温度推定方法。単電池の内部では、放熱性等の理由により、温度分布にバラツキが生じている。一般的には、単電池の中心部における温度は、単電池の外面における温度よりも高くなりやすい。このような温度分布を有する単電池において、単電池の外面に取り付けた温度センサの出力だけでは、単電池の内部における温度を取得することができない。本願第1の発明である蓄電システムは、充放電を行う蓄電素子と、蓄電素子の外部における温度を取得するための温度センサと、蓄電素子の外部における温度と、熱の移動を表す式とを用いて、蓄電素子の内部抵抗に対応した温度を示す基準点の温度を算出するコントローラと、を有する。コントローラは、蓄電素子の充放電を再開するとき、蓄電素子の内部における温度バラツキが解消しているか否かを判別する。コントローラは、温度バラツキが解消しているときには、蓄電素子の表面における表面温度を基準点の温度として用い、温度バラツキが解消していないときには、基準点の温度を算出する。蓄電素子の充放電を再開させるとき、充放電が停止しているときの蓄電素子の状態によっては、蓄電素子の内部における温度バラツキが解消していることがある。そこで、本願第1の発明では、温度バラツキが解消しているか否かに応じて、基準点の温度を算出する方法を変更することにより、基準点の温度の推定精度を向上させることができる。温度バラツキが解消しているときには、基準点の温度として、表面温度を用いればよい。温度バラツキが解消していないときには、蓄電素子の外部における温度と、熱の移動を表す式とを用いて、基準点の温度を算出することができる。これにより、内部において温度バラツキが発生しやすい蓄電素子において、蓄電素子の内部抵抗に対応した温度を推定することができる。基準点の温度を算出するとき、蓄電素子の外面における温度および熱伝導方程式を用いることができる。温度バラツキが解消しているか否かを判別する第1の方法として、まず、蓄電素子の充放電を行っていない間の経過時間を計測する。そして、経過時間が所定時間よりも長いときには、温度バラツキが解消していると判別することができる。所定時間は、温度バラツキを解消させる時間として予め設定することができる。温度バラツキが解消しているか否かを判別する第2の方法として、まず、蓄電素子の充放電を行っていない間において、所定周期で表面温度を取得する。そして、所定周期における表面温度の変化量が閾値よりも小さいときには、温度バラツキが解消していると判別することができる。温度バラツキが解消し始めると、表面温度の変化量が減少し始める。したがって、表面温度の変化量を閾値と比較することにより、温度バラツキが解消しているか否かを判別することができる。閾値は、温度バラツキの解消を許容する範囲内において、適宜設定することができる。基準点の温度を算出するときには、蓄電素子の充放電を停止したときから、蓄電素子の充放電を再開させるときまでの間における温度変化を考慮することができる。具体的には、温度センサによる検出温度の経時変化に基づいて、基準点の温度を算出することができる。基準点の温度は、蓄電素子の内部状態を推定するために用いることができる。内部状態としては、SOCStateOfCharge又はSOHStateOfHealthがある。SOCは、蓄電素子の満充電容量Cfに対する、現在の充電容量Ccの割合Cc/Cfを示す値である。SOHは、初期状態の満充電容量Cfiniと劣化後の満充電容量Cfcとの比率Cfc/Cfiniであり、蓄電素子の劣化状態を示す指標となる。蓄電素子から出力されたエネルギは、車両を走行させるために用いることができる。蓄電素子を車両に搭載するときには、複数の蓄電素子を直列に接続することができる。本願第2の発明である蓄電素子の温度推定方法は、充放電を行う蓄電素子の外部における温度を取得するステップと、蓄電素子の充放電を再開するとき、蓄電素子の内部における温度バラツキが解消しているか否かを判別するステップと、を有する。温度バラツキが解消しているときには、蓄電素子の内部抵抗に対応した温度を示す基準点の温度として、蓄電素子の表面における表面温度を用いる。温度バラツキが解消していないときには、蓄電素子の外部における温度と、熱の移動を表す式とを用いて、基準点の温度を算出する。本願第2の発明においても、本願第1の発明と同様の効果を得ることができる。以下、本発明の実施例について説明する。本発明の実施例1である電池システム蓄電システムに相当するについて、図1を用いて説明する。図1は、電池システムの構成を示す図である。本実施例の電池システムは、車両に搭載することができる。車両としては、ハイブリッド自動車や電気自動車がある。ハイブリッド自動車は、車両を走行させる動力源として、後述する組電池に加えて、エンジン又は燃料電池を備えている。電気自動車は、車両を走行させる動力源として、後述する組電池だけを備えている。本実施例の電池システムは、組電池10を有する。組電池10は、直列に接続された複数の単電池11を有する。単電池蓄電素子に相当する11としては、ニッケル水素電池やリチウムイオン電池といった二次電池を用いることができる。また、二次電池の代わりに、電気二重層キャパシタコンデンサを用いることができる。単電池11の数は、組電池10の要求出力などに基づいて、適宜設定することができる。本実施例では、組電池10を構成する、すべての単電池11が直列に接続されているが、並列に接続された複数の単電池11が、組電池10に含まれていてもよい。電圧センサ21は、組電池10の端子間電圧総電圧を検出し、検出結果をコントローラ30に出力する。ここで、電圧センサ21を用いて、各単電池11の電圧を検出したり、少なくとも2つの単電池11を含む電池ブロックの電圧を検出したりすることができる。電流センサ22は、組電池10に流れる電流値を検出し、検出結果をコントローラ30に出力する。温度センサ23は、組電池10単電池11の温度を検出し、検出結果をコントローラ30に出力する。温度センサ23の数は、適宜設定することができる。複数の温度センサ23を用いるときには、互いに異なる位置に配置された単電池11に対して、温度センサ23を配置することができる。コントローラ30は、メモリ30aを有しており、メモリ30aは、コントローラ30が所定処理を行うための各種の情報を格納している。本実施例では、メモリ30aが、コントローラ30に内蔵されているが、コントローラ30の外部にメモリ30aを設けることもできる。組電池10の正極ラインPLには、システムメインリレーSMR-Bが設けられている。システムメインリレーSMR-Bは、コントローラ30からの制御信号を受けることにより、オンおよびオフの間で切り替わる。組電池10の負極ラインNLには、システムメインリレーSMR-Gが設けられている。システムメインリレーSMR-Gは、コントローラ30からの制御信号を受けることにより、オンおよびオフの間で切り替わる。システムメインリレーSMR-Gに対しては、システムメインリレーSMR-Pおよび電流制限抵抗Rが並列に接続されている。システムメインリレーSMR-Pおよび電流制限抵抗Rは、直列に接続されている。システムメインリレーSMR-Pは、コントローラ30からの制御信号を受けることにより、オンおよびオフの間で切り替わる。電流制限抵抗Rは、組電池10を負荷後述するインバータ31と接続するときに、突入電流が流れるのを抑制するために用いられる。イグニッションスイッチがオフからオンに切り替わったときに、組電池10は、負荷と接続される。組電池10を負荷と接続するとき、コントローラ30は、システムメインリレーSMR-B,SMR-Pをオフからオンに切り替える。これにより、電流制限抵抗Rに電流が流れることになる。次に、コントローラ30は、システムメインリレーSMR-Gをオフからオンに切り替えた後に、システムメインリレーSMR-Pをオンからオフに切り替える。これにより、組電池10および負荷の接続が完了する。一方、イグニッションスイッチがオンからオフに切り替わったとき、組電池10および負荷の接続が遮断される。組電池10および負荷の接続を遮断するとき、コントローラ30は、システムメインリレーSMR-B,SMR-Gをオンからオフに切り替える。インバータ31は、組電池10からの直流電力を交流電力に変換し、交流電力をモータジェネレータ32に出力する。モータジェネレータ32としては、例えば、三相交流モータを用いることができる。モータジェネレータ32は、インバータ31からの交流電力を受けて、車両を走行させるための運動エネルギを生成する。モータジェネレータ32によって生成された運動エネルギは、車輪に伝達される。車両を減速させたり、停止させたりするとき、モータジェネレータ32は、車両の制動時に発生する運動エネルギを電気エネルギ交流電力に変換する。インバータ31は、モータジェネレータ32が生成した交流電力を直流電力に変換し、直流電力を組電池10に出力する。これにより、組電池10は、回生電力を蓄えることができる。本実施例では、組電池10をインバータ31に接続しているが、これに限るものではない。具体的には、組電池10およびインバータ31の間の電流経路に、昇圧回路を配置することができる。昇圧回路を用いることにより、組電池10の出力電圧を昇圧することができる。また、昇圧回路は、インバータ31から組電池10への出力電圧を降圧することができる。次に、単電池11の構造について、図2を用いて説明する。図2は、単電池11の構造を示す概略図である。図1において、X軸、Y軸およびZ軸は、互いに直交する軸であり、X軸、Y軸およびZ軸の関係は、他の図面においても同様である。単電池11は、発電要素111と、発電要素111を収容する電池ケース112とを有する。発電要素111は、充放電を行う要素であり、図3に示すように、正極板111aと、負極板111bと、正極板111aおよび負極板111bの間に配置されるセパレータ111cとを有する。セパレータ111cは、電解液を含んでいる。発電要素111は、正極板111a、セパレータ111cおよび負極板111bを重ねた積層体図3に示す構成を、Y軸図2参照の周りで捲くことによって構成されている。本実施例では、正極板111a、セパレータ111cおよび負極板111bを重ねた積層体を捲くことにより、発電要素111を構成しているが、これに限るものではない。例えば、正極板111a、セパレータ111cおよび負極板111bを重ねただけで、発電要素111を構成することができる。また、本実施例では、電解液を含むセパレータ111cを用いているが、セパレータ111cの代わりに、固体電解質を用いることもできる。固体電解質としては、高分子固体電解質や無機固体電解質を用いることができる。正極板111aは、図4に示すように、集電板111a1の表面に正極活物質層111a2を形成したものであり、正極活物質層111a2は、正極活物質や導電剤等を含んでいる。負極板111bは、集電板111b1の表面に負極活物質層111b2を形成したものであり、負極活物質層111b2は、負極活物質や導電剤等を含んでいる。単電池11として、リチウムイオン二次電池を用いるとき、例えば、正極板111aの集電板111a1をアルミニウムで形成し、負極板111bの集電板111b1を銅で形成することができる。また、正極活物質として、LiO2を用い、負極活物質として、カーボンを用いることができる。なお、正極板111aおよび負極板111bは、図4に示す構成に限るものではない。例えば、集電板の一方の面に正極活物質層を形成し、集電板の他方の面に負極活物質層を形成した電極板いわゆるバイポーラ電極を用いることができる。電池ケース112は、例えば、金属で形成することができる。図2に示すように、電池ケース112の上面には、正極端子113および負極端子114が設けられている。正極端子113は、発電要素111の正極板111aと電気的に接続されており、負極端子114は、発電要素111の負極板111bと電気的に接続されている。図2に示す構成では、電池ケース112の上面に、温度センサ23が設けられている。温度センサ23は単電池11電池ケース112の外面に取り付けられているため、温度センサ23による検出温度は、単電池11の外面における温度となる。温度センサ23としては、例えば、熱電対を用いることができる。また、電池ケース112に対する温度センサ23の取り付け位置は、適宜設定することができる。ここで、複数の単電池11をX方向に並べて配置する場合には、電池ケース112の上面に温度センサ23を配置することができる。次に、単電池11の内部における温度特性について、図5を用いて説明する。図5では、縦軸を温度とし、横軸を単電池11の厚さとした座標系と、単電池11の内部構造とを重ねて示している。単電池11の厚さとは、X方向における単電池11の長さである。図5に示す横軸の方向は、正極板111a、セパレータ111cおよび負極板111bが重なっている方向である。図5には、単電池11の内部における温度分布一例を示している。図5に示す中心点Oは、単電池11の厚さ方向における発電要素111の中心に相当する位置を示している。単電池11発電要素111は、充放電によって発熱するが、放熱性等によって、図5に示す温度分布が生じる。発電要素111の外面は、最も熱が逃げやすくなっており、最も温度が低くなりやすい。発電要素111の外面は、電池ケース112の内壁面と隣り合っている。一方、中心点Oに近づくにつれて、熱が逃げにくくなり、温度が高くなりやすい。図5に示すように、単電池11の内部では、単電池11の厚さ方向の位置に応じて温度が異なる。本実施例では、発電要素111の内部抵抗に対応した温度を、単電池11の温度以下、性能温度というとして用い、以下に説明するように、単電池11の性能温度を推定している。本実施例では、単電池11の性能温度を推定するために、下記式1に示す熱伝導方程式を用いる。式1において、Tは温度、tは時間、λは熱伝導率、ρは密度、cは比熱、xは熱拡散距離、qは単位体積当たりの発熱量を示している。式1の右辺において、第1項は熱拡散項を示し、第2項は発熱項を示している。本実施例では、一次元の熱伝導方程式を用いているが、二次元や三次元の熱伝導方程式を用いることもできる。一次元の熱伝導方程式を用いれば、単電池11の性能温度を推定するための演算処理を簡素化することができる。式1は、下記式2のように差分化することができる。式2において、iは、単電池11の厚さ方向における格子点を示している。格子点とは、図6および図7に示すように、単電池11の厚さ方向において、中心点Oおよび点Sの間の領域を複数に分割したときの、各領域内の点を示す。点Sは、単電池11の厚さ方向において、中心点Oから最も離れた位置にあり、電池ケース112の外面上に位置している。格子点の数は、適宜設定することができる。格子点の数を増やせば、単電池11の厚さ方向の位置に応じた温度の推定精度を向上させることができる。また、格子点の数を減らせば、単電池11の厚さ方向の位置に応じた温度を推定するときの演算処理を簡素化することができる。式2に示すように、格子点iの温度は、格子点iと隣り合う2つの格子点i-1,i1における温度の影響を受ける。点Sの温度は、温度センサ23によって検出される温度とみなしている。すなわち、点Sの温度と、温度センサ23が取り付けられた部分の温度とは、略等しいものとみなしている。熱伝導性に優れた金属によって電池ケース112を形成した場合には、点Sの温度と、温度センサ23が取り付けられた部分の温度とは、略等しくなる。本実施例では、単電池11の外面に温度センサ23を取り付けることにより、点Sの温度を、温度センサ23によって検出される温度と見なしているが、これに限るものではない。すなわち、単電池11の外面から離れた位置に、温度センサ23を配置することもできる。この場合には、単電池11および温度センサ23の間に存在する空気層における熱伝達を考慮する必要がある。単電池11および温度センサ23の間における熱伝達を、実験等に基づいて予め特定しておけば、温度センサ23によって検出される温度から、単電池11の外面における温度点Sの温度を推定することができる。本実施例では、図5に示すように、単電池11の厚さ方向X方向に着目しているが、これに限るものではない。発電要素111は、3次元の直方体であるため、X方向の位置だけでなく、Z方向やY方向における位置を考慮することができる。ここで、考慮すべき位置は、単電池11の伝熱経路に応じて変化する。なお、伝熱経路は、寸法、伝熱面積、単電池11の材質、空気層の有無などによっても影響を受ける。本実施例では、一例として、単電池11の内部における伝熱経路として、X方向に沿った伝熱経路を最も支配的としている。次に、単電池11の性能温度を示す格子点iを特定する方法について説明する。図8は、性能温度を示す格子点iを特定する方法を説明するフローチャートである。図8に示すフローチャートに沿って、格子点iの特定方法について説明する。ステップS101において、単電池11の抵抗および温度の関係を示すマップを作成する。具体的には、単電池11における抵抗および温度の関係を取得する。ここでは、単電池11として、単電池11の内部における温度のバラツキが緩和されたものを用いる。すなわち、X方向の位置にかかわらず、単電池11の内部における温度を略均一な温度とした後に、単電池11の抵抗を測定する。単電池11の全体を略均一な温度にするためには、例えば、特定温度の下で単電池11を十分な時間、放置しておけばよい。単電池11の温度を変えながら、抵抗を測定すれば、例えば、図9に示すマップが得られる。後述する図10のパターンでSOCが変動する場合には、SOCの依存性を考慮したマップとすることが望ましい。図9に示すマップは、単電池11の抵抗内部抵抗および性能温度が対応関係にあることを示している。温度のバラツキが十分に抑えられた状態の単電池11を用いて、抵抗を測定すれば、単電池11の内部抵抗および性能温度の対応関係が分かる。そして、図9に示すマップを用いれば、単電池11の抵抗を測定することにより、性能温度を特定することができる。図8のステップS102において、単電池11の性能温度を特定する。この性能温度は、格子点iを特定するために用いられる。まず、図10および図11に示すパターンに基づいて、単電池11の充放電を行うことにより、温度センサ23の出力に基づいて単電池11の温度を検出するとともに、図9に示すマップを用いて性能温度を特定する。図10および図11において、横軸は、時間を示し、縦軸は、電流値を示す。ここで、単電池11を充電するとき、電流値を正の値とし、単電池11を放電するとき、電流値を負の値としている。図10に示すパターンの充放電は、単電池11を発熱させる期間発熱期間で行われる。図11に示すパターンの充放電は、単電池11の温度バラツキを緩和させる期間温度緩和期間で行われる。図10に示すパターンの充放電を行った後に、図11に示すパターンの充放電を行う。発熱期間では、図10に示すように、第1パターンPcおよび第2パターンPhの充放電を1サイクルとし、このサイクルを繰り返して行う。第1パターンPcは、単電池11の抵抗を測定するために用いられる。第2パターンPhは、単電池11発電要素111を発熱させるために用いられる。図10に示す充放電サイクルの回数は、適宜設定することができる。具体的には、単電池11を発熱させて、単電池11の温度が変化しにくくなるまで、充放電サイクルを繰り返すことができる。本実施例では、第1パターンPcの充放電を開始してから2秒後における抵抗を測定している。なお、抵抗の測定は、第1パターンPcの充放電を開始してから2秒後に限るものではなく、他の時間に設定することもできる。抵抗を測定するタイミングは、各サイクルにおいて同一のタイミングであればよい。例えば、第1パターンPcの充放電を開始してから、1秒又は10秒における抵抗を測定することができる。また、第2パターンPhの充放電を行うことにより、単電池11を発熱させるとともに、所定時間後の抵抗を測定することができる。この場合には、第1パターンPcの充放電を省略することができる。さらに、第1パターンPcの充放電を行うことにより、単電池11を発熱させることもできる。単電池11の抵抗を測定するために用いられるパターンは、図10に示すパターンに限るものではない。図10に示す第1パターンPcでは、充電および放電のパルスを発生させているが、充電又は放電のパルスを発生させるだけでもよい。また、図10に示す第2パターンPhについても、充電又は放電のパルスを発生させるだけでもよい。第2パターンPhについては、単電池11を発熱させることができればよい。なお、図10に示す第1パターンPcや第2パターンPhのように、充電および放電をクーロン量が等しくなるように交互に行えば、単電池11のSOCStateOfChargeを略一定に維持することができる。この場合において、図9に示すマップに関して、SOCの依存性を取得することを省略することもできる。図10に示す充放電サイクルが終了した後、言い換えれば、発熱期間が終了した後は、図11に示す充放電サイクルを繰り返して行う。図11に示す充放電では、第1パターンPcの充放電だけを1サイクルとして行い、このサイクルを繰り返して行う。第1パターンPcの充放電を行った後、次の充放電を行うまでの間に、十分な休止時間充放電を行わない時間を確保しているため、単電池11発電要素111の温度変化は非常に小さい。図11に示す充放電サイクルの回数は、適宜設定することができる。具体的には、単電池11の発熱を停止させた後に、単電池11の温度が変化しにくくなるまで、図11に示す充放電サイクルを繰り返すことができる。図12には、発熱期間および温度緩和期間において測定された抵抗を示している。図12において、縦軸は単電池11の抵抗を示し、横軸は充放電サイクルの数言い換えれば、時間を示している。図12に示すように、発熱期間から温度緩和期間に切り替わると、単電池11の抵抗が上昇している。図12に示す抵抗と、図9に示すマップとに基づいて、性能温度を特定することができる。図13には、温度センサ23による検出温度Tsと、図9に示すマップを用いて特定された性能温度Tpとの関係を示している。図13において、縦軸は温度を示し、横軸は充放電サイクルの数言い換えれば、時間を示している。また、図13の一点鎖線で示す分布は、温度センサ23による検出温度Tsを示し、図13の実線で示す分布は、性能温度Tpを示している。図13に示すように、検出温度Tsおよび性能温度Tpは、互いに似た挙動を示しているが、発熱期間において、性能温度Tpは、検出温度Tsよりも高くなっている。また、温度緩和期間における性能温度Tpおよび検出温度Tsの差は、発熱期間における性能温度Tpおよび検出温度Tsの差よりも小さくなっている。図8のステップS103において、性能温度Tpの温度変化に最も近い温度変化を示す格子点iを特定する。温度センサ23による検出温度Tsと、式2に示す熱伝導方程式とに基づいて、各格子点の温度を算出することができる。具体的には、図6および図7において、点Sの温度が検出温度Tsとなるため、式2に示す熱伝導方程式を用いれば、点Sと隣り合う格子点の温度を算出することができる。また、この格子点と隣り合う格子点の温度は、点Sと隣り合う格子点の温度に基づいて算出することができる。この算出方法により、複数の格子点における温度を算出することができる。複数の格子点における温度のうち、性能温度Tpに最も近い温度を特定すれば、性能温度を示す格子点を特定できる。特定された格子点に関する情報は、メモリ30aに格納しておくことができる。単電池11の温度を推定する場合には、温度センサ23による検出温度と、式2に示す熱伝導方程式とに基づいて、性能温度に対応した格子点iの温度を算出する。この算出処理は、コントローラ30図1,2参照によって行われる。これにより、単電池11の内部抵抗に対応した温度を推定することができる。格子点iの温度は、単電池11の温度として、単電池11の各種制御に用いられる。本実施例では、温度センサ23による検出温度および熱伝導方程式に基づいて、性能温度に対応した格子点iの温度を算出しているが、これに限るものではない。性能温度に対応した格子点iの温度は、単電池11の外面における温度と、熱伝導方程式とに基づいて算出できればよい。例えば、温度センサ23を用いなくても、単電池11の外面における温度を特定又は推定できるのであれば、特定又は推定される温度と、熱伝導方程式とに基づいて、性能温度に対応した格子点iの温度を算出できる。次に、組電池10の充放電を制御する処理について、図14に示すフローチャートを用いて説明する。図14に示す処理は、コントローラ30によって実行される。組電池10の充放電を行うときには、例えば、組電池10の入力電力および出力電力が、予め定められた上限値を超えないように、充放電が制御される。上限値は、組電池10単電池11を保護するために予め定められた値であり、適宜設定される。すなわち、組電池10の入力電力や出力電力が上限値に到達するまでは、組電池10の充放電が許容される。例えば、組電池10を放電するとき、組電池10の出力電力が上限値出力に対応する上限値に到達すると、組電池10の放電が制限される。組電池10の放電を制限する場合には、上限値を低下させる。上限値を低下させれば、組電池10の出力として許容される電力が抑制されることになる。上限値を0[kW]まで低下させれば、組電池10の放電が行われなくなる。また、組電池10を充電するとき、組電池10の入力電力が上限値入力に対応する上限値に到達すると、組電池10の充電が制限される。組電池10の充電を制限する場合には、上限値を低下させる。上限値を低下させれば、組電池10の入力として許容される電力が抑制されることになる。上限値を0[kW]まで低下させれば、組電池10の充電が行われなくなる。組電池10の充放電制御では、例えば、組電池10の状態を監視して充放電を制御する場合や、単電池11の状態を監視して充放電を制御する場合がある。組電池10や単電池11の状態としては、電圧、電流および温度が含まれる。図14に示す処理は、組電池10の充放電制御で用いられる上限値を設定する処理である。ステップS201において、コントローラ30は、上述したように、性能温度に対応した格子点iの温度を算出推定する。ステップS202において、コントローラ30は、単電池11のSOCStateofChargeを推定する。SOCは、単電池11の満充電容量に対する、現在の単電池11の充電容量の割合を示す。単電池11のSOCは、公知の様々な方法のうち、いずれかの方法を用いて推定することができる。例えば、コントローラ30は、電流センサ22の出力から、単電池11を充放電したときの電流値を取得する。コントローラ30は、取得した電流値を積算することによって、単電池11のSOCを推定することができる。ここで、例えば、電流センサ22によって検出された値に関して、充電電流を正の値とし、放電電流を負の値とすることができる。一方、単電池11のOCVOpenCircuitVoltageを測定し、測定されたOCVからSOCを推定することができる。OCVおよびSOCは、対応関係にあるため、この対応関係を示すマップを予め用意しておけば、OCVを測定することによって、SOCを特定することができる。ステップS203において、コントローラ30は、ステップS201で算出した温度と、ステップS202で推定したSOCとに基づいて、単電池11を放電出力および充電入力するときの上限値を設定する。図15は、単電池11の温度およびSOCと、単電池11を放電するときの上限値出力上限値というとの関係を示すマップである。図15の矢印で示す高および低は、出力上限値の高低関係を示している。図15に示すマップは、予め用意しておき、メモリ30aに格納しておくことができる。図15では、SOCおよび温度が変化したときの出力上限値の変化を示している。具体的には、SOCが変化しないとき、温度が低下することに応じて、出力上限値が低下する。温度が低下するほど、放電時の化学反応が抑制されるため、単電池11の出力を抑制することが好ましい。したがって、温度が低下することに応じて、出力上限値を低下させている。図15に示すマップには、温度が変化しても、出力上限値が変化しない場合が含まれていてもよい。一方、温度が変化しないとき、SOCが低下することに応じて、出力上限値が低下する。単電池11のSOCが低下するほど、単電池11を放電し難くなるため、SOCが低下することに応じて、出力上限値を低下させている。図15に示すマップには、SOCが変化しても、出力上限値が変化しない場合が含まれていてもよい。図15に示すマップにおける温度は、ステップS201で算出した温度である。図15のマップにおけるSOCは、ステップS202で推定したSOCである。単電池11の温度およびSOCを特定できれば、図15に示すマップから、出力上限値を特定することができる。そして、設定された出力上限値を基準として、単電池11組電池10の放電制御が行われる。図15は、単電池11を放電するときの上限値を示しているが、単電池11を充電するときの上限値についても、図15に対応したマップを用意しておけばよい。そして、出力上限値を特定する方法と同様の方法によって、単電池11を充電するときの上限値入力上限値というを特定することができる。単電池11のSOCが変化しないとき、単電池11の温度が低下することに応じて、入力上限値を低下させることができる。ここで、温度が変化しても、入力上限値を変化させないこともできる。また、温度が変化しないとき、SOCが上昇することに応じて、入力上限値を低下させることができる。単電池11のSOCが上昇するほど、単電池11を充電し難くなるため、SOCが上昇することに応じて、入力上限値を低下させている。設定された入力上限値を基準として、単電池11組電池10の充電制御が行われる。本実施例では、図15に示すマップを用いて、出力上限値を特定しているが、これに限るものではない。具体的には、出力上限値が、単電池11の温度およびSOCを変数とした関数で表されるとき、この関数を用いて、出力上限値を算出することができる。この場合には、温度およびSOCを変数とした関数に関する情報を、メモリ30aに格納しておけばよい。また、本実施例では、単電池11の温度およびSOCに基づいて、充放電制御の上限値を設定しているが、これに限るものではない。具体的には、単電池11の温度だけに基づいて、充放電制御の上限値を設定することができる。本実施例によれば、単電池11の性能温度に対応した上限値出力上限値や入力上限値を設定しているため、単電池11の入出力特性に応じた充放電を行うことができる。単電池11の性能温度は、単電池11の内部抵抗に対応した温度であり、単電池11の入出力特性は、単電池11の内部抵抗に依存する。したがって、単電池11の性能温度に対応した上限値を設定すれば、単電池11の性能に応じた充放電を行うことができる。図5で説明したように、単電池11の表面温度点Sの温度は、単電池11の性能温度よりも低い傾向がある。このため、単電池11の表面温度に基づいて、出力上限値や入力上限値を設定すれば、単電池11の充放電を過度に制限してしまうことがある。一方、単電池11発電要素111の中心における温度は、性能温度よりも高い傾向がある。このため、単電池11の中心温度最も高い温度に基づいて、出力上限値や入力上限値を設定すれば、単電池11の充放電を過度に行ってしまうことがある。一方、格子点iの温度性能温度に基づいて、単電池11の温度調節を行うことができる。格子点iの温度が上昇していれば、単電池11に対して冷却用の熱交換媒体を供給することにより、単電池11の温度上昇を抑制することができる。また、格子点iの温度性能温度に基づいて、単電池11のSOHを推定することができる。特開2004-257781号公報には、開放電圧および積算電気量に基づいて、SOHを推定しており、開放電圧を温度で補正している。この場合には、開放電圧を補正するための温度情報として、本実施例で説明した性能温度を用いることができる。本実施例によれば、単電池11の内部抵抗に対応した温度を推定することができる。また、推定した温度性能温度に基づいて、単電池11の内部状態SOC等を推定すれば、内部状態の推定精度を向上させることができる。特に、温度が低下するとともに、SOCが低下した状態において、SOCの推定精度を向上させることができる。上述した実施例では、単電池11の厚さ方向において設けられた複数の格子点のうち、性能温度に対応した格子点を特定しておき、特定された格子点の温度を、単電池11の温度として推定している。ここで、3つの格子点だけを用いて、単電池11の性能温度を算出することができる。図16には、格子点を3つに設定した場合において、単電池11の厚さ方向における位置と、温度との関係を示している。図16に示す例では、3つの格子点のうち、単電池11の内部における格子点の温度を性能温度Tpとし、単電池11の表面における格子点の温度を温度センサ23による検出温度Tsとしている。ここで、検出温度Tsの代わりに、他の格子点の温度を用いることもできる。ここで、格子点を3つに設定したときの熱伝導方程式は、下記式3のように簡易化することができる。また、式3は、下記式4で表すことができる。ここで、αおよびβは、下記式5、6によって表される。式3や式5、6におけるk1,k2は、補正係数を示しており、例えば、以下に説明する方法によって決定することができる。図13を用いて説明したように、単電池11の性能温度を算出しておく。次に、補正係数k1,k2のそれぞれを変更しながら、式3や式4を用いて、性能温度と推定される温度推定温度を算出する。そして、性能温度および推定温度の差が最小となる補正係数k1,k2を特定する。図17には、性能温度および推定温度の関係一例を示している。発熱期間では、第1パターンPcおよび第2パターンPhの充放電図10参照を交互に行うことにより、単電池11を発熱させている。温度緩和期間では、第1パターンPcだけの充放電図11参照を行うことにより、単電池11を発熱させずに、単電池11の温度を環境に応じた温度に到達させている。なお、第2パターンPhの充放電を行うことにより、単電池11を発熱させるとともに、所定時間後の抵抗を測定することができる。この場合には、第1パターンPcの充放電を省略することができる。また、第1パターンPcの充放電を行うことにより、単電池11を発熱させることもできる。発熱期間および温度緩和期間において、単電池11の抵抗を測定し、測定された抵抗と図9に示すマップを用いて、性能温度を特定することができる。性能温度の分布を、図17の実線で示している。一方、発熱期間および温度緩和期間において取得した温度センサ23の検出温度を、式3に示す熱伝導方程式に代入し、補正係数k1,k2を適宜設定することにより、推定温度を特定する。推定温度の分布一例を図17の点線で示している。図17に示すように、推定温度が性能温度よりも高いときには、推定温度を低下させて性能温度に近づけるように、補正係数k1,k2を変更する。また、推定温度が性能温度よりも低いときには、推定温度を上昇させて性能温度に近づけるように、補正係数k1,k2を変更する。すなわち、推定温度および性能温度の差ΔTがゼロに近づくように、補正係数k1,k2を決定する。ここで、補正係数k1,k2を決定する処理は、発熱期間における推定温度および性能温度の差ΔTに基づいて行ったり、温度緩和期間における推定温度および性能温度の差ΔTに基づいて行ったりすることができる。得られた補正係数k1,k2又は、αおよびβは、メモリ30aに格納しておくことができる。これにより、温度センサ23による検出温度Tsを取得すれば、式3に基づいて、性能温度Tpを算出することができる。本実施例の電池システムにおいて、組電池10の充放電を停止した後では、単電池11の内部における温度バラツキが緩和される。温度バラツキが緩和されることにより、単電池11の外面における温度、言い換えれば、温度センサ23による検出温度と、単電池11の性能温度とが略等しくなる。組電池10の充放電を停止した直後では、単電池11の内部において、温度のバラツキが発生している。すなわち、図5を用いて説明したように、単電池11の内部における温度は、単電池11の外面における温度よりも高くなっている。一方、組電池10の充放電を停止する時間が長くなるほど、単電池11の内部における温度のバラツキが緩和される。すなわち、単電池11の内部における温度は、単電池11の外面における温度と略等しくなる。そこで、組電池10の充放電を再開させるときには、組電池10の充放電を停止する時間に応じて、単電池11の性能温度を推定する方法を変更することができる。図18は、単電池11の性能温度を推定する処理を説明するフローチャートである。図18に示す処理は、コントローラ30によって実行される。ステップS301において、コントローラ30は、イグニッションスイッチがオンからオフに切り替わったことを確認する。イグニッションスイッチのオンおよびオフに関する情報は、コントローラ30に入力される。コントローラ30は、イグニッションスイッチのオフに関する情報を受けると、図1に示す電池システムにおいて、システムメインリレーSMR-B,SMR-Gをオンからオフに切り替える。これにより、組電池10および負荷の接続が遮断され、組電池10の充放電が停止する。ステップS302において、コントローラ30は、イグニッションスイッチがオンからオフに切り替わったタイミングから、タイマを用いて時間toffの計測を開始する。また、コントローラ30は、イグニッションスイッチがオンからオフに切り替わったタイミングにおいて、温度センサ23による検出温度Ts1を取得する。検出温度Tsに関する情報は、メモリ30aに記憶することができる。ステップS303において、コントローラ30は、イグニッションスイッチがオフからオンに切り替わったことを確認する。イグニッションスイッチのオンに関する情報がコントローラ30に入力されることにより、コントローラ30は、イグニッションスイッチがオンであることを確認できる。ステップS304において、コントローラ30は、時間toffの計測を終了するとともに、温度センサ23による検出温度Ts2を取得する。検出温度Ts2は、組電池10単電池11の充放電が停止された後の単電池11の温度となる。検出温度Ts2に関する情報は、メモリ30aに記憶することができる。ステップS305において、コントローラ30は、計測時間toffが所定時間tthよりも長いか否かを判別する。所定時間tthは、単電池11の内部における温度バラツキを解消させる観点から、予め定められた時間である。単電池11の充放電を停止した状態において、単電池11を放熱させる実験を行うことにより、単電池11の内部における温度バラツキが解消する時間を求めることができる。所定時間tthに関する情報は、メモリ30aに記憶することができる。所定時間tthは、固定値であってもよいし、変更してもよい。所定時間tthを変更するときには、例えば、ステップS302で取得した検出温度Ts1に応じて、所定時間tthを変更することができる。具体的には、検出温度Ts1および所定時間tthの関係を示すマップを予め求めておき、このマップを用いて、ステップS302で取得した検出温度Ts1に対応した所定時間tthを特定することができる。検出温度Ts1および所定時間tthの対応関係を示すマップは、メモリ30aに記憶することができる。ステップS305において、計測時間toffが所定時間tthよりも長いときには、ステップS306の処理に進む。計測時間toffが所定時間tthよりも短いときには、ステップS307の処理に進む。ステップS306において、コントローラ30は、単電池11の性能温度Tpが、温度センサ23による検出温度Ts2と等しいと判別する。組電池10の充放電を再開させたとき、コントローラ30は、性能温度Tpが検出温度Ts2であるとして、単電池11の内部状態を推定したり、単電池11組電池10の入出力を制御したりすることができる。ステップS307において、コントローラ30は、単電池11の性能温度Tpが、温度センサ23による検出温度Ts2とは異なっていると判別する。この場合において、コントローラ30は、検出温度Ts2を用いて、単電池11の性能温度Tpを推定することができる。図18に示す処理によれば、図19に示すように、イグニッションスイッチがオンからオフに切り替わったときの温度センサ23による検出温度Ts1と、イグニッションスイッチがオフからオンに切り替わったときの温度センサ23による検出温度Ts2とを取得できる。検出温度Ts1,Ts2を取得できれば、イグニッションスイッチがオフからオンに切り替わったときに、単電池11の性能温度Tpを推定することができる。例えば、イグニッションスイッチがオフの間の温度変化が時間に対して線形や指数関数などとなる関係を予め実験により求めておけば、イグニッションスイッチがオフの間検出温度Ts1から検出温度Ts2になる間の時間刻みのみを規定し、このときの温度センサ23の検出温度を補完することで、式4を用いることにより、性能温度Tpの変化量を算出することができる。本実施例では、計測時間toffが所定時間tthよりも長いか否かに基づいて、単電池11の内部における温度バラツキが解消されているか否かを判別しているが、これに限るものではない。具体的には、イグニッションスイッチがオフとなっている間において、コントローラ30は、所定周期Δtが経過するたびに、検出温度Tsを取得する。所定周期Δtで取得される検出温度Tsは、メモリ30aに記憶することができる。コントローラ30は、所定周期Δtが経過したときの検出温度Tsの変化量ΔTsを算出し、所定周期Δtに対する変化量ΔTsの比ΔTs/Δtを算出する。コントローラ30は、比ΔTs/Δtが閾値よりも小さいか否かに基づいて、単電池11の内部における温度バラツキが解消されているか否かを判別することができる。所定周期Δtは、適宜設定することができる。所定周期Δtを短くするほど、検出温度Tsが変化する挙動を追跡しやすくなるが、メモリ30aに記憶される検出温度Tsのデータ量が増加してしまう。一方、所定周期Δtを長くするほど、メモリ30aに記憶される検出温度Tsのデータ量を低減できるが、検出温度Tsが変化する挙動を追跡し難くなる。この点を考慮して、所定周期Δtを設定することができる。単電池11の内部における温度バラツキが解消し始めると、比ΔTs/Δtが低下し始める。このため、比ΔTs/Δtが閾値よりも小さいか否かを判別することにより、単電池11の内部における温度バラツキが解消されているか否かを判別することができる。閾値は、検出温度Tsのバラツキを考慮して適宜設定することができる。温度バラツキが完全に解消され、温度差ΔTsが変化しないときには、比ΔTs/Δtが0となる。したがって、閾値を0に設定しておけば、温度バラツキが完全に解消されたことを判別することができる。比ΔTs/Δtが閾値よりも小さいとき、コントローラ30は、単電池11の内部における温度バラツキが解消されていると判別することができる。温度バラツキが解消されているときには、ステップS306の処理を行うことができる。また、比ΔTs/Δtが閾値よりも大きいとき、コントローラ30は、単電池11の内部における温度バラツキが解消されていないと判別することができる。温度バラツキが解消されていないときには、ステップS307の処理を行うことができる。例えば、周期的に計測した時間と、温度センサ23の検出温度の履歴から、式4を用いて、性能温度Tpを推定することができる。性能温度Tpの推定計算は、イグニッションスイッチのオンと同時に行ってもよい。本実施例によれば、イグニッションスイッチがオフからオンに切り替わったとき、言い換えれば、組電池10単電池11の充放電が再開されたとき、単電池11の内部における温度バラツキの解消状態に応じて、単電池11の性能温度を特定することができる。本実施例では、単電池11の外面における温度と、熱伝導方程式とに基づいて、性能温度に対応した格子点iの温度を算出しているが、これに限るものではない。具体的には、単電池11の外部における温度を取得しておき、熱の移動を表す式を用いて、性能温度に対応した格子点iの温度を算出することができる。単電池11の外部における温度とは、単電池11の外面における温度だけでなく、単電池11の外面から離れた位置における温度も含まれる。単電池11の外面から離れた位置における温度は、温度センサ23を用いて検出することができる。単電池11の外面から離れた位置としては、例えば、単電池11に熱交換媒体気体や液体を供給して、単電池11の温度を調節するときには、熱交換媒体が単電池11に供給される箇所例えば、熱交換媒体の供給口とすることができる。熱の移動を表す式とは、単電池11の外部から、性能温度に対応した格子点iまでの間における熱の移動を表す式である。熱の移動を表す式には、本実施例で説明した熱伝導方程式の他にも、例えば、熱伝達を表す式や、熱等価回路のモデルに基づく式が含まれる。単電池11の外面から離れた位置における温度を取得するときには、単電池11の外面から離れた位置から単電池11の外面までの間の熱伝達と、単電池11の内部における熱伝導とを考慮して、性能温度に対応した格子点iの温度を算出することができる。ここで、熱伝達および熱伝導を個別に評価して、格子点iの温度を算出することもできるし、熱伝達および熱伝導の両者をまとめて評価して、格子点iの温度を算出することもできる。熱伝達を表す式を用いるときには、単電池11の外面から離れた位置から単電池11の外面までの間における熱の移動を考慮して、熱伝達を表す式を適宜決めることができる。単電池11の内部における熱伝導や、単電池11の外面から離れた位置から単電池11の外面までの間の熱伝達を評価するときには、本実施例で説明した熱伝導方程式の他にも、熱等価回路のモデルを用いることができる。熱等価回路のモデルでは、熱の伝わり方を考慮して、熱容量Cおよび熱抵抗Rを用いて表すことができる。そして、熱等価回路のモデルに基づいて、熱伝導や熱伝達を表す式を適宜決めることができる。"