(57)【要約】【課題】打具のスイングに含まれる個々の挙動に基づいて、新たな挙動を推定する。【解決手段】スイング解析装置は、プレイヤーによる打具のスイング時の複数の注目点の挙動を時系列に表す時系列挙動データを取得するデータ取得部と、前記時系列挙動データを特異値分解することにより、前記挙動を分解するモード展開部とを備える。前記モード展開部は、前記時系列挙動データを特異値分解することにより、前記挙動から複数のモードにそれぞれ対応する複数の特徴的な挙動を抽出し、前記複数のモードのうち、少なくとも１つの注目モードに対応する前記特徴的な挙動を構成する、前記複数の注目点の挙動の励起パターンを表す基底である第１基底を修正し、前記修正後の第１基底に基づいて、前記注目モードに対応する前記特徴的な挙動を再構成する。【選択図】図３
2020185370.tif 000002
【特許請求の範囲】
【請求項１】
プレイヤーによる打具のスイング時の複数の注目点の挙動を時系列に表す時系列挙動データを取得するデータ取得部と、
前記時系列挙動データを特異値分解することにより、前記挙動を分解するモード展開部と
を備え、
前記モード展開部は、
前記時系列挙動データを特異値分解することにより、前記挙動から複数のモードにそれぞれ対応する複数の特徴的な挙動を抽出し、
前記複数のモードのうち、少なくとも１つの注目モードに対応する前記特徴的な挙動を構成する、前記複数の注目点の挙動の励起パターンを表す基底である第１基底を修正し、
前記修正後の第１基底に基づいて、前記注目モードに対応する前記特徴的な挙動を再構成する、
スイング解析装置。
【請求項２】
前記モード展開部は、前記第１基底に含まれる前記複数の注目点の挙動にそれぞれ対応する複数の成分のうち、少なくとも１つの成分の値を修正する、
請求項１に記載のスイング解析装置。
【請求項３】
前記モード展開部は、ユーザからの入力に従って、前記複数のモードのうち修正の対象となる少なくとも１つの前記注目モード、前記複数の成分のうち値を修正する少なくとも１つの注目成分、及び、前記注目成分の修正量の少なくとも１つを決定する、
請求項２に記載のスイング解析装置。
【請求項４】
前記モード展開部は、
前記注目モードの前記第１基底を、前記複数のモードのうち前記注目モードとは別のペアモードの前記第１基底とともに、前記注目モードの前記第１基底及び前記ペアモードの前記第１基底に直交する軸周りで回転させることにより修正する、
請求項１に記載のスイング解析装置。
【請求項５】
前記モード展開部は、
前記複数のモードのうち、前記注目モードとは別の複数の候補モードの各々に対し、前記注目モードの前記第１基底を、前記注目モードの前記第１基底が目標値に近づくように、前記注目モードの前記第１基底及び前記候補モードの前記第１基底に直交する軸周りで回転させ、
前記回転後の注目モードの前記第１基底と前記目標値との相関係数が最大となる回転を与える前記候補モードを前記ペアモードとして選択する、
請求項４に記載のスイング解析装置。
【請求項６】
前記モード展開部は、
回転後の前記注目モードの前記第１基底とその目標値との誤差を表す項を含む誤差関数を最小化する回転角度を算出し、
前記注目モードの前記第１基底を、前記ペアモードの前記第１基底とともに、前記軸周りで前記回転角度だけ回転させることにより修正する、
請求項４に記載のスイング解析装置。
【請求項７】
前記誤差関数は、回転後の前記ペアモードの前記第１基底とその目標値との誤差を表す項をさらに含む、
請求項６に記載のスイング解析装置。
【請求項８】
前記時系列挙動データは、前記複数の注目点の位置、姿勢、速度、加速度、角速度、角加速度、力及びトルクの少なくとも１つを時系列に表すデータである、
請求項１から７のいずれかに記載のスイング解析装置。
【請求項９】
前記打具は、ゴルフクラブである、
請求項１から８のいずれかに記載のスイング解析装置。
【請求項１０】
プレイヤーによる打具のスイング時の複数の注目点の挙動を時系列に表す時系列挙動データを取得することと、
前記時系列挙動データを特異値分解することにより、前記挙動を分解することと
を含み、
前記挙動を分解することは、
前記時系列挙動データを特異値分解することにより、前記挙動から複数のモードにそれぞれ対応する複数の特徴的な挙動を抽出することと、
前記複数のモードのうち、少なくとも１つの注目モードに対応する前記特徴的な挙動を構成する、前記複数の注目点の挙動の励起パターンを表す基底である第１基底を修正することと、
前記修正後の第１基底に基づいて、前記注目モードに対応する前記特徴的な挙動を再構成することと
を含む、
スイング解析方法。
【請求項１１】
プレイヤーによる打具のスイング時の複数の注目点の挙動を時系列に表す時系列挙動データを取得することと、
前記時系列挙動データを特異値分解することにより、前記挙動を分解することと
をコンピュータに実行させ、
前記挙動を分解することは、
前記時系列挙動データを特異値分解することにより、前記挙動から複数のモードにそれぞれ対応する複数の特徴的な挙動を抽出することと、
前記複数のモードのうち、少なくとも１つの注目モードに対応する前記特徴的な挙動を構成する、前記複数の注目点の挙動の励起パターンを表す基底である第１基底を修正することと、
前記修正後の第１基底に基づいて、前記注目モードに対応する前記特徴的な挙動を再構成することと
を含む、
スイング解析プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ゴルフスイング等の打具のスイングを解析するスイング解析装置、方法及びプログラムに関する。
【背景技術】
【０００２】
プレイヤーによる打具のスイングは、様々な挙動が含まれる複合動作である。例えば、ゴルフスイングは、腕の振り上げ及び振り下げ、腰の回転、身体の平行移動等の様々な挙動を組み合わせた動作である。非特許文献１及び２は、ゴルフスイングを計測した計測データから得られる観測行列を特異値分解することにより、ゴルフスイングを複数の挙動（モード）に分解する手法を開示している。
【先行技術文献】
【非特許文献】
【０００３】
【非特許文献１】松本賢太，他５名，「クラブ設計を目的とした特異値分解によるゴルフスイングの動作分析」，設計工学，Ｖｏｌ．５３，Ｎｏ．６，ｐｐ．４４７−４６２，２０１８年６月
【非特許文献２】畑中崚志，他４名，「ゴルフスイングの動作分解に基づくクラブ特性影響評価」，Ｎｏ．１８−１５  日本機械学会  シンポジウム：スポーツ工学・ヒューマンダイナミクス２０１８講演論文集，Ｂ−１０，２０１８年１１月
【発明の概要】
【発明が解決しようとする課題】
【０００４】
非特許文献１及び２の手法によれば、ゴルフスイングを構成する様々な挙動を個別に評価することができる。しかしながら、このような試みは始まったばかりであり、解析手法のさらなる改良が望まれる。なお、ゴルフクラブ以外の打具、例えば、テニスラケットやベースボールバット等のスイングの解析にも、同様のことが言える。
【０００５】
本発明は、打具のスイングに含まれる個々の挙動に基づいて、新たな挙動を推定することができるスイング解析装置、方法及びプログラムを提供することを目的とする。
【課題を解決するための手段】
【０００６】
第１観点に係るスイング解析装置は、プレイヤーによる打具のスイング時の複数の注目点の挙動を時系列に表す時系列挙動データを取得するデータ取得部と、前記時系列挙動データを特異値分解することにより、前記挙動を分解するモード展開部とを備える。前記モード展開部は、前記時系列挙動データを特異値分解することにより、前記挙動から複数のモードにそれぞれ対応する複数の特徴的な挙動を抽出し、前記複数のモードのうち、少なくとも１つの注目モードに対応する前記特徴的な挙動を構成する、前記複数の注目点の挙動の励起パターンを表す基底である第１基底を修正し、前記修正後の第１基底に基づいて、前記注目モードに対応する前記特徴的な挙動を再構成する。
【０００７】
第２観点に係るスイング解析装置は、第１観点に係るスイング解析装置であって、前記モード展開部は、前記第１基底に含まれる前記複数の注目点の挙動にそれぞれ対応する複数の成分のうち、少なくとも１つの成分の値を修正する。
【０００８】
第３観点に係るスイング解析装置は、第２観点に係るスイング解析装置であって、前記モード展開部は、ユーザからの入力に従って、前記複数のモードのうち修正の対象となる少なくとも１つの前記注目モード、前記複数の成分のうち値を修正する少なくとも１つの注目成分、及び、前記注目成分の修正量の少なくとも１つを決定する。
【０００９】
第４観点に係るスイング解析装置は、第１観点に係るスイング解析装置であって、前記モード展開部は、前記注目モードの前記第１基底を、前記複数のモードのうち前記注目モードとは別のペアモードの前記第１基底とともに、前記注目モードの前記第１基底及び前記ペアモードの前記第１基底に直交する軸周りで回転させることにより修正する。
【００１０】
第５観点に係るスイング解析装置は、第４観点に係るスイング解析装置であって、前記モード展開部は、前記複数のモードのうち、前記注目モードとは別の複数の候補モードの各々に対し、前記注目モードの前記第１基底を、前記注目モードの前記第１基底が目標値に近づくように、前記注目モードの前記第１基底及び前記候補モードの前記第１基底に直交する軸周りで回転させ、前記回転後の注目モードの前記第１基底と前記目標値との相関係数が最大となる回転を与える前記候補モードを前記ペアモードとして選択する。
【００１１】
第６観点に係るスイング解析装置は、第４観点に係るスイング解析装置であって、前記モード展開部は、回転後の前記注目モードの前記第１基底とその目標値との誤差を表す項を含む誤差関数を最小化する回転角度を算出し、前記注目モードの前記第１基底を、前記ペアモードの前記第１基底とともに、前記軸周りで前記回転角度だけ回転させることにより修正する。
【００１２】
第７観点に係るスイング解析装置は、第６観点に係るスイング解析装置であって、前記誤差関数は、回転後の前記ペアモードの前記第１基底とその目標値との誤差を表す項をさらに含む。
【００１３】
第８観点に係るスイング解析装置は、第１観点から第７観点のいずれかに係るスイング解析装置であって、前記時系列挙動データは、前記複数の注目点の位置、姿勢、速度、加速度、角速度、角加速度、力及びトルクの少なくとも１つを時系列に表すデータである。
【００１４】
第９観点に係るスイング解析装置は、第１観点から第８観点のいずれかに係るスイング解析装置であって、前記打具は、ゴルフクラブである。
【００１５】
第１０観点に係るスイング解析方法は、以下のことを含む。また、第１１観点に係るスイング解析プログラムは、以下のことをコンピュータに実行させる。
・プレイヤーによる打具のスイング時の複数の注目点の挙動を時系列に表す時系列挙動データを取得すること
・前記時系列挙動データを特異値分解することにより、前記挙動を分解すること
【００１６】
前記挙動を分解することは、以下のことを含む。
・前記時系列挙動データを特異値分解することにより、前記挙動から複数のモードにそれぞれ対応する複数の特徴的な挙動を抽出すること
・前記複数のモードのうち、少なくとも１つの注目モードに対応する前記特徴的な挙動を構成する、前記複数の注目点の挙動の励起パターンを表す基底である第１基底を修正すること
・前記修正後の第１基底に基づいて、前記注目モードに対応する前記特徴的な挙動を再構成すること
【発明の効果】
【００１７】
以上の観点によれば、打具のスイングに含まれる個々の挙動に基づいて、新たな挙動を推定することができる。
【図面の簡単な説明】
【００１８】
【図１】本発明の第１実施形態に係るスイング解析装置を含むスイング解析システムの全体構成を示す機能ブロック図。
【図２】ゴルファーの身体に取り付けられたマーカーの位置を示す図。
【図３】第１実施形態に係るスイング解析方法の流れを示すフローチャート。
【図４】（Ａ）時系列挙動データの時間平均を表すスティック線図。（Ｂ）疑似拡張データの時間平均を表すスティック線図。
【図５】時系列挙動データと疑似拡張データとの関係を表す概念図。
【図６】２５個の観測点の第１モードのスティック線図。
【図７】同じゴルファーによる異なるゴルフスイング間でインパクトのタイミングにおける各モードのスティック線図を比較した比較図。
【図８】異なるゴルファーによるゴルフスイング間でインパクトのタイミングにおける各モードのスティック線図を比較した比較図。
【図９】第１実施形態に係る図３のフローチャートに含まれる新たな挙動を推定するステップＳ８のサブステップを示すフローチャート。
【図１０】右特異ベクトル（第１基底）の成分のグラフ。
【図１１】左特異ベクトル（第２基底）の成分のグラフ。
【図１２】第１実施形態に係る右特異ベクトルの修正後の挙動を示すスティック線図。
【図１３】第１実施形態に係る右特異ベクトルの修正後の挙動を示す別のスティック線図。
【図１４】本発明の第２実施形態に係る図３のフローチャートに含まれる新たな挙動を推定するステップＳ８のサブステップを示すフローチャート。
【図１５Ａ】第１実施形態に係る右特異ベクトルの変換を示すイメージ図。
【図１５Ｂ】第２実施形態に係る右特異ベクトルの変換を示すイメージ図。
【図１６】第２実施形態に係る右特異ベクトルの修正後の挙動を示すスティック線図。
【図１７】図６と同じ計測データから作成された、５３点の観測点の第１モードのスティック線図（参考図）。
【図１８】変形例に係る時系列挙動データと疑似拡張データとの関係を表す概念図。
【図１９】別の変形例に係る平均挙動をインパクト時の挙動に近付ける場合の時系列挙動データと疑似拡張データとの関係を表す概念図。
【図２０】本発明の第３実施形態に係る図３のフローチャートに含まれる新たな挙動を推定するステップＳ８のサブステップを示すフローチャート。
【発明を実施するための形態】
【００１９】
以下、図面を参照しつつ、本発明の幾つかの実施形態に係るスイング解析装置、方法及びプログラムについて説明する。
【００２０】
＜１．第１実施形態＞
＜１−１．スイング解析システムの概要＞
図１に、第１実施形態に係るスイング解析装置１を含むスイング解析システム１００の全体構成図を示す。スイング解析システム１００は、ゴルフスイングを解析するためのシステムである。ゴルフスイングは、様々な挙動が含まれる複合動作であり、例えば、腕の振り上げ及び振り下げ、腰の回転、身体の平行移動等の様々な挙動を組み合わせることにより構成される。スイング解析システム１００は、ゴルフスイングの特性をより正確に捉えることができるよう、ゴルフスイングに含まれるこのような特徴的な挙動を個別に評価する。より具体的には、ゴルファーがゴルフクラブをスイングする時の挙動を時系列に表す時系列挙動データを特異値分解によりモード展開することにより、ゴルフスイングを特徴的な挙動に分解する。ユーザは、以上の特異値分解の結果に基づいて、ゴルフスイングに含まれる特徴的な挙動を理解し、ひいてはゴルファーの特性やゴルフクラブの特性等を理解することができる。なお、ここでいうユーザとは、ゴルファー自身又はそのインストラクター、ゴルファーに適したゴルフクラブのフィッティングを行うフィッター、或いはゴルフクラブの開発者等、ゴルフスイングの解析の結果を必要とする者の総称である。
【００２１】
特異値分解の対象となる時系列挙動データは、計測機器２により計測される計測データに基づいて取得される。計測機器２は、スイング解析装置１とともに、スイング解析システム１００を構成する。以下、スイング解析システム１００の各部の構成を説明した後、スイング解析システム１００によるスイング解析方法について説明する。
【００２２】
＜１−２．各部の詳細＞
＜１−２−１．計測機器＞
本実施形態に係る計測機器２は、複数台のカメラ２１，２１，・・・を備えるモーションキャプチャシステムである。モーションキャプチャシステムとしては、例えば、ＶＩＣＯＮ社製の三次元動作分析システムを好ましく使用することができる。複数台のカメラ２１，２１，・・・は、ゴルファーの身体動作の三次元計測が可能なように、ゴルファーがゴルフクラブをスイングする様子を様々な方向から撮影することができる位置に配置される。
【００２３】
本実施形態では、カメラ２１，２１，・・・によりゴルファーの身体の挙動を捉え易いように、身体における観測点の位置にマーカーが取り付けられる。より具体的には、図２に示すように、ゴルファー７の頭、手首、手の指先、肘、肩、腰、膝、踝、脚の指先等のＩ0個（本実施形態では、Ｉ0≧３であり、特にＩ0＝５３）の観測点の位置に、光反射性の球体のマーカー２０，２０，・・・が取り付けられる。この例では、マーカー２０，２０，・・・は、ゴルファー７が着用するボディスーツに取り付けられる。
【００２４】
カメラ２１，２１，・・・は、ゴルファー７がゴルフクラブを使用してゴルフスイングを行う間、その様子を所定のサンプリング周波数で連続撮影する。サンプリング周波数は、例えば、５００Ｈｚとすることができる。本実施形態では、少なくともアドレスから、バック９時、トップ、ダウン９時、インパクト、フォロー３時を順に経て、フィニッシュまでの期間において、ゴルフスイングが計測される。なお、バック９時とは、ゴルファー７を正面から見て、バックスイング中にゴルフクラブが時計の９時の方向を指すタイミングであり、ダウン９時とは、同方向から見て、ダウンスイング中にゴルフクラブが時計の９時の方向を指すタイミングであり、フォロー３時とは、同方向から見て、インパクト後のフォロースイング中にゴルフクラブが時計の３時の方向を指すタイミングである。本実施形態では、計測機器２による計測データとして、少なくともアドレスからフィニッシュまでの時系列の画像データが撮影される。カメラ２１，２１，・・・による撮影は、同期が取られており、同じ時刻に様々な方向から見たゴルフスイングが撮影される。
【００２５】
特に図示されないが、計測機器２には、カメラ２１，２１，・・・の他、カメラ２１，２１，・・・により撮影された画像データを、通信線１７を介して外部のデバイスであるスイング解析装置１に送信するための通信装置も搭載されている。通信装置は、スイング動作の妨げにならないように無線式とすることもできるし、ケーブルを介して有線式にスイング解析装置１に接続することもできる。本実施形態では、カメラ２１，２１，・・・により撮影された画像データは、通信装置を介してリアルタイムにスイング解析装置１に送信される。しかしながら、例えば、カメラ２１，２１，・・・内の記憶装置に画像データを格納しておき、ゴルフスイングの終了後に当該記憶装置から画像データを取り出して、スイング解析装置１に受け渡すようにしてもよい。
【００２６】
＜１−２−２．スイング解析装置＞
スイング解析装置１は、ハードウェアとしては汎用のコンピュータであり、例えば、デスクトップ型コンピュータ、ノート型コンピュータ、タブレットコンピュータ、スマートフォンとして実現される。図１に示すとおり、スイング解析装置１は、コンピュータで読み取り可能なＣＤ−ＲＯＭ等の記録媒体３０から、或いはインターネット等の通信回線を介して、スイング解析プログラム６を汎用のコンピュータにインストールすることにより製造される。スイング解析プログラム６は、計測機器２から送られてくる計測データに基づいて、ゴルフスイングを解析するソフトウェアであり、スイング解析装置１に後述する動作を実行させる。
【００２７】
スイング解析装置１は、表示部１１、入力部１２、記憶部１３、制御部１４及び通信部１５を備える。これらの部１１〜１５は、互いにバス線１６を介して接続されており、相互に通信可能である。表示部１１は、液晶ディスプレイ等で構成することができ、ゴルフスイングの解析結果等をユーザに対し表示する。入力部１２は、マウス、キーボード、タッチパネル等で構成することができ、スイング解析装置１に対するユーザからの操作を受け付ける。
【００２８】
記憶部１３は、ハードディスク等で構成することができる。記憶部１３内には、スイング解析プログラム６が格納されている他、計測機器２から送られてくる計測データが保存される。制御部１４は、ＣＰＵ、ＲＯＭおよびＲＡＭ等から構成することができる。制御部１４は、記憶部１３内のスイング解析プログラム６を読み出して実行することにより、仮想的にデータ取得部１４ａ、モード展開部１４ｂ及び画面作成部１４ｃとして動作する。各部１４ａ〜１４ｃの動作の詳細については、後述する。通信部１５は、計測機器２等の外部のデバイスとの間でデータを送受信する通信インターフェースとして機能する。
【００２９】
＜１−３．スイング解析方法＞
以下、スイング解析システム１００によるスイング解析方法について説明する。本方法では、ゴルファー７によるゴルフスイングが計測され、これが特異値分解され、様々な特徴的な挙動に分解される。本方法は、ゴルフのレッスンや、ゴルファーに対するゴルフクラブのフィッティング、ゴルフクラブの開発等の場面において、ゴルフスイングの特性を理解するのに利用される。以下、詳細に説明する。
【００３０】
図３は、本実施形態に係るスイング解析方法の流れを示すフローチャートである。同図に示すように、最初のステップＳ１では、ゴルフクラブが用意され、ゴルファー７によりゴルフクラブがスイングされる。このとき、計測機器２による計測が行われ、ゴルファー７がゴルフクラブをスイングする様子を表す計測データが収集される。本実施形態に係る計測データとは、上記のとおり、ゴルファー７がゴルフクラブをスイングする時のＩ0個（Ｉ0＝５３）の観測点の挙動を時系列に計測したデータである。また、本実施形態に係る計測データとは、少なくともアドレスからフィニッシュまでの期間において、カメラ２１，２１，・・・により撮影される時系列の画像データである。なお、ここでの計測データには、複数台のカメラ２１，２１，・・・の撮影方向にそれぞれ対応する複数の系統の時系列の画像データが含まれる。計測機器２により計測された計測データは、通信装置からスイング解析装置１に送信される。スイング解析装置１側では、データ取得部１４ａが通信部１５を介してこの計測データを取得し、記憶部１３内に格納する。
【００３１】
続くステップＳ２〜Ｓ４では、データ取得部１４ａは、計測機器２からの計測データに基づいて、ゴルファー７がゴルフクラブをスイングする時の挙動を時系列に表す時系列挙動データを取得する。まず、ステップＳ２では、データ取得部１４ａは、Ｉ0個（本実施形態では、Ｉ0＝５３）の観測点を、Ｉ0個よりも少ないＩ1個（本実施形態では、Ｉ0＞Ｉ1≧２であり、特にＩ1＝２５）のグループに分ける。本実施形態では、グループの分け方は予め設定されており、物理的に近傍にある観測点どうしが同じグループに分けられる。本実施形態では、同じ観測点が重複して複数のグループに分けられることはない。また、ここでのグループは、複数個の観測点を要素とする少なくとも１つの複数要素グループを含む。本実施形態では、Ｉ1個のグループは、全て複数要素グループであり、２つ以上の複数要素グループを含む。
【００３２】
続くステップＳ３では、データ取得部１４ａは、計測データに含まれる画像データを画像処理することにより、画像データに写るマーカー２０，２０，・・・の像を検出し、その位置を特定する。マーカー２０，２０，・・・の像は、複数の撮影方向から検出され、その位置は三次元座標として特定される。マーカー２０，２０，・・・の三次元座標、すなわち、ゴルファー７の身体のＩ0個の観測点の三次元座標は、アドレスからフィニッシュまでの期間におけるサンプリング周波数に対応する時間間隔刻みの各時刻に対して導出される。
【００３３】
言い換えると、ステップＳ３では、ｎ＝１，２，・・・，Ｎ（Ｎは、２以上の整数）、かつ、ｉ＝１，２，・・・，Ｉ0（Ｉ0＝５３）に対し、第ｎ時刻における第ｉ観測点の三次元座標ｒi（ｎ）が導出される。なお、第１時刻がアドレスであり、第Ｎ時刻がフィニッシュである。
【数１】
2020185370.tif 000003
【００３４】
上式のｘi（ｎ）、ｙi（ｎ）及びｚi（ｎ）は、それぞれｒi（ｎ）のＸ、Ｙ及びＺ方向成分である。Ｘ、Ｙ及びＺ方向は、図２に示す通り定義される。すなわち、Ｘ方向は、ゴルファー７の背から腹に向かう方向であり、Ｙ方向は、飛球線方向であり、Ｚ軸は、鉛直下から上に向かう方向である。ｒi（ｎ）は、第ｎ時刻における第ｉ観測点の位置ベクトルである。
【００３５】
今、第ｉ観測点に対し、以下の行列［ｒi］を定義する。行列［ｒi］は、Ｎ行３列の行列であり、第１行から第Ｎ行に、時系列に沿って順に、第ｉ観測点の三次元座標ｒi（１），ｒi（２），・・・，ｒi（Ｎ）を並べた行列である。ステップＳ３では、Ｉ0個の観測点の各々について、その三次元位置を時系列に表すデータ［ｒi］が導出される。
【数２】
2020185370.tif 000004
【００３６】
続くステップＳ４では、データ取得部１４ａは、Ｉ1個（Ｉ1＝２５）のグループに含まれる各グループについて、観測点データを導出する。第ｉグループについての観測点データとは、第ｉグループに属する複数個の観測点を代表する代表点である第ｉ代表点の挙動を時系列に表すデータである（ｉ＝１，２，・・・，Ｉ1）。すなわち、ステップＳ４では、ｉ＝１，２，・・・，Ｉ1、かつ、ｎ＝１，２，・・・，Ｎ（Ｎは、２以上の整数）に対し、第ｎ時刻における第ｉ代表点の挙動を表すｓi（ｎ）が導出される。本実施形態では、ｓi（ｎ）は、第ｎ時刻における第ｉ代表点の三次元座標を表す１行３列の行列である。第ｎ時刻における第ｉ代表点の挙動（三次元座標）は、第ｉ代表点が代表する複数個の観測点の第ｎ時刻における挙動（三次元座標）を平均化することにより導出される。本実施形態では、単純平均により、代表点の挙動が算出される。しかし、これに限らず、相乗平均や重み付け平均等により、代表点の挙動を算出してもよい。
【００３７】
第ｉグループについての観測点データは、以下の行列［ｓi］により表される。行列［ｓi］は、Ｎ行３列の行列であり、第１行から第Ｎ行に、時系列に沿って順に、第ｉ代表点の三次元座標ｓi（１），ｓi（２），・・・，ｓi（Ｎ）を並べた行列である。ステップＳ４では、Ｉ1個の代表点の各々について、その三次元位置を時系列に表すデータ［ｓi］が導出される。
【数３】
2020185370.tif 000005
【００３８】
さらに、全ての代表点に対するＩ1個（Ｉ1＝２５）の行列［ｓ1］，［ｓ2］，・・・，［ｓ25］を、この順に列方向（横方向）に並べた行列［Ｒ］を定義する。行列［Ｒ］は、Ｎ行７５（＝３×２５）列の行列である。
【数４】
2020185370.tif 000006
【００３９】
行列［Ｒ］は、第１時刻から第Ｎ時刻までの２５個の代表点の挙動（位置）を時系列に表す時系列挙動データである。行列［Ｒ］には、Ｉ1個のグループについての観測点データ［ｓ1］，［ｓ2］，・・・，［ｓ25］が含まれる。ステップＳ４では、このような時系列挙動データである行列［Ｒ］が導出される。なお、［Ｒ］の第ｎ行には、第ｎ時刻でのゴルファー７の身体における２５個の代表点の三次元座標が含まれるため、［Ｒ］の第ｎ行は、第ｎ時刻のゴルファー７の姿勢を表している。
【００４０】
続くステップＳ５及びＳ６では、ステップＳ４の時系列挙動データ［Ｒ］を特異値分解することにより、時系列挙動データ［Ｒ］により表されるゴルフスイングの挙動が特徴的な挙動に分解される。より正確には、時系列挙動データ［Ｒ］を疑似的に拡張した疑似拡張データ［Ｒa］が作成され、これが特異値分解される。ゴルフスイングは、非周期的な運動である。そのため、本実施形態では、ゴルフスイングを複製及び合成した疑似拡張データ［Ｒa］を特異値分解することにより、ゴルフスイングを疑似的に周期的な運動とみなして解析が実行される。
【００４１】
ステップＳ５では、データ取得部１４ａが、時系列挙動データ［Ｒ］に基づいて、その時間平均が基準挙動に近付くように、疑似拡張データ［Ｒa］を作成する。疑似拡張データ［Ｒa］は、特異値分解の基準点を基準挙動に近付けるために作成される。基準挙動とは、本実施形態では、アドレス時の挙動である。なお、図４（Ａ）は、実際にあるゴルファーにゴルフスイングを行わせた時の、時系列挙動データ［Ｒ］に含まれる２５個の代表点の位置の時間平均を結んでできたスティック線図であり、図４（Ｂ）は、同ゴルフスイングに基づく疑似拡張データ［Ｒa］に含まれる２５個の代表点の位置の時間平均を結んでできたスティック線図である。同図から分かるとおり、前者の時間平均は、アドレスの姿勢にはなっておらず、少し腕をテイクバックした姿勢になっている。よって、時系列挙動データ［Ｒ］をそのまま特異値分解すると、第１モードをはじめとする各モードの意味を、アドレスを基準として理解できない虞がある。そこで、本実施形態では、疑似拡張データ［Ｒa］が作成される。
【００４２】
より具体的には、データ取得部１４ａは、時系列挙動データ［Ｒ］に基づいて、以下の追加データ［Ｒ（１）］及び［Ｒt］を作成する。今、アドレス時の第ｉ代表点の三次元座標を表す１行３列の行列ｓi（１）を、行方向（縦方向）に２Ｎ個並べた２Ｎ行３列の行列［ｓi（１）］を定義する。このとき、［Ｒ（１）］は、全ての代表点に対する２５個の行列［ｓ1（１）］，［ｓ2（１）］，・・・，［ｓ25（１）］を、この順に列方向（横方向）に並べた２Ｎ行７５（＝３×２５）列の行列として導出される。
【数５】
2020185370.tif 000007
【００４３】
一方、追加データ［Ｒt］は、下式のとおり導出される。すなわち、追加データ［Ｒt］は、時系列挙動データ［Ｒ］を時系列方向に反転した行列である。
【数６】
2020185370.tif 000008
【００４４】
続いて、データ取得部１４ａは、時系列挙動データ［Ｒ］に下式のとおり追加データ［Ｒ（１）］及び［Ｒt］を結合することにより、疑似拡張データ［Ｒa］を作成する。
【数７】
2020185370.tif 000009
【００４５】
図５は、本実施形態に係る時系列挙動データ［Ｒ］と疑似拡張データ［Ｒa］との関係を表す概念図である。すなわち、アドレス時の挙動を時間長２Ｎだけ繰り返し複製した追加データ［Ｒ（１）］を、時系列挙動データ［Ｒ］のアドレスの時刻ｎ＝１の前に結合する。さらに、時系列挙動データ［Ｒ］を時系列方向に反転した時間長Ｎの追加データ［Ｒt］を、時系列挙動データ［Ｒ］のフィニッシュの時刻ｎ＝Ｎの後に結合し、さらにその後に、追加データ［Ｒ（１）］を追加する。このように、追加データ［Ｒ（１）］及び［Ｒt］は、時系列挙動データ［Ｒ］に連続するように結合され、こうして作成される疑似拡張データ［Ｒa］は、連続的に値が変化するデータとなる。疑似拡張データ［Ｒa］は、６Ｎ行７５（＝３×２５）列の行列である。以上のとおり作成される疑似拡張データ［Ｒa］も、ゴルファー７がゴルフクラブをスイングする時の挙動を時系列に表す時系列挙動データである。
【００４６】
続くステップＳ６では、モード展開部１４ｂは、下式に従って、ゴルフスイングの挙動を表す疑似拡張データ［Ｒa］を特異値分解する。より具体的には、モード展開部１４ｂは、疑似拡張データ［Ｒa］に基づき、初期挙動［Ｒ0］を導出する。そして、モード展開部１４ｂは、疑似拡張データ［Ｒa］を初期挙動［Ｒ0］を基準として特異値分解することにより、疑似拡張データ［Ｒa］により表されるゴルフスイングの挙動から、特徴的な挙動λjｖjｚjTを抽出する。なお、ｊ＝１，２，・・・，Ｊ（Ｊは、７５以下の整数）である。
【数８】
2020185370.tif 000010
【００４７】
上式中の初期挙動［Ｒ0］は、［Ｒa］の時間平均である。［Ｒ0］は、［Ｒa］を時間方向（行方向）に平均化した平均行列を作成し、その平均行列を６Ｎ個行方向に並べた６Ｎ行７５（＝３×２５）列の行列である。すなわち、［Ｒ0］の第ｌ列は、［Ｒa］の第ｌ列の６Ｎ個の成分の平均値を、行方向に６Ｎ個並べた値列である。ところで、疑似拡張データ［Ｒa］には、上記のとおり、時間長４Ｎ分のアドレス時の挙動（本実施形態では、ゴルファー７の姿勢）を表すデータが追加されている。そのため、疑似拡張データ［Ｒa］の時間平均［Ｒ0］により表される挙動は、時系列挙動データ［Ｒ］の同様の時間平均により表される挙動よりも、アドレス時の挙動に近付いている。なお、本実施形態では、単純平均により平均行列が算出されるが、これに限らず、相乗平均や重み付け平均等により、平均行列を算出してもよい。また、初期挙動［Ｒ0］は、［Ｒa］の全体の時間平均ではなく、［Ｒa］に含まれるアドレスの前後の期間の行列を時間方向（行方向）に平均化した平均行列を作成し、その平均行列を６Ｎ個行方向に並べたものとしてもよい。また、初期挙動［Ｒ0］は、［Ｒa］に含まれるアドレスの時刻に対応する行を、６Ｎ個行方向に並べたものとしてもよい。
【００４８】
以上のとおり、本実施形態では、疑似拡張データ［Ｒa］とその時間平均［Ｒ0］との差分である［Ｒa］−［Ｒ0］が特異値分解され、ｊ＝１，２，・・・，Ｊ（Ｊは、７５以下の整数）に対し、λj，ｚj及びｖjが導出される。つまり、疑似拡張データ［Ｒa］が、その時間平均［Ｒ0］を基準として特異値分解される。λjは、第ｊモードの特異値であり、［Ｒa］−［Ｒ0］に対する第ｊモードの割合を表す。ｖjは、［Ｒa］−［Ｒ0］の第ｊモードの左特異ベクトルであり、Ｎ次元である。ｚjは、［Ｒa］−［Ｒ0］の第ｊモードの右特異ベクトルであり、７５次元である。ｚjは、第ｊモードでのゴルファー７の身体における各代表点の挙動の励起パターンを表す基底であり、ｖjは、第ｊモードでのｚjの時間に対する励起パターンを表す基底である。
【００４９】
また、モード展開部１４ｂは、ｊ＝１，２，・・・，Ｊに対し、特異値λjに基づいて、［Ｒa］に対する第ｊモードの貢献度である寄与率Ｃjを下式に従って導出する。
【数９】
2020185370.tif 000011
【００５０】
以上の特異値分解により、時系列挙動データ［Ｒ］と同じくゴルフスイングの挙動を時系列に表す疑似拡張データ［Ｒa］から、ゴルフスイングに含まれる第１モードから第Ｊモードの各モードに対応する特徴的な挙動λjｖjｚjTが抽出される。この結果を利用し、モード展開部１４ｂは、ｊ＝１，２，・・・，Ｊに対し、下式に従って、第ｊモードの特徴的な挙動λjｖjｚjTに初期挙動［Ｒ0］を合成した第ｊモードの挙動を表す［Ｒj］を導出する。［Ｒj］は、ゴルフスイングの挙動から分解された第ｊモードの特徴的な挙動λjｖjｚjTに基づくため、以下、［Ｒj］の表す挙動を、分解挙動ということがある。
【数１０】
2020185370.tif 000012
【００５１】
なお、Ｊは、最大で７５であるが、何次モードまでの挙動を導出するかは、解析の目的等に応じて適宜選択され得る。なお、本発明者らの検討によると、一般に、第１モードの寄与率だけで９０％程度であり、第１モードから第２モードまでの累積寄与率は９７％程度であり、第１モードから第３モードまでの累積寄与率は９９％を超え、第１モードから第５モードまでの累積寄与率は約９９．８％に達する。よって、第５モードまでの挙動を導出することにより、ゴルフスイングを概ね評価し得る。
【００５２】
また、モード展開部１４ｂは、ｋ＝１，２，・・・，Ｊに対し、第１モードから第ｋモードまでの挙動を累積した挙動（以下、累積挙動ということがある）を表す［Ｒ1〜k］を下式に従って導出する。［Ｒ1〜k］は、ゴルフスイングの挙動から分解された第１モードから第ｋモードまでの特徴的な挙動λjｖjｚjTに基づくため、以下、［Ｒ1〜k］の表す挙動についても、分解挙動ということがある。
【数１１】
2020185370.tif 000013
【００５３】
［Ｒ1〜k］は、各代表点の各時刻における第１モードから第ｋモードまでの累積挙動（三次元座標）を成分として含む。よって、ｋ＝１，２，・・・，Ｊに対し、［Ｒ1〜k］を導出することは、各モードの各時刻におけるゴルファーの姿勢を導出することを意味する。
【００５４】
以上に説明したステップＳ１〜Ｓ６は、同じ又は異なるゴルファーが同じ又は異なるゴルフクラブを用いて行う複数回のゴルフスイングに対し繰り返し実行することができる。
【００５５】
続くステップＳ７では、画面作成部１４ｃが、以上の特異値分解の結果を表示する画面（以下、結果画面という）を作成し、表示部１１上に表示させる。例えば、結果画面には、図６に示すようなスティック線図を表示することができる。図６は、あるゴルファーがゴルフスイングを行った時の第１モードの挙動を表す［Ｒ1］（＝［Ｒ1〜1］）を算出し、これに含まれる２５個の代表点の三次元座標を結んでできたスティック線図である。同図には、アドレスからフィニッシュまでの７つの時刻におけるスティック線図が示されている。なお、［Ｒ1］に代えて又は加えて、ｋ≧２に対しても、［Ｒ1〜k］を算出し、これに基づく第ｋモードの様々な時刻におけるスティック線図を結果画面に表示してもよい。このようなスティック線図を見たユーザは、ゴルフスイングに含まれる個々の挙動を適切に評価し、ゴルフスイングの特性を理解することができる。
【００５６】
また、結果画面には、同じ又は異なるゴルファーが同じ又は異なるゴルフクラブをスイングした時のスティック線図の比較図を表示することもできる。比較図とは、例えば、複数のスティック線図を重ねて又は並べて表示する図である。このとき、比較される複数のスティック線図としては、典型的には、ステップＳ１〜Ｓ６の繰り返しにより導出された異なるゴルフスイングについてのスティック線図が想定される。しかしながら、ステップＳ１〜Ｓ６により導出されたスティック線図の比較の対象として、予め記憶部１３内に格納されていた、又は通信部１５を介して外部から取得された挙動[Ｒ1-k]に基づくスティック線図が選択されてもよい。例えば、上級者の理想的なゴルフスイングに関する挙動[Ｒ1-k]に基づくスティック線図が、比較の対象とされてもよい。
【００５７】
一例として、図７は、同じゴルファーが同じゴルフクラブを用いてドロースイングとフェードスイングを行った時の、インパクトのタイミングにおける様々なモードのスティック線図を重ねた図である。また、別の例を挙げると、図８は、異なるゴルファーＧ１及びＧ２がそれぞれのゴルフクラブを用いてスイングを行った時の、インパクトのタイミングにおける様々なモードのスティック線図を重ねた図である。このような図を見たユーザは、異なるゴルフスイングをゴルフスイングに含まれる個々の挙動のレベルで比較することができ、ゴルフスイング間の特性の違いを適切に理解することができる。
【００５８】
なお、非特許文献１及び２にも示されるとおり、第１〜第５モードの挙動は、以下のような挙動であることが知られている。このような知見と合わせて、スティック線図に表れるゴルフスイングの特性についての理解をより深めることもできる。
【表１】
2020185370.tif 000014
【００５９】
結果画面には、ｊ＝１，２,・・・,Ｊに対し、特異値λj及び寄与率Ｃjを表示することもできる。また、ｋ＝１，２,・・・,Ｊに対し、累積寄与率Ｃ1＋Ｃ2＋・・・＋Ｃkを表示することもできる。
【００６０】
また、モード展開部１４ｂは、スティック線図のようにモード別の挙動を図として再現するだけでなく、ｊ＝１，２，・・・，Ｊ、ｋ＝１，２,・・・,Ｊに対し、モード別の挙動［Ｒj］、［Ｒj］−［Ｒ0］、［Ｒ1〜k］及び［Ｒ1〜k］−［Ｒ0］を様々に分析することができる。例えば、モード展開部１４ｂは、第２モードがスイング全体における肩の回転を表すことに鑑みて（表１参照）、［Ｒ2］、［Ｒ2］−［Ｒ0］、［Ｒ1〜2］及び［Ｒ1〜2］−［Ｒ0］に基づき、肩の回転量を定性的又は定量的に評価することができる。また、モード展開部１４ｂは、比較されるゴルフスイングについて、挙動の差が顕著な部分を検出するようにしてもよい。画面作成部１４ｃは、このようなさらなる分析の結果についても結果画面に表示することができる。
【００６１】
続くステップＳ８では、モード展開部１４ｂは、ステップＳ１〜Ｓ６により導出されたあるゴルフスイングについての右特異ベクトルｚjを修正し、ゴルファー７の新たな挙動を推定する。本実施形態では、ステップＳ８は、ユーザが望む場合に実行され、詳細には、図９のとおりに実行される。
【００６２】
まず、ステップＳ８のサブステップＳ１１及びＳ１２として、モード展開部１４ｂは、第１〜第Ｊモードのうち、少なくとも１つのモードｊに対応する特徴的な挙動λjｖjｚjTを構成する右特異ベクトルｚjを修正する。以下、修正の対象となる右特異ベクトルｚjに対応するモードを、注目モードと呼ぶことがある。ここで、右特異ベクトルｚjとは、上記のとおり、第ｊモードでのゴルファー７の身体における各代表点の挙動の励起パターンを表す基底である。
【００６３】
図１０は、図７の例に対応する第４モードの右特異ベクトルｚ4に含まれる７５個の成分の値を示している。既に述べたことから明らかなとおり、右特異ベクトルｚjに含まれる７５個の成分は、２５個の代表点のＸ、Ｙ及びＺ方向の挙動（位置）の励起の大きさをそれぞれ表す。図１０では、横軸に沿って７５個の成分が並べられており、縦軸は各成分の値（励起の大きさ）を示している。一方、図１１は、図７の例に対応する第４モードの左特異ベクトルｖ4に含まれる多数の成分の値を示している。既に述べたとおり、左特異ベクトルｖjは、第ｊモードでのｚjの時間に対する励起パターンを表す基底である。図１１では、横軸が時間軸を、縦軸が各時刻の成分の値（励起の大きさ）を示している。
【００６４】
サブステップＳ１１では、モード展開部１４ｂは、ユーザからの入力に従って、第１〜第Ｊモードのうち、修正の対象となる１又は複数の注目モードを決定する。また、モード展開部１４ｂは、ユーザからの入力に従って、各注目モードの右特異ベクトルｚjに含まれる複数の成分のうち、値を修正する１又は複数の成分（以下、注目成分ということがある）と、各注目成分の修正量とを決定する。
【００６５】
ユーザが、注目モード、注目成分及びその修正量を設定する際には、ステップＳ７での結果画面を参照することができる。例えば、ユーザが、インパクト時におけるドロースイングの姿勢をフェードスイングの姿勢に近づけるような修正をしようとする場合を考える。この場合、図７のインパクト時のスティック線図を参照すると、第４モードから両ゴルフスイングの違いが顕著になり、特に右肘の動きに差が見られる。よって、図７のスティック線図から（特に点線の囲みを参照）、注目モードとして第４モードを、注目成分として右肘の挙動（三次元位置）に対応する３つの成分を選択することができる。また、ユーザは、図７のスティック線図を参照しながら、図１０の右肘のＸ、Ｙ及びＺ方向の挙動にそれぞれ対応する３つの成分の値に対し、それぞれ任意の修正量を設定する。なお、ユーザによる修正量の設定を補助するために、結果画面には、ｊ＝１，２，・・・Ｊに対し、図１０及び図１１のような右特異ベクトルｚj及び左特異ベクトルｖjの成分のグラフを表示してもよい。このとき、ユーザは、図１１の左特異ベクトルｖjのグラフから、現在注目しているインパクト時の励起の大きさを読み取り、図７のスティック線図から第４モードでの目標とのずれを読み取り、さらに図１０の右特異ベクトルｚjのグラフから、第４モードでの注目成分の値を読み取り、これらを総合して修正量を判断することができる。なお、修正量の設定のために、１．１倍、０．９８倍といった修正倍率を入力してもよい。
【００６６】
ユーザが、入力部１２を介して注目モード、注目成分及びその修正量を設定すると、モード展開部１４ｂは、これらの情報に基づいて、注目モードの右特異ベクトルｚjの修正を行う（サブステップＳ１２）。より具体的には、モード展開部１４ｂは、各注目モードの右特異ベクトルｚjに含まれる各注目成分の値を、それぞれの修正量に従って修正する。以下では、注目成分の修正後の右特異ベクトルｚjを、ｚj'と表す。
【００６７】
続くサブステップＳ１３では、モード展開部１４ｂは、修正後の右特異ベクトルｚj'に基づいて、注目モードに対応する特徴的な挙動λjｖjｚjTを再構成し、特徴的な挙動λjｖjｚj'Tを導出する。また、モード展開部１４ｂは、注目モード以降の各モードに対し、修正後の特徴的な挙動λjｖjｚj'Tに基づいて、分解挙動［Ｒj］及び［Ｒ1〜k］を再構成する。このとき、数１０及び１１の計算式において、注目モードのλjｖjｚjTをλjｖjｚj' Tに置換することにより、修正後の分解挙動［Ｒj］'及び［Ｒ1〜k］'が導出される。
【００６８】
注目モードから第Ｊモードまでの挙動［Ｒj］'及び［Ｒ1〜k］'は、元の挙動［Ｒj］及び［Ｒ1〜k］を修正した新たな挙動である。図１２は、図７のドロースイングをフェードスイングに近づけるように、第４モードでの右肘のＸ、Ｙ及びＺ方向の挙動（位置）を修正したときの、［Ｒ1〜8］'に基づくインパクト時のスティック線図である。なお、図１２では、右肘より先端側の右手等の部位の挙動（位置）については、右肘の修正に伴って平行移動されている。このように、現実の挙動の変化をイメージし易いように、修正される部位より先端側の拘束されていない部位の挙動については、修正される部位と同様に修正してもよい。また、図１３は、上体が前に進みすぎとの指摘を受けていた初級者の［Ｒ1〜8］'に基づくインパクト時のスティック線図を、ステップＳ１１〜Ｓ１３の方法で修正した図である。図１３の例では、第８モードの右特異ベクトルｚjに含まれる上体に対応する注目成分について、上体が後ろに下がるような修正が加えられた。図１３でも、図１２と同様に、上体の挙動の修正に伴って、上体より先端側の拘束されていない頭、右腕及び左腕等の部位の挙動が平行移動されている。
【００６９】
画面作成部１４ｃは、結果画面上に、ステップＳ８の分析結果を表示する。ここでいうステップＳ８の分析結果には、図１２及び図１３に示すようなスティック線図が含まれる。なお、図１２及び図１３は、修正後の挙動と、修正前の挙動及び／又は修正の目標となる挙動との比較図が示されているが、比較図ではなく、修正後の挙動のみを示すスティック線図を表示してもよい。また、モード展開部１４ｂは、修正後の挙動をスティック線図のような図として再現するだけでなく、ｊ＝１，２，・・・，Ｊ、ｋ＝１，２,・・・,Ｊに対し、修正前の挙動の分析と同様に、モード別の修正後の挙動［Ｒj］'、［Ｒj］'−［Ｒ0］、［Ｒ1〜k］'及び［Ｒ1〜k］'−［Ｒ0］を様々に分析することができる。画面作成部１４ｃは、このようなさらなる分析の結果についても結果画面に表示することができる。
【００７０】
＜２．第２実施形態＞
次に、本発明の第２実施形態について説明する。第２実施形態は、第１実施形態と比較して、図３に示すスイング解析方法のステップＳ８の内容のみが異なる。よって、以下では、第２実施形態に係るステップ８の内容についてのみ説明する。第２実施形態に係るステップＳ８でも、モード展開部１４ｂは、ステップＳ１〜Ｓ６により導出されたあるゴルフスイングについての右特異ベクトルｚjを修正し、ゴルファー７の新たな挙動を推定する。
【００７１】
第２実施形態に係るステップＳ８は、図１４のとおりに実行される。まず、ステップＳ８のサブステップＳ２１〜Ｓ２３として、モード展開部１４ｂは、第１〜第Ｊモードのうち、少なくとも１つのモードｊに対応する特徴的な挙動λjｖjｚjTを構成する右特異ベクトルｚjを修正する。第２実施形態でも、修正の対象となる右特異ベクトルｚjに対応するモードを、注目モードと呼ぶことがある。第２実施形態では、あるゴルフスイングをある別のゴルフスイング（以下、目標スイングという）に近づける修正が行われ、注目モードの右特異ベクトルｚjが、目標スイングの挙動を構成する目標基底ｗjに近づくように修正される。
【００７２】
サブステップＳ２１では、モード展開部１４ｂは、ユーザからの入力に従って、第１〜第Ｊモードのうち、修正の対象となる１又は複数の注目モードを決定する。また、モード展開部１４ｂは、ユーザからの入力に従って、各注目モードに対し、目標基底ｗjを決定する。ユーザが注目モードを設定する方法は、第１実施形態と同様である。一方、目標基底ｗjの設定のためには、ユーザは、目標スイングを指定する。モード展開部１４ｂは、ユーザにより指定された目標スイングの挙動を構成する注目モードの右特異ベクトルｚjを取得し、これを目標基底ｗjとして決定する。なお、目標基底ｗjは、ステップＳ１〜Ｓ６の繰り返しにより取得されてもよいし、予め記憶部１３内に格納されていてもよいし、通信部１５を介して外部から取得されてもよい。目標スイングは、例えば、上級者の理想的なゴルフスイングとすることができる。
【００７３】
第２実施形態では、注目モードの右特異ベクトルｚjが目標基底ｗjに近づくように修正されるが、このとき、注目モードの右特異ベクトルｚjだけでなく、注目モードとは別のモード（以下、ペアモードという）の右特異ベクトルｚjも修正される。これは、注目モードの右特異ベクトルｚjの修正後も、全てのモードに対応する右特異ベクトルｚ1，ｚ2，・・・，ｚJ間の直交性を確保するためである。すなわち、第１実施形態では、注目モードの右特異ベクトルｚjの修正により、注目モードの右特異ベクトルｚjと、他のモードの右特異ベクトルｚjとの間で直交性が損なわれる可能性が高い（図１５Ａ参照）。直交性が損なわれたとしても、ゴルフスイングのモデル化は可能であるが、直交性の損なわれる程度が大きくなると、モデルとして不適切となり得る。第２実施形態の修正のアルゴリズムは、この問題を解消することができる。
【００７４】
続くサブステップＳ２２では、モード展開部１４ｂは、第１〜第Ｊモードの中から、各注目モードとペアリングするペアモードを選択する。本実施形態では、第１〜第Ｊモードのうち、注目モード以外のモード（以下、候補モードという）は全て、ペアモードとしての候補となる。モード展開部１４ｂは、全ての候補モードに対し、順次、注目モードとの仮のペアリングを行う。そして、各仮のペアリングに対し、注目モードの右特異ベクトルｚjを、目標基底ｗjに近づくように、仮のペアリングに対応する軸周りで回転させる（図１５Ｂ参照）。仮のペアリングに対応する軸とは、注目モードの右特異ベクトルｚjと、これと仮のペアリングがされている候補モードの右特異ベクトルｚjとの両方に直交する軸である。
【００７５】
このとき、本実施形態では、注目モードの右特異ベクトルｚjは、目標基底ｗjに最も近づくように回転させられる。注目モードの右特異ベクトルｚjと目標基底ｗjとの距離ｄは、両ベクトルの残差二乗和（ｚj−ｗj）2により表すことができる。なお、ｚj2＝ｗj2＝１であるため、残差二乗和（ｚj−ｗj）2＝ｚj2＋ｗj2−２ｚj・ｗj＝２−２ｚj・ｗjである。一方で、注目モードの右特異ベクトルｚjと目標基底ｗjとの相関係数ｃは、下式で表され、その分母はモードのノルムであるため１である。よって、注目モードの右特異ベクトルｚjが目標基底ｗjに最も近づくときは、距離ｄが最小となる、すなわち、相関係数ｃが最大となる場合として判断される。
【数１２】
2020185370.tif 000015
【００７６】
モード展開部１４ｂは、注目モードの右特異ベクトルｚjを、一定の範囲内で小刻みに回転させながら、回転後の注目モードの右特異ベクトルｚjと目標基底ｗjとの相関係数ｃを算出し、相関係数ｃが最大となるときの回転角度θ及び当該最大の相関係数ｃmaxを特定する。このとき、例えば、−π／２からπ／２の範囲内で、π／１０００刻みで回転させることができる。なお、距離ｄ＝（ｚj−ｗj）2は、回転後の注目モードの右特異ベクトルｚjと目標基底ｗjとの誤差を表す誤差関数であり、ここでは、このような誤差関数を最小化する回転角度θが導出される。そして、各候補ベクトルに対し、以上のとおりに相関係数ｃmaxを導出した後、全ての候補ベクトルの中から、それぞれに対応する相関係数ｃmaxが最大となる回転を与える候補モードをペアモードとして選択する。
【００７７】
続くサブステップＳ２３では、モード展開部１４ｂは、注目モードの右特異ベクトルｚjを、ペアモードの右特異ベクトルｚjとともに、これらの両ベクトルに直交する軸周りで回転させる（図１５Ｂ参照）。このときの回転角度は、注目モードの右特異ベクトルｚjに最大の相関係数ｃmaxを与えた（すなわち、誤差関数を最小化する）回転角度θであり、既にサブステップＳ２２で特定されている。以上により、注目モードの右特異ベクトルｚj及びペアモードの右特異ベクトルｚjが、同じ方向に同じ量だけ回転させられることにより修正される。第２実施形態でも、修正後の右特異ベクトルｚjを、ｚj'と表す。このような回転後の注目モードの右特異ベクトルｚj'及びペアモードの右特異ベクトルｚj'は、残りのモードの右特異ベクトルｚjと直交性を維持する。
【００７８】
サブステップＳ２３の後は、第１実施形態と同様のサブステップＳ１３が実行される。図１６は、あるゴルフスイングを目標スイングに近づけるように、サブステップＳ２１〜Ｓ２３により右特異ベクトルｚjを修正したときの［Ｒ1〜8］'に基づくスティック線図である。この例では、第１モード及び第３モードが注目モードとして選択され、それぞれのペアモードとして第２モード及び第４モードが選択された。同図に示すとおり、第２実施形態によっても、元のゴルフスイングを目標スイングに近づけることができる。
【００７９】
＜３．第３実施形態＞
次に、本発明の第３実施形態について説明する。第３実施形態は、第１及び第２実施形態と比較して、図３に示すスイング解析方法のステップＳ８の内容のみが異なる。よって、以下では、第３実施形態に係るステップＳ８の内容についてのみ説明する。第３実施形態に係るステップＳ８でも、モード展開部１４ｂは、ステップＳ１〜Ｓ６により導出されたあるゴルフスイングについての右特異ベクトルｚjを修正し、ゴルファー７の新たな挙動を推定する。
【００８０】
第３実施形態に係るステップＳ８は、図２０のとおりに実行される。まず、ステップＳ８のサブステップＳ３１〜Ｓ３４として、モード展開部１４ｂは、第１〜第Ｊモードのうち、少なくとも１つのモードｊに対応する特徴的な挙動λjｖjｚjTを構成する右特異ベクトルｚjを修正する。第３実施形態でも、修正の対象となる右特異ベクトルｚjに対応するモードを、注目モードと呼ぶことがある。
【００８１】
サブステップＳ３１では、モード展開部１４ｂは、ユーザからの入力に従って、第１〜第Ｊモードのうち、修正の対象となる１又は複数の注目モードを決定する。ユーザが注目モードを設定する方法は、第１及び第２実施形態と同様である。また、モード展開部１４ｂは、ユーザからの入力に従って、各注目モードに対し、目標基底ｗjを決定する。ユーザが目標基底ｗjを設定する方法は、第２実施形態と同様であり、目標スイングを指定する。
【００８２】
第３実施形態でも、第２実施形態と同様に、注目モードの右特異ベクトルｚjが、目標スイングの挙動を構成する目標基底ｗjに近づくように修正されるとともに、注目モードとペアリングされるペアモードの右特異ベクトルｚjも修正される。これは、注目モードの右特異ベクトルｚjの修正後も、全てのモードに対応する右特異ベクトルｚ1，ｚ2，・・・，ｚJ間の直交性を確保するためである。以下、区別のために、修正の対象となる注目モードの右特異ベクトルｚjをｚtjと呼び、修正の対象となるペアモードの右特異ベクトルｚjをｚpjと呼ぶことがある。なお、第３実施形態では、注目モードの右特異ベクトルｚtjが目標基底ｗjに近づくように修正されるだけでなく、ペアモードの右特異ベクトルｚpjも、所定の目標基底ｗpjに近づくように修正される。
【００８３】
続くサブステップＳ３２では、モード展開部１４ｂは、第１〜第Ｊモードの中から、各注目モードとペアリングするペアモードを選択する。また、モード展開部１４ｂは、各ペアモードに対し、目標基底ｗpjを決定する。本実施形態では、ペアモードは、予め定められているルールに従って、注目モードがどのモードであるかに応じて選択される。例えば、注目モードが第１モードであれば、第２モードをペアモードとし、注目モードが第２モードであれば、第１モードをペアモードとし、注目モードが第３モードであれば、第４モードをペアモードとする、といったルールが予め定められている。また、モード展開部１４ｂは、ステップＳ３１においてユーザにより指定された目標スイングの挙動を構成するペアモードの右特異ベクトルｚjを取得し、これを目標基底ｗpjとして決定する。
【００８４】
続くサブステップＳ３３では、モード展開部１４ｂは、注目モードの右特異ベクトルｚtj及びペアモードの右特異ベクトルｚpjを回転させる回転角度θrを算出する。θrは、以下のとおりに算出することができる。
【００８５】
まず、注目モードの右特異ベクトルｚtj及びペアモードの右特異ベクトルｚpjを、これらのベクトルｚtj及びｚpjに直交する軸周りで角度θだけ回転させることを考える。このとき、この回転は、下式のとおり表すことができる。なお、ｚtj'及びｚpj'は、それぞれ、回転後の右特異ベクトルｚtj及びｚpjである。
【数１３】
2020185370.tif 000016
【００８６】
また、下式のとおり、誤差関数Ｅを定義する。第１項は、回転後の注目モードの右特異ベクトルｚtj'と、これを近づける目標値である目標基底ｗjとの誤差を表している。第２項は、回転後のペアモードの右特異ベクトルｚpj'と、これを近づける目標値である目標基底ｗpjとの誤差を表している。
【数１４】
2020185370.tif 000017
【００８７】
上式は、以下のとおり展開される。なお、［Ｉ］は、単位行列である。
【数１５】
2020185370.tif 000018
【００８８】
さらに、数１５の式に数１３の式を代入すると、下式のとおりとなる。誤差関数Ｅは、θの関数である。
【数１６】
2020185370.tif 000019
【００８９】
右特異ベクトルｚtj及びｚpjをそれぞれの目標基底ｗj及びｗpjに近づけるということは、誤差関数Ｅを最小化することである。また、誤差関数Ｅが最小化されるとき、誤差関数Ｅは極値をとる。よって、下式のとおり、数１６の誤差関数Ｅを微分し、これがゼロになるθを算出する。
【数１７】
2020185370.tif 000020
【００９０】
上式は、下式のとおり変形される。モード展開部１４ｂは、下式に従って、誤差関数Ｅを最小化するθを算出し、これを回転角度θrとする。
【数１８】
2020185370.tif 000021
【００９１】
続くステップＳ３４では、モード展開部１４ｂは、注目モードの右特異ベクトルｚtjを、ペアモードの右特異ベクトルｚpjとともに、これらの両ベクトルｚtj及びｚpjに直交する軸周りで回転角度θrだけ回転させる。以上により、注目モードの右特異ベクトルｚtj及びペアモードの右特異ベクトルｚpjが、同じ方向に同じ量だけ回転させられ、修正される。
【００９２】
続くステップＳ３５では、上述したステップＳ１３と同様に、モード展開部１４ｂは、修正後の右特異ベクトルｚtj'に基づいて、注目モードに対応する特徴的な挙動λjｖjｚtj'Tを再構成する。また、モード展開部１４ｂは、修正後の右特異ベクトルｚpj'に基づいて、ペアモードに対応する特徴的な挙動λjｖjｚpj'Tを再構成する。さらに、モード展開部１４ｂは、上述したステップＳ１３と同様に、修正後の特徴的な挙動λjｖjｚtj'T及びλjｖjｚpj'T基づいて、分解挙動［Ｒj］及び［Ｒ1〜k］を再構成する。このとき、数１０及び１１の計算式において、注目モード及びペアモードのλjｖjｚjTを、それぞれλjｖjｚtj'T及びλjｖjｚpj'Tに置換することにより、修正後の分解挙動［Ｒj］'及び［Ｒ1〜k］'が導出される。ステップＳ８の分析結果は、第１及び第２実施形態と同様に、適宜、結果画面上に表示される。
【００９３】
＜４．特徴＞
＜４−１＞
第１〜第３実施形態のステップＳ８によれば、計測された元の挙動に対し、これに含まれる一部のモードの挙動のみを修正することができる。よって、修正後の新たな挙動には、元の挙動の特性が残る。例えば、初級者は、理想的なゴルフスイングを示されても、自身のゴルフスイングとの違いが大きすぎると、多くの場合、どのように自身の動作を修正すればよいのか分からない。しかしながら、自分らしさの残る動作例を示されると、どのように修正すればよいかをイメージし易い。例えば、ｋ＝１０等、高次のｋに対する［Ｒ1〜k］'から再現される姿勢は、ゴルファーにゴルフスイングの改善のための有用な指針を与えることができる。以上のようなステップＳ８による修正は、新たなゴルフスイングのシミュレーションを可能にする。また、修正後の挙動は、ゴルフのレッスンの場面の他、ゴルフクラブのフィッティングの場面、ゴルフクラブの開発の場面等で様々な知見をもたらし得る。
【００９４】
＜４−２＞
以上の実施形態によれば、Ｉ0個の観測点を、これよりも少ないＩ1個の代表点にまとめることにより、その後の特異値分解によるモード数（次元数）が削減される。すなわち、特異値分解の対象となる注目点（観測点又は代表点）の数が多いと、スイングの自由度は高くなるが、その分、モード数が増えるため、１つのモードが有する情報が少なくなる。この場合、特異値分解の結果から、ゴルフスイングに含まれる個々の挙動を適切に評価することがむしろ困難になる場合もある。かといって、モード数を減らすべく、やみくもに観測点の数を減らしてしまうと、スイングの解析の精度が低下し得る。この点、以上の実施形態では、複数の観測点をグループ化し、複数の観測点の情報を１つの代表点の情報に圧縮することが行われる。よって、スイングの解析の精度を維持したまま、各モードの情報量を維持することができる。よって、ゴルフスイングに含まれる個々の挙動（各モードの挙動）をより適切に評価することができる。
【００９５】
既に説明した図６は、２５個の代表点の第１モードのスティック線図である。これに対し、図１７は、図６と同じ計測データに基づき、代表点を作成せずに、第１モードでの５３個の観測点の三次元座標を結んでできたスティック線図である。図６及び図１７を比較すれば分かるが、注目点を減らしても、スイングの挙動全体を表現することに成功している。一方で、注目点の多い図１７のスティック線図からは、上体の伸縮が見える。これに対し、注目点が少なく、シンプルな図６のスティック線図では、上体の伸縮というよりはむしろ、回転運動が見える。このように、注目点を圧縮し、挙動をよりシンプルに表現することにより、新たな情報を得ることができる。
【００９６】
また、以上のモード数の削減によりシンプルに挙動が表現される場合、異なるスイングどうしの比較がより容易になる。すなわち、特異値分解の対象となる注目点が多すぎると、ゴルファーの体格や姿勢等が正確に再現されすぎるため、むしろ比較が難しくなり得る。このような場合、シンプルな挙動でスイングを表現することにより、問題を解消することができる。
【００９７】
また、複数の観測点の挙動を１つの代表点の挙動に平均化することにより、回転の情報をカットすることもできる。例えば、飛球線方向周りの腰の回転は、ゴルフスイングを分析する上で、余り必要とされないことがある。このような場合に、腰の近傍の観測点を１つの代表点にまとめることにより、不要な情報を減らし、より必要な情報に注目することができる。
【００９８】
＜５．変形例＞
以上、本発明の幾つかの実施形態について説明したが、本発明は上記実施形態に限定されるものではなく、その趣旨を逸脱しない限りにおいて、種々の変更が可能である。例えば、以下の変更が可能である。また、以下の変形例の要旨は、適宜組み合わせることができる。
【００９９】
＜５−１＞
上述した計測機器２の構成は例示であり、例えば、計測機器２は、慣性センサ、距離画像センサ等であってもよいし、慣性センサ、距離画像センサ及びモーションキャプチャシステムからなる群から選択される２つ以上の機器を任意に組み合わせたものであってもよい。
【０１００】
＜５−２＞
上記実施形態では、ゴルフスイングが解析されたが、ゴルフクラブ以外の打具、例えば、テニスラケットやベースボールバット等のスイングを解析の対象とすることもできる。
【０１０１】
＜５−３＞
上記実施形態では、基準挙動は、アドレス時の挙動とされた。しかしながら、基準挙動は、他の時刻の挙動、例えば、インパクト又はトップの挙動とすることもできる。例えば、トップ時の挙動を基準挙動とする場合には、図１８に示すように疑似拡張データを作成することができ、インパクト時の挙動を基準挙動とする場合には、図１９に示すように疑似拡張データを作成することができる。この場合、初期挙動［Ｒ0］は、インパクト又はトップ時の挙動に近づく。
【０１０２】
＜５−４＞
上記実施形態では、特異値分解の対象となる観測点データは、代表点の位置を時系列に表すデータとされたが、代表点の姿勢、速度、加速度、角速度、角加速度、力又はトルク等を時系列に表すデータであってもよいし、これらと位置とからなる群から選択される２以上の要素を時系列に表すデータを組み合わせたものであってもよい。なお、ここでいう姿勢とは、クォータニオンのような角度情報を有するパラメータである。
【０１０３】
＜５−５＞
上記実施形態では、ゴルファーの身体においてのみ観測点が設定されたが、観測点は、ゴルフクラブ上に設定してもよい。
【０１０４】
＜５−６＞
Ｉ0個の観測点を分けるＩ1個のグループには、１個の観測点を要素とする単一要素グループが含まれてもよい。この場合、単一要素グループに対応する観測点データは、当該単一要素グループに属する１個の観測点の挙動を時系列に表すデータとすることができる。
【０１０５】
＜５−７＞
上記実施形態では、同じ観測点が重複して複数のグループに分けられることがないように、観測点がグループ化された。しかしながら、１又は複数の観測点を複数のグループに重複して分けることも可能である。
【０１０６】
＜５−８＞
上記実施形態において、Ｉ0個の観測点を、Ｉ0個よりも少ないＩ1個のグループに分ける処理を省略してもよい。この場合、ステップＳ４において、全ての観測点に対するＩ0個の行列［ｒ1］，［ｒ2］，・・・，［ｒ25］をこの順に列方向（横方向）に並べた行列を、時系列挙動データ［Ｒ］として導出するようにすれば、その後の処理は、上記実施形態と同様に実行することができる。
【０１０７】
＜５−９＞
上記実施形態では、疑似拡張データ［Ｒa］が特異値分解されたが、時系列挙動データ［Ｒ］を特異値分解してもよい。
【０１０８】
＜５−１０＞
代表点の挙動は、代表点が代表する複数個の観測点の挙動を平均化することにより導出するのではなく、その他の方法で導出してもよい。代表点の挙動として、例えば、代表点が代表する複数個の観測点の挙動のうち最も大きいものを選択してもよいし、最も大きいものと最も小さいものを平均化する等してもよい。
【０１０９】
＜５−１１＞
第１実施形態のサブステップＳ１１では、ユーザからの入力に従って、注目モード、注目成分及び修正量が決定された。しかしながら、これらの少なくとも１つは、修正後の挙動が目標となる挙動に近づくように最適化により自動的に決定されてもよい。同様に、第２実施形態のサブステップＳ２１及び第３実施形態のサブステップＳ３１でも、ユーザからの入力に従って、注目モードが決定されたが、これも、修正後の挙動が目標となる挙動に近づくように最適化により自動的に決定されてもよい。
【０１１０】
＜５−１２＞
第２実施形態において、ペアモードとしては、相関係数ｃmaxを最大にするモードが選択されたが、これに限らず、注目モード以外のモードの中から任意に選択することができる。
【０１１１】
＜５−１３＞
第２実施形態のサブステップＳ２３において、注目モードの右特異ベクトルｚj及びペアモードの右特異ベクトルｚjを回転させるとき、前者が目標基底ｗjに近づくように回転させる必要はない。例えば、ユーザにより指定される任意の回転量だけ回転させてもよい。第３実施形態についても同様に、注目モードの右特異ベクトルｚj及びペアモードの右特異ベクトルｚjを、任意の回転量だけ回転させることができる。
【０１１２】
＜５−１４＞
第３実施形態のステップＳ３２では、ペアモードは、予め定められているルールに従って、注目モードがどのモードであるかに応じて選択された。しかしながら、第２実施形態と同様に、第１〜第Ｊモードのうち、注目モード以外の全てのモードを、ペアモードの候補となる候補モードとしてもよい。この場合、モード展開部１４ｂは、全ての候補モードに対し、順次、注目モードとの仮のペアリングを行う。そして、各仮のペアリングに対し、誤差関数Ｅを最小にする回転角度θrを算出し、このθrを数１６の式に代入し、誤差関数Ｅの値を算出する。そして、全ての候補モードのうち、最小の誤差関数Ｅの値を与える候補モードを、最終的にペアモードとして決定してもよい。
【０１１３】
＜５−１５＞
誤差関数Ｅの定義式は、数１４に示すものに限られず、例えば、下式のように、回転後の注目モードの右特異ベクトルｚtj'と目標基底ｗjとの誤差を表す項のみで構成してもよい。
【数１９】
2020185370.tif 000022
【０１１４】
なお、この場合、誤差関数Ｅを最小化するθは、下式により算出される。
【数２０】
2020185370.tif 000023
【符号の説明】
【０１１５】
１  スイング解析装置（コンピュータ）
１００  スイング解析システム
１４ａ  データ取得部
１４ｂ  モード展開部
１４ｃ  画面作成部
２  計測装置
６  スイング解析プログラム
７  ゴルファー
【図１】
2020185370.tif 000024
【図２】
2020185370.tif 000025
【図３】
2020185370.tif 000026
【図４】
2020185370.tif 000027
【図５】
2020185370.tif 000028
【図６】
2020185370.tif 000029
【図７】
2020185370.tif 000030
【図８】
2020185370.tif 000031
【図９】
2020185370.tif 000032
【図１０】
2020185370.tif 000033
【図１１】
2020185370.tif 000034
【図１２】
2020185370.tif 000035
【図１３】
2020185370.tif 000036
【図１４】
2020185370.tif 000037
【図１５Ａ】
2020185370.tif 000038
【図１５Ｂ】
2020185370.tif 000039
【図１６】
2020185370.tif 000040
【図１７】
2020185370.tif 000041
【図１８】
2020185370.tif 000042
【図１９】
2020185370.tif 000043
【図２０】
2020185370.tif 000044