"固体燃料を燃焼して推力を得るロケットエンジンを複数段搭載した多段式ロケットを誘導する多段式ロケット誘導装置において、前記多段式ロケットの位置ベクトルと速度ベクトルと推力加速度ベクトルとを計測する航法装置から現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルとを入力し、入力した現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルと、現在時刻と、推力を得ている現段のロケットエンジンが現段の固体燃料を燃焼し終わる予定時刻として予め記憶する現段燃焼終了時刻と、前記現段燃焼終了時刻までの推力加速度の予測値を示す現段推力加速度プロファイルと、に基づいて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、最終段のロケットエンジンが最終段の固体燃料を燃焼し終わる予定時刻として予め記憶する最終段燃焼終了時刻と、前記最終段燃焼終了時刻までの推力加速度の予測値を示す最終段推力加速度プロファイルとに基づいて、前記最終段燃焼終了時刻に前記多段式ロケットが飛行する予測の楕円軌道の軌道長半径と近地点高度とを積分計算によって算出する予測演算部と、前記予測演算部によって算出された前記予測の楕円軌道の軌道長半径と、前記予測演算部によって算出された前記予測の楕円軌道の近地点高度と目標の近地点高度との差とに基づいて前記多段式ロケットの新たな推力方向の角度を算出し、算出した前記多段式ロケットの新たな推力方向の角度を前記多段式ロケットの推力方向を制御するロケット制御装置に出力することによって前記目標の近地点高度に位置する近地点を通って前記予測の楕円軌道より遠地点高度が大きい楕円軌道へ前記多段式ロケットを誘導する誘導演算部とを備えることを特徴とする多段式ロケット誘導装置。前記予測演算部は、前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、次段のロケットエンジンが固体燃料を燃焼し始める予定時刻として予め記憶する次段燃焼開始時刻と、に基づいて前記次段燃焼開始時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記次段燃焼開始時刻の予測位置ベクトルと予測速度ベクトルと、前記次段燃焼終了時刻と、前記最終段燃焼終了時刻と、前記次段燃焼終了時刻までの推力加速度の予測値を示す次段推力加速度プロファイルと、前記最終段推力加速度プロファイルと、に基づいて前記予測の楕円軌道の前記軌道長半径と前記近地点高度とを積分計算によって算出することを特徴とする請求項1記載の多段式ロケット誘導装置。前記予測演算部は、前記ロケット制御装置が前記多段式ロケットの推力方向の角度を前記新たな推力方向の角度に変更した後、前記航法装置から新たな現在の位置ベクトルと新たな現在の速度ベクトルと新たな現在の推力加速度ベクトルとを入力し、前記新たな現在の位置ベクトルと前記新たな現在の速度ベクトルと新たな現在時刻と前記現段燃焼終了時刻と前記現段推力加速度プロファイルとを用いて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを新たに算出することを特徴とする請求項1または請求項2記載の多段式ロケット誘導装置。固体燃料を燃焼して推力を得るロケットエンジンを複数段搭載した多段式ロケットを誘導する多段式ロケット誘導装置において、前記多段式ロケットの位置ベクトルと速度ベクトルと推力加速度ベクトルとを計測する航法装置から現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルとを入力し、入力した現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルと、現在時刻と、推力を得ている現段のロケットエンジンが現段の固体燃料を燃焼し終わる予定時刻として予め記憶する現段燃焼終了時刻と、前記現段燃焼終了時刻までの推力加速度の予測値を示す現段推力加速度プロファイルと、に基づいて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、最終段のロケットエンジンが最終段の固体燃料を燃焼し終わる予定時刻として予め記憶する最終段燃焼終了時刻と、前記最終段燃料終了時刻までの推力加速度の予測値を示す最終段推力加速度プロファイルとに基づいて、前記最終段燃焼終了時刻に前記多段式ロケットが飛行する予測の楕円軌道の軌道長半径と遠地点高度とを積分計算によって算出する予測演算部と、前記予測演算部によって算出された前記予測の楕円軌道の軌道長半径と、前記予測演算部によって算出された前記予測の楕円軌道の遠地点高度と目標の遠地点高度との差とに基づいて前記多段式ロケットの新たな推力方向の角度を算出し、算出した前記多段式ロケットの新たな推力方向の角度を前記多段式ロケットの推力方向を制御するロケット制御装置に出力することによって前記目標の遠地点高度に位置する遠地点を通って前記予測の楕円軌道より近地点高度が大きい楕円軌道へ前記多段式ロケットを誘導する誘導演算部とを備えることを特徴とする多段式ロケット誘導装置。固体燃料を燃焼して推力を得るロケットエンジンを複数段搭載した多段式ロケットを誘導する多段式ロケット誘導装置を機能させる多段式ロケット誘導プログラムであって、前記多段式ロケットの位置ベクトルと速度ベクトルと推力加速度ベクトルとを計測する航法装置から現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルとを入力し、入力した現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルと、現在時刻と、推力を得ている現段のロケットエンジンが現段の固体燃料を燃焼し終わる予定時刻として予め記憶する現段燃焼終了時刻と、前記現段燃焼終了時刻までの推力加速度の予測値を示す現段推力加速度プロファイルと、に基づいて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、最終段のロケットエンジンが最終段の固体燃料を燃焼し終わる予定時刻として予め記憶する最終段燃焼終了時刻と、前記最終段燃焼終了時刻までの推力加速度の予測値を示す最終段推力加速度プロファイルとに基づいて、前記最終段燃焼終了時刻に前記多段式ロケットが飛行する予測の楕円軌道の軌道長半径と近地点高度とを積分計算によって算出する予測演算部と、前記予測演算部によって算出された前記予測の楕円軌道の軌道長半径と、前記予測演算部によって算出された前記予測の楕円軌道の近地点高度と目標の近地点高度との差とに基づいて前記多段式ロケットの新たな推力方向の角度を算出し、算出した前記多段式ロケットの新たな推力方向の角度を前記多段式ロケットの推力方向を制御するロケット制御装置に出力することによって前記目標の近地点高度に位置する近地点を通って前記予測の楕円軌道より遠地点高度が大きい楕円軌道へ前記多段式ロケットを誘導する誘導演算部として前記多段式ロケット誘導装置を機能させることを特徴とする多段式ロケット誘導プログラム。前記予測演算部は、前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、次段のロケットエンジンが固体燃料を燃焼し始める予定時刻として予め記憶する次段燃焼開始時刻と、に基づいて前記次段燃焼開始時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記次段燃焼開始時刻の予測位置ベクトルと予測速度ベクトルと、前記次段燃焼終了時刻と、前記最終段燃焼終了時刻と、前記次段燃焼終了時刻までの推力加速度の予測値を示す次段推力加速度プロファイルと、前記最終段推力加速度プロファイルと、に基づいて前記予測の楕円軌道の前記軌道長半径と前記近地点高度とを積分計算によって算出することを特徴とする請求項5記載の多段式ロケット誘導プログラム。前記予測演算部は、前記ロケット制御装置が前記多段式ロケットの推力方向の角度を前記新たな推力方向の角度に変更した後、前記航法装置から新たな現在の位置ベクトルと新たな現在の速度ベクトルと新たな現在の推力加速度ベクトルとを入力し、前記新たな現在の位置ベクトルと前記新たな現在の速度ベクトルと新たな現在時刻と前記現段燃焼終了時刻と前記現段推力加速度プロファイルとを用いて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを新たに算出することを特徴とする請求項5または請求項6記載の多段式ロケット誘導プログラム。固体燃料を燃焼して推力を得るロケットエンジンを複数段搭載した多段式ロケットを誘導する多段式ロケット誘導装置を機能させる多段式ロケット誘導プログラムであって、前記多段式ロケットの位置ベクトルと速度ベクトルと推力加速度ベクトルとを計測する航法装置から現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルとを入力し、入力した現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルと、現在時刻と、推力を得ている現段のロケットエンジンが現段の固体燃料を燃焼し終わる予定時刻として予め記憶する現段燃焼終了時刻と、前記現段燃焼終了時刻までの推力加速度の予測値を示す現段推力加速度プロファイルと、に基づいて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、最終段のロケットエンジンが最終段の固体燃料を燃焼し終わる予定時刻として予め記憶する最終段燃焼終了時刻と、前記最終段燃焼終了時刻までの推力加速度の予測値を示す最終段推力加速度プロファイルとに基づいて、前記最終段燃焼終了時刻に前記多段式ロケットが飛行する予測の楕円軌道の軌道長半径と遠地点高度とを積分計算によって算出する予測演算部と、前記予測演算部によって算出された前記予測の楕円軌道の軌道長半径と、前記予測演算部によって算出された前記予測の楕円軌道の遠地点高度と目標の遠地点高度との差とに基づいて前記多段式ロケットの新たな推力方向の角度を算出し、算出した前記多段式ロケットの新たな推力方向の角度を前記多段式ロケットの推力方向を制御するロケット制御装置に出力することによって前記目標の遠地点高度に位置する遠地点を通って前記予測の楕円軌道より近地点高度が大きい楕円軌道へ前記多段式ロケットを誘導する誘導演算部として前記多段式ロケット誘導装置を機能させることを特徴とする多段式ロケット誘導プログラム。固体燃料を燃焼して推力を得るロケットエンジンを複数段搭載した多段式ロケットを誘導する多段式ロケット誘導装置によって実行する多段式ロケット誘導方法において、予測演算部が、前記多段式ロケットの位置ベクトルと速度ベクトルと推力加速度ベクトルとを計測する航法装置から現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルとを入力し、入力した現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルと、現在時刻と、推力を得ている現段のロケットエンジンが現段の固体燃料を燃焼し終わる予定時刻として予め記憶する現段燃焼終了時刻と、前記現段燃焼終了時刻までの推力加速度の予測値を示す現段推力加速度プロファイルと、に基づいて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、最終段のロケットエンジンが最終段の固体燃料を燃焼し終わる予定時刻として予め記憶する最終段燃焼終了時刻と、前記最終段燃焼終了時刻までの推力加速度の予測値を示す最終段推力加速度プロファイルとに基づいて、前記最終段燃焼終了時刻に前記多段式ロケットが飛行する予測の楕円軌道の軌道長半径と近地点高度とを積分計算によって算出し、誘導演算部が、前記予測演算部によって算出された前記予測の楕円軌道の軌道長半径と、前記予測演算部によって算出された前記予測の楕円軌道の近地点高度と目標の近地点高度との差とに基づいて前記多段式ロケットの新たな推力方向の角度を算出し、算出した前記多段式ロケットの新たな推力方向の角度を前記多段式ロケットの推力方向を制御するロケット制御装置に出力することによって前記目標の近地点高度に位置する近地点を通って前記予測の楕円軌道より遠地点高度が大きい楕円軌道へ前記多段式ロケットを誘導することを特徴とする多段式ロケット誘導方法。固体燃料を燃焼して推力を得るロケットエンジンを複数段搭載した多段式ロケットを誘導する多段式ロケット誘導装置によって実行する多段式ロケット誘導方法おいて、予測演算部が、前記多段式ロケットの位置ベクトルと速度ベクトルと推力加速度ベクトルとを計測する航法装置から現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルとを入力し、入力した現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルと、現在時刻と、推力を得ている現段のロケットエンジンが現段の固体燃料を燃焼し終わる予定時刻として予め記憶する現段燃焼終了時刻と、前記現段燃焼終了時刻までの推力加速度の予測値を示す現段推力加速度プロファイルとに基づいて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、最終段のロケットエンジンが最終段の固体燃料を燃焼し終わる予定時刻として予め記憶する最終段燃焼終了時刻と、前記最終段燃焼終了時刻までの最終段推力加速度プロファイルとに基づいて、前記最終段燃焼終了時刻に前記多段式ロケットが飛行する予測の楕円軌道の軌道長半径と遠地点高度とを積分計算によって算出し、誘導演算部が、前記予測演算部によって算出された前記予測の楕円軌道の軌道長半径と、前記予測演算部によって算出された前記予測の楕円軌道の遠地点高度と目標の遠地点高度との差とに基づいて前記多段式ロケットの新たな推力方向の角度を算出し、算出した前記多段式ロケットの新たな推力方向の角度を前記多段式ロケットの推力方向を制御するロケット制御装置に出力することによって前記目標の遠地点高度に位置する遠地点を通って前記予測の楕円軌道より近地点高度が大きい楕円軌道へ前記多段式ロケットを誘導することを特徴とする多段式ロケット誘導方法。固体燃料を燃焼して推力を得るロケットエンジンを複数段搭載した多段式ロケットの位置ベクトルと速度ベクトルとを計測する航法装置と、前記多段式ロケットの推力方向を制御するロケット制御装置と、前記多段式ロケットを誘導する多段式ロケット誘導装置と、を備える多段式ロケット誘導システムにおいて、前記多段式ロケット誘導装置は、前記航法装置から現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルとを入力し、入力した現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルと、現在時刻と、推力を得ている現段のロケットエンジンが現段の固体燃料を燃焼し終わる予定時刻として予め記憶する現段燃焼終了時刻と、前記現段燃焼終了時刻までの推力加速度の予測値を示す現段推力加速度プロファイルに基づいて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、最終段のロケットエンジンが最終段の固体燃料を燃焼し終わる予定時刻として予め記憶する最終段燃焼終了時刻と、前記最終段燃焼終了時刻までの推力加速度の予測値を示す最終段推力加速度プロファイルとに基づいて、前記最終段燃焼終了時刻に前記多段式ロケットが飛行する予測の楕円軌道の軌道長半径と近地点高度とを積分計算によって算出する予測演算部と、前記予測演算部によって算出された前記予測の楕円軌道の軌道長半径と、前記予測演算部によって算出された前記予測の楕円軌道の近地点高度と目標の近地点高度との差とに基づいて前記多段式ロケットの新たな推力方向の角度を算出し、算出した前記多段式ロケットの新たな推力方向の角度を前記ロケット制御装置に出力することによって前記目標の近地点高度に位置する近地点を通って前記予測の楕円軌道より遠地点高度が大きい楕円軌道へ前記多段式ロケットを誘導する誘導演算部とを備えることを特徴とする多段式ロケット誘導システム。固体燃料を燃焼して推力を得るロケットエンジンを複数段搭載した多段式ロケットの位置ベクトルと速度ベクトルとを計測する航法装置と、前記多段式ロケットの推力方向を制御するロケット制御装置と、前記多段式ロケットを誘導する多段式ロケット誘導装置と、を備える多段式ロケット誘導システムにおいて、前記多段式ロケット誘導装置は、前記航法装置から現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルとを入力し、入力した現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルと、現在時刻と、推力を得ている現段のロケットエンジンが現段の固体燃料を燃焼し終わる予定時刻として予め記憶する現段燃焼終了時刻と、前記現段燃焼終了時刻までの推力加速度の予測値を示す現段推力加速度プロファイルとに基づいて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、最終段のロケットエンジンが最終段の固体燃料を燃焼し終わる予定時刻として予め記憶する最終段燃焼終了時刻と、前記最終段燃焼終了時刻までの推力加速度を示す最終段推力加速度プロファイルとに基づいて、前記最終段燃焼終了時刻に前記多段式ロケットが飛行する予測の楕円軌道の軌道長半径と遠地点高度とを積分計算によって算出する予測演算部と、前記予測演算部によって算出された前記予測の楕円軌道の軌道長半径と、前記予測演算部によって算出された前記予測の楕円軌道の遠地点高度と目標の遠地点高度との差とに基づいて前記多段式ロケットの新たな推力方向の角度を算出し、算出した前記多段式ロケットの新たな推力方向の角度を前記ロケット制御装置に出力することによって前記目標の遠地点高度に位置する遠地点を通って前記予測の楕円軌道より近地点高度が大きい楕円軌道へ前記多段式ロケットを誘導する誘導演算部とを備えることを特徴とする多段式ロケット誘導システム。本発明は、例えば、空中発射システムまたは中間の目標地点を要しない場合の地上発射システムで発射する固体ロケットをできるだけ遠地点高度の大きな楕円軌道、もしくは、できるだけ近地点高度の大きな楕円軌道に誘導できるようにすることを目的とする。本発明の多段式ロケット誘導装置は、固体燃料を燃焼して推力を得るロケットエンジンを複数段搭載した多段式ロケットを誘導する。前記多段式ロケット誘導装置は、前記多段式ロケットの位置ベクトルと速度ベクトルと推力加速度ベクトルを計測する航法装置から現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルを入力し、入力した現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルと、現在時刻と、推力を得ている現段のロケットエンジンが現段の固体燃料を燃焼し終わる予定時刻として予め記憶する現段燃焼終了時刻と、前記現段燃焼終了時刻までの推力加速度の予測値を示す現段推力加速度プロファイルとに基づいて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、最終段のロケットエンジンが最終段の固体燃料を燃焼し終わる予定時刻として予め記憶する最終段燃焼終了時刻と、前記最終段燃焼終了時刻までの推力加速度の予測値を示す最終段推力加速度プロファイルとに基づいて、前記最終段燃焼終了時刻に前記多段式ロケットが飛行する予測の楕円軌道の軌道長半径と近地点高度とを積分計算によって算出する予測演算部と、前記予測演算部によって算出された前記予測の楕円軌道の軌道長半径と、前記予測演算部によって算出された前記予測の楕円軌道の近地点高度と目標の近地点高度との差とに基づいて前記多段式ロケットの新たな推力方向の角度を算出し、算出した前記多段式ロケットの新たな推力方向の角度を前記多段式ロケットの推力方向を制御するロケット制御装置に出力することによって前記目標の近地点高度に位置する近地点を通って前記予測の楕円軌道より遠地点高度が大きい楕円軌道へ前記多段式ロケットを誘導する誘導演算部とを備える。前記予測演算部は、前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、次段のロケットエンジンが固体燃料を燃焼し始める予定時刻として予め記憶する次段燃焼開始時刻と、に基づいて前記次段燃焼開始時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記次段燃焼開始時刻の予測位置ベクトルと予測速度ベクトルと、前記次段燃焼終了時刻と、前記最終段燃焼終了時刻と、前記次段燃焼終了時刻までの推力加速度の予測値を示す次段推力加速度プロファイルと、前記最終段推力加速度プロファイルと、に基づいて前記予測の楕円軌道の前記軌道長半径と前記近地点高度とを積分計算によって算出する。前記予測演算部は、前記ロケット制御装置が前記多段式ロケットの推力方向の角度を前記新たな推力方向の角度に変更した後、前記航法装置から新たな現在の位置ベクトルと新たな現在の速度ベクトルと新たな現在の推力加速度ベクトルとを入力し、前記新たな現在の位置ベクトルと前記新たな現在の速度ベクトルと新たな現在時刻と前記現段燃焼終了時刻と前記現段推力加速度プロファイルとを用いて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを新たに算出する。本発明の多段式ロケット誘導装置は、固体燃料を燃焼して推力を得るロケットエンジンを複数段搭載した多段式ロケットを誘導する。前記多段式ロケット誘導装置は、前記多段式ロケットの位置ベクトルと速度ベクトルと推力加速度ベクトルを計測する航法装置から現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルを入力し、入力した現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルと、現在時刻と、推力を得ている現段のロケットエンジンが現段の固体燃料を燃焼し終わる予定時刻として予め記憶する現段燃焼終了時刻と、前記現段燃焼終了時刻までの推力加速度の予測値を示す現段推力加速度プロファイルとに基づいて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、最終段のロケットエンジンが最終段の固体燃料を燃焼し終わる予定時刻として予め記憶する最終段燃焼終了時刻と、前記最終段燃焼終了時刻までの推力加速度の予測値を示す最終段推力加速度プロファイルとに基づいて、前記最終段燃焼終了時刻に前記多段式ロケットが飛行する予測の楕円軌道の軌道長半径と遠地点高度とを積分計算によって算出する予測演算部と、前記予測演算部によって算出された前記予測の楕円軌道の軌道長半径と、前記予測演算部によって算出された前記予測の楕円軌道の遠地点高度と目標の遠地点高度との差とに基づいて前記多段式ロケットの新たな推力方向の角度を算出し、算出した前記多段式ロケットの新たな推力方向の角度を前記多段式ロケットの推力方向を制御するロケット制御装置に出力することによって前記目標の遠地点高度に位置する遠地点を通って前記予測の楕円軌道より近地点高度が大きい楕円軌道へ前記多段式ロケットを誘導する誘導演算部とを備える。本発明の多段式ロケット誘導プログラムは、固体燃料を燃焼して推力を得るロケットエンジンを複数段搭載した多段式ロケットを誘導する多段式ロケット誘導装置を機能させる。前記多段式ロケット誘導プログラムは、前記多段式ロケットの位置ベクトルと速度ベクトルと推力加速度ベクトルとを計測する航法装置から現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルとを入力し、入力した現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルと、現在時刻と、推力を得ている現段のロケットエンジンが現段の固体燃料を燃焼し終わる予定時刻として予め記憶する現段燃焼終了時刻と、前記現段燃焼終了時刻までの推力加速度の予測値を示す現段推力加速度プロファイルとに基づいて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、最終段のロケットエンジンが最終段の固体燃料を燃焼し終わる予定時刻として予め記憶する最終段燃焼終了時刻と、前記最終段燃焼終了時刻までの推力加速度の予測値を示す最終段推力加速度プロファイルとに基づいて、前記最終段燃焼終了時刻に前記多段式ロケットが飛行する予測の楕円軌道の軌道長半径と近地点高度とを積分計算によって算出する予測演算部と、前記予測演算部によって算出された前記予測の楕円軌道の軌道長半径と、前記予測演算部によって算出された前記予測の楕円軌道の近地点高度と目標の近地点高度との差とに基づいて前記多段式ロケットの新たな推力方向の角度を算出し、算出した前記多段式ロケットの新たな推力方向の角度を前記多段式ロケットの推力方向を制御するロケット制御装置に出力することによって前記目標の近地点高度に位置する近地点を通って前記予測の楕円軌道より遠地点高度が大きい楕円軌道へ前記多段式ロケットを誘導する誘導演算部として前記多段式ロケット誘導装置を機能させる。前記予測演算部は、前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、次段のロケットエンジンが固体燃料を燃焼し始める予定時刻として予め記憶する次段燃焼開始時刻と、に基づいて前記次段燃焼開始時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記次段燃焼開始時刻の予測位置ベクトルと予測速度ベクトルと、前記次段燃焼終了時刻と、前記最終段燃焼終了時刻と、前記次段燃焼終了時刻までの推力加速度の予測値を示す次段推力加速度プロファイルと、前記最終段推力加速度プロファイルと、に基づいて前記予測の楕円軌道の前記軌道長半径と前記近地点高度とを積分計算によって算出する。前記予測演算部は、前記ロケット制御装置が前記多段式ロケットの推力方向の角度を前記新たな推力方向の角度に変更した後、前記航法装置から新たな現在の位置ベクトルと新たな現在の速度ベクトルと新たな現在の推力加速度ベクトルとを入力し、前記新たな現在の位置ベクトルと前記新たな現在の速度ベクトルと新たな現在時刻と前記現段燃焼終了時刻と前記現段推力加速度プロファイルとを用いて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを新たに算出する。本発明の多段式ロケット誘導プログラムは、固体燃料を燃焼して推力を得るロケットエンジンを複数段搭載した多段式ロケットを誘導する多段式ロケット誘導装置を機能させる。前記多段式ロケット誘導プログラムは、前記多段式ロケットの位置ベクトルと速度ベクトルと推力加速度ベクトルとを計測する航法装置から現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルとを入力し、入力した現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルと、現在時刻と、推力を得ている現段のロケットエンジンが現段の固体燃料を燃焼し終わる予定時刻として予め記憶する現段燃焼終了時刻と、前記現段燃焼終了時刻までの推力加速度の予測値を示す現段推力加速度プロファイルとに基づいて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、最終段のロケットエンジンが最終段の固体燃料を燃焼し終わる予定時刻として予め記憶する最終段燃焼終了時刻と、前記最終段燃焼終了時刻までの推力加速度の予測値を示す最終段推力加速度プロファイルとに基づいて、前記最終段燃焼終了時刻に前記多段式ロケットが飛行する予測の楕円軌道の軌道長半径と遠地点高度とを積分計算によって算出する予測演算部と、前記予測演算部によって算出された前記予測の楕円軌道の軌道長半径と、前記予測演算部によって算出された前記予測の楕円軌道の遠地点高度と目標の遠地点高度との差とに基づいて前記多段式ロケットの新たな推力方向の角度を算出し、算出した前記多段式ロケットの新たな推力方向の角度を前記多段式ロケットの推力方向を制御するロケット制御装置に出力することによって前記目標の遠地点高度に位置する遠地点を通って前記予測の楕円軌道より近地点高度が大きい楕円軌道へ前記多段式ロケットを誘導する誘導演算部として前記多段式ロケット誘導装置を機能させる。本発明の多段式ロケット誘導方法は、固体燃料を燃焼して推力を得るロケットエンジンを複数段搭載した多段式ロケットを誘導する多段式ロケット誘導装置によって実行する。前記多段式ロケット誘導方法は、予測演算部が、前記多段式ロケットの位置ベクトルと速度ベクトルと推力加速度ベクトルとを計測する航法装置から現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルとを入力し、入力した現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルと、現在時刻と、推力を得ている現段のロケットエンジンが現段の固体燃料を燃焼し終わる予定時刻として予め記憶する現段燃焼終了時刻と、前記現段燃焼終了時刻までの推力加速度の予測値を示す現段推力加速度プロファイルとに基づいて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、最終段のロケットエンジンが最終段の固体燃料を燃焼し終わる予定時刻として予め記憶する最終段燃焼終了時刻と、前記最終段燃焼終了時刻までの推力加速度の予測値を示す最終段推力加速度プロファイルとに基づいて、前記最終段燃焼終了時刻に前記多段式ロケットが飛行する予測の楕円軌道の軌道長半径と近地点高度とを積分計算によって算出し、誘導演算部が、前記予測演算部によって算出された前記予測の楕円軌道の軌道長半径と、前記予測演算部によって算出された前記予測の楕円軌道の近地点高度と目標の近地点高度との差とに基づいて前記多段式ロケットの新たな推力方向の角度を算出し、算出した前記多段式ロケットの新たな推力方向の角度を前記多段式ロケットの推力方向を制御するロケット制御装置に出力することによって前記目標の近地点高度に位置する近地点を通って前記予測の楕円軌道より遠地点高度が大きい楕円軌道へ前記多段式ロケットを誘導する。本発明の多段式ロケット誘導方法は、固体燃料を燃焼して推力を得るロケットエンジンを複数段搭載した多段式ロケットを誘導する多段式ロケット誘導装置によって実行する。前記多段式ロケット誘導方法は、予測演算部が、前記多段式ロケットの位置ベクトルと速度ベクトルと推力加速度ベクトルとを計測する航法装置から現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルとを入力し、入力した現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルと、現在時刻と、推力を得ている現段のロケットエンジンが現段の固体燃料を燃焼し終わる予定時刻として予め記憶する現段燃焼終了時刻と、前記現段燃焼終了時刻までの推力加速度の予測値を示す現段推力加速度プロファイルとに基づいて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、最終段のロケットエンジンが最終段の固体燃料を燃焼し終わる予定時刻として予め記憶する最終段燃焼終了時刻と、前記最終段燃焼終了時刻までの推力加速度の予測値を示す最終段推力加速度プロファイルとに基づいて、前記最終段燃焼終了時刻に前記多段式ロケットが飛行する予測の楕円軌道の軌道長半径と遠地点高度とを積分計算によって算出し、誘導演算部が、前記予測演算部によって算出された前記予測の楕円軌道の軌道長半径と、前記予測演算部によって算出された前記予測の楕円軌道の遠地点高度と目標の遠地点高度との差とに基づいて前記多段式ロケットの新たな推力方向の角度を算出し、算出した前記多段式ロケットの新たな推力方向の角度を前記多段式ロケットの推力方向を制御するロケット制御装置に出力することによって前記目標の遠地点高度に位置する遠地点を通って前記予測の楕円軌道より近地点高度が大きい楕円軌道へ前記多段式ロケットを誘導する。本発明の多段式ロケット誘導システムは、固体燃料を燃焼して推力を得るロケットエンジンを複数段搭載した多段式ロケットの位置ベクトルと速度ベクトルと推力加速度ベクトルとを計測する航法装置と、前記多段式ロケットの推力方向を制御するロケット制御装置と、前記多段式ロケットを誘導する多段式ロケット誘導装置と、を備える。前記多段式ロケット誘導装置は、前記航法装置から現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルとを入力し、入力した現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルと、現在時刻と、推力を得ている現段のロケットエンジンが現段の固体燃料を燃焼し終わる予定時刻として予め記憶する現段燃焼終了時刻と、前記現段燃焼終了時刻までの推力加速度の予測値を示す現段推力加速度プロファイルとに基づいて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、最終段のロケットエンジンが最終段の固体燃料を燃焼し終わる予定時刻として予め記憶する最終段燃焼終了時刻と、前記最終段燃焼終了時刻までの推力加速度の予測値を示す最終段推力加速度プロファイルとに基づいて、前記最終段燃焼終了時刻に前記多段式ロケットが飛行する予測の楕円軌道の軌道長半径と近地点高度とを積分計算によって算出する予測演算部と、前記予測演算部によって算出された前記予測の楕円軌道の軌道長半径と、前記予測演算部によって算出された前記予測の楕円軌道の近地点高度と目標の近地点高度との差とに基づいて前記多段式ロケットの新たな推力方向の角度を算出し、算出した前記多段式ロケットの新たな推力方向の角度を前記ロケット制御装置に出力することによって前記目標の近地点高度に位置する近地点を通って前記予測の楕円軌道より遠地点高度が大きい楕円軌道へ前記多段式ロケットを誘導する誘導演算部とを備える。本発明の多段式ロケット誘導システムは、固体燃料を燃焼して推力を得るロケットエンジンを複数段搭載した多段式ロケットの位置ベクトルと速度ベクトルと推力加速度ベクトルとを計測する航法装置と、前記多段式ロケットの推力方向を制御するロケット制御装置と、前記多段式ロケットを誘導する多段式ロケット誘導装置と、を備える。前記多段式ロケット誘導装置は、前記航法装置から現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルとを入力し、入力した現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルと、現在時刻と、推力を得ている現段のロケットエンジンが現段の固体燃料を燃焼し終わる予定時刻として予め記憶する現段燃焼終了時刻と前記現段燃焼終了時刻までの推力加速度の予測値を示す現段推力加速度プロファイルと、に基づいて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、最終段のロケットエンジンが最終段の固体燃料を燃焼し終わる予定時刻として予め記憶する最終段燃焼終了時刻と、前記最終段燃焼終了時刻までの推力加速度の予測値を示す最終段推力加速度プロファイルとに基づいて、前記最終段燃焼終了時刻に前記多段式ロケットが飛行する予測の楕円軌道の軌道長半径と遠地点高度とを積分計算によって算出する予測演算部と、前記予測演算部によって算出された前記予測の楕円軌道の軌道長半径と、前記予測演算部によって算出された前記予測の楕円軌道の遠地点高度と目標の遠地点高度との差とに基づいて前記多段式ロケットの新たな推力方向の角度を算出し、算出した前記多段式ロケットの新たな推力方向の角度を前記ロケット制御装置に出力することによって前記目標の遠地点高度に位置する遠地点を通って前記予測の楕円軌道より近地点高度が大きい楕円軌道へ前記多段式ロケットを誘導する誘導演算部とを備える。本発明によれば、例えば、空中発射システムまたは中間の目標地点を要しない場合の地上発射システムで発射する固体ロケットをできるだけ遠地点高度の大きな楕円軌道、もしくは、できるだけ近地点高度の大きな楕円軌道に誘導することができる。実施の形態1.多段式ロケットをできるだけ遠地点高度の大きな楕円軌道、もしくは、できるだけ近地点高度の大きな楕円軌道に誘導するシステムの形態について説明する。図1は、実施の形態1における多段式ロケット誘導システム100の構成図である。実施の形態1における多段式ロケット誘導システム100の構成について、図1に基づいて説明する。多段式ロケット誘導システム100は、固体燃料を燃焼して推力を得るロケットエンジンを複数段搭載した多段式ロケットを誘導するシステムである。例えば、多段式ロケット誘導システム100は多段式ロケットに搭載される。但し、多段式ロケット誘導システム100の構成の少なくとも一部を多段式ロケットの外部例えば、地上の施設に設けても構わない。多段式ロケット誘導システム100は、航法装置110と、ロケット制御装置120と、ロケット誘導装置200とを備える。航法装置110は、多段式ロケットの位置ベクトルと速度ベクトルと加速度ベクトル推力加速度ベクトルとを航法演算によって計測する装置である。例えば、航法装置110は慣性計測装置IMUとも呼ばれ、ジャイロと加速度計とを備える。航法装置110は、ジャイロによって計測される3軸方向の角速度と加速度計によって計測される3軸方向の加速度加速度ベクトルとを積分して3軸方向の速度速度ベクトルの変化量と3軸方向の位置位置ベクトルの変化量とを算出する。航法装置110は、速度ベクトルの変化量を前回の速度ベクトルに加算して現在の速度ベクトルを算出し、位置ベクトルの変化量を前回の位置ベクトルに加算して現在の位置ベクトルを算出する。ロケット制御装置120は、多段式ロケットの推力方向ピッチ姿勢角を制御する装置である。例えば、ロケット制御装置120は、多段式ロケットの機体やロケットエンジンの噴射向きを制御して多段式ロケットの推力方向を制御する。ロケット誘導装置200多段式ロケット誘導装置の一例は、多段式ロケットを誘導するための演算を行う装置コンピュータである。ロケット誘導装置200は、推力加速度予測部210と、予測演算部220と、誘導演算部230と、誘導装置記憶部290とを備える。推力加速度予測部210は、航法装置110によって計測された推力加速度ベクトルに基づいて推力加速度プロファイルを更新する。推力加速度プロファイルは、ロケットエンジンの各時刻経過時間の推力加速度の予測値を示すデータである。例えば、推力加速度プロファイルは、各段のロケットエンジンが固体燃料を燃焼し始める予定時刻推力加速度がゼロから上昇する時刻や各段のロケットエンジンが固体燃料を燃焼し終わる予定時刻推力加速度がゼロになる時刻を示す。予測演算部220は、航法装置110から多段式ロケットの現在の位置ベクトルと現在の速度ベクトルとを入力し、誘導装置記憶部290から推力加速度プロファイルを取得する。予測演算部220は、現在の位置ベクトルと現在の速度ベクトルと、現在時刻と、推力を得ている現段のロケットエンジンが現段の固体燃料を燃焼し終わる現段燃焼終了時刻と、現段燃焼終了時刻までの推力加速度の予測値を示す現段推力加速度プロファイルと、に基づいて現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出する。予測演算部220は、算出した現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、現段燃焼終了時刻と、最終段のロケットエンジンが最終段の固体燃料を燃焼し終わる最終段燃焼終了時刻と、前記最終段燃焼終了時刻までの推力加速度の予測値を示す最終段推力加速度プロファイルと、に基づいて最終段燃焼終了時刻に多段式ロケットが飛行する予測の楕円軌道の軌道長半径軌道要素予測値の一例と近地点高度または遠地点高度軌道要素予測値の一例とを積分計算によって算出する。例えば、予測演算部220は、以下のように予測の楕円軌道の軌道長半径と近地点高度または遠地点高度とを算出する。予測演算部220は、現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、現段燃焼終了時刻と、次段のロケットエンジンが固体燃料を燃焼し始める次段燃焼開始時刻と、に基づいて次段燃焼開始時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出する。予測演算部220は、算出した次段燃焼開始時刻の予測位置ベクトルと予測速度ベクトルと、次段燃焼終了時刻と、最終段燃焼終了時刻と、次段燃焼終了時刻までの推力加速度の予測値を示す次段推力加速度プロファイルと、最終段推力加速度プロファイルと、に基づいて予測の楕円軌道の軌道長半径と近地点高度または遠地点高度とを積分計算によって算出する。予測演算部220は、ロケット制御装置120が多段式ロケットの推力方向の角度を新たな推力方向の角度に変更した後、航法装置110から新たな現在の位置ベクトルと新たな現在の速度ベクトルとを入力する。予測演算部220は、新たな現在の位置ベクトルと新たな現在の速度ベクトルと新たな現在時刻と現段燃焼終了時刻とを用いて現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを新たに算出し、最終段燃焼終了時刻に多段式ロケットが飛行する予測の楕円軌道の軌道長半径と近地点高度または遠地点高度とを新たに算出する。誘導演算部230は、予測演算部220によって算出された予測の楕円軌道の軌道長半径と、予測演算部220によって算出された予測の楕円軌道の近地点高度と目標の近地点高度軌道要素目標値の一例との差とに基づいて多段式ロケットの新たな推力方向の角度を算出する。誘導演算部230は、算出した多段式ロケットの新たな推力方向の角度をロケット制御装置120に出力することによって、目標の近地点高度に位置する近地点を通って予測の楕円軌道より遠地点高度が大きい楕円軌道へ多段式ロケットを誘導する。または、誘導演算部230は、予測演算部220によって算出された予測の楕円軌道の軌道長半径と、予測演算部220によって算出された予測の楕円軌道の遠地点高度と目標の遠地点高度軌道要素目標値の一例との差とに基づいて多段式ロケットの新たな推力方向の角度を算出する。誘導演算部230は、算出した多段式ロケットの新たな推力方向の角度をロケット制御装置120に出力することによって、目標の遠地点高度に位置する遠地点を通って予測の楕円軌道より近地点高度が大きい楕円軌道へ多段式ロケットを誘導する。誘導装置記憶部290は、ロケット誘導装置200で使用するデータを記憶する。推力加速度プロファイル、位置ベクトル、速度ベクトル、推力加速度ベクトル、軌道要素予測値、軌道要素目標値およびピッチ姿勢角は、誘導装置記憶部290に記憶するデータの一例である。以下、多段式ロケット誘導システム100が3段式ロケット多段式ロケットの一例を誘導する場合について説明する。但し、多段式ロケット誘導システム100が誘導する多段式ロケットは2段式ロケットであっても4段以上の多段式ロケットであっても構わない。図2は、実施の形態1における多段式ロケット誘導システム100による3段式ロケット誘導処理を示すフローチャートである。実施の形態1における多段式ロケット誘導システム100による3段式ロケット誘導処理多段式ロケット誘導方法、多段式ロケット誘導プログラムの一例について、図2に基づいて説明する。以下、3段式ロケット多段式ロケットの一例を単にロケットという。以下、第n段n:1以上の整数のロケットエンジンが固体燃料を燃焼しているときに行うロケットの誘導をn段誘導という。また、前段のロケットエンジンが固体燃料を燃焼し終わって切り離されてから次段のロケットエンジンが固体燃料を燃焼し始めるまでのロケットの誘導、特に、2段誘導の終了時から3段誘導の開始時までのロケットの誘導をコースト誘導という。図2に示す3段式ロケット誘導処理では2段誘導以降の処理を説明する。図示省略する1段誘導1段誘導の終了時から2段誘導の開始時までのコースト誘導を含むにおいて、多段式ロケット誘導システム100は、所定の誘導処理例えば、従来の誘導処理によってロケットのヨー姿勢角方位角を調整してロケットを目標の軌道面に誘導する。目標の軌道面は、目標の楕円軌道遠地点高度または近地点高度が大きい楕円軌道を含んだ平面である。多段式ロケット誘導システム100は、図2に示す2段誘導以降の処理によってロケットのピッチ姿勢角を調整してロケットを目標の楕円軌道に誘導する。多段式ロケット誘導システム100は、図2に示す3段式ロケット誘導処理を3段の燃焼が終了するまで周期的例えば、1または2秒毎に実行する。つまり、多段式ロケット誘導システム100は、ロケットのピッチ姿勢角推力方向を周期的に正しい角度に調整する。S110において、予測演算部220は、現在の誘導が2段誘導とコースト誘導と3段誘導とのいずれであるかを判定する。現在の誘導がいずれの誘導であるかを示すデータはロケットエンジンの切り離し時や固体燃料の燃焼開始時に誘導装置記憶部290に記憶されるものとする。現在の誘導が2段誘導である場合、S120に進む。現在の誘導がコースト誘導である場合、ロケットのピッチ姿勢角を変更せず、今回の周期の3段式ロケット誘導処理を終了する。現在の誘導が3段誘導である場合、S140に進む。S120において、予測演算部220は、後述する2段誘導予測演算処理によって最終段3段の燃焼終了時の軌道要素予測値軌道長半径、近地点高度、遠地点高度、近地点引数などの予測値を算出する。S120の後、S150に進む。S140において、予測演算部220は、後述する3段誘導予測演算処理によって最終段3段の燃焼終了時の軌道要素予測値を算出する。S140の後、S150に進む。S150において、誘導演算部230は誘導オプションの種類を判定する。誘導オプションの種類は誘導装置記憶部290に予め記憶する。誘導オプションの種類は、オプションA投入誤差の最小化と、オプションB遠地点高度の最大化と、オプションC近地点高度の最大化との3種類である。オプションA投入誤差の最小化は、目標の遠地点と目標の近地点とを通る楕円軌道目標軌道へロケットを誘導するための誘導オプションである。オプションB遠地点高度の最大化は、目標の近地点を通って遠地点高度ができるだけ大きい楕円軌道へロケットを誘導するための誘導オプションである。オプションC近地点高度の最大化は、目標の遠地点を通って近地点高度ができるだけ大きい楕円軌道へロケットを誘導するための誘導オプションである。誘導オプションの種類がオプションA投入誤差の最小化である場合、S151に進む。誘導オプションの種類がオプションB遠地点高度の最小化である場合、S152に進む。誘導オプションの種類がオプションC近地点高度の最小化である場合、S153に進む。S151において、誘導演算部230は、オプションA投入誤差の最小化のΔyを算出する。オプションAのΔyは、遠地点高度の目標値軌道要素目標値の一例と予測値との値差を用いた成分と、近地点高度の目標値軌道要素目標値の一例と予測値との値差を用いた成分と、を含んだベクトルである。オプションAのΔyの詳細については別途説明する。S151の後、S154に進む。S152において、誘導演算部230は、オプションB遠地点高度の最大化のΔyを算出する。オプションBのΔyは、軌道長半径の予測値の逆数を用いた成分と、近地点高度の目標値と予測値との値差を用いた成分と、を含んだベクトルである。オプションBのΔyの詳細については別途説明する。S152の後、S154に進む。S153において、誘導演算部230は、オプションC近地点高度の最大化のΔyを算出する。オプションCのΔyは、軌道長半径の予測値の逆数を用いた成分と、遠地点高度の目標値と予測値との値差を用いた成分と、を含んだベクトルである。オプションCのΔyの詳細については別途説明する。S153の後、S154に進む。S151からS153で用いる軌道要素目標値は誘導装置記憶部290に予め記憶する。図3は、実施の形態1における誘導オプションとΔyとの関係を示す表である。誘導演算部230は、誘導オプションの種類S150と現在の誘導S110との組み合わせに基づいて、図3に示すΔy更新ベクトルの計算式を計算してΔyを算出するS151からS153。図3に示す記号の意味は以下の通りである。Reは赤道半径である。aは軌道長半径である。rpは近地点高度の予測値であり、rpfは近地点高度の目標値である。raは遠地点高度の予測値であり、rafは遠地点高度の目標値である。ωは近地点引数の予測値であり、ωfは近地点引数の目標値である。例えば、近地点高度と遠地点高度との少なくともいずれかの目標値に対する予測値の誤差を最小化したい場合、誘導演算部230は以下のようにΔyを算出する。誘導オプションがオプションA投入誤差の最小化である場合S151、誘導演算部230は図3の項番1に示すΔyを算出する。2段誘導と3段誘導とでΔyの計算式は同じである。誘導オプションがオプションB遠地点高度の最大化である場合S152、誘導演算部230は図3の項番3に示すΔyを算出する。2段誘導と3段誘導とでΔyの計算式は同じである。誘導オプションがオプションC近地点高度の最大化である場合S153、誘導演算部230は図3の項番5に示すΔyを算出する。2段誘導と3段誘導とでΔyの計算式は異なる。但し、誘導演算部230は、近地点高度または遠地点高度以外の軌道要素の誤差を最小化しても良い。例えば、近地点高度または遠地点高度以外に近地点引数軌道要素の一例の誤差を最小化したい場合、誘導演算部230は以下のようにΔyを算出する。誘導オプションがオプションA投入誤差の最小化である場合S151、誘導演算部230は図3の項番2に示すΔyを算出する。2段誘導と3段誘導とでΔyの計算式は同じである。誘導オプションがオプションB遠地点高度の最大化である場合S152、誘導演算部230は図3の項番4に示すΔyを算出する。2段誘導と3段誘導とでΔyの計算式は同じである。誘導オプションがオプションC近地点高度の最大化である場合S153、誘導演算部230は図3の項番6に示すΔyを算出する。2段誘導と3段誘導とでΔyの計算式は同じである。図2に戻り、S154から説明を続ける。S154において、誘導演算部230は、S151からS153のいずれかで算出したΔyの二乗和Δy2を算出する。二乗和Δy2はΔyに含まれる各成分を二乗して合計した値である。誘導演算部230は、二乗和Δy2と誘導装置記憶部290に予め記憶する誘導閾値とを比較する。二乗和Δy2が誘導閾値未満である場合、ロケットは目標の軌道に向けて正しく誘導されているため、今回の周期の3段式ロケット誘導処理を終了する。二乗和Δy2が誘導閾値以上である場合、S160に進む。S160において、誘導演算部230は、S151からS153のいずれかで算出したΔyに基づいてロケットの新たな独立変数を算出する。独立変数とは、例えばピッチ姿勢角、後述するステアリング係数、最終段燃焼開始時刻である。誘導演算部230は、算出したピッチ姿勢角をロケット制御装置120に出力する。例えば、誘導演算部230は、以下のようにロケットの新たな独立変数を算出する。誘導演算部230は、以下の関係式1を用いて式2に含まれる独立変数の更新量を表す更新量ベクトルΔCを算出する。関係式1において感度行列Jは、軌道要素予測値を独立変数で偏微分して得られるヤコビ行列である。係数λは、繰り返し計算においてΔy2が前回値より小さい場合に小さく、Δy2が前回値より大きい場合に大きく設定する値である。重み付け行列Wは、誘導オプションの種類に応じて予め定義する行列である。誘導演算部230は、関係式1を用いて算出した更新量ベクトルΔCが表す独立変数の更新量を現在の独立変数に加算して新たな独立変数を算出する。誘導演算部230は、新たな独立変数のうちステアリング係数を式2に設定し、時間変数tに現在時刻から現段燃焼終了時刻までの時間および3段燃焼開始時刻から3段燃焼終了時刻までの時間を代入して式2を計算する。式2は、ロケットの新たなピッチ姿勢角θのタンジェント値tanθを計算するための時間変数tの2次式である。式2を2次タンジェント則という。誘導演算部230は、式2の代わりに式3を用いても構わない。式3は、ロケットの新たなピッチ姿勢角θのタンジェント値tanθを計算するための時間変数tの1次式である。式3をリニアタンジェント則という。例えば、誘導演算部230は、2段誘導と3段誘導との両方で2次タンジェント則を用いる。但し、誘導演算部230は、2段誘導で2次タンジェント則を用いて、3段誘導でリニアタンジェント則を用いてもよい。また、誘導演算部230は、各段の誘導時に2次タンジェント則とリニアタンジェント則とのいずれを選択しても構わない。S160の後、S110に戻り、軌道要素の目標値と予測値との二乗和Δy2が誘導閾値未満に収束するまでS110からS160の処理を繰り返す最小二乗法、レーベンバーグマルカート法。図2に示した3段式ロケット誘導処理を3段の燃焼が終了するまで周期的に繰り返すことにより、ロケットを目標の軌道に誘導することができる。図4は、実施の形態1における2段誘導演算処理S120を示すフローチャートである。実施の形態1における2段誘導演算処理S120について、図4に基づいて説明する。S121において、推力加速度予測部210は、航法装置110からロケットの現在の加速度ベクトル推力加速度ベクトルを取得し、取得した現在の加速度ベクトルに基づいて所定の更新処理によって2段推力加速度プロファイルを更新する。また、予測演算部220は、航法装置110からロケットの現在の位置ベクトルと速度ベクトルとを入力する。さらに、予測演算部220は、2段の固体燃料の燃焼が終了する2段燃焼終了時刻と、現在時刻から2段燃焼終了時刻までの各時刻の推力加速度とを推力加速度プロファイルから取得する。予測演算部220は、現在の位置ベクトルと、現在の速度ベクトルと、現在時刻と、2段燃焼終了時刻と、現在時刻から2段燃焼終了時刻までの各時刻の推力加速度と、を入力にして所定の予測演算関数を演算し、2段燃焼終了時刻の位置ベクトルと速度ベクトルとを算出する。予測演算関数の演算において、予測演算部220は、現在時刻から2段燃焼終了時刻まで推力加速度の積分計算を行い、現在時刻から2段燃焼終了時刻までの速度ベクトルと位置ベクトルとの変位量を算出する。さらに、予測演算部220は、速度ベクトルの変位量を現在の速度ベクトルに加算して2段燃焼終了時刻の速度ベクトルを算出し、位置ベクトルの変位量を現在の位置ベクトルに加算して2段燃焼終了時刻の位置ベクトルを算出する。S121の後、S122に進む。S122において、2段燃焼終了時刻から3段燃焼開始時刻までの間はいずれのロケットエンジンも稼働していないため、各時刻の推力加速度はゼロである。予測演算部220は、2段燃焼終了時刻の位置ベクトルと、2段燃焼終了時刻の速度ベクトルと、現在時刻と、3段燃焼開始時刻と、を入力にして予測演算関数を演算する。予測演算関数の演算において、予測演算部220は、S121と同様に積分計算によって3段燃焼開始時刻の位置ベクトルと速度ベクトルとを算出する。S122の後、S123に進む。S123において、予測演算部220は、3段の固体燃料の燃焼が終了する3段燃焼終了時刻と、3段燃焼開始時刻から3段燃焼終了時刻までの各時刻の推力加速度とを推力加速度プロファイルから取得する。予測演算部220は、3段燃焼開始時刻の位置ベクトルと、3段燃焼開始時刻の速度ベクトルと、現在時刻と、3段燃焼終了時刻と、3段燃焼開始時刻から3段燃焼終了時刻までの各時刻の推力加速度と、を入力にして予測演算関数を演算する。予測演算関数の演算において、予測演算部220は、S121と同様に積分計算によって3段燃焼終了時刻の位置ベクトルと速度ベクトルとを算出する。さらに、予測演算部220は、3段燃焼終了時刻の位置ベクトルと速度ベクトルとに基づいて軌道要素の予測値軌道要素予測値を算出する。軌道要素および軌道要素予測値の算出方法については別途説明する。S123により、2段誘導演算処理S120は終了する。図5は、実施の形態1における3段誘導演算処理S140を示すフローチャートである。実施の形態1における3段誘導演算処理S140について、図5に基づいて説明する。S141において、推力加速度予測部210は、航法装置110からロケットの現在の加速度ベクトル推力加速度ベクトルを取得し、取得した現在の加速度ベクトルに基づいて所定の更新処理によって3段推力加速度プロファイルを更新する。また、予測演算部220は、航法装置110からロケットの現在の位置ベクトルと速度ベクトルとを入力する。さらに、予測演算部220は、3段の固体燃料の燃焼が終了する3段燃焼終了時刻と、現在時刻から3段燃焼終了時刻までの各時刻の推力加速度とを推力加速度プロファイルから取得する。予測演算部220は、現在の位置ベクトルと、現在の速度ベクトルと、現在時刻と、3段燃焼終了時刻と、現在時刻から3段燃焼終了時刻までの各時刻の推力加速度と、を入力にして予測演算関数を演算する。予測演算関数の演算において、予測演算部220は、図3のS121と同様に積分計算によって3段燃焼終了時刻の位置ベクトルと速度ベクトルとを算出する。さらに、予測演算部220は、3段燃焼終了時刻の位置ベクトルと速度ベクトルとに基づいて軌道要素予測値を算出する。軌道要素予測値の算出方法については別途説明する。S141により、3段誘導演算処理S140は終了する。図6、図7は、実施の形態1における軌道要素の説明図である。実施の形態1における軌道要素について、図6、図7に基づいて説明する。軌道要素は、軌道の大きさを表す項目である。軌道要素として、近地点高度、遠地点高度、軌道長半径、軌道傾斜角、近地点引数、真近点離角、などが挙げられる。これらの軌道要素について以下に説明する。目標軌道は、地球の重心を焦点とする楕円軌道である。目標軌道内でのロケットの進行方向を矢印で記す。近地点は目標軌道内で地球重心から最も近い点であり、地球重心から近地点までの距離を近地点半径という。近地点高度rpは近地点半径から地球の赤道半径を引いた距離である。遠地点は目標軌道内で地球重心から最も遠い点であり、地球重心から遠地点までの距離を遠地点半径という。遠地点高度raは遠地点半径から地球の赤道半径を引いた距離である。軌道長半径aは目標軌道の中心oから遠地点までの距離である目標軌道の中心oから近地点までの距離と等しい。軌道傾斜角iは目標軌道と地球の赤道面とが成す角度である。昇交点は目標軌道と地球の赤道面とが交わる点地球の南半球から北半球への通過点である。近地点引数ωは昇交点から近地点までの角度ロケットの進行方向側の角度である。真近点離角fは近地点からロケットの位置までの角度ロケットの進行方向側の角度である。また、目標軌道の中心oから地球重心までの距離aeの軌道長半径aに対する比率を離心率eといい、昇交点からロケットの位置までの角度ロケットの進行方向側の角度を緯度引数lという。以下に、位置ベクトルrと速度ベクトルvとに基づいて各軌道要素の予測値軌道要素予測値を算出する算出式を以下に示す。図8は、実施の形態1における多段式ロケット誘導システム100で用いる座標系の一例を示す図である。図8に示すように、例えば、地球中心または地球重心を座標系の原点oとし、原点oからグリニッジ子午線と赤道との交点を通る方向をX軸方向とし、北極方向北方向をZ軸方向とし、Z軸とX軸とに直交し右手系を成す方向をY軸方向とする。図9は、実施の形態1におけるロケットのピッチ姿勢角θを定義する図である。図9に示すように、ロケットのピッチ角θは接線方向と角運動量方向とを含んだ平面に対して推力方向が成す角度である。ロケットのピッチ姿勢角θが表す方向は、ロケットの推力方向とみなすことができる。接線方向は位置ベクトル×速度ベクトル×位置ベクトル×は外積を表すで定義される方向である。軌道半径方向は位置ベクトルの向きに相当する方向である。角運動量方向は接線方向と軌道半径方向とに対して垂直な方向であり、位置ベクトル×速度ベクトル×は外積を表すで定義される。次に、図2から図5で説明した多段式ロケットの誘導アルゴリズムをシミュレーションした結果について説明する。図10は、実施の形態1における2段誘導をシミュレーションした際のパラメータを示す表である。図11は、実施の形態1における3段誘導をシミュレーションした際のパラメータを示す表である。図12は、実施の形態1におけるシミュレーションパターンを示す表である。図13は、実施の形態1における誘導アルゴリズム遠地点高度の最大化のシミュレーション結果を示すグラフである。図14は、実施の形態1における誘導アルゴリズム近地点高度の最大化のシミュレーション結果を示すグラフである。図10、図11においてシミュレーションで使用したパラメータを○で記す。図10、図11に示す各記号の意味は以下の通りである。Aは実施の形態1における誘導アルゴリズム中間目標値を設定せずに全段を統合したロジックを意味する。Qは2次タンジェント則を意味する。Lはリニアタンジェント則を意味する。CXはステアリング係数を意味する。θは3段燃焼中時刻のロケットのピッチ姿勢角を示す。Tは2段燃焼終了時刻から3段燃焼開始時刻までの時間を示す。rpは近地点高度を示す。raは遠地点高度を示す。aは軌道長半径を示す。ωは近地点引数を示す。図10、図11において、Aの従属変数rp,ra,a,ωは3段燃焼終了時刻の値である。図12に示す9パターンのシミュレーション結果のRSSRootSummedSquare値を図13、14に示す。500×250キロメートルの楕円軌道遠地点高度:500km、近地点高度:250kmを当初の予定軌道として誘導オプションB遠地点高度の最大化を実行した結果、標準偏差1σと標準偏差ノミナルと標準偏差-1σのRSS値はいずれも遠地点高度を200キロメートル以上大きくすることができることを示した図13参照。同様に、誘導オプションC近地点高度の最大化を実行した結果、標準偏差1σと標準偏差ノミナルと標準偏差-1σのRSS値はいずれも近地点高度を75キロメートル以上大きくすることができることを示した図14参照。図15は、実施の形態1における誘導オプションの一例を示す表である。実施の形態1において、例えば、図15に示す誘導オプションに従ってロケットを目標軌道へ誘導する誘導アルゴリズムについて説明した。図15において、○は軌道要素近地点高度、遠地点高度、近地点引数の目標値に対する予測値の誤差を最小にすることを意味し、maxは軌道要素近地点高度、遠地点高度をできるだけ大きくすることを意味する。図16は、実施の形態1におけるロケット誘導装置200のハードウェア資源の一例を示す図である。図16において、ロケット誘導装置200は、CPU911CentralProcessingUnitを備えている。CPU911は、バス912を介してROM913、RAM914、通信ボード915と接続され、これらのハードウェアデバイスを制御する。通信ボード915は、有線または無線でネットワークに接続している。ROMまたはRAMには、OSオペレーティングシステム、プログラム群、ファイル群が記憶されている。プログラム群には、実施の形態において~部として説明する機能を実行するプログラムが含まれる。プログラム例えば、多段式ロケット誘導プログラムは、CPU911により読み出され実行される。すなわち、プログラムは、~部としてコンピュータを機能させるものであり、また~部の手順や方法をコンピュータに実行させるものである。ファイル群には、実施の形態において説明する~部で使用される各種データ入力、出力、判定結果、計算結果、処理結果などが含まれる。実施の形態において構成図およびフローチャートに含まれている矢印は主としてデータや信号の入出力を示す。フローチャートなどに基づいて説明する実施の形態の処理はCPU911やその他のハードウェアを用いて実行される。実施の形態において~部として説明するものは~回路、~装置、~機器であってもよく、また~ステップ、~手順、~処理であってもよい。すなわち、~部として説明するものは、ファームウェア、ソフトウェア、ハードウェアまたはこれらの組み合わせのいずれで実装されても構わない。実施の形態1において、多段式ロケットを目標軌道近地点高度または遠地点高度ができるだけ大きい軌道に誘導するシステムについて説明した。実施の形態1により、固体燃料を用いる多段式ロケットを中間の目標地点を経由せずに高い精度で目標軌道に誘導することができる。また、固体燃料を用いることにより発射整備作業の時間を短縮し、中間の目標地点を設定しないことによりミッション解析期間を短縮することができる。つまり、実施の形態1により、緊急の打ち上げ要請にもタイムリーに応えられるという即応性を有した上で、多段式ロケットを目標軌道に高い精度で誘導することができる。また、近地点高度または遠地点高度を大きくすることができることにより、ユーザの要求投入条件に柔軟に対応することができる。例えば、ペイロード人工衛星、探査機などを搭載したロケットをペイロードの目標軌道またはそれに近い軌道へ投入することができる。これにより、ペイロードの燃料を節約することができ、余った燃料をペイロードの軌道修正などに利用することもできる。実施の形態1では3段式ロケットを例に挙げて説明したが、2段式ロケットまたは4段以上の多段式ロケットを誘導しても構わない。誘導方法は3段式ロケットと同様である。実施の形態1では2段誘導以降でロケットのピッチ姿勢角推力方向を調整するアルゴリズムを例に挙げて説明したが、1段誘導においてもロケットのピッチ姿勢角を調整しても構わない。調整方法は2段誘導と同様である。"