pinescript版日銀の為替介入検知フーリエ逆変換ソースコード

承知いたしました。提示されたMQL4のコードを参考に、同様の機能を持つPine Script（TradingViewのスクリプト言語）のコードを記述します。

  ```pinescript //@version=5 indicator(title="Fourier Analysis (Pine)", shorttitle="Fourier Pine", overlay=false) // 入力パラメータ sampleSize = input.integer(256, title="Sample Size (Power of 2)", minval=64, maxval=4096, step=64) dominantComponents = input.integer(10, title="Dominant Components to Show", minval=1, maxval=20, step=1) noiseThreshold = input.float(0.05, title="Noise Threshold", minval=0.0, maxval=1.0, step=0.01) // バッファ（Pine Scriptでは直接的なバッファ操作は少ないですが、配列で代替します） var float originalBuffer = array.new_float(0) var float realPart = array.new_float(0) var float imagPart = array.new_float(0) var float magnitudes = array.new_float(0) var int componentIndices = array.new_int(0) var float componentWaves = array.new_float(dominantComponents, 0.0) // フーリエ変換関数 (簡易実装 - TradingViewにはFFTの組み込み関数がないため) // **注意:** これは教育目的の簡易的な実装であり、MQL4版のFFTと同等の精度や効率を持つわけではありません。 // TradingViewの制限により、完全なFFTアルゴリズムをPine Scriptで効率的に実装するのは困難です。 // より高度な分析には、外部ライブラリやAPIの利用を検討する必要があります。 function simpleFFT(data, n, inverse) if array.size(data) != n runtime.error("Data size must be equal to n") real = array.copy(data) imag = array.new_float(n, 0.0) for m = 1 to math.log2(n) step = pow(2, m) halfStep = step / 2 wPhaseStep = -2 * math.pi / step for k = 0 to halfStep - 1 wPhase = k * wPhaseStep wReal = math.cos(wPhase) wImag = math.sin(wPhase) for i = k to n - 1 step step j = i + halfStep tempReal = wReal * real.get(j) - wImag * imag.get(j) tempImag = wReal * imag.get(j) + wImag * real.get(j) real.set(j, real.get(i) - tempReal) imag.set(j, imag.get(i) - tempImag) real.set(i, real.get(i) + tempReal) imag.set(i, imag.get(i) + tempImag) endfor endfor endfor if inverse for i = 0 to n - 1 real.set(i, real.get(i) / n) imag.set(i, imag.get(i) / n) endfor endif [real, imag] // メイン処理 if bar_index >= sampleSize - 1 // 最新のサンプルデータを取得 array.clear(originalBuffer) for i = 0 to sampleSize - 1 array.push(originalBuffer, close[bar_index - sampleSize + 1 + i]) // フーリエ変換を実行 fftResult = simpleFFT(originalBuffer, sampleSize, false) realPart := array.get(fftResult, 0) imagPart := array.get(fftResult, 1) // 振幅を計算 array.clear(magnitudes) for i = 0 to sampleSize / 2 - 1 magnitude = math.sqrt(pow(array.get(realPart, i), 2) + pow(array.get(imagPart, i), 2)) array.push(magnitudes, magnitude) // 振幅でソートされたインデックスを取得 array.clear(componentIndices) for i = 0 to array.size(magnitudes) - 1 array.push(componentIndices, i) array.sort(componentIndices, (a, b) => array.get(magnitudes, b) - array.get(magnitudes, a)) // 上位成分のインデックスを保存 topIndices = array.slice(componentIndices, 0, math.min(dominantComponents, array.size(componentIndices))) // 各主要成分の波形を計算 array.clear(componentWaves) array.new_float(dominantComponents, 0.0) // サイズをリセット for i = 0 to sampleSize - 1 for j = 0 to array.size(topIndices) - 1 index = array.get(topIndices, j) if array.get(magnitudes, index) > noiseThreshold waveValue = array.get(realPart, index) * math.cos(2 * math.pi * index * i / sampleSize) - array.get(imagPart, index) * math.sin(2 * math.pi * index * i / sampleSize) array.set(componentWaves, j, array.get(componentWaves, j) + waveValue / sampleSize) endif endfor endfor // 波形をプロット plot(array.get(originalBuffer, array.size(originalBuffer) - 1), title="Original Price", color=color.red) for i = 0 to array.size(componentWaves) - 1 plot(array.get(componentWaves, i), title="Component " + (i + 1), color=color.blue, linewidth=1, offset=-(sampleSize - 1)) endfor // 色の定義 (TradingViewのplot関数内で直接指定することも可能です) // color componentColors[] = ... ```  

 Pine Script版の注意点:

 * FFT関数の簡易実装: TradingViewには組み込みのFFT関数が存在しないため、上記のコードでは非常に基本的な Cooley-Tukeyアルゴリズムのバタフライ演算の一部を簡略化した関数 simpleFFT を実装しています。この実装は、MQL4版のFFTと同等の精度や効率を持つわけではない ことを強く理解しておいてください。特に、大きな sampleSize では計算に時間がかかる可能性があります。より高度なフーリエ解析を行う場合は、TradingViewの制限を超えるため、外部ツールやAPIの利用を検討する必要があります。 

 * リアルタイム性の制約: Pine Scriptは、MQL4と同様に、複雑な計算をリアルタイムに行うには制約があります。特に、上記の簡易FFT関数は計算負荷が高いため、チャートの更新頻度に影響を与える可能性があります。 

 * バッファの扱い: Pine Scriptでは、MQL4のような直接的なバッファ操作は一般的ではありません。代わりに、配列を使用してデータを管理します。 

 * 上位成分の抽出: MQL4版と同様に、振幅の大きい順に上位の成分を選択しています。 

 * 成分波形の計算とプロット: 各上位成分に対応する正弦波を計算し、元の価格と合わせてプロットしています。offset パラメータを使用して、過去のデータに対応するようにプロット位置を調整しています。 

 * 色の指定: plot() 関数内で直接色を指定しているため、MQL4版のような色の配列はここでは直接的には使用していません。必要に応じて、色の配列を定義して使用することも可能です。 

改善の方向性 (Pine Script): 

 * より効率的なFFTの実装の検討: もし高度なフーリエ解析が必要な場合は、Pine Scriptの範囲内で可能な限り効率的なFFTアルゴリズムを調査・実装する必要があります。ただし、TradingViewの制限により限界があることを理解しておく必要があります。 

 * 外部連携の検討: TradingViewのWebhook機能などを利用して、外部の高性能な数値計算ライブラリ（PythonのSciPyなど）でフーリエ解析を行い、その結果をTradingViewに送ってプロットすることを検討するのも一つの方法です。ただし、リアルタイム性は損なわれます。 

 * TradingViewの組み込み関数の活用: TradingViewには、移動平均や標準偏差など、様々な組み込みのテクニカル指標関数があります。これらの関数を組み合わせることで、フーリエ解析の代替となるような分析アプローチを検討することもできます。 このPine Scriptコードは、MQL4版の基本的なアイデアをTradingView上で実現しようとしたものですが、特にFFTの実装においては簡略化されている点に注意してください。実際の利用においては、パフォーマンスや精度を確認し、必要に応じて調整や代替手段の検討を行ってください。