※重いのでTurboWarp推奨 https://turbowarp.org/1287783974/editor?fps=240&hqpen&turbomode 参考にしたもの： リアルタイムの加速度記録から長周期地震動階級を計算する https://qiita.com/NoneType1/items/cfd0f9ea2802c9c383a5 使い方： 0. 加速度データをもらう 気象庁｜強震観測データ https://www.data.jma.go.jp/eqev/data/kyoshin/jishin/index.html で加速度のCSVファイルがダウンロードできます 1. 加速度データを入れる 1-1. 中を見る 1-2. 緑の旗を押す 1-3. リスト「元データ(csv)」を右クリック→読み込み→CSVファイルを選択 2. 長周期地震動階級計算 2-1. ターボモードにする 2-2. 緑の旗を押してリアルタイム震度が計算される その他の操作 1キー: N-S成分を波形を表示 2キー: E-W成分の波形を表示 3キー: U-D成分の波形を表示 スペースキー: 1画面に収まるように表示倍率を自動調整 Dキー: すべてリセット K-NET ASCII形式のデータも一応読み込めます(未検証) 0. もしプロジェクトが実行中なら停止する 1. Aキーを押してリストが3つ表示される 2. .EWと.NSと.UDの元データをリストに読み込む 3. もう一度Aキーを押して加速度が自動変換される 4. 緑の旗を押して長周期地震動階級が計算される 内部計算の手順 1. 台形積分を使って現在の地面の速度を計算 2. 現在の地面加速度の値と一ターン前の地面加速度の値と一ターン前の地面との相対変位と一ターン前の地面との相対速度と各定数(A,B)を使って現在の相対速度を計算 3. 現在の地面との相対速度と現在の地面の速度を足して現在の絶対速度を取得 4. 1~3を周期秒を変更(1.6~7.8秒まで0.2秒毎)して計算する 5. 1~4を毎ターン繰り返しその中の最大値を対応表に入れて階級を求める 多分こんな感じです で本来は地震ごとに一つの値を示す長周期地震動階級ですがここではリアルタイムで求めることを目的に時間軸ごとに各周期の最大を取得し、過去n秒間のデータの中の最大を使用しています
