Microsoft PowerPoint 集い横田.ppt [互換モード]

Size: px
Start display at page:

Download "Microsoft PowerPoint 集い横田.ppt [互換モード]"

Transcription

1 波浪推算モデルによる海面抵抗係数の推定に関する研究 九州大学大学院工学研究院助教新キャンパス計画推進室勤務横田雅紀

2 第 世代波浪推算モデル (WAM) 波浪推算 : 風情報を入力値として波高を予測 有義波法 : 風速, 吹送時間, 吹送距離と波高 周期の関係 スペクトル法 : 時空間的に変化する風場での方向別周波数別成分波 エネルギー平衡方程式 E( f t, θ ) + Cg E( f, θ ) = S 発達 S = Sin + Sds + 減衰相互作用 S nl 既往の擾乱事例に対して 概ね観測された波高を再現可能 温暖化により 台風が強大化すると

3 風から波へのエネルギー輸送 Ocean Sin = BE Miles の不安定機構 Janssen B = ρa β ρ W U* max cos C 強風時に波高が過大 [ ( )], θ θ σ W 海面からの高度 z における風速分布 U(z) U U + * e o ( z) = ln κ ze z e z z zo = U τ W τ =. g z * z o U * : 摩擦速度 κ: カルマン定数 z e : 有効水面粗度 z : 粗度長

4 海面抵抗係数 Ocean C D U Wu 98 * = CDU Drag Coefficient C (* - ) D Honda & Mitsuyasu 98 Andreas (spray) Wu Honda & Mitsuyasu Wu Andreas Mitsuyasu & Honda 5 5 Wind Speed U Wind Speed U (m/s) (m/s) Aerodynamic Observation Powell et al,

5 波浪推算結果の比較 Ocean Wave Height Hs(m) Janssen Wu Honda & Mitsuyasu Andreas Wilson 風向 風速一定 フェッチ 5km Wind Speed U (m/s) Wilson gh / = + gf.. U U / 5

6 強風条件における海面抵抗係数の推定 海面上の風速分布から推定 風速 m/s 以上の条件を狙った観測暴風 暴浪を再現しての室内実験 費用面で困難 理論も成り立っているのか疑問 C D (U ) データ同化が有効 海面抵抗係数 波浪推算モデルを間接的な観測装置とし波浪観測値を満たす最適な係数を推定 風速 U 6

7 データ同化手法 ADWAM (Hersbach 998) C D (U ) 離散的一定値関数 C D st 初期値 nd ~ 修正 Run Backward Descent Method 最適値 WAM 推算波高 評価関数 Adjoint WAM J (x) 海上風 x x n- Run U 観測波高 観測誤差 J ( x) = T t= ( H ( x) y ) t + W t N Τ n= R t ( x t x n ( H ( x) y ) n x t x N Τ n ) Βt ( xn xn ) 背景誤差 ( 先験条件 )

8 双子実験 5 t=6h OBS point Wind Speed(m/s) OBS point t=h typhoon course 度 度格子間隔 :.5 度 85hPa : 中心気圧 km : 最大風速半径 5km/h : 移動速度 5 6 Time(h) 観測地点で風速 5~5m/s が発生この範囲の CD はデータ同化で推定出来ると期待される 8

9 実験の手順 Drag Coefficient C D ( - ) 変数の目標値を設定 目標値を変数とする波浪推算値 ( 観測 ) 任意の初期値を仮定 データ同化により推算値に合う変数を探索 TargetValue InitialValue DeducedValue 5 WindSpeed(m/s) Height Hs(m) Wave WAM ADWAM ObservedData BeforeAssim. AfterAssim. 5 6 Time(hr) 9

10 妥当性の検証 Division Number of Wind Speed = Number of Unknown Parameter Drag Coefficient C D (* - ) Target Value Initial Value Deduced Value 風速間隔.5m/s ( パラメータ数 ) Wave Data : 66 目標値周辺の値を推定不安定

11 精度検証 Drag Coefficient C D (* - ) Drag Coefficient C D (* - ) 5 Wind Speed U (m/s) Target Value Initial Value Deduced Value Target Value Initial Value Deduce Value 風速間隔.5m/s パラメータ数 5 Wind Speed U (m/s) 先験条件なし W= 風速間隔.m/s ( パラメータ数 5) 先験条件あり W=

12 数値実験 強風による波浪の発達 風域を外れたうねりの時系列データ

13 計算条件 無風域 Wind Speed (m/s) 5 a b c d.8m/s.5m/s 台風モデル 85hPa : 中心気圧 5km/h : 移動速度 wave height H s (m) d d d d d5 Time (hr) 5 6 time (h) 波高時系列

14 台風経路から外れた海域の波高による推定 Ocean Drag Coefficient Drag Coefficient Drag Coefficient Drag Coefficient Drag Coefficient ficient C D (* - ) C D (* - ) C D (* - ) C D (* - ) C D (* - ) Target Value Initial Value Deduced Value a Maximum Wind Speed Target Value Initial Value Deduced Value b Maximum Wind Speed Target Value Initial Value Deduced Value c Maximum Wind Speed Target Value Initial Value Deduced Value e Maximum Wind Speed Target Value Initial Value Deduced Value f a b c e f Maximum Wind Speed 5 各海域の最大風速まで高精度

15 無風域の波高データによる推定 Drag Coefficient Drag Coefficient Drag Coefficient Drag Coefficient C (* - ) C (* - ) C (* - ) D D D C (* - ) D.5.5 Target Value Initial Value Deduced Value d Target Value Initial Value Deduced Value d Target Value Initial Value Deduced Value d 5 Target Value Initial Value Deduced Value d 5 波高減衰後の地点でも可 5

16 海域 Oceana の波高修正状況 Drag Coefficient wave height H s (m) C D (* - ) Target Value Initial Value Deduced Value a Maximum Wind Speed 5 station a with Target Value (observed data) with Initial Value (before assim.) with Deduced Value (after assim.) 5 6 time (h) a 観測情報が不足 6

17 現地適用 Ocean Honda&Mitsuyasu NOWPHAS 月 /9 7 月 / 7 月 / 7 月 / 8 月 / 8 月 / 8 月 / 8 月 / 8 月 / 8 月 /5.. NOWPHAS. - ) 海面抵抗係数 ( Honda& Mitsuyasu Honda&Mitsuyasu

18 波浪観測値を用いたデータ同化 対象擾乱台風 号台風 号台風 号台風 5 号台風 7 号台風 75 号 項目 計算諸元 領域 緯度.5 ~5 経度 5 ~5 格子幅.5 計算格子数 5 6 TimeStep 6s パラメータ数 5 個 (.m/s 間隔 ) 同化繰り返し回数 回

19 逆推定結果 観測波高観測波高初期波高初期波高同化後波高 月 /8.5 5 月 /9 5 月 / 5 月 / 5 月 / 5 月 / 5 月 / 5 月 / 観測波高観測波高初期波高初期波高同化後波高.5 5 月 /8 5 月 /9 5 月 / 5 月 / 5 月 / 5 月 / 5 月 / 5 月 / 観測波高初期波高観測波高同化後波高初期波高 月 /8 5 月 /9 5 月 / 5 月 / 5 月 / 5 月 / 5 月 / 5 月 / A B C - ) 海面抵抗係数 ( - ) 海面抵抗係数 ( - ) 海面抵抗係数 ( 海面抵抗係数 ( - ) 初期値逆推定値 初期値逆推定値 初期値逆推定値 A A B C Honda&Mitsuyasu B C

20 結論 同化後の波高は観測波高に近い値に修正され, 海面抵抗係数は初期値によらず本多 光易の提案式に近い値に修正された ADWAM の有用性が確認できたといえる 今後,m/s 以上の強風速を受けて発達した波浪の観測データが収集できれば,ADWAM により強風速範囲の海面抵抗係数が解明できるものと期待される.

21 現地適用に向けての課題 Ocean パラメータ数の設定 ( 風速範囲 風速間隔 ) 現地風データ 地形を用いた双子実験により 推定可能な風速範囲 風速間隔は把握出来そう 先験条件の重み 修正回数 観測誤差に対する脆弱性 観測情報 ( 同化期間 地点数 波高以外の情報 ) 入力風データの精度 ( 時間間隔 格子間隔 台風コース 到達時刻 )

22 Ocean 観測波高初期波高同化後波高 海面抵抗係数 ( - ).... 初期値 逆推定値 T T. 9 月 /5 9 月 /6 9 月 /7 9 月 /8 9 月 /9 9 月 / 月 / 月 / 観測波高初期波高同化後波高 T 海面抵抗係数 ( - ).... 初期値 逆推定値 T. 9 月 /5 9 月 /6 9 月 /7 9 月 /8 9 月 /9 9 月 / 月 / 月 /

23 Ocean 観測波高初期波高同化後波高 海面抵抗係数 ( - ).... 初期値 逆推定値 T. 月 /6 月 /7 月 /8 月 /9 月 / 月 / 月 / 月 / 観測波高初期波高同化後波高 海面抵抗係数 ( - ).... 初期値 逆推定値 T. 月 /6 月 /7 月 /8 月 /9 月 / 月 / 月 / 月 /

24 Ocean 観測波高初期波高同化後波高 海面抵抗係数 ( - ).... 初期値 逆推定値 T5. 8 月 / 8 月 / 8 月 / 8 月 / 8 月 / 8 月 /5 8 月 /6 8 月 /7 8 月 / 観測波高初期波高同化後波高 海面抵抗係数 ( - ).... 初期値 逆推定値 T5. 8 月 / 8 月 / 8 月 / 8 月 / 8 月 / 8 月 /5 8 月 /6 8 月 /7 8 月 /

25 その他のケース T 観測波高初期波高同化後波高 T. 7 月 / 7 月 / 7 月 / 7 月 / 7 月 / 7 月 /5 7 月 /6 7 月 /7 7 海面抵抗係数 ( - ).... 初期値 逆推定値 T7 T 観測波高初期波高同化後波高. 7 月 / 7 月 / 7 月 / 7 月 / 7 月 / 7 月 /5 7 月 /6 7 月 /7 7 海面抵抗係数 ( - ).... 初期値 逆推定値 T

26 その他のケース T 観測波高初期波高同化後波高 海面抵抗係数 ( - ).... 初期値 逆推定値 T 月 /9 7 月 / 7 月 / 8 月 / 8 月 / 8 月 / 8 月 / 8 月 / 観測波高初期波高同化号波高 海面抵抗係数 ( - ).... 初期値 逆推定値 T75. 7 月 /9 7 月 / 7 月 / 8 月 / 8 月 / 8 月 / 8 月 / 8 月 /

Microsoft Word _資料2(最終版).doc

Microsoft Word _資料2(最終版).doc 資料 -2 今回の潮位の特性 ( 追加 ) 2.1 被災時の潮位特性 (1) 概要被災時の潮位特性について 観測潮位データを基に調べた 先ず 28 年 2 月 24 日を含む 日本海沿岸の検潮データを収集して整理した 先ず 小樽 深浦 佐渡 富山 能登 舞鶴 境 浜田 ( いずれも気象庁 ) の8 地点の2 月の観測潮位と潮位偏差を収集して整理した 8 地点の位置を下図に示した 次に 伏木富山 新湊

More information

<4D F736F F D208E9197BF ED8DD08E9E82CC8B438FDB8A438FDB F4390B3816A>

<4D F736F F D208E9197BF ED8DD08E9E82CC8B438FDB8A438FDB F4390B3816A> 被災時の気象 海象状況について 1. 台風 23 号の概況 1.1 台風 23 号の経路と気圧変化 1.2 既往の台風との比較 423 号台風 (TOKAGE) の来襲により ナウファス * の室津波浪観測計で波高 H1/3=13.m 周期 T1/3=.8sec( 速報値 ) の有義波が記録された この値は 1 年の観測開始以来最大の有義波である / : /21 9 時 : 温帯低気圧に変わる /

More information

風力発電インデックスの算出方法について 1. 風力発電インデックスについて風力発電インデックスは 気象庁 GPV(RSM) 1 局地気象モデル 2 (ANEMOS:LAWEPS-1 次領域モデル ) マスコンモデル 3 により 1km メッシュの地上高 70m における 24 時間の毎時風速を予測し

風力発電インデックスの算出方法について 1. 風力発電インデックスについて風力発電インデックスは 気象庁 GPV(RSM) 1 局地気象モデル 2 (ANEMOS:LAWEPS-1 次領域モデル ) マスコンモデル 3 により 1km メッシュの地上高 70m における 24 時間の毎時風速を予測し 風力発電インデックスの算出方法について 1. 風力発電インデックスについて風力発電インデックスは 気象庁 GPV(RSM) 1 局地気象モデル 2 (ANEMOS:LAWEPS-1 次領域モデル ) マスコンモデル 3 により 1km メッシュの地上高 70m における 24 時間の毎時風速を予測し 2000kW 定格風車の設備利用率として表示させたものです 数値は風車の定格出力 (2000kW)

More information

85 4

85 4 85 4 86 Copright c 005 Kumanekosha 4.1 ( ) ( t ) t, t 4.1.1 t Step! (Step 1) (, 0) (Step ) ±V t (, t) I Check! P P V t π 54 t = 0 + V (, t) π θ : = θ : π ) θ = π ± sin ± cos t = 0 (, 0) = sin π V + t +V

More information

A

A A04-164 2008 2 13 1 4 1.1.......................................... 4 1.2..................................... 4 1.3..................................... 4 1.4..................................... 5 2

More information

2009 年 11 月 16 日版 ( 久家 ) 遠地 P 波の変位波形の作成 遠地 P 波の変位波形 ( 変位の時間関数 ) は 波線理論をもとに P U () t = S()* t E()* t P() t で近似的に計算できる * は畳み込み積分 (convolution) を表す ( 付録

2009 年 11 月 16 日版 ( 久家 ) 遠地 P 波の変位波形の作成 遠地 P 波の変位波形 ( 変位の時間関数 ) は 波線理論をもとに P U () t = S()* t E()* t P() t で近似的に計算できる * は畳み込み積分 (convolution) を表す ( 付録 遠地 波の変位波形の作成 遠地 波の変位波形 ( 変位の時間関数 ) は 波線理論をもとに U () t S() t E() t () t で近似的に計算できる は畳み込み積分 (convolution) を表す ( 付録 参照 ) ここで St () は地震の断層運動によって決まる時間関数 1 E() t は地下構造によって生じる種々の波の到着を与える時間関数 ( ここでは 直達 波とともに 震源そばの地表での反射波や変換波を与える時間関数

More information

Freak wave estimation by WAVEWATCH III and HOSM

Freak wave estimation by WAVEWATCH III and HOSM 2017 年 9 月 25 日第 2 回理研データ同化ワークショップ アンサンブル 4 次元変分法を活用した巨大波浪の再現 藤本航 早稲田卓爾 東京大学新領域創成科学研究科海洋技術環境学専攻 1 導入 : フリーク波とは 海洋に突発的に生じる巨大波 波浪場の平均的な振幅 (H s : 有義波高 ) を大きく超える波 H > 2 or η crest > 1.3 H s H s (Dysthe et

More information

スライド 1

スライド 1 非線形数理秋の学校 パターン形成の数理とその周辺 - 反応拡散方程式理論による時 空間パターンの解析を中心に - 2007 年 9 月 25 日 -27 日 モデル方程式を通してみるパターン解析ー進行波からヘリカル波の分岐を例としてー 池田勉 ( 龍谷大学理工学部 ) 講義概要, 講義資料, 講義中に使用する C 言語プログラムと初期値データ, ヘリカル波のアニメーションをウェブで公開しています :

More information

Microsoft Word - 0_0_表紙.doc

Microsoft Word - 0_0_表紙.doc 2km Local Forecast Model; LFM Local Analysis; LA 2010 11 2.1.1 2010a LFM 2.1.1 2011 3 11 2.1.1 2011 5 2010 6 1 8 3 1 LFM LFM MSM LFM FT=2 2009; 2010 MSM RMSE RMSE MSM RMSE 2010 1 8 3 2010 6 2010 6 8 2010

More information

Microsoft PowerPoint - 海の波.ppt [互換モード]

Microsoft PowerPoint - 海の波.ppt [互換モード] 008/5/6 風波の発生 ( 深海域 波浪について ( 初級編 波の浅水変形と砕波 波による質量輸送と海浜流系の発達 離岸流の発生 ( 浅海域 3.6m の高波が襲来 水粒子の軌道 ( スト - クスドリフト 離岸流 沿岸流 (1 波 浅海域の波浪変形の研究 1. 海岸保全に関して (A 季節風時 台風時の波と流れによる地形変化 ( 侵食 (B 漁港 港湾 海岸保全構造物による波 流れ 地形変化

More information

An ensemble downscaling prediction experiment of summertime cool weather caused by Yamase

An ensemble downscaling prediction experiment of summertime cool weather caused by Yamase ヤマセによる冷夏をターゲットにした アンサンブルダウンスケール予報実験 東北大学 福井真 1. Introduction 1.1 ヤマセ 海洋性極気団を起源とした冷湿な北東風 => 水平規模 ~1000 kmの現象 (Kodama et al. 2009) 冷夏となり 農作物に大きな被害 ( 冷害 ) をもたらすことも => 重要な中期予報の対象 背が低く 複雑な地形の影響を大きく受ける ( 工藤

More information

Microsoft PowerPoint - 2_6_shibata.ppt [互換モード]

Microsoft PowerPoint - 2_6_shibata.ppt [互換モード] 圧密問題への逆問題の適用 一次元圧密と神戸空港の沈下予測 1. 一次元圧密の解析 2. 二次元圧密問題への適用 3. 神戸空港の沈下予測 1. 一次元圧密の解析 一次元圧密の実験 試験システムの概要 分割型圧密試験 逆解析の条件 未知量 ( 同定パラメータ ) 圧縮指数 :, 透水係数 :k 初期体積ひずみ速度 : 二次圧密係数 : 観測量沈下量 ( 計 4 点 ) 逆解析手法 粒子フィルタ (SIS)

More information

Microsoft Word - ME-webGIS使用マニュアル_160615

Microsoft Word - ME-webGIS使用マニュアル_160615 Marine Energy-webGIS(ME-webGIS) 使用マニュアル 2016 年 6 月 15 日版 目次 1 ME-webGIS へのアクセス方法... 1 2 波浪時系列情報... 3 2.1 NOWPAHS 観測地点波浪時系列表示... 3 2.2 波浪エネルギー発電の主要候補地波浪時系列表示... 12 3 波浪情報マップ... 19 3.1 波浪情報マップ ( 気候値 )...

More information

B2 Vol. B2-65 No Construction of a Data Set of 45-Year Sea Wind Distribution on the Inner Bay and Inland Sea of Japan Using SDP Winds M

B2 Vol. B2-65 No Construction of a Data Set of 45-Year Sea Wind Distribution on the Inner Bay and Inland Sea of Japan Using SDP Winds M B2 Vol. B2-65 No.1 2009 186-190 Construction of a Data Set of 45-Year Sea Wind Distribution on the Inner Bay and Inland Sea of Japan Using SDP Winds Measurement Data Masataka YAMAGUCHI, Manabu OHFUKU,

More information

No.18SP1-4 Reports of RIAM Symposium No.18SP1-4 Study on features and generation mechanisms of freak waves II Proceedings of a symposium held at Resea

No.18SP1-4 Reports of RIAM Symposium No.18SP1-4 Study on features and generation mechanisms of freak waves II Proceedings of a symposium held at Resea No.8SP-4 Reports of RIAM Symposium No.8SP-4 Study on features and generation mechanisms of freak waves II Proceedings of a symposium held at Research Institute for Applied Mechanics, Kyushu Universiy,

More information

スライド 1

スライド 1 移動体観測を活用した交通 NW の リアルタイムマネジメントに向けて : プローブカーデータを用いた動的 OD 交通量のリアルタイム推定 名古屋大学山本俊行 背景 : マルチモード経路案内システム PRONAVI 2 プローブカーデータの概要 プローブカー : タクシー 157 台 蓄積用データ収集期間 : 22 年 1 月 ~3 月,1 月 ~23 年 3 月 データ送信はイベントベース : 車両発進

More information

内湾流動に及ぼす大気の影響 名古屋大学村上智一

内湾流動に及ぼす大気の影響 名古屋大学村上智一 内湾流動に及ぼす大気の影響 名古屋大学村上智一 研究背景 沿岸域の海水流動の海水流動はは, 風による吹送流, 日射による成層化, 降水 蒸発などの気象場からの影響気象場からの影響を強く受ける. そのため, 沿岸域の海水流動計算を高精度で行うには, 気象場からの影響を適切に評価することが必須となる. 水収支 熱収支 降水蒸発 短波放射 運動量収支 風 対流顕熱潜熱長波放射 海流 気象場からの影響を正しく評価するためには,

More information

フィードバック ~ 様々な電子回路の性質 ~ 実験 (1) 目的実験 (1) では 非反転増幅器の増幅率や位相差が 回路を構成する抵抗値や入力信号の周波数によってどのように変わるのかを調べる 実験方法 図 1 のような自由振動回路を組み オペアンプの + 入力端子を接地したときの出力電圧 が 0 と

フィードバック ~ 様々な電子回路の性質 ~ 実験 (1) 目的実験 (1) では 非反転増幅器の増幅率や位相差が 回路を構成する抵抗値や入力信号の周波数によってどのように変わるのかを調べる 実験方法 図 1 のような自由振動回路を組み オペアンプの + 入力端子を接地したときの出力電圧 が 0 と フィードバック ~ 様々な電子回路の性質 ~ 実験 (1) 目的実験 (1) では 非反転増幅器の増幅率や位相差が 回路を構成する抵抗値や入力信号の周波数によってどのように変わるのかを調べる 実験方法 図 1 のような自由振動回路を組み オペアンプの + 入力端子を接地したときの出力電圧 が 0 となるように半固定抵抗器を調整する ( ゼロ点調整のため ) 図 1 非反転増幅器 2010 年度版物理工学実験法

More information

<4D F736F F F696E74202D208A438ADD8BDF82AD82CC97AC82EA82C697A48B4E8CB A82E B8CDD8AB B83685D>

<4D F736F F F696E74202D208A438ADD8BDF82AD82CC97AC82EA82C697A48B4E8CB A82E B8CDD8AB B83685D> 海岸近くの流れと陸起源漂着物 の移動メカニズム 西隆一郎 ( 鹿児島大学 水産学部 環境情報科学講座 ) 専門 ; 沿岸域の波浪 流れ 砂移動 航路保全 自然災害 環境アセス 沿岸域の海底地形と流れ 漂着物の漂着空間分布 砂浜はどのようにしてできる? 天然材料でできた海岸以外に漂着物 ( ゴミ ) でできた海岸が出現するのでは? カリフォルニア州には廃棄ガラスでできた海岸がある! 海岸の砂浜の構成材料

More information

領域シンポ発表

領域シンポ発表 1 次元の減衰運動の中の強制振動 ) ( f d d d d d e f e ce ) ( si ) ( 1 ) ( cos ω =ω -γ とおくと 一般解は 外力 f()=f siω の場合 f d d d d si f ce f ce si ) cos( cos si ) cos( この一般解は 1 φ は外力と変位との間の位相差で a 時間が経つと 第 1 項は無視できる この場合の振幅を

More information

(00)顕著速報(表紙).xls

(00)顕著速報(表紙).xls 対象地域 : 八重山地方 ( 速報 ) 平成 24 年 台風第 14 号について 目 次 資料 1 気象概況 資料 2 台風経路図 資料 3 地上天気図 沿岸波浪実況図 資料 4 気象衛星画像 気象レーダー 資料 5 暴風域 強風域入出日時 資料 6 台風期間における観測値 資料 7 気象警報 注意報発表状況 資料 8 時系列図 ( 降水量 風向 風速 現地気圧 ) 資料 9 観測値表 ( 降水量

More information

Microsoft Word - 01.doc

Microsoft Word - 01.doc 国立研究開発法人海洋研究開発機構国立大学法人京都大学 エルニーニョ予測の新展開 ~ 春先からの予測精度向上に新たな可能性 ~ 1. 概要国立研究開発法人海洋研究開発機構 ( 理事長平朝彦 以下 JAMSTEC という) 地球環境観測研究開発センター海洋循環研究グループの増田周平グループリーダーらは 京都大学と共同で力学解析に基づいたエルニーニョ現象の新しい予測手法を考案しました 猛暑や旱魃 豪雨など社会的に影響の大きな異常気象を各地で引き起こすエルニーニョ現象の精度の高い予測方法の開発は

More information

橡Ⅰ.企業の天候リスクマネジメントと中長期気象情

橡Ⅰ.企業の天候リスクマネジメントと中長期気象情 1 1 2 1 2 2 3 4 4 3 4 3 5 1400 53 8.8 11 35 6 5 6 20012Q 926 1,438 15.032.2 4 ART 7 8 9 7 8 9 5 19712000 30 33 60 10 33 10 60 70 30 40 6 12 3000 2000 7 沈降した後 付近の流れに乗って海中を水平に漂流するように設計されている その後 予め設定した時間間隔

More information

A Precise Calculation Method of the Gradient Operator in Numerical Computation with the MPS Tsunakiyo IRIBE and Eizo NAKAZA A highly precise numerical

A Precise Calculation Method of the Gradient Operator in Numerical Computation with the MPS Tsunakiyo IRIBE and Eizo NAKAZA A highly precise numerical A Precise Calculation Method of the Gradient Operator in Numerical Computation with the MPS Tsunakiyo IRIBE and Eizo NAKAZA A highly precise numerical calculation method of the gradient as a differential

More information

Microsoft Word doc

Microsoft Word doc . 正規線形モデルのベイズ推定翠川 大竹距離減衰式 (PGA(Midorikawa, S., and Ohtake, Y. (, Attenuation relationships of peak ground acceleration and velocity considering attenuation characteristics for shallow and deeper earthquakes,

More information

1校_本扉.indd

1校_本扉.indd 第 回風工学シンポジウム () 台風時の風向特性と複雑地形の増速特性を考慮したの評価手法の提案 Prediction of Wind Speedup Ratio Considering Wind Direction of yphoon and Speed-up of Complex errain 菊地由佳 石原孟 Yuka KIKUCHI ) akeshi ISHIHARA ) ABSRAC It

More information

0 スペクトル 時系列データの前処理 法 平滑化 ( スムージング ) と微分 明治大学理 学部応用化学科 データ化学 学研究室 弘昌

0 スペクトル 時系列データの前処理 法 平滑化 ( スムージング ) と微分 明治大学理 学部応用化学科 データ化学 学研究室 弘昌 0 スペクトル 時系列データの前処理 法 平滑化 ( スムージング ) と微分 明治大学理 学部応用化学科 データ化学 学研究室 弘昌 スペクトルデータの特徴 1 波 ( 波数 ) が近いと 吸光度 ( 強度 ) の値も似ている ノイズが含まれる 吸光度 ( 強度 ) の極大値 ( ピーク ) 以外のデータも重要 時系列データの特徴 2 時刻が近いと プロセス変数の値も似ている ノイズが含まれる プロセス変数の極大値

More information

Microsoft PowerPoint - 資料04 重回帰分析.ppt

Microsoft PowerPoint - 資料04 重回帰分析.ppt 04. 重回帰分析 京都大学 加納学 Division of Process Control & Process Sstems Engineering Department of Chemical Engineering, Koto Universit manabu@cheme.koto-u.ac.jp http://www-pse.cheme.koto-u.ac.jp/~kano/ Outline

More information

重力方向に基づくコントローラの向き決定方法

重力方向に基づくコントローラの向き決定方法 ( ) 2/Sep 09 1 ( ) ( ) 3 2 X w, Y w, Z w +X w = +Y w = +Z w = 1 X c, Y c, Z c X c, Y c, Z c X w, Y w, Z w Y c Z c X c 1: X c, Y c, Z c Kentaro Yamaguchi@bandainamcogames.co.jp 1 M M v 0, v 1, v 2 v 0 v

More information

ポリトロープ、対流と輻射、時間尺度

ポリトロープ、対流と輻射、時間尺度 宇宙物理学 ( 概論 ) 6/6/ 大阪大学大学院理学研究科林田清 ポリトロープ関係式 1+(1/) 圧力と密度の間にP=Kρ という関係が成り立っていると仮定する K とは定数でをポリトロープ指数と呼ぶ 5 = : 非相対論的ガス dlnp 3 断熱変化の場合 断熱指数 γ, と dlnρ 4 = : 相対論的ガス 3 1 = の関係にある γ 1 等温変化の場合は= に相当 一様密度の球は=に相当

More information

11-1 最近の地震観測の精度 ~気象庁における地震観測業務~

11-1 最近の地震観測の精度 ~気象庁における地震観測業務~ 11-1 最近の地震観測の精度 ~ 気象庁における地震観測業務 ~ The precision of recent seismic observation. - Observation and data processing by JMA - 原田智史 ( 気象庁地震予知情報課 ) Satoshi HARADA (Earthquake Prediction Information Division,

More information

資料 2-3 平成 28 年 11 月 21 日火力部会資料 相馬港天然ガス発電所 ( 仮称 ) 設置計画 環境影響評価準備書 補足説明資料 平成 28 年 11 月 福島ガス発電株式会社 1

資料 2-3 平成 28 年 11 月 21 日火力部会資料 相馬港天然ガス発電所 ( 仮称 ) 設置計画 環境影響評価準備書 補足説明資料 平成 28 年 11 月 福島ガス発電株式会社 1 資料 2-3 平成 28 年 11 月 21 日火力部会資料 相馬港天然ガス発電所 ( 仮称 ) 設置計画 環境影響評価準備書 補足説明資料 平成 28 年 11 月 福島ガス発電株式会社 1 火力部会補足説明資料目次 1. 煙突高さを変更した経緯 理由について 3 2. ベキ法則を用いた上層気象の設定について 9 3. 準備書の記載の誤りについて 14 2 1. 煙突高さを変更した経緯 理由について

More information

5. 数値解析 5.2. サンゴ浮遊幼生ネットワークモデルの検討

5. 数値解析 5.2. サンゴ浮遊幼生ネットワークモデルの検討 5. 数値解析 5.2. サンゴ浮遊幼生ネットワークモデルの検討 目 次 5.2. サンゴ浮遊幼生ネットワークモデルの検討 V-5-2-15 1) 追跡計算の概要 V-5-2-15 2) 追跡手法 V-5-2-16 3) 追跡モデル 追跡手法の妥当性確認 V-5-2-17 (1) 慶良間列島 ~ 沖縄本島西岸 V-5-2-17 (2) 石西礁湖 ~ 沖縄本島 V-5-2-18 4) 追跡計算結果

More information

Microsoft PowerPoint - 01_内田 先生.pptx

Microsoft PowerPoint - 01_内田 先生.pptx 平成 24 年度 SCOPE 研究開発助成成果報告会 ( 平成 22 年度採択 ) 塩害劣化した RC スラブの一例 非破壊評価を援用した港湾コンクリート構造物の塩害劣化予測手法の開発 かぶりコンクリートのはく落 大阪大学大学院鎌田敏郎佐賀大学大学院 内田慎哉 の腐食によりコンクリート表面に発生したひび割れ ( 腐食ひび割れ ) コンクリート構造物の合理的な維持管理 ( 理想 ) 開発した手法 点検

More information

Microsoft PowerPoint _HARU_Keisoku_LETKF.ppt [互換モード]

Microsoft PowerPoint _HARU_Keisoku_LETKF.ppt [互換モード] を用いた大阪平野南部で 発達した雷雨の再現実験 ( のネストシステムを目指して ) 瀬古弘 露木義 斉藤和雄 ( 気象研究所 ) 黒田徹 ( 海洋研究開発機構 ) 藤田匡 ( 気象庁 ) 三好建正 ( メリーランド大 ) を用いたアンサンブル予報 観測やに誤差はつきもの大気の初期状態はある存在確率で把握する方が望ましい ( 特に局地豪雨は初期値に敏感で 決定論的な予報は困難 ) 単独予報値摂動予報値

More information

気象庁 札幌管区気象台 資料 -6 Sapporo Regional Headquarters Japan Meteorological Agency 平成 29 年度防災気象情報の改善 5 日先までの 警報級の可能性 について 危険度を色分けした時系列で分かりやすく提供 大雨警報 ( 浸水害 )

気象庁 札幌管区気象台 資料 -6 Sapporo Regional Headquarters Japan Meteorological Agency 平成 29 年度防災気象情報の改善 5 日先までの 警報級の可能性 について 危険度を色分けした時系列で分かりやすく提供 大雨警報 ( 浸水害 ) 気象庁 札幌管区気象台 資料 -6 平成 29 年度防災気象情報の改善 5 日先までの 警報級の可能性 について 危険度を色分けした時系列で分かりやすく提供 大雨警報 ( 浸水害 ) を改善するための表面雨量指数の導入及び大雨警報 ( 浸水害 ) の危険度分布の提供 洪水警報を改善するための流域雨量指数の精緻化及び洪水警報の危険度分布の提供 メッシュ情報 ( 危険度分布 ) の技術を活用した大雨特別警報の発表対象区域の改善

More information

Kumamoto University Center for Multimedia and Information Technologies Lab. 熊本大学アプリケーション実験 ~ 実環境における無線 LAN 受信電波強度を用いた位置推定手法の検討 ~ InKIAI 宮崎県美郷

Kumamoto University Center for Multimedia and Information Technologies Lab. 熊本大学アプリケーション実験 ~ 実環境における無線 LAN 受信電波強度を用いた位置推定手法の検討 ~ InKIAI 宮崎県美郷 熊本大学アプリケーション実験 ~ 実環境における無線 LAN 受信電波強度を用いた位置推定手法の検討 ~ InKIAI プロジェクト @ 宮崎県美郷町 熊本大学副島慶人川村諒 1 実験の目的 従来 信号の受信電波強度 (RSSI:RecevedSgnal StrengthIndcator) により 対象の位置を推定する手法として 無線 LAN の AP(AccessPont) から受信する信号の減衰量をもとに位置を推定する手法が多く検討されている

More information

PowerPoint プレゼンテーション

PowerPoint プレゼンテーション 非線形カルマンフィルタ ~a. 問題設定 ~ 離散時間非線形状態空間表現 x k + 1 = f x k y k = h x k + bv k + w k f : ベクトル値をとるx k の非線形関数 h : スカラ値をとるx k の非線形関数 v k システム雑音 ( 平均値 0, 分散 σ v 2 k ) x k + 1 = f x k,v k w k 観測雑音 ( 平均値 0, 分散 σ w

More information

19 σ = P/A o σ B Maximum tensile strength σ % 0.2% proof stress σ EL Elastic limit Work hardening coefficient failure necking σ PL Proportional

19 σ = P/A o σ B Maximum tensile strength σ % 0.2% proof stress σ EL Elastic limit Work hardening coefficient failure necking σ PL Proportional 19 σ = P/A o σ B Maximum tensile strength σ 0. 0.% 0.% proof stress σ EL Elastic limit Work hardening coefficient failure necking σ PL Proportional limit ε p = 0.% ε e = σ 0. /E plastic strain ε = ε e

More information

Microsoft PowerPoint - 時系列解析(10)_講義用.pptx

Microsoft PowerPoint - 時系列解析(10)_講義用.pptx 時系列解析 () 季節調整モデルと成分分解 信号抽出 東京 学数理 情報教育研究センター 北川源四郎 東京 学北川源四郎数理 法 VII ( 時系列解析 ) 季節調整とは.5 WHARD 月次データ.3..9.7 5 49 73 97 45 何らかの原因で特定の周期で繰り返す成分を除去して本質的な現象を抽出する方法 東京 学北川源四郎数理 法 VII ( 時系列解析 ) 季節調整モデル 観測モデル

More information

SAP11_03

SAP11_03 第 3 回 音声音響信号処理 ( 線形予測分析と自己回帰モデル ) 亀岡弘和 東京大学大学院情報理工学系研究科日本電信電話株式会社 NTT コミュニケーション科学基礎研究所 講義内容 ( キーワード ) 信号処理 符号化 標準化の実用システム例の紹介情報通信の基本 ( 誤り検出 訂正符号 変調 IP) 符号化技術の基本 ( 量子化 予測 変換 圧縮 ) 音声分析 合成 認識 強調 音楽信号処理統計的信号処理の基礎

More information

Microsoft PowerPoint nakagawa.ppt [互換モード]

Microsoft PowerPoint nakagawa.ppt [互換モード] 気象庁現業全球モデルによる 台風予報の現状と課題 2013 年 3 月 6 日 第 6 回気象庁数値モデル研究会 数値モデルによる台風予報の課題と展望 気象庁予報部数値予報課中川雅之 檜垣将和 氏家将志 1 内容 気象庁全球数値予報システムの概要 台風進路予報の現状と課題 台風強度予報の現状と課題 今後の開発計画とまとめ 2 気象庁全球数値予報システムの概要 3 気象庁の全球数値予報システムの概要

More information

<4D F736F F F696E74202D E738A5889BB8BE688E68A4F82CC926E89BF908492E882C98AD682B782E98CA48B862E707074>

<4D F736F F F696E74202D E738A5889BB8BE688E68A4F82CC926E89BF908492E882C98AD682B782E98CA48B862E707074> 市街化区域外の地価推定に関する研究 不動産 空間計量研究室 筑波大学第三学群社会工学類都市計画主専攻宮下将尚筑波大学大学院システム情報工学研究科社会システム工学専攻高野哲司 背景 日本の国土の区域区分 都市計画区域 市街化区域 市街化を促進する区域 市街化調整区域 市街化を抑制する区域 非線引都市計画区域 上記に属さない区域 非線引き市街化調整区域市街化区域 都市計画区域 本研究での対象区域 都市計画区域外

More information

スライド タイトルなし

スライド タイトルなし アンテナ狭小化に伴う方位分解能劣化 の改善と東京湾での評価結果 - 民需等の利活用拡大を目指して - 直線 4 アレイ ( 八木 ) 菱形 4 アレイ ( ダイポール ) 伊藤浩之, 千葉修, 小海尊宏, 大西喬之 *1 山田寛喜 *2 長野日本無線 ( 株 ) *1 新潟大学 *2 08 年 12 月 17 日 08 年海洋レーダ研究集会 No.1 目次 1. はじめに : 海洋レーダの課題 2.

More information

CAEシミュレーションツールを用いた統計の基礎教育 | (株)日科技研

CAEシミュレーションツールを用いた統計の基礎教育 | (株)日科技研 CAE シミュレーションツール を用いた統計の基礎教育 ( 株 ) 日本科学技術研修所数理事業部 1 現在の統計教育の課題 2009 年から統計教育が中等 高等教育の必須科目となり, 大学でも問題解決ができるような人材 ( 学生 ) を育てたい. 大学ではコンピューター ( 統計ソフトの利用 ) を重視した教育をより積極的におこなうのと同時に, 理論面もきちんと教育すべきである. ( 報告 数理科学分野における統計科学教育

More information

I-2 (100 ) (1) y(x) y dy dx y d2 y dx 2 (a) y + 2y 3y = 9e 2x (b) x 2 y 6y = 5x 4 (2) Bernoulli B n (n = 0, 1, 2,...) x e x 1 = n=0 B 0 B 1 B 2 (3) co

I-2 (100 ) (1) y(x) y dy dx y d2 y dx 2 (a) y + 2y 3y = 9e 2x (b) x 2 y 6y = 5x 4 (2) Bernoulli B n (n = 0, 1, 2,...) x e x 1 = n=0 B 0 B 1 B 2 (3) co 16 I ( ) (1) I-1 I-2 I-3 (2) I-1 ( ) (100 ) 2l x x = 0 y t y(x, t) y(±l, t) = 0 m T g y(x, t) l y(x, t) c = 2 y(x, t) c 2 2 y(x, t) = g (A) t 2 x 2 T/m (1) y 0 (x) y 0 (x) = g c 2 (l2 x 2 ) (B) (2) (1)

More information

Probit , Mixed logit

Probit , Mixed logit Probit, Mixed logit 2016/5/16 スタートアップゼミ #5 B4 後藤祥孝 1 0. 目次 Probit モデルについて 1. モデル概要 2. 定式化と理解 3. 推定 Mixed logit モデルについて 4. モデル概要 5. 定式化と理解 6. 推定 2 1.Probit 概要 プロビットモデルとは. 効用関数の誤差項に多変量正規分布を仮定したもの. 誤差項には様々な要因が存在するため,

More information

大気環境シミュレーション

大気環境シミュレーション 第 3 回 (Q) 各自 eelを用いて 次の漸化式 + = の解の初期値依存性を調べよ.は50まで () 0 =.0 () 0 =.5 (3) 0 =.0 締切 04 年 月 6 日 ( 月 ) 夕方まで 提出先 347 室 オーバーフロー失敗ゴメンなさい (Q) 各自 eelを用いて 次の漸化式 + = の解の初期値依存性を調べよ.は50まで () 0 =.330 () 0 =.33 (3) 0

More information

25 3 4

25 3 4 25 3 4 1 µ e + ν e +ν µ µ + e + +ν e + ν µ e e + TAC START STOP START veto START (2.04 ± 0.18)µs 1/2 STOP (2.09 ± 0.11)µs 1/8 G F /( c) 3 (1.21±0.09) 5 /GeV 2 (1.19±0.05) 5 /GeV 2 Weinberg θ W sin θ W

More information

数値計算で学ぶ物理学 4 放物運動と惑星運動 地上のように下向きに重力がはたらいているような場においては 物体を投げると放物運動をする 一方 中心星のまわりの重力場中では 惑星は 円 だ円 放物線または双曲線を描きながら運動する ここでは 放物運動と惑星運動を 運動方程式を導出したうえで 数値シミュ

数値計算で学ぶ物理学 4 放物運動と惑星運動 地上のように下向きに重力がはたらいているような場においては 物体を投げると放物運動をする 一方 中心星のまわりの重力場中では 惑星は 円 だ円 放物線または双曲線を描きながら運動する ここでは 放物運動と惑星運動を 運動方程式を導出したうえで 数値シミュ 数値計算で学ぶ物理学 4 放物運動と惑星運動 地上のように下向きに重力がはたらいているような場においては 物体を投げると放物運動をする 一方 中心星のまわりの重力場中では 惑星は 円 だ円 放物線または双曲線を描きながら運動する ここでは 放物運動と惑星運動を 運動方程式を導出したうえで 数値シミュレーションによって計算してみる 4.1 放物運動一様な重力場における放物運動を考える 一般に質量の物体に作用する力をとすると運動方程式は

More information

PowerPoint プレゼンテーション

PowerPoint プレゼンテーション 反応工学 Raction Enginring 講義時間 ( 場所 : 火曜 限 (8-A 木曜 限 (S-A 担当 : 山村 火 限 8-A 期末試験中間試験以降 /7( 木 まで持ち込みなし要電卓 /4( 木 質問受付日講義なし 授業アンケート (li campus の入力をお願いします 晶析 (crystallization ( 教科書 p. 濃度 溶解度曲線 C C s A 安定 液 ( 気

More information

Microsoft Word - 予稿集表紙.doc

Microsoft Word - 予稿集表紙.doc ミクロ組織に基づくフェライト セメンタイト鋼の脆性破壊発生予測 柴沼一樹東京大学大学院工学系研究科 ミクロ組織に基づくフェライト セメンタイト鋼の脆性破壊発生予測 柴沼一樹 東京大学 大学院工学系研究科システム創成学専攻 113-8656 東京都文京区本郷 7-3-1 shibanuma@struct.t.-u-tokyo.ac.jp 近年, 構造物に使用される鋼材の高張力化や使用環境の過酷化が進み,

More information

気象庁の現業数値予報システム一覧 数値予報システム ( 略称 ) 局地モデル (LFM) メソモデル (MSM) 全球モデル (GSM) 全球アンサンブル予報システム 全球アンサンブル予報システム 季節アンサンブル予報システム 水平分解能 2km 5km 約 20km 約 40km 約 40km(1

気象庁の現業数値予報システム一覧 数値予報システム ( 略称 ) 局地モデル (LFM) メソモデル (MSM) 全球モデル (GSM) 全球アンサンブル予報システム 全球アンサンブル予報システム 季節アンサンブル予報システム 水平分解能 2km 5km 約 20km 約 40km 約 40km(1 気象庁データを利用した気象研究の現状と展望 気象庁現業メソモデルの 最近の開発とその成果 2017 年 5 月 25 日 気象庁予報部数値予報課 原旅人 1 気象庁の現業数値予報システム一覧 数値予報システム ( 略称 ) 局地モデル (LFM) メソモデル (MSM) 全球モデル (GSM) 全球アンサンブル予報システム 全球アンサンブル予報システム 季節アンサンブル予報システム 水平分解能 2km

More information

図 1 COBE-SST のオリジナル格子から JCDAS の格子に変換を行う際に用いられている海陸マスク 緑色は陸域 青色は海域 赤色は内海を表す 内海では気候値 (COBE-SST 作成時に用いられている 1951~2 年の平均値 ) が利用されている (a) (b) SST (K) SST a

図 1 COBE-SST のオリジナル格子から JCDAS の格子に変換を行う際に用いられている海陸マスク 緑色は陸域 青色は海域 赤色は内海を表す 内海では気候値 (COBE-SST 作成時に用いられている 1951~2 年の平均値 ) が利用されている (a) (b) SST (K) SST a 平成 22 年 2 月 JCDAS における 内海の海面水温の取り扱いの不具合について 気象庁地球環境 海洋部気候情報課 気候データ同化システム (JCDAS) では COBE-SST 累年値データを境界条件とする 6 時間予報及び客観解析を行っておりますが 25 年 1 月の JCDAS のルーチン運用開始以降 一部の内海において SST 観測値ではなく気候値が適用されていることが判明しました 原因

More information

スライド 1

スライド 1 ディスカッション ( データ同化相談会 ) 縮約モデルと粒子フィルタを用いた 低層風擾乱の流れ場予測 菊地亮太 ( 東北大学流体科学研究所 ) 1 研究背景低層風擾乱の運航影響 空港周辺での低層風擾乱による運航影響が問題となっている 着陸復行 代替地着陸 ハードランディング 地形や ビルなどの影響を強く受ける 小規模な気象現象 地形性乱気流 2 航空機運航において重要な課題となっている 研究背景低層風擾乱の解析

More information

平成 27 年度 新潟県たにはま海水浴場流況調査 報告書 平成 27 年 6 月調査 第九管区海上保安本部

平成 27 年度 新潟県たにはま海水浴場流況調査 報告書 平成 27 年 6 月調査 第九管区海上保安本部 平成 27 年度 新潟県たにはま海水浴場流況調査 報告書 平成 27 年 6 月調査 第九管区海上保安本部 1 目的平成 27 年度海洋情報業務計画に基づき 新潟県たにはま海水浴場において流況調査を実施し 海浜事故の防止に資する基礎資料を得るとともに 管内海上保安部署と連携してマリンレジャー安全推進に係る啓発活動を実施する 2 調査区域 新潟県上越市大字長浜たにはま海水浴場 ( 図 1 参照 ) 3

More information

3 数値解の特性 3.1 CFL 条件 を 前の章では 波動方程式 f x= x0 = f x= x0 t f c x f =0 [1] c f 0 x= x 0 x 0 f x= x0 x 2 x 2 t [2] のように差分化して数値解を求めた ここでは このようにして得られた数値解の性質を 考

3 数値解の特性 3.1 CFL 条件 を 前の章では 波動方程式 f x= x0 = f x= x0 t f c x f =0 [1] c f 0 x= x 0 x 0 f x= x0 x 2 x 2 t [2] のように差分化して数値解を求めた ここでは このようにして得られた数値解の性質を 考 3 数値解の特性 3.1 CFL 条件 を 前の章では 波動方程式 f x= x = f x= x t f c x f = [1] c f x= x f x= x 2 2 t [2] のように差分化して数値解を求めた ここでは このようにして得られた数値解の性質を 考える まず 初期時刻 t=t に f =R f exp [ik x ] [3] のような波動を与えたとき どのように時間変化するか調べる

More information

様々なミクロ計量モデル†

様々なミクロ計量モデル† 担当 : 長倉大輔 ( ながくらだいすけ ) この資料は私の講義において使用するために作成した資料です WEB ページ上で公開しており 自由に参照して頂いて構いません ただし 内容について 一応検証してありますが もし間違いがあった場合でもそれによって生じるいかなる損害 不利益について責任を負いかねますのでご了承ください 間違いは発見次第 継続的に直していますが まだ存在する可能性があります 1 カウントデータモデル

More information

予報時間を39時間に延長したMSMの初期時刻別統計検証

予報時間を39時間に延長したMSMの初期時刻別統計検証 第 1 章領域拡張 予報時間 39 時間化されたメソモデルの特性 1.1 メソモデルの領域拡張 予報時間 39 時間化の概 1 要メソモデル (MSM) は 2013 年 3 月に予報領域が拡張された また 2013 年 5 月に全初期時刻における予報時間が39 時間に延長された 表 1.1.1に今回の変更前後の主な仕様を また 図 1.1.1に領域拡張前後の予報領域を示す 本節では 仕様拡張の目的及び概要を説明する

More information

データ解析

データ解析 データ解析 ( 前期 ) 最小二乗法 向井厚志 005 年度テキスト 0 データ解析 - 最小二乗法 - 目次 第 回 Σ の計算 第 回ヒストグラム 第 3 回平均と標準偏差 6 第 回誤差の伝播 8 第 5 回正規分布 0 第 6 回最尤性原理 第 7 回正規分布の 分布の幅 第 8 回最小二乗法 6 第 9 回最小二乗法の練習 8 第 0 回最小二乗法の推定誤差 0 第 回推定誤差の計算 第

More information

.. 9 (NAPS9) NAPS km, km, 3 km, 8 (GSM) 64 : / 64 : / 64 : / (UTC) (UTC) (, UTC) 3 : 3 / 3 : 3 / 3 : / (, 6, 8UTC) (, 6, 8UTC) (6, 8UTC) 4 km,

.. 9 (NAPS9) NAPS km, km, 3 km, 8 (GSM) 64 : / 64 : / 64 : / (UTC) (UTC) (, UTC) 3 : 3 / 3 : 3 / 3 : / (, 6, 8UTC) (, 6, 8UTC) (6, 8UTC) 4 km, . 8 6 99 Cray XC Cray SR6/M (NAPS) 3 3 4, 6, 8UTC 84 3 6 6 (LFM)..3 NAPS (7) XC 8 6 TOP (https://www.top.org), 6 3 (NAPS) 4 8 UTC 64 6 8 6. NAPS.. NAPS (7) (7).. NAPS 8 6 NAPS.. () () NAPS 8 6 6, 6, 8UTC

More information

#A A A F, F d F P + F P = d P F, F y P F F x A.1 ( α, 0), (α, 0) α > 0) (x, y) (x + α) 2 + y 2, (x α) 2 + y 2 d (x + α)2 + y 2 + (x α) 2 + y 2 =

#A A A F, F d F P + F P = d P F, F y P F F x A.1 ( α, 0), (α, 0) α > 0) (x, y) (x + α) 2 + y 2, (x α) 2 + y 2 d (x + α)2 + y 2 + (x α) 2 + y 2 = #A A A. F, F d F P + F P = d P F, F P F F A. α, 0, α, 0 α > 0, + α +, α + d + α + + α + = d d F, F 0 < α < d + α + = d α + + α + = d d α + + α + d α + = d 4 4d α + = d 4 8d + 6 http://mth.cs.kitmi-it.c.jp/

More information

<4D F736F F F696E74202D CD90EC8B5A8F708DA7926B89EF816993BF938791E BA816A8F4390B E707074>

<4D F736F F F696E74202D CD90EC8B5A8F708DA7926B89EF816993BF938791E BA816A8F4390B E707074> 安価で効率的な水位 - 流量曲線の作成と流出モデルパラメータ同定法, その応用 徳島大学准教授徳島大学名誉教授 田村隆雄端野道夫 概略概略 流出モデルと水位 - 流量曲線とを連動させることにより, 実測流量データがない河川でも, 洪水解析 ( 流量ハイドログラフの推定 ) と H-Q 曲線の作成を可能とします. 以下の特徴があります. 雨量データと水位データがあれば適用できる 河道横断形状が不明でも適用できる

More information

untitled

untitled MRR Physical Basis( 1.8.4) METEK MRR 1 MRR 1.1 MRR 24GHz FM-CW(frequency module continuous wave) 30 r+ r f+ f 1.2 1 4 MRR 24GHz 1.3 50mW 1 rf- (waveguide) (horn) 60cm ( monostatic radar) (continuous wave)

More information

DVIOUT

DVIOUT 第 章 離散フーリエ変換 離散フーリエ変換 これまで 私たちは連続関数に対するフーリエ変換およびフーリエ積分 ( 逆フーリエ変換 ) について学んできました この節では フーリエ変換を離散化した離散フーリエ変換について学びましょう 自然現象 ( 音声 ) などを観測して得られる波 ( 信号値 ; 観測値 ) は 通常 電気信号による連続的な波として観測機器から出力されます しかしながら コンピュータはこの様な連続的な波を直接扱うことができないため

More information

Microsoft PowerPoint SIGAL.ppt

Microsoft PowerPoint SIGAL.ppt アメリカン アジアンオプションの 価格の近似に対する 計算幾何的アプローチ 渋谷彰信, 塩浦昭義, 徳山豪 ( 東北大学大学院情報科学研究科 ) 発表の概要 アメリカン アジアンオプション金融派生商品の一つ価格付け ( 価格の計算 ) は重要な問題 二項モデルにおける価格付けは計算困難な問題 目的 : 近似精度保証をもつ近似アルゴリズムの提案 アイディア : 区分線形関数を計算幾何手法により近似 問題の説明

More information

, COMPUTATION OF SHALLOW WATER EQUATION WITH HIERARCHICAL QUADTREE GRID SYSTEM 1 2 Hiroyasu YASUDA and Tsuyoshi HOSHINO

, COMPUTATION OF SHALLOW WATER EQUATION WITH HIERARCHICAL QUADTREE GRID SYSTEM 1 2 Hiroyasu YASUDA and Tsuyoshi HOSHINO , 2 11 8 COMPUTATION OF SHALLOW WATER EQUATION WITH HIERARCHICAL QUADTREE GRID SYSTEM 1 2 Hiroyasu YASUDA and Tsuyoshi HOSHINO 1 9-2181 2 8 2 9-2181 2 8 Numerical computation of river flows have been employed

More information

Microsoft PowerPoint - e-stat(OLS).pptx

Microsoft PowerPoint - e-stat(OLS).pptx 経済統計学 ( 補足 ) 最小二乗法について 担当 : 小塚匡文 2015 年 11 月 19 日 ( 改訂版 ) 神戸大学経済学部 2015 年度後期開講授業 補足 : 最小二乗法 ( 単回帰分析 ) 1.( 単純 ) 回帰分析とは? 標本サイズTの2 変数 ( ここではXとY) のデータが存在 YをXで説明する回帰方程式を推定するための方法 Y: 被説明変数 ( または従属変数 ) X: 説明変数

More information

スライド 1

スライド 1 本資料について 本資料は下記論文を基にして作成されたものです. 文書の内容の正確さは保障できないため, 正確な知識を求める方は原文を参照してください. 著者 : 伊藤誠吾吉田廣志河口信夫 論文名 : 無線 LANを用いた広域位置情報システム構築に関する検討 出展 : 情報処理学会論文誌 Vol.47 No.42 発表日 :2005 年 12 月 著者 : 伊藤誠悟河口信夫 論文名 : アクセスポイントの選択を考慮したベイズ推定による無線

More information

Chap2.key

Chap2.key . f( ) V (V V ) V e + V e V V V V ( ) V V ( ) E. - () V (0 ) () V (0 ) () V (0 ) (4) V ( ) E. - () V (0 ) () V (0 ) O r θ ( ) ( ) : (r θ) : { r cos θ r sn θ { r + () V (0 ) (4) V ( ) θ θ arg( ) : π π

More information

PowerPoint プレゼンテーション

PowerPoint プレゼンテーション 第 6 回基礎ゼミ資料 Practice NL&MXL from R 平成 30 年 5 月 18 日 ( 金 ) 朝倉研究室修士 1 年小池卓武 使用データ 1 ~ 横浜プローブパーソンデータ ~ 主なデータの中身 トリップ ID 目的 出発, 到着時刻 総所要時間 移動距離 交通機関別の時間, 距離 アクセス, イグレス時間, 距離 費用 代表交通手段 代替手段生成可否 性別, 年齢等の個人属性

More information

台風21号越波等検証委員会 -海象状況と浸水状況の再現-

台風21号越波等検証委員会 -海象状況と浸水状況の再現- 台風 21 号越波等検証委員会 - 海象状況と浸水状況の再現 - 2018.12.11 Committee on Assessing Overtopping Wave Caused by Typhoon Jebi -Third Committee Meeting- Reproducing marine phenomenon and flooding situation 目次 1. 関西国際空港における台風

More information

解析センターを知っていただく キャンペーン

解析センターを知っていただく キャンペーン 005..5 SAS 問題設定 目的 PKパラメータ (AUC,Cmax,Tmaxなど) の推定 PKパラメータの群間比較 PKパラメータのバラツキの評価! データの特徴 非反復測定値 個体につき 個の測定値しか得られない plasma concentration 非反復測定値のイメージ図 測定時点間で個体の対応がない 着目する状況 plasma concentration 経時反復測定値のイメージ図

More information

1 Jensen et al.[6] GRT S&P500 GRT RT GRT Kiriu and Hibiki[8] Jensen et al.[6] GRT 3 GRT Generalized Recovery Theorem (Jensen et al.[6])

1 Jensen et al.[6] GRT S&P500 GRT RT GRT Kiriu and Hibiki[8] Jensen et al.[6] GRT 3 GRT Generalized Recovery Theorem (Jensen et al.[6]) Generalized Recovery Theorem Ross[11] Recovery Theorem(RT) RT forward looking Kiriu and Hibiki[8] Generalized Recovery Theorem(GRT) Jensen et al.[6] GRT RT Kiriu and Hibiki[8] 1 backward looking forward

More information

ニュートン重力理論.pptx

ニュートン重力理論.pptx 3 ニュートン重力理論 1. ニュートン重力理論の基本 : 慣性系とガリレイ変換不変性 2. ニュートン重力理論の定式化 3. 等価原理 4. 流体力学方程式とその基礎 3.1 ニュートン重力理論の基本 u ニュートンの第一法則 = 力がかからなければ 等速直線運動を続ける u 等速直線運動に見える系を 慣性系 と呼ぶ ² 直線とはどんな空間の直線か? ニュートン理論では 3 次元ユークリッド空間

More information

第122号.indd

第122号.indd -1- -2- -3- 0852-36-5150 0852-36-5163-4- -5- -6- -7- 1st 1-1 1-2 1-3 1-4 1-5 -8- 2nd M2 E2 D2 J2 C2-9- 3rd M3 E3 D3 J3 C3-10- 4th M4 E4 D4 J4 C4-11- -12- M5 E5 J5 D5 C5 5th -13- -14- NEWS NEWS -15- NEWS

More information

技術資料 台風モデルによる波浪の再現計算と経路変更による感度実験 2018 年 9 月の台風 18 号を例に 岩﨑慎介 大塚淳一 1. はじめに気候変動に関する政府間パネル (IPCC) の第 5 次報告書 1) ( 以下 AR5) の公表以降 気候変動による沿岸部への影響は今後中長期的に避けること

技術資料 台風モデルによる波浪の再現計算と経路変更による感度実験 2018 年 9 月の台風 18 号を例に 岩﨑慎介 大塚淳一 1. はじめに気候変動に関する政府間パネル (IPCC) の第 5 次報告書 1) ( 以下 AR5) の公表以降 気候変動による沿岸部への影響は今後中長期的に避けること 技術資料 台風モデルによる波浪の再現計算と経路変更による感度実験 2018 年 9 月の台風 18 号を例に 岩﨑慎介 大塚淳一 1. はじめに気候変動に関する政府間パネル (IPCC) の第 5 次報告書 1) ( 以下 AR5) の公表以降 気候変動による沿岸部への影響は今後中長期的に避けることができないものと認識することが必要とされ 平成 27 年 2 月に変更された 海岸保全区域等に係る海岸の保全に関する基本的な方針

More information

最新耐震構造解析 ( 第 3 版 ) サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. このサンプルページの内容は, 第 3 版 1 刷発行時のものです.

最新耐震構造解析 ( 第 3 版 ) サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます.   このサンプルページの内容は, 第 3 版 1 刷発行時のものです. 最新耐震構造解析 ( 第 3 版 ) サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/052093 このサンプルページの内容は, 第 3 版 1 刷発行時のものです. i 3 10 3 2000 2007 26 8 2 SI SI 20 1996 2000 SI 15 3 ii 1 56 6

More information

1. 設計手順 ディープウェル工事の設計は 下記に示す手順で実施する 掘削区域内への排水量の検討 ディープウェル仕様の仮定 ( 径 深さ ) ディープウェル 1 本当たりの揚水能力の検討 ディープウェル本数 配置の設定 井戸配置で最も不利な点を所要水位低下させるのに必要な各井戸の合計排水量の検討 -

1. 設計手順 ディープウェル工事の設計は 下記に示す手順で実施する 掘削区域内への排水量の検討 ディープウェル仕様の仮定 ( 径 深さ ) ディープウェル 1 本当たりの揚水能力の検討 ディープウェル本数 配置の設定 井戸配置で最も不利な点を所要水位低下させるのに必要な各井戸の合計排水量の検討 - 管理記号 : 0001 作成年月日 : 2018/6/18 工事名称 : 仮設計画ガイドブック ( 全日本建設技術協会 ) 工区名称 : page209~page214 設計条件 設計結果 ディープウェル工事設計計算書 1. 掘削寸法 ( 幅 )40.0m ( 長さ )40.0m ( 深さ )12.0m 2. 滞水層厚 D=19.0m 3. 地下水位 GL-3.0m 4. 計画水位 GL-13.0m

More information

密集市街地における換気・通風性能簡易評価ツールの開発 (その2 流体計算部分の開発)」

密集市街地における換気・通風性能簡易評価ツールの開発 (その2 流体計算部分の開発)」 OpenCAE ワークショップ 2013 2013.6.21 密集市街地における換気 通風性能簡易評価ツールの開発その 2 : 流体計算部分の開発 福本雅彦 ( 株式会社森村設計 ) 小縣信也 ( 株式会社森村設計 ) 勝又済 ( 国土交通省国土技術政策総合研究所 ) 西澤繁毅 ( 国土交通省国土技術政策総合研究所 ) 岩見達也 ( 国土交通省国土技術政策総合研究所 ) 概要 換気 通風性能簡易評価ツール

More information

29

29 9 .,,, 3 () C k k C k C + C + C + + C 8 + C 9 + C k C + C + C + C 3 + C 4 + C 5 + + 45 + + + 5 + + 9 + 4 + 4 + 5 4 C k k k ( + ) 4 C k k ( k) 3 n( ) n n n ( ) n ( ) n 3 ( ) 3 3 3 n 4 ( ) 4 4 4 ( ) n n

More information

r d 2r d l d (a) (b) (c) 1: I(x,t) I(x+ x,t) I(0,t) I(l,t) V in V(x,t) V(x+ x,t) V(0,t) l V(l,t) 2: 0 x x+ x 3: V in 3 V in x V (x, t) I(x, t

r d 2r d l d (a) (b) (c) 1: I(x,t) I(x+ x,t) I(0,t) I(l,t) V in V(x,t) V(x+ x,t) V(0,t) l V(l,t) 2: 0 x x+ x 3: V in 3 V in x V (x, t) I(x, t 1 1 2 2 2r d 2r d l d (a) (b) (c) 1: I(x,t) I(x+ x,t) I(0,t) I(l,t) V in V(x,t) V(x+ x,t) V(0,t) l V(l,t) 2: 0 x x+ x 3: V in 3 V in x V (x, t) I(x, t) V (x, t) I(x, t) V in x t 3 4 1 L R 2 C G L 0 R 0

More information

untitled

untitled ( 9:: 3:6: (k 3 45 k F m tan 45 k 45 k F m tan S S F m tan( 6.8k tan k F m ( + k tan 373 S S + Σ Σ 3 + Σ os( sin( + Σ sin( os( + sin( os( p z ( γ z + K pzdz γ + K γ K + γ + 9 ( 9 (+ sin( sin { 9 ( } 4

More information

ダイポールアンテナ標準:校正の実際と不確かさ

ダイポールアンテナ標準:校正の実際と不確かさ ダイポールアンテナ標準 校正の実際と不確かさ ( 独 ) 産業技術総合研究所 森岡健浩 概要 アンテナ係数 3アンテナ法 ( 半自由空間と自由空間 ) 置換法 不確かさ積算 異なるアンテナ校正によるアンテナ係数の一意性 まとめ アンテナ係数の定義 z 波源 V 付属回路 受信アンテナ図 アンテナ係数の定義 V 測定量 : アンテナ係数 ( 水平偏波.0 m 高 または自由空間 ) 校正方法 : 3アンテナ法

More information

背景 ヤマセと海洋の関係 図 1: 親潮の流れ ( 気象庁 HP より ) 図 2:02 年 7 月上旬の深さ 100m の水温図 ( )( 気象庁 HP より ) 黒潮続流域 親潮の貫入 ヤマセは混合域の影響を強く受ける現象 ヤマセの気温や鉛直構造に沿岸の海面水温 (SST) や親潮フロントの影響

背景 ヤマセと海洋の関係 図 1: 親潮の流れ ( 気象庁 HP より ) 図 2:02 年 7 月上旬の深さ 100m の水温図 ( )( 気象庁 HP より ) 黒潮続流域 親潮の貫入 ヤマセは混合域の影響を強く受ける現象 ヤマセの気温や鉛直構造に沿岸の海面水温 (SST) や親潮フロントの影響 大気再解析データで表現されるヤマセ - モデルによる SST の違いと解析された気温への影響 - 弘前大学大学院理工学専攻佐々木実紀 背景 ヤマセと海洋の関係 図 1: 親潮の流れ ( 気象庁 HP より ) 図 2:02 年 7 月上旬の深さ 100m の水温図 ( )( 気象庁 HP より ) 黒潮続流域 親潮の貫入 ヤマセは混合域の影響を強く受ける現象 ヤマセの気温や鉛直構造に沿岸の海面水温

More information

技術資料 JARI Research Journal OpenFOAM を用いた沿道大気質モデルの開発 Development of a Roadside Air Quality Model with OpenFOAM 木村真 *1 Shin KIMURA 伊藤晃佳 *2 Akiy

技術資料 JARI Research Journal OpenFOAM を用いた沿道大気質モデルの開発 Development of a Roadside Air Quality Model with OpenFOAM 木村真 *1 Shin KIMURA 伊藤晃佳 *2 Akiy 技術資料 176 OpenFOAM を用いた沿道大気質モデルの開発 Development of a Roadside Air Quality Model with OpenFOAM 木村真 *1 Shin KIMURA 伊藤晃佳 *2 Akiyoshi ITO 1. はじめに自動車排出ガスの環境影響は, 道路沿道で大きく, 建物など構造物が複雑な気流を形成するため, 沿道大気中の自動車排出ガス濃度分布も複雑になる.

More information

lim lim lim lim 0 0 d lim 5. d 0 d d d d d d 0 0 lim lim 0 d

lim lim lim lim 0 0 d lim 5. d 0 d d d d d d 0 0 lim lim 0 d lim 5. 0 A B 5-5- A B lim 0 A B A 5. 5- 0 5-5- 0 0 lim lim 0 0 0 lim lim 0 0 d lim 5. d 0 d d d d d d 0 0 lim lim 0 d 0 0 5- 5-3 0 5-3 5-3b 5-3c lim lim d 0 0 5-3b 5-3c lim lim lim d 0 0 0 3 3 3 3 3 3

More information

1 3 1.1.......................... 3 1............................... 3 1.3....................... 5 1.4.......................... 6 1.5........................ 7 8.1......................... 8..............................

More information

GFS-WRF-SWAN システムによる 3 シーズンの波浪予測とシステムの検証 VERIFICATION OF GFS-WRF-SWAN WAVE PREDICTION SYSTEM BY THREE SEASONS' COMPARISON 間瀬肇 1 勝井伸悟 2 安田誠宏 3 Tracey H

GFS-WRF-SWAN システムによる 3 シーズンの波浪予測とシステムの検証 VERIFICATION OF GFS-WRF-SWAN WAVE PREDICTION SYSTEM BY THREE SEASONS' COMPARISON 間瀬肇 1 勝井伸悟 2 安田誠宏 3 Tracey H GFS--SWAN システムによる シーズンの波浪予測とシステムの検証 VERIFICATION OF GFS--SWAN WAVE PREDICTION SYSTEM BY THREE SEASONS' COMPARISON 間瀬肇 勝井伸悟 安田誠宏 Tracey H. Tom 小川和幸 Hajime MASE, Shingo KATSUI, Tomohiro YASUDA, Tracey H.

More information

Microsoft PowerPoint - H22制御工学I-10回.ppt

Microsoft PowerPoint - H22制御工学I-10回.ppt 制御工学 I 第 回 安定性 ラウス, フルビッツの安定判別 平成 年 6 月 日 /6/ 授業の予定 制御工学概論 ( 回 ) 制御技術は現在様々な工学分野において重要な基本技術となっている 工学における制御工学の位置づけと歴史について説明する さらに 制御システムの基本構成と種類を紹介する ラプラス変換 ( 回 ) 制御工学 特に古典制御ではラプラス変換が重要な役割を果たしている ラプラス変換と逆ラプラス変換の定義を紹介し

More information

微分方程式による現象記述と解きかた

微分方程式による現象記述と解きかた 微分方程式による現象記述と解きかた 土木工学 : 公共諸施設 構造物の有用目的にむけた合理的な実現をはかる方法 ( 技術 ) に関する学 橋梁 トンネル ダム 道路 港湾 治水利水施設 安全化 利便化 快適化 合法則的 経済的 自然および人口素材によって作られた 質量保存則 構造物の自然的な性質 作用 ( 外力による応答 ) エネルギー則 の解明 社会的諸現象のうち マスとしての移動 流通 運動量則

More information

kouenyoushi_kyoshida

kouenyoushi_kyoshida 大規模アンサンブルシミュレーションによる熱帯低気圧の将来変化 吉田康平, 杉正人, 水田亮, 石井正好 ( 気象研究所 ) 村上裕之 ( プリンストン大学, 米国地球流体力学研究所 ) 1. はじめに地球温暖化の熱帯低気圧 ( 以下 台風と表記 ) への影響は 科学的重要性に加え その社会的な影響から大きな関心を集める話題である 気候変動に関する政府間パネル (IPCC) 第五次評価報告書では 温暖化の進行とともに地球全体での台風の発生数が減少または実質的に変化しないことと

More information

T2K 実験 南野彰宏 ( 京都大学 ) 他 T2Kコラボレーション平成 25 年度宇宙線研究所共同利用成果発表会 2013 年 12 月 20 日 1

T2K 実験 南野彰宏 ( 京都大学 ) 他 T2Kコラボレーション平成 25 年度宇宙線研究所共同利用成果発表会 2013 年 12 月 20 日 1 T2K 実験 南野彰宏 ( 京都大学 ) 他 T2Kコラボレーション平成 25 年度宇宙線研究所共同利用成果発表会 2013 年 12 月 20 日 1 T2K 実験 J- PARC でほぼ純粋な ν µμ ビームを生成 生成点直後の前置検出器と 295km 離れたスーパーカミオカンデでニュートリノを観測 ニュートリノ振動の精密測定 T2K 実験における振動モード 1. ν µμ ν e (ν e

More information

台風解析の技術 平成 21 年 10 月 29 日気象庁予報部

台風解析の技術 平成 21 年 10 月 29 日気象庁予報部 台風解析の技術 平成 21 年 10 月 29 日気象庁予報部 内 容 1 台風解析の概要とポイント 2 台風解析の実例 台風第 18 号の事例 3 まとめ 2 1 台風解析の概要とポイント 1-1 諸元の決定 中心位置 最大風速 中心気圧 暴風域半径 強風域半径などをまとめて 台風諸元と呼んでいる 台風が日本の責任海域に存在する時は 3 時間毎に諸元決定を行っている 台風が日本に接近した時は 1

More information

<4D F736F F D2091E6358FCD31328B438FDB A5182F08ADC82DE816A2E646F6378>

<4D F736F F D2091E6358FCD31328B438FDB A5182F08ADC82DE816A2E646F6378> 注 : 大阪管区気象台における平成 15 年から平成 24 年の観測データを元に作成 図 5-12-3 日最大平均風速の風向出現頻度 346 5.12. 2 施設の存在に伴う影響の予測 評価 (1) 予測内容施設の存在に伴う影響として 建築物の出現が事業計画地周辺の風環境に及ぼす影響について 風洞実験により予測した 予測内容は表 5-12-2 に示すとおりである 表 5-12-2 予測内容 予測項目予測範囲

More information

津波警報等の留意事項津波警報等の利用にあたっては 以下の点に留意する必要があります 沿岸に近い海域で大きな地震が発生した場合 津波警報等の発表が津波の襲来に間に合わない場合があります 沿岸部で大きな揺れを感じた場合は 津波警報等の発表を待たず 直ちに避難行動を起こす必要があります 津波警報等は 最新

津波警報等の留意事項津波警報等の利用にあたっては 以下の点に留意する必要があります 沿岸に近い海域で大きな地震が発生した場合 津波警報等の発表が津波の襲来に間に合わない場合があります 沿岸部で大きな揺れを感じた場合は 津波警報等の発表を待たず 直ちに避難行動を起こす必要があります 津波警報等は 最新 2.3 津波に関する防災気象情報 (1) 大津波警報 津波警報 津波注意報 津波による災害の発生が予想される場合には 地震が発生してから約 3 分を目標に大津波警報 津波警報または津波注意報を発表 地震が発生した時は地震の規模や位置を即時に推定し これらをもとに沿岸で予想 される津波の高さを求め 津波による災害の発生が予想される場合には 地震が発生 してから約 3 分を目標に津波予報区ごとに大津波警報

More information