【パワーメーター2019】VBAでトルク波形解析Pgm作った＜ＷorksheetFunction教訓得た＞

左クランクの基板は、３．７ＶＰＯ電池用電源をＴＩのTPS63000にすることにして実装をこつこつとやってます。７月末から８月お盆過ぎまで、家事等が多く活動のＤＵＹＴが低くなります。更に、９月８日の信越五高原ロングライドに参加する予定で練習に集中しないと完走できそうもないので、練習を最優先させます。

●波形解析方法の検討
トルク波形とフレームに各部ひずみ波形の相関係数をとって
相関係数が高いひずみ波形がトルクと比例関係が強いので
それを採用していくという考え方で一旦、解析してみます。
ＶＢＡで解析用プログラムを作ってみました。

●ＶＢＡ波形解析Ｐｇｍ
やり方：ＣＰＬＴでログした波形データを各ＣＨ間で相関係数を
求めます。何故ＶＢかとういうと、ＥＸＣＥＬの関数を
ＶＢＡプログラムで使えるので、自分で苦労して関数を
作らなくて済むのと、ＥＸＣＥＬの関数は、最速の計算速度
を提供してくれますので、自分でやるよりすべて良いです。
ＶＢＡでなくてもＶＢ．ＮＥＴでも対応できてますが、私が未だＶＢ．ＮＥＴに慣れないので、しばらくは、ＶＢＡメインで使います。
ＭＳの説明
https://docs.microsoft.com/ja-jp/office/vba/api/excel.worksheetfunction

ＥＸＣＥＬ関数をＶＢＡで使う使い方はこちらのＨＰが有名ですプロの方が
長年にわたって気づき上げたＨＰで、網羅的にＶＢＡ解説されてます。
判り易さを一貫しているので、ＶＢＡプログラムを作成するときは、こちらを
何回も利用させていただいております。

http://officetanaka.net/excel/vba/speed/s9.htm

ＴＮＡＫＡ様のＨＰで不足した場合は、詳細な解説のＨＰを探します。
https://www.moug.net/tech/exvba/0100035.html

WorksheetFunction.correl(range(cells(1,A),cells(x,A)),range(cellls(1,B),cells(x,B)))

ＶＢＡのソース　VBA-CPLT_Correl

ＵＩ：ＣＰＬＴで波形を見ながら、範囲値をＶＢＡに記入して、その範囲内の
相関係数を計算させます。
ＣＰＬＴのグラフのカーソルで比較したい波形のデータ範囲を決めてから
ＣＰＬＴのデータ番号をＶＢＡプログラムのフォームへ手入力します。
データ範囲をいれたら計算ボタンを押せば、Max,Min、相関係数がでます

おかしい？？　SeatTubeねじり波形とトルク波形が相関１となってます。
確かにピークの波形はぎざぎざ部分も一致してますが、谷側の形が違うのでこの相関係数１はおかしいですので、
＝＞どうもworksheet function で関数はトラブルが多いようで、ググると大量にでてきます。　関数パラメータが仕様通りになっているかチェックしてみます。
▼はまった事例
https://www.vba-ie.net/qanda2/qanda.cgi?mode=viewthread&id=24

●correlの仕様
ＭＳサイトから確認
https://docs.microsoft.com/ja-jp/office/vba/api/excel.worksheetfunction.correl

https://support.office.com/ja-jp/article/correl-%E9%96%A2%E6%95%B0-995dcef7-0c0a-4bed-a3fb-239d7b68ca92

型定義がＤＯＵＢＬＥにしてなかったのが原因でした。ＤＯＵＢＬＥに直したらＥＸＣＥＬの手計算と同じ値がでたので、バグＦＩＸしたと判断します。WorksheetFuncitonを使うときは、使うFunctionの仕様を入念にチェックしないと関係ないところでいろいろなエラーが発生してしまうという教訓です。ということでデバッグ後の計算結果は

●計算結果　　踏み始め３波形のサンプリング右クランクのみ

		相関係数
1-2ch	SeatTube 曲げ波形　vs 　SeatTubeねじり波形	-0.9672
1-3ch	SeatTube 曲げ波形　vs 　クランクトルク波形	0.9255
2-3ch	SeatTube ねじり波形　vs 　クランクトルク波形	-0.8333