開発が欠かせません。そういった中で技術波及効果が高く、産業的に広い裾野を持つ航空 機産業の成長を促し、今後の我が国の基幹産業の一つにしていくことを目指して国際競争 力のある技術を産業界が受け止められるレベルまで持ち上げることがJAXA航空の役割 であり目標であると考えています。
文部科学省 科学技術・学術審議会研究計画・評価分科会航空科学技術委員会において
「航空科学技術に関する研究開発の推進のためのロードマップ(2012)」(平成24年 8月)として了承されたところですが、航空機産業を研究開発面で促進していくためには、
技術開発の面でも産官学の連携をよりシステマティックにして、お互いが相補的な役割を 担いながら、基礎的な研究が応用研究、実用化研究へとうまくつながって行くことが重要 であることは、言うまでもありません。
JAXA航空部門では、成果を効果的に創出するとともに、外部機関との協力を拡大し、我 が国の航空基盤を強化する等の観点から、従来の個別に進めてきた共同研究や委託研究の 枠を広げて、大きなテーマと課題を設定し、そのための解決策を外部の知的リソースに幅 広く求めていく仕組みとして、平成22年度から公募型研究制度を導入しています。
平成25年3月21日に各研究間の情報交流、意見交換の場を提供するために、平成2 4年度に採択されました「次世代運航システム(DREAMS)の研究開発」に係る課題 のうち研究継続中のものについて、成果報告会を開催しました。本報告書は、その際の成 果をまとめたものです。
JAXAの目指している産学官の相互交流による人材発掘、人材交流の促進及び若手人 材への機会提供への貢献の一端を理解していただくための資料として活用していただけれ ば幸いです。
宇宙研究開発機構航空プログラムグループ 航空プログラムディレクタ 岩宮敏幸
1.低層風擾乱のデータ同化シミュレーションに関する研究...1
○菊地 亮太、三坂 孝志、大林 茂(東北大学流体科学研究所)
2.不確実環境下でのコンフリクトフリー4次元軌道生成アルゴリズム...11
○松野 賀宣、土屋 武司(東京大学)
3.マルチGNSS高精度測位およびその信頼性評価アルゴリズムに関する研究 ...23
久保 幸弘(立命館大学)
4. Cognitive aspects of curved approaches ...37
ヨルグ エントジンガー(東京大学)
5.指向性を考慮した騒音伝搬予測手法に関する研究...45
坂本 慎一(東京大学生産技術研究所)
6.動的計画法によるコンフリクトフリー最適軌道の生成に関する研究...59
宮沢 与和(九州大学)
7.混雑空港周辺およびコンフリクト発生空域における
最適交通流制御アルゴリズムの開発 ...67
○請川 克之、矢野 夏子、宇田川 佑介(株式会社 構造計画研究所)
(所属は発表時)
పᒙ㢼ᨐ䛾䝕 䝍ྠ䝅䝭䝳䝺 䝅䝵䞁䛻㛵䛩䜛◊✲
పᒙ㢼ᨐ䛾䝕䞊䝍ྠ䝅䝭䝳䝺䞊䝅䝵䞁䛻㛵䛩䜛◊✲
ŨᓎעʰٽᲦɤקܑ࣓Ღٻᒔ Ტி҅ٻܖ ්˳ᅹܖᄂᆮᲣ
1 24࠰ࡇ ഏɭˊᢃᑋǷǹȆȠπѪᄂᆮɶ᧓إԓ˟ 2013࠰3உ21ଐ
◊✲⫼ᬒ పᒙ㢼ᨐ䛾㐠⯟ᙳ㡪
◊✲⫼ᬒ పᒙ㢼ᨐ䛾㐠⯟ᙳ㡪
ȷᑋᆰೞƷܲᡛƸفьͼӼ ȷᢃᑋƷјྙ҄ƕ൭NJǒǕƯƍǔ ȷјྙ҄ƴݣƢǔൢᝋƷࢨ᪪ƸٻƖƍ јྙ҄ƴݣƢǔൢᝋƷࢨ᪪ƸٻƖƍ
ଐஜᑋᆰೞႆң˟
ȷᆰลԗᡀưƷ˯ޖʏƴǑǔࢨ᪪ƸཎƴٻƖƍ ȷბᨕࣄᘍ
ˊஆעბᨕ ȷˊஆעბᨕ
ȷȏȸȉȩȳȇǣȳǰ
࠽ϋᆰลưƸ˯ޖʏƴǑǓ
ܤμƔ јྙႎƳᢃᑋ ƨNJ
࠽ϋᆰลưƸ˯ޖʏƴǑǓ ݼᑋྙƴࢨ᪪ƕưƯƍǔ
ܤμƔƭјྙႎƳᢃᑋƷƨNJƴ
˯ޖʏǛʖยƢǔƜƱƸ᩼ࠝƴᙲ
࠽ϋᆰลƷܭ̝ƷஉКݼᑋ 18~20◊✲⫼ᬒ పᒙ㢼ᨐ䛾ண
◊✲⫼ᬒ పᒙ㢼ᨐ䛾ண
˯ޖʏᲴᆰ᧓ႎƴݱᙹƳൢᝋྵᝋ
ࡽဇᲴJAXA
ൢᝋʖยȢȇȫᲴȡǽǹDZȸȫˌɥƷྵᝋǛݣᝋ
ൢᝋ ย ྵᝋ ᝋ
ൢᝋʖยȢȇȫƩƚưƷʖยƸᩊ
3
ሗ࿌ෆᐜ ሗ࿌ෆᐜ
˯ޖʏƷᚐௌȷʖยƴӼƚƨ
˯ޖʏƷᚐௌ ʖยƴӼƚƨ
ൢᝋᚐௌඥƷನሰ
(1) ൢᝋȢȇȫƱȩȸǸǨȇǣǷȥȟȬȸ Ƿ ǛኵLjӳǘƤƨ᭗ᚐࡇൢᝋᚐௌ
ǷȧȳǛኵLjӳǘƤƨ᭗ᚐࡇൢᝋᚐௌ
(2)
◊✲ᡭἲ 㧗ゎീᗘẼ㇟ゎᯒᡭἲ
◊✲ᡭἲ 㧗ゎീᗘẼ㇟ゎᯒᡭἲ
ٻ᪸؏ƷᚘምኽௐǛݱ᪸؏ƷᚘምƴМဇ ٻ᪸؏ƷᚘምኽௐǛݱ᪸؏ƷᚘምƴМဇ
ൢᝋʖย ȇ
ൢᝋʖยȢȇȫ ȩȸǸǨȇǣ
ǷȟȥȬȸǷȧȳ
5
˯ޖʏƷᚐௌǛӧᏡƱƢǔ
◊✲ᡭἲ ᗉෆ✵ ࿘㎶ᆅᙧ
◊✲ᡭἲ ᗉෆ✵ ࿘㎶ᆅᙧ
ɯᲢ᭗ƞኖ 60m)
ᆰลȓȫሁ
࠼࠘؏ȬȸȀ
๖ឥែᲴ RWY09/27 2000m
⤖ᯝ (1) ୣ䛾ᙳ㡪 - 䝺䞊䝎ほ ್䛸ẚ㍑
⤖ᯝ (1) ୣ䛾ᙳ㡪 䝺 䝎ほ ್䛸ẚ㍑
࠼࠘؏ȬȸȀᚇยᲢˏᚌ4.5ࡇᲣ RWY27Ʒᡶλኺែɥᡮ
ɯ
᭗ᚐࡇൢᝋᚐௌኽௐᲢˏᚌ4.5ࡇᲣ
ɯ
LESኽௐǛȬ ȀᚇยƱ
ȀᚇยƱӷಮ ͼӼǛ 1.5[(m/s)/m] ƷǦǤȳȉǷǢ
LESኽௐǛȬȸȀᚇยƱ
ӷಮƷ࢟ࡸƴƯЈщ
ȬȸȀᚇยƱӷಮƷͼӼǛ ᚐௌưਵƑǔƜƱưƖƨ
7
⤖ᯝ (2) ✵ 䝡䝹䛾ᙳ㡪
⤖ᯝ (2) ✵ 䝡䝹䛾ᙳ㡪
(1)
ᆰลȓȫ
(2)(1)
RWY27 ᡶλ
3[(m/s)/m] Ʒ
䜎䛸䜑 㧗ゎീᗘẼ㇟ゎᯒ
䜎䛸䜑 㧗ゎീᗘẼ㇟ゎᯒ
᭗ᚐࡇൢᝋᚐௌǛဇƍƯŴ˯ޖʏʙ̊ Ŵ ƷᚐௌǛᘍƬƨ
ȷ˯ޖʏǛਵƑǔƜƱƕưƖƨ
ႆဃᙲ׆Ჴɯ ᆰลȓȫ
ႆဃᙲ׆ᲴɯŴᆰลȓȫ ȷȬ Ȁ ᚇยƱǑƍɟᐲǛᅆƠƨ ȷȬȸȀȸᚇยƱǑƍɟᐲǛᅆƠƨ
ǦǤȳȉǷǢƷႆဃ᭗ࡇŴ
٭ѣƷࠢǛਵƑƨ
9
◊✲ᡭἲ 䝸䜰䝹䝍䜲䝮ண
◊✲ᡭἲ 䝸䜰䝹䝍䜲䝮ண
᭗ᚐࡇൢᝋᚐௌᲴᚘምdzǹȈƕ᩼ࠝƴٻƖƍ
ൢᝋžʖยſǛᘍƏƜƱƸᩊ
˯ޖʏǛʖยƢǔƜƱƸưƖƳ ƍƩǖƏƔᲹ
ƍƩǖƏƔᲹ
(1) ȇ ǿӷ҄
(1) ȇȸǿӷ҄
ĺ ᚇยȇȸǿƱǷȟȥȬȸǷȧȳǛኵLjӳǘƤǔඥ (2) ஊႺʩЎᚐ (Proper Orthogonal Decomposition: POD) ĺ ٻᙹȇȸǿƔǒɼЎǛࡽƖЈƢወᚘܖႎඥ
˯ޖʏƷȪǢȫǿǤȠʖยƴӼƚƯ
◊✲ᡭἲ 䝕䞊䝍ྠ䛸䛿
◊✲ᡭἲ 䝕 䝍ྠ䛸䛿
ᚇยƱǷȟȥȬȸǷȧȳǛᗡӳƢǔඥ ᚇยƱǷȟȥȬ ǷȧȳǛᗡӳƢǔඥ
ǷȟȥȬȸǷȧȳ
ȇȸǿӷ҄ƋǓ
ᚇย͌
ᚇย͌ ȇȸǿ
ӷ҄
ᚇย͌Ǜ̅ƬƯǷȟȥȬ ǷȧȳǛ̲ദƢǔ ӷ҄
11
ᚇย͌Ǜ̅ƬƯǷȟȥȬȸǷȧȳǛ̲ദƢǔ
◊✲ᡭἲ ᅛ᭷┤ศゎ (POD)
P O th l D iti◊✲ᡭἲ ᅛ᭷┤ศゎ (POD)
ٶഏΨȇȸǿƔǒ˯ഏΨȇȸǿǛਁЈƢǔඥ
Proper Orthogonal Decomposition
ĺ ᙐᩃƳ්ǕئƔǒཎࣉႎƳನᡯǛǨȍȫǮȸƝ ƱƴਁЈƢǔƜƱƕưƖǔ
▐ሙ 1
▐ሙ 2 ▐ሙ k
▐ሙ 1
▐ሙ 2 ▐ሙ k
◊✲ 䝕䞊䝍ྠᡭἲ䜢⏝䛔䛯ゎᯒ
◊✲ 䝕 䝍ྠᡭἲ䜢⏝䛔䛯ゎᯒ
ჺ᧓έƷ˯ޖʏǛʖยƢǔǷǹȆȠ ჺ᧓έƷ˯ޖʏǛʖยƢǔǷǹȆȠ
Ƿȟ Ȭ Ƿ ȳ ǷȟȥȬȸǷȧȳ
ஊႺʩЎᚐƴǑǔ
ஊႺʩЎᚐ Ǒǔ ቇତȢȇȫ
ȇ ǿӷ҄ඥᲢቩ܇ȕǣȫǿᲣ
ൢᝋʖยȢȇȫƩƚưƸʖإưƖƳƍ
ȇȸǿӷ҄ඥᲢቩ܇ȕǣȫǿᲣ
13
ൢᝋʖยȢȇȫƩƚưƸʖإưƖƳƍ
˯ޖʏǛʖยǛႸਦƢ
◊✲ᡭἲ Ꮚᐇ㦂䛾タᐃ
◊✲ᡭἲ Ꮚᐇ㦂䛾タᐃ
᭗ᚐࡇൢᝋᚐௌ POD (2010/2/6) POD ؕࡁșǯȈȫ ᭗ᚐࡇൢᝋᚐௌ
ȷJMA-NHM ȷLES
ȷ POD ؕࡁșǯȈȫ
ȷ Temporal Coeficient
ȬȸȀƱ ӷಮƷ࢟ࡸ
2010/1/25
Ⴧ͌ ᚇย͌ ቩ܇ȕǣȫǿ
Ⴧ͌ ᚇ
ൔ᠋ ਖ਼ܭኽௐ
᭗ᚐࡇൢᝋᚐௌƷኽௐǛžჇ͌ſƱƢǔ
ൔ᠋ ਖ਼ܭኽௐ
ȷ᭗ᚐࡇൢᝋᚐௌƷኽௐǛžჇ͌ſƱƢǔ ȷჇ͌ƔǒŴᚇย͌Ǜွ˩ႎƴ˺ǓЈƢ
ቩ܇ȕǣȫǿǛဇƍƯ Ⴧ͌ƴᡈƮƘƔǛ౨ᚰ
ȷቩ܇ȕǣȫǿǛဇƍƯŴჇ͌ƴᡈƮƘƔǛ౨ᚰ
⤖ᯝ (3) ᩘ್ᐇ㦂䛾⤖ᯝ
⤖ᯝ (3) ᩘ್ᐇ㦂䛾⤖ᯝ
ȇȸǿӷ҄Ǜ10ׅᘍƬƨᨥƷ᭗ࡇ40mƷ൦ᡮЎࠋ
Ⴧ͌ POD
JMA ࢽǒǕƨ්ǕئƷڲࣱ࢘Ǜᚸ̖Ƣ ǔƨNJƴ R tM S ᚘ JMA-NHM ǔƨNJƴŴRootMeanSquareƷᚘ
ምǛᘍƏŵ
¦
(u u )2 ᱁ᏊⅬᩘU
LES = 4.84g10-2 POD = 4.98g10-2
¦
(uuave) ᱁ᏊⅬᩘURMS
15
ૺ᩿ϋƷʏ්ࢍࡇƸಒƶɟᐲƠƯƍǔŵ
⤖ᯝ (2) ᩘ್ᐇ㦂䛾⤖ᯝ
⤖ᯝ (2) ᩘ್ᐇ㦂䛾⤖ᯝ
ᑋᆰೞƷᡶλኺែɥƷദݣ
ǦǤȳȉǷǢ
ȷ 4m/s Ɣǒ 12m/s ƴ
ദݣƕ࣯٭
䜎䛸䜑 䜎䛸䜑
(1) ᭗ᚐࡇൢᝋᚐௌඥƷನሰ
ȷ˯ޖʏƷžᚐௌſǛӧᏡƱƠƨ ˯ޖʏ ᚐௌſǛӧᏡ Ơ ȷȬȸȀᚇยƱಒƶɟᐲƠƨ
(2) ȇȸǿӷ҄ඥǛဇƍƨʖยඥƷನሰ
˯ޖʏƷžʖยſǛႸਦƢ ȷ˯ޖʏƷžʖยſǛႸਦƢ
ȷૠ᬴ܱ͌ưƸʏ්ࢍࡇŴǦǤȳȉǷǢƷ ˮፗǛʖยƢǔƜƱƕưƖƨ
17
☜ᐇ⎔ቃୗ䛷䛾䝁䞁䝣䝸䜽䝖䝣䝸䞊
☜ᐇ⎔ቃୗ䛷䛾 䞁 䝸䜽䝖 䝸 4 ḟඖ㌶㐨⏕ᡂ䜰䝹䝂䝸䝈䝮
ᮾிᏛᏛ㝔
⯟✵ᏱᐂᕤᏛᑓᨷ ᯇ㔝㈡ᐉ䠈ᅵᒇṊྖ
⫼ᬒ
2
୰㛫ሗ࿌
2013/3/21
⫼ᬒ
ᅾ ⯟✵⟶ไ
• ⌧ᅾ䛾⯟✵⟶ไ䝅䝇䝔䝮
¾
⯟✵⟶ไᐁ䛜㐠⏝
¾
✵ᇦ䛾Ᏻ䜢┘ど
¾
⯟✵ᶵ㛫䛾䝁䞁䝣䝸䜽䝖䜢ᅇ㑊
⫼ᬒ
⫼ᬒ
⟶ 䜢ᑟ ḟ ㌶㐨䝧
• 㛫⟶⌮䜢ᑟධ䛧䛯 4 ḟඖ㌶㐨䝧䞊䝇㐠⏝
¾
䛩䜉䛶䛾⯟✵ᶵ䛜䛭䜜䛮䜜䛾 4 ḟඖ㌶㐨䛻ᚑ䛖
䝁䞁䝣䝸䜽䝖䝣䝸䞊
¾ ⯟✵ᶵ㛫䛾Ᏻ䜢☜ಖ
䜴䜵䜲䝫䜲䞁䝖
¾ ⯟✵ᶵ㛫䛾Ᏻ䜢☜ಖ
¾ 䝁䞁䝣䝸䜽䝖䝣䝸䞊㌶㐨䜢㣕⾜
Æ䝁䞁䝣䝸䜽䝖䝣䝸䞊 4 ḟඖ㌶㐨᭱㐺
┠ⓗ
4
୰㛫ሗ࿌
2013/3/21
┠ⓗ
ḟ ㌶㐨 㐺
• 䝁䞁䝣䝸䜽䝖䝣䝸䞊 4 ḟඖ㌶㐨᭱㐺
¾
㣕⾜୰䛾☜ᐇ⎔ቃ䜢⪃៖
Ẽ㇟ሗ䛾☜䛛䛥
⯟ἲㄗᕪ
¾
☜⋡ⓗ䛺⯟✵ᶵ䛾㐠ື䝰䝕䝹
Æ ☜⋡ㄽⓗ㌶㐨᭱㐺
¾
☜ᐇせ⣲䜢⪃៖䛧䛯䝻䝞䝇䝖䛺㣕⾜㌶㐨䜢⏕ᡂ
☜ᐇ⎔ቃ
5
୰㛫ሗ࿌
2013/3/21
☜ᐇ⎔ቃ
• 2 ḟඖᖹ㠃ෆ • ☜ᐇせ⣲
1.
㢼㏿䛾ண ㄗᕪ
㢼㏿
2.
ᶵయ㏿ᗘ䛾 ᐃㄗᕪ
ᶵయ㏿ᗘ
㢼㏿
ᶵ
3.
⌧ᅾ⨨䛾᥎ᐃㄗᕪ ⌧ᅾ⨨ ᥎ᐃㄗᕪ
Æ ☜⋡ኚᩘ䛸䛧䛶ᢅ䛖 Æ ☜⋡ኚᩘ䛸䛧䛶ᢅ䛖
⌧ᅾ⨨
☜ᐇ⎔ቃ
6
୰㛫ሗ࿌
2013/3/21
☜ᐇ ⣲ ☜⋡ኚᩘ
☜ᐇ⎔ቃ
• ☜ᐇせ⣲䠄☜⋡ኚᩘ䠅
¾
㢼㏿䛾ண ㄗᕪ
• 5.17 m/s
¾
ᶵయ㏿ᗘ䛾 ᐃㄗᕪ
• 5 knots (2.57 m/s)
⌧ᅾ⨨䛾᥎ᐃㄗᕪ
☜⋡ㄽⓗᡭἲ
☜⋡ㄽⓗᡭἲ
ᑗ᮶ 䜢䜺䜴 ศᕸ ⾲
1. ᑗ᮶䛾⨨䜢䜺䜴䝇ศᕸ䛷⾲⌧
¾
ㄗᕪయ䜢䜺䜴䝇ศᕸ䛾ᶆ‽೫ᕪ䛸䛧䛶䝰䝕䝹
ඹ ቑ
䠄㛫䛸ඹ䛻ቑ䠅
㢼㏿䛾ண ㄗᕪ
㐍⾜᪉ྥ䠖
୍ᐃ್䛻᮰䛥䛫䜛
┤᪉ྥ䠖
୍ᐃ್䛻᮰䛥䛫䜛 䠄ᐇ㝿䛾᧯⦪䜢ᶍᨃ䠅
ᶵయ㏿ᗘ䛾 ᐃㄗᕪ
⌧ᅾ⨨䛾᥎ᐃㄗᕪ
ᶵయ㏿ᗘ䛾 ᐃㄗᕪ
☜⋡ㄽⓗᡭἲ
2013/3/21 ୰㛫ሗ࿌ 8
☜⋡ㄽⓗᡭἲ
ᑗ᮶ 䜢䜺䜴 ศᕸ ⾲
1. ᑗ᮶䛾⨨䜢䜺䜴䝇ศᕸ䛷⾲⌧
¾
ᶵయ㏿ᗘ䠖 600 knots
10 䢢
¾
10 ศᚋ䛻 45 ᗘ᪉䜢ኚ᭦
-10 0
M
t = 10 min.
t = 5 min.
-30
NM -20
15 i
๓䛻ィ⏬䛥䜜䛯㌶㐨
0 20 40 60 80 100 120 140 160
-50 -40
NM
䢢
t = 15 min.
NM
☜⋡ㄽⓗᡭἲ
2013/3/21 ୰㛫ሗ࿌ 9
☜⋡ㄽⓗᡭἲ
ᑗ᮶ 䜢䜺䜴 ศᕸ ⾲
1. ᑗ᮶䛾⨨䜢䜺䜴䝇ศᕸ䛷⾲⌧
¾
ᶵయ㏿ᗘ䠖 600 knots
¾
10 ศᚋ䛻 45 ᗘ᪉䜢ኚ᭦
๓䛻ィ⏬䛥䜜䛯㌶㐨
☜⋡ㄽⓗᡭἲ
2013/3/21 ୰㛫ሗ࿌ 10
☜⋡ㄽⓗᡭἲ
ᑗ᮶ 䜢䜺䜴 ศᕸ ⾲
1. ᑗ᮶䛾⨨䜢䜺䜴䝇ศᕸ䛷⾲⌧
¾
」ᩘᶵ䛾ศᕸ䜢ồ䜑䜛䛣䛸䛷䝁䞁䝣䝸䜽䝖䜢⪃៖
䝁䞁䝣䝸䜽䝖䝣䝸䞊 䝁䞁䝣䝸䜽䝖䝣䝸
☜⋡ㄽⓗᡭἲ
☜⋡ㄽⓗᡭἲ
ᑗ᮶ 䜢䜺䜴 ศᕸ ⾲
1. ᑗ᮶䛾⨨䜢䜺䜴䝇ศᕸ䛷⾲⌧
¾
⯟✵ᶵ䛾㐠ື䜢༢⣧䛧䛯䝰䝕䝹
¾
ไᚚධຊ䜢⪃៖䛧䛶䛔䛺䛔
•
䛘䜀䠈᪕ᅇ䛸䛔䛳䛯⯟✵ᶵ䛾㐠ື䜢 ヲ⣽䛻ᢅ䛖䛣䛸䛜ᅔ㞴
Æ☜⋡㐣⛬䜢ྵ䜐㐠ື᪉⛬ᘧ䜢⪃䛘䜛
☜⋡ㄽⓗᡭἲ
12
୰㛫ሗ࿌
2013/3/21
☜⋡ㄽⓗᡭἲ
2. ☜⋡ⓗ䛺䝎䜲䝘䝭䜽䝇
䠄☜⋡ᚤศ᪉⛬ᘧ䠅 • ☜⋡㐣⛬䜢ྵ䜐᪉⛬ᘧ ☜⋡㐣⛬䜢ྵ䜐᪉⛬ᘧ
㢼㏿䛾ண ㄗᕪ
• ☜⋡ㄽⓗゎ
•
⤫ィⓗሗ
¾ᮇᚅ್
ᶵయ㏿ᗘ䛾 ᐃㄗᕪ
¾ศᩓ
ᶵయ㏿ᗘ䛾 ᐃㄗᕪ
☜⋡ㄽⓗᡭἲ
13
୰㛫ሗ࿌
2013/3/21
☜⋡ㄽⓗᡭἲ
• 䝰䞁䝔䜹䝹䝻ἲ ἲ
¾
☜⋡ኚᩘ䛾䝷䞁䝎䝮䝃䞁䝥䝸䞁䜾
¾
ྛ䝃䞁䝥䝹Ⅼ䛷Ỵᐃㄽⓗၥ㢟䛻᭩䛝䛘䜛
¾
䝃䞁䝥䝹Ⅼ䛷䛾ゎ୍ᘧ䛛䜙
☜⋡ㄽⓗゎ䠄ᮇᚅ್䞉ศᩓ䠅䜢ồ䜑䜛
¾
㐺⏝䛧䜔䛩䛔䛜䠈ィ⟬㈇Ⲵ䛜㧗䛔
☜⋡ㄽⓗᡭἲ
14
୰㛫ሗ࿌
2013/3/21
☜⋡ㄽⓗᡭἲ
⯡ከ㡯ᘧ䜹䜸䝇
• ୍⯡ከ㡯ᘧ䜹䜸䝇
¾
䝰䞁䝔䜹䝹䝻ἲ䛸ྠᵝ䠈㐺⏝䛧䜔䛩䛔
¾
ィ⟬㈇Ⲵ䜢ῶᑡ䛥䛫䜛≉ᚩ
1ゎ䛾ከ㡯ᘧ㏆ఝ
1.
ゎ䛾ከ㡯ᘧ㏆ఝ
¾ゎ䜢☜⋡ኚᩘ䛾┤ከ㡯ᘧ䛸䛧䛶⾲䛩
䝅 䝅 ၥ㢟タᐃ 䝅䝭䝳䝺䞊䝅䝵䞁 - ၥ㢟タᐃ
• 4 ḟඖ㌶㐨᭱㐺
¾
2 ḟඖᖹ㠃ෆ䠈㧗ᗘ 34000 ft
¾≧ែኚᩘ䠖 ⤒ᗘ䠈⦋ᗘ䠈 ᶵ㤳᪉ゅ
⦋ᗘ 䜴䜵䜲䝫䜲䞁䝖
(139.875, 35.1)
¾
ไᚚኚᩘ䠖 䝞䞁䜽ゅ
ᶵ㤳᪉ゅ ᗘ
⤊➃㛫ᅛᐃ
¾
┠ⓗ㛵ᩘ䠖
⌧ᅾ⨨ ⛊ (139.375, 34.7)
⤊➃㛫ᅛᐃ
᭱ᑠ
⤒ᗘ
= 500 knots
䝅 䝅 ၥ㢟タᐃ
2013/3/21 ୰㛫ሗ࿌ 16
䝅䝭䝳䝺䞊䝅䝵䞁 - ၥ㢟タᐃ
¾
㢼㏿䛾ண
• 4 ḟඖ㌶㐨᭱㐺
¾
2 ḟඖᖹ㠃ෆ䠈㧗ᗘ 34000 ft
⦋ᗘ 䜴䜵䜲䝫䜲䞁䝖
(139.875, 35.1)
ᶵ㤳᪉ゅ ᗘ
⤊➃㛫ᅛᐃ
⌧ᅾ⨨ ⛊ (139.375, 34.7)
⤊➃㛫ᅛᐃ
•
Ẽ㇟ᗇ䛻䜘䜚බ㛤
•
⦋ᗘ䞉⤒ᗘ䛾䜾䝸䝑䝗䝕䞊䝍
•
⿵㛫䛧䛶㐺⏝
= 500 knots
䝅 䝅 ⤖ᯝ
2013/3/21 ୰㛫ሗ࿌ 17
䝅䝭䝳䝺䞊䝅䝵䞁 - ⤖ᯝ
•
୍⯡ከ㡯ᘧ䜹䜸䝇䛸䝰䞁䝔䜹䝹䝻ἲ䛸䛾ẚ㍑
¾
㛫 - ⤒ᗘ䠄ᮇᚅ್䠅
¾㛫 - ⦋ᗘ䠄ᮇᚅ್䠅
䝅 䝅 ⤖ᯝ
2013/3/21 ୰㛫ሗ࿌ 18
䝅䝭䝳䝺䞊䝅䝵䞁 - ⤖ᯝ
•
୍⯡ከ㡯ᘧ䜹䜸䝇䛸䝰䞁䝔䜹䝹䝻ἲ䛸䛾ẚ㍑
¾
㛫 - ⤒ᗘ䠄ศᩓ䠅
¾㛫 - ⦋ᗘ䠄ศᩓ䠅
䝅 䝅 ⤖ᯝ 䝅䝭䝳䝺䞊䝅䝵䞁 - ⤖ᯝ
•
୍⯡ከ㡯ᘧ䜹䜸䝇䛸䝰䞁䝔䜹䝹䝻ἲ䛸䛾ẚ㍑
ᡭἲ 䝃䞁䝥䝹Ⅼ ィ⟬㛫
୍⯡ከ㡯ᘧ䜹䜸䝇 401 14.5ศ 䝰䞁䝔䜹䝹䝻ἲ 20000 10.5㛫
䝅 䝅 ⤖ᯝ
2013/3/21 ୰㛫ሗ࿌ 20
䝅䝭䝳䝺䞊䝅䝵䞁 - ⤖ᯝ
•
୍⯡ከ㡯ᘧ䜹䜸䝇䛸䝰䞁䝔䜹䝹䝻ἲ䛸䛾ẚ㍑
¾
⤒ᗘ - ⦋ᗘ䠄ᮇᚅ್䠅
┤⥺㣕⾜
㎽ᅇ
䝅 䝅 ⤖ᯝ
2013/3/21 ୰㛫ሗ࿌ 21
䝅䝭䝳䝺䞊䝅䝵䞁 - ⤖ᯝ
¾
⤒ᗘ - ⦋ᗘ 䠄⤒ᗘ䛾ㄗᕪ䠅
¾
㛫 - ⤒ᗘ䠄ศᩓ䠅
䝅 䝅 ⤖ᯝ
2013/3/21 ୰㛫ሗ࿌ 22
䝅䝭䝳䝺䞊䝅䝵䞁 - ⤖ᯝ
¾
⤒ᗘ - ⦋ᗘ 䠄⦋ᗘ䛾ㄗᕪ䠅
¾
㛫 - ⦋ᗘ䠄ศᩓ䠅
䜎䛸䜑 䜎䛸䜑
☜⋡ㄽⓗ ḟ ㌶㐨 㐺
• ☜⋡ㄽⓗ 4 ḟඖ㌶㐨᭱㐺
¾
☜⋡ㄽⓗᡭἲ
•
㣕⾜୰䛾☜ᐇせ⣲䜢⪃៖
•
☜⋡㐣⛬䜢ྵ䜐㐠ື᪉⛬ᘧ
¾
୍⯡ከ㡯ᘧ䜹䜸䝇
䜹 ἲ䛸 ᵝ 㐺 䜔䛩
•
䝰䞁䝔䜹䝹䝻ἲ䛸ྠᵝ䠈㐺⏝䛧䜔䛩䛔
•
ィ⟬㈇Ⲵ䜢ᖜ䛻ῶᑡ
¾
☜⋡ㄽⓗゎ
㌶㐨䛾ᮇᚅ್䛸ศᩓ䛾㛫ᒚṔ
•
㌶㐨䛾ᮇᚅ್䛸ศᩓ䛾㛫ᒚṔ
ᚋ ணᐃ
24
୰㛫ሗ࿌
2013/3/21
ᚋ䛾ணᐃ
ḟ ✵ ື
• 3 ḟඖ✵㛫䛷䛾㐠ື
• 」ᩘᶵ䛻䛴䛔䛶⪃៖
¾
䝍䞊䝭䝘䝹✵ᇦ䛻䛚䛡䜛䝬䞊䝆䞁䜾
¾
䝁䞁䝣䝸䜽䝖䝣䝸䞊㌶㐨䜢⏕ᡂ
• 㛫䛻౫Ꮡ䛩䜛☜⋡ኚᩘ䛾ᑟධ 㛫䛻౫Ꮡ䛩䜛☜⋡ኚᩘ䛾ᑟධ
ᖹᡂ䠎䠐ᖺᗘḟୡ௦㐠⯟䝅䝇䝔䝮බເᆺ◊✲
ᖹᡂ䠎䠐ᖺᗘḟୡ௦㐠⯟䝅䝇䝔䝮බເᆺ◊✲
◊✲ㄢ㢟ྡ䠖 䝬䝹䝏㻳㻺㻿㻿㧗⢭ᗘ 䛚䜘䜃䛭䛾ಙ㢗ᛶホ౯䜰䝹䝂 䝸䝈䝮䛻㛵䛩䜛◊✲
◊✲௦⾲⪅䠖 ஂಖᖾᘯ䠄❧㤋Ꮫ ⌮ᕤᏛ㒊㟁Ẽ㟁ᏊᕤᏛ⛉䠅
ᖹᡂ㻞㻠ᖺᗘḟୡ௦㐠⯟䝅䝇䝔䝮බເᆺ◊✲୰㛫ሗ࿌
㻞㻜㻝㻟ᖺ㻟᭶㻞㻝᪥ 䠖㻶㻭㼄㻭ㄪᕸ⯟✵Ᏹᐂ䝉䞁䝍䞊
1/26
㻯㻻㻺㼀㻱㻺㼀㻿
◊✲䛾⫼ᬒ䞉┠ⓗ䞉ᴫせ
◊✲ ⫼ᬒ ┠ⓗ ᴫせ 䝬䝹䝏㻳㻺㻿㻿⢭ᐦ༢⊂
䛾ಙ㢗ᛶホ౯䜰䝹䝂䝸䝈䝮䛾㛤Ⓨ
㻳㻺㻿㻿㻛㻵㻺㻿」ྜ⯟ἲ䜰䝹䝂䝸䝈䝮䛾㛤Ⓨ
㻳㻺㻿㻿㻛㻵㻺㻿」ྜ⯟ἲ䜰䝹䝂䝸䝈䝮䛾㛤Ⓨ
䜎䛸䜑
GNSS⾨ᫍᩘண
2010ᖺ௦ ๓༙
2010ᖺ௦ ᚋ༙
2020ᖺ௦ ๓༙
GPS 32 32 32
GPS 32 32 32
GLONASS 24(+3) 24(+3) 24(+3)
GALILEO 4 18 27(+3)
COMPASS 12 30 32
QZS 1 4 7
ྜィ 73 108 122
zྛᅜ䛻䛚䛡䜛 ⾨ᫍ䛾㛤Ⓨ䛜㐍䜑䜙䜜䜜䜀 ከ䝅䝇䝔䝮 ከ zྛᅜ䛻䛚䛡䜛 ⾨ᫍ䛾㛤Ⓨ䛜㐍䜑䜙䜜䜜䜀䠈ከ䝅䝇䝔䝮䠈ከ
࿘Ἴ䠄䝬䝹䝏GNSS䠅䛾⎔ቃ䛜⏝ྍ⬟
z⯟✵ᶵ䛻䛚䛡䜛㧗⢭ᗘ⾨ᫍ⯟ἲᢏ⾡䛾◊✲㛤Ⓨ䛻䛚䛔䛶䜒 䝬
3/26
z⯟✵ᶵ䛻䛚䛡䜛㧗⢭ᗘ⾨ᫍ⯟ἲᢏ⾡䛾◊✲㛤Ⓨ䛻䛚䛔䛶䜒䠈䝬 䝹䝏GNSS⎔ቃ䜢ぢᤣ䛘䛯ᑐᛂ䛜ྍḞ
◊✲䛾ᴫせ䛸┠ⓗ
z䝬䝹䝏GNSS⎔ቃୗ䛻䛚䛡䜛 䐟 㧗⢭ᗘ 䜰䝹䝂䝸䝈䝮 䐟 㧗⢭ᗘ 䜰䝹䝂䝸䝈䝮
䐠 䠄⯟ἲ䠅䛾ಙ㢗ᛶホ౯䜰䝹䝂䝸䝈䝮 䐡GNSS/INS」ྜ䜰䝹䝂䝸䝈䝮
䐡GNSS/INS」ྜ䜰䝹䝂䝸䝈䝮
䠄䝔䞊䝬䐟䠖㧗⢭ᗘ 䜰䝹䝂䝸䝈䝮䠅
ᚑ᮶䛾䝁䞊䝗ᨃఝ㊥㞳䛻ຍ䛘䛶 ᦙ㏦Ἴ┦ほ 㔞䜒✚ᴟⓗ䛻⏝䛧䛯⢭ᐦ༢⊂ ᚑ᮶䛾䝁 䝗ᨃఝ㊥㞳䛻ຍ䛘䛶䠈ᦙ㏦Ἴ┦ほ 㔞䜒✚ᴟⓗ䛻⏝䛧䛯⢭ᐦ༢⊂ 䠄PPP䠅䜢ά⏝䛧䛯䝬䝹䝏GNSS-PPP䛾㛤Ⓨ
䠄䝔䞊䝬䐠䠖 䛾ಙ㢗ᛶホ౯䜰䝹䝂䝸䝈䝮䠅
ྛ䝅䝇䝔䝮䛾≉ᛶ䠄ㄗᕪ䜔䝞䜲䜰䝇䛾⤫ィⓗᛶ㉁䠅䜢⪃៖䛷䛝䜛䝣䜱䝹䝍䠄᥎ᐃჾ䠅䜢ά⏝䛧
ྛ䝅䝇䝔䝮䛾≉ᛶ䠄ㄗᕪ䜔䝞䜲䜰䝇䛾⤫ィⓗᛶ㉁䠅䜢⪃៖䛷䛝䜛䝣䜱䝹䝍䠄᥎ᐃჾ䠅䜢ά⏝䛧䠈 ㄗᕪ䛾ศᩓ䠄ඹศᩓ䠅䜢ṇ☜䛻ぢ✚䜒䜛䛣䛸䛾ฟ᮶䜛䜰䝹䝂䝸䝈䝮㛤Ⓨ
䠄䝔䞊䝬䐡䠖GNSS/INS」ྜ䜰䝹䝂䝸䝈䝮䠅
✀䚻䛾㠀⥺ᙧ䝣䜱䝹䝍䜢ẚ㍑᳨ウ䛧 ⯟ἲ䛾⢭ᗘ䛾䜏䛺䜙䛪 ಙ㢗ᛶホ౯䛾ほⅬ䛛䜙䜒᭷
✀䚻䛾㠀⥺ᙧ䝣䜱䝹䝍䜢ẚ㍑᳨ウ䛧䠈⯟ἲ䛾⢭ᗘ䛾䜏䛺䜙䛪䠈ಙ㢗ᛶホ౯䛾ほⅬ䛛䜙䜒᭷
䛸䛺䜛」ྜ䜰䝹䝂䝸䝈䝮䛾㛤Ⓨ
zḟୡ௦㐠⯟䝅䝇䝔䝮䛾㧗⢭ᗘ⾨ᫍ⯟ἲᢏ⾡㛤Ⓨ䛻䛚䛡䜛 䝬䝹䝏GNSS⏝䛾ᢏ zḟୡ௦㐠⯟䝅䝇䝔䝮䛾㧗⢭ᗘ⾨ᫍ⯟ἲᢏ⾡㛤Ⓨ䛻䛚䛡䜛䠈䝬䝹䝏GNSS⏝䛾ᢏ
⾡☜❧䛻㈉⊩䛩䜛䛣䛸䜢┠ⓗ䛸䛩䜛
㧗⢭ᗘ 䜰䝹䝂䝸䝈䝮䛾㛤Ⓨ
ᚑ᮶䛾༢⊂ 䠄Aug. 11, 2005, 15:00—16:00 (UTC)䠅
z⢭ᐦ༢⊂ 䛻䛴䛔䛶
䞉S.Sugimoto:GNSSRegressiveModelsandPrecisePointPositioning,Proc.
ofthe36thISCIEInt.Symp.onStochasticSystemsTheoryandIts Applications,pp.159Ͳ164,2004.
ᨃఝ㊥㞳䛻ຍ䛘䛶 ᦙ㏦Ἴ┦ほ 䜢⏝䛩䜛
PPP ⤖ᯝ
ᨃఝ㊥㞳䛻ຍ䛘䛶䠈ᦙ㏦Ἴ┦ほ 䜢⏝䛩䜛 䛣䛸䛷䠈RTK᪉ᘧ䛻༉ᩛ䛩䜛⢭ᗘ䜢ᐇ⌧ྍ⬟
䝬䝹䝏
GNSS
ᅗ: PPP 䛾䜲䝯䞊䝆䠄❧㤋ᏛBKCW.W.7㝵ᒇୖ䠅
䝬䝹䝏
GNSS
䞉GLONASS䛾ά⏝
GPS䛸䛾┦㐪Ⅼ 䞉࿘Ἴᩘ
䞉้⣔้⣔
䞉ᗙᶆ⣔
5/26
㻳㻼㻿㻛㻳㻸㻻㻺㻭㻿㻿」ྜほ ᪉⛬ᘧ
䠘ホ౯᪉ἲ䠚 䠘ほ 䞉 ᮲௳䠚 䞉GPS/GLONASS」ྜ⯟ἲ
䞉GPS༢⊂ GLONASS༢⊂ 䞉GLONASS༢⊂
7/26
㻳㻼㻿㻛㻳㻸㻻㻺㻭㻿㻿」ྜ ⤖ᯝ䛾
㻳㻼㻿㻛㻳㻸㻻㻺㻭㻿㻿」ྜ 䛾Ⅼ
ྍど⾨ᫍ(㟷䠖GPS ㉥䠖GLONASS)
ྍど⾨ᫍᩘ䛾ኚ
9/26
ಙ㢗ᛶ䛾⿵ᙉ䝅䝇䝔䝮
Ẹ㛫⯟✵ᶵ䛾⯟ἲ䛻䛚䛔䛶ಙ㢗ᛶ䜢⿵ᙉ䛩䜛᪉ἲ
ABAS䠄AircraftBasedAugmentationSystem䠅
ᶵୖ⨨䛻䜘䜛⿵ᙉ
Ẽᅽ㧗ᗘィ䠈INS䞉IRS 䠈RAIMཷಙᶵ
㟼Ṇ⾨ᫍ䛛䜙䛾⿵ᙉሗ䜢ᥦ౪
SBAS䠄SatelliteBasedAugmentationSystem䠅㟼Ṇ⾨ᫍ䛛䜙䛾⿵ᙉሗ䜢ᥦ౪
RAIM: Receiver Autonomous Integrity Monitoring
ཷಙᶵ䛻䜘䜛
⮬ᚊⓗ䛺䜲䞁䝔䜾䝸䝔䜱┘ど
RAIM: Receiver Autonomous Integrity Monitoring
⮬ᚊⓗ䛺䜲䞁䝔䜾䝸䝔䜱┘ど
⯟ἲ䝅䝇䝔䝮䛾ಖド⢭ᗘ䛜つᐃ䛥䜜䛯⢭ᗘ䛻‶䛯䛥䛺䛔ሙྜ䛻䝅 䝇䝔䝮⮬㌟䛜䛭䛾≧ἣ䜢᳨ฟ䛧䝴䞊䝄䛻㐺㆙࿌䜢Ⓨ䛩䜛⬟ຊ
つᐃ䛥䜜䛯⢭ᗘ
⯟ἲ䝅䝇䝔䝮䛾ಖド⢭ᗘ
㆙࿌Ⓨ⏕
11/26
䜲䞁䝔䜾䝸䝔䜱䜢☜ಖ䛩䜛䛯䜑䛻
䜰䝷䞊䝖䝸䝭䝑䝖 䠄 AL 䠅
⯟ἲ䝰䞊䝗䛻ᛂ䛨䛶
ICAO䛻䜘䛳䛶䛒䜙䛛䛨䜑ᐃ䜑䜙䜜䛯್
䜲䞁䝔䜾䝸䝔䜱 䜲䞁䝔䜾䝸䝔䜱
☜ಖ䛾せ௳
䝥䝻䝔䜽䝅䝵䞁䝺䝧䝹䠄 PL 䠅
ㄗᕪ䛜AL್௨ୗ䛷䛒䜛☜⋡䛜99.99999%௨ୖ䛸䛔䛖䛣䛸䜢ಖド䛩䜛
) 10 1 ( 7
ಙ㢗Ỉ‽䜢 䛸䛧䛯ሙྜ䛾 ⤖ᯝ䛾ಙ㢗㝈⏺䛾ᣦᶆ
㞀᳨ ᕪ䛜 ☜
ᨾ㞀᳨ฟ䜢⾜䛳䛯ୖ䛷䠈ᐃ䛥䜜䜛᭱ᝏ䛾 ㄗᕪ䛜PL䜢㉸䛘䜛☜⋡
䜰䝷䞊䝖䝸䝭䝑䝖䛸䝥䝻䝔䜽䝅䝵䞁䝺䝧䝹
₇⟬䛻䛚䛔䛶䛿䠈䝥䝻䝔䜽䝅䝵䞁䝺䝧䝹䜢ồ䜑䛶
AL<PL 䛾ሙྜ ㆙࿌Ⓨ⏕
AL>PL AL<PL
Alert
AL AL
AL PL AL
PL
HPLVPL
HPL
13/26
ほ 㞧㡢䛾ศᕸ
ᨃఝ㊥㞳ほ 㞧㡢䛾୍
ほ 㞧㡢䛾ศᕸ 䜺 䜺䜴䝇ศᕸ䜢௬ᐃ
䜹 䝍
䜹䝹䝬䞁䝣䜱䝹䝍
䜹䝹䝬䞁䝣䜱䝹䝍 䜺䜴䝇䝃䝮䝣䜱䝹䝍 ℐἼ䞉ண ᥎ᐃ್
䜹䝹䝬䞁䝀䜲䞁
᥎ ᩓ
᥎ᐃㄗᕪඹศᩓ⾜ิ
15/26
ᐇ㦂⎔ቃ
ほ ሙᡤ ඵᓥ㏆㎶
ཷಙᶵ/
䜰䞁䝔䝘 ASHTECHG12/ S67Ͳ1575Ͳ96 ほ 㛤ጞ㛫 2008/2/7 0䠖40䠖50䠄GPST䠅 ほ 㛤ጞ㛫 / / 䠄 䠅 ほ ⤊㛫 2008/2/7 2䠖08䠖14䠄GPST䠅
⏝䝕䞊䝍 CA䝁䞊䝗ᨃఝ㊥㞳, L1ᦙ㏦Ἴ┦
L1ᦙ㏦Ἴ┦
᪉ἲ 㧗⢭ᗘ༢⊂ 䠄䜻䝛䝬䝔䜱䝑䜽䠅
ᐇ㦂⎔ቃ
ほ 㞧㡢タᐃ
ほ 㞧㡢䛾ศᩓ
CA
䝁䞊䝗
0.54054䜹䝹䝬䞁䝣䜱䝹䝍
L1ᦙ㏦Ἴ┦ 0.000054054
ほ 㞧㡢䛾ศᩓ ほ 㞧㡢䛾ศᩓ ศᕸ䛾㔜䜏
0.595730 0.404270CA䝁䞊䝗 0.2525 0.9650
䜺䜴䝇䝣䜱䝹䝍
L1ᦙ㏦Ἴ┦ 0.00002525 0.00009650
䝟䝷䝯 䝍タᐃ 䝟䝷䝯䞊䝍タᐃ
ㄗ㆙ሗ⋡
.17/26
᳨ฟ⋡
.ᐇ㦂⤖ᯝ䠄 ⤖ᯝ䠅
㉥䠖ᆶ┤᪉ྥ䛾PL
㟷䠖Ỉᖹ᪉ྥ䛾
PL19/26
㉥䠖ᆶ┤᪉ྥ䛾PL
㟷䠖Ỉᖹ᪉ྥ䛾
PLᐇ⥺䠖䜺䜴䝇䝣䜱䝹䝍 ◚⥺䠖䜹䝹䝬䞁䝣䜱䝹䝍
㻳㻼㻿㻛㻵㻺㻿」ྜ⯟ἲ
Closed loop䠉Loosely coupled᪉ᘧ
9GPSሗ䜢⏝䛔䛶䠈INS䛾⯟ἲ ㄗᕪ䜢᥎ᐃ䠈⿵ṇ䛩䜛 9 ㄗᕪ䛿㠀⥺ 䛺䝎䜲䝘 9INSㄗᕪ䛿㠀⥺ᙧ䛺䝎䜲䝘䝭
䜽䝇䛾䛯䜑䠈㠀⥺ᙧ䝣䜱䝹䝍 䛜ᚲせ
䞉
Extended Kalman Filter (EKF)䞉
Quasi Linear Filter (QLF)G i S Filt (GSF)
䞉
Gaussian Sum Filter (GSF)䞉
Gaussian sum-Quasi Linear Filter (GQLF)Jorg ENTZINGER (Ph.D.)
ᮾிᏛᏛ㝔ᕤᏛ⣔◊✲⛉㻌 ຓᩍ ᅜ㝿ᕤᏛᩍ⫱᥎㐍ᶵᵓ㻌
䝞䜲䝸䞁䜺䝹䜻䝱䞁䝟䝇᥎㐍䝉䞁䝍䞊
[email protected]䝶䝹䜾㻌 㻌 䜶䞁䝖䝆䞁䜺䞊㻌 㻌 䠄ᕤᏛ༤ኈ䠅
What is the goal?
● Understand differences
in mental models and cognitive processes between curved and straight approaches
● Support the pilot
(through training or interfaces)
in his supervision of automation
and in decision making
Why curved approaches?
● Limit noise impact on the ground
Tokyo Haneda airport ILS approach (left) versus RNP-AR curved approach (right) 3
Why curved approaches?
● Limit noise impact on the ground
● Enabling approaches with challenging terrain
Radar-vectoring approach with ILS RNP-AR curved approach with
manual (FD) landing RNP-AR curved
approach with ILS Radar-vec
approach RNP-AR c approach manual (F RNP-AR c approach
5
What are we looking at?
RNP-AR curved approach with manual (FD) landing
Straight final:
FD (+ Auto Throttle) Curved part:
Auto Pilot + Auto Throttle
Situational Awareness
Decision Making
Cognitive Workload
What is the problem?
● Verification of navigation database during flight is practically impossible
● Workload reduction under normal conditions, but likely to increase with abnormalities
● Supervision of automation is more difficult
● Less time to confirm alignment with runway
○ Decision height may be in a curve
○ Go-around may be initiated in a curve
○ Location of runway in visual field uncertain
7What is the problem?
A curved approach is less safe than a straight approach
Mental model of curves
is more complex Worse decisions, less
stabilized, larger errors
Briefings and training of students by a
veteran airline pilot
Combining research and actual experience
9
How to measure cognitive workload?
● Subjective measures (e.g. Task Load Index)
How to measure cognitive workload?
● Subjective measures (e.g. Task Load Index)
● Pupil size
11
How to measure cognitive workload?
● Subjective measures (e.g. Task Load Index)
● Pupil size
● Heart-rate-variability y
● Research: mostly hardware or trajectory
● Pilot training: systems and procedures
● Missing: human factors
● Irregular cases are not yet considered
(challenging weather, one engine out, etc.)
● Late changes of plan due to mixed-mode airport operations occur
● Small differences between procedures,
certifications, and systems confuse crew
13Where do we go from here?
● What are the most demanding cases?
● What support/training is needed?
Challenges:
● Current simulator has no RNP-AR capability
● Still few pilots with RNP-AR experience
(especially in Japan)
指向性を考慮した騒音伝搬予測手法に関する研究
•騒音は航空機運航ルートの決定における大きな要因
•長距離伝搬に関する風や温度分布の影響と 航空機
•長距離伝搬に関する風や温度分布の影響と、航空機 騒音の指向特性の影響
•本研究では、音源の指向性を考慮した騒音の伝搬予 測計算法を波動数値解析により開発することを目的
•計算と実験によって音源指向性が騒音レベル分布に 及ぼす影響を調べる。
?
? ? ? ? ? ? ?
及ぼす影響を調べる。
東京大学生産技術研究所 坂本 慎一 東京大学生産技術研究所 坂本 慎一
(専門分野:応用音響工学)
指向性を考慮した騒音伝搬予測手法に関する研究
④全体研究計画
H24(2012) H25(2013) H24(2012) H25(2013)
(各年度目標) ①波動数値解析に対
する指向性の導入手 法 確立
②フィールドにおける 検証実験
③音源指向性デ タ 法の確立 ③音源指向性データ の取得と騒音伝搬特 性に及ぼすケースス タデ
(課題と計画)
①指向性を考慮した波動数値解析手法
タディ
①指向性を考慮した波動数値解析手法 の確立と整備
デ 向
using Fourier expansion (in 2‐D)
(アイデア):指向性の分解⇒合成をFDTDに取り込む
g p ( )
( ) ∑ { ( ) ( ) } ∑
∞( )
=
∞
=
= +
= Φ
0 0
cos sin
cos
m m
m Am mϕ Bm mϕ A mϕ
ϕ
h l h ( )
( ) 2 n + 1 ( n − m ) ( ) !
jusing spherical harmonic expansion (in 3‐D)
( ) ( )
( ) ( θ )
ϕϕ π
θ
nm jmnm
P e
m n
m n
Y n cos
!
! 4
1 , 2
+
= +
Express by initial distribution of pressure p y p
3
y i d fi d
i d fi d
1 FDTDM
課題:指向性を考慮した波動数値解析手法
( ) ( )
( ) ( )
{ }
∑
+ ++
−
− + Δ +
−
+
= +
M n n
n n
k j m i p k j m i p t D
k j i u k j i u
2 / 1 2
/ 1 1
, , ,
,1 , ,2 / 1 ,
,2 / 1
y pis defined uis defined
1. FDTDM
( ) ( )
{ }
∑
=+
Δ m Dm p i+m jk p i m jk h 0
, , ,
ρ ,1
( ) ( )
+ + = n +
n i j k v i j k
v 1 , 1/2, , 1/2, Δh
( ) ( )
{ }
∑
=+
+ + + − −
Δ
− Δ M
m
n
m pn i j m k p i j mk
h D t
0
2 / 1 2
/
1 , ,1 , ,
ρ x
( ) ( )
( ) ( )
{ }
∑
+ ++
−
− + Δ +
− Δ
+
= +
M n n
m n n
m k j i p m k j i p h D t
k j i w k
j i w
0
2 / 1 2
/ 1 1
, , 1
, ,
2 / 1 , , 2 / 1 , , ρΔhm=0
ρ
( ) ( )
( ) ( )
{ }
⎡ Δ
=
∑
− +
M n n
n n
k j i
k j i
t D
k j i p k j i
p 1/2 1/2
2 / 1 2
/ 1
, , ,
, κ
vis defined
( ) ( )
{ }
( ) ( )
{
+ + − − −}
+
⎢⎣
⎡ + + − − −
−Δ
∑
∑
=M n n
m m
n m n
k m j i v k m j i v D
k j m i u k j m i u h 0D
, 2 / 1 , , 2 / 1 ,
, , 2 / 1 ,
, 2 / 1
2
( ) ( )
{ }
⎥⎦
− ⎤
−
− + + +
∑
=
= M m
n m n
m
m k j i w m
k j i w D
0 0
2 / 1 , , 2 / 1 ,
,
Initial-value problem
課題:指向性を考慮した波動数値解析手法
課題:指向性を考慮した波動数値解析手法
No.
(0, 0) Input spatial
Distribution [0, 0] Spectral
correction filter [0, 0]
Spatial Spectral Transient
Inverse filter ORDER‐0 (1, 1)
(1 0)
[ , ] Input spatial Distribution [1, 1]
Input spatial
[ , ] Spectral correction filter [1, 1]
Spectral Inverse
filter ORDER‐1 (1, 0)
(1,‐1)
Distribution [1, 0]
Input spatial Distribution [1,‐1]
i l
correction filter [1, 0]
Spectral correction filter [1,‐1]
S l
d nt)
(2, 2) (2, 1)
Input spatial Distribution [2, 2]
Input spatial
Distribution [2 1] Spectral
correction filter [2 1]
Spectral correction filter [2, 2]
Inverse
nd fiel d ce poi n
(2, 0)
Distribution [2, 1]
Input spatial
correction filter [2, 1]
Spectral filter ORDER‐2
Sou n c
characteristics for (n,m)=(3,3)
( , ) (3,3)
z
y x
7
x-z plane
航空機の指向性を入手し、解析に取り込むことが課題
指向性を考慮した騒音伝搬予測手法に関する研究
④全体研究計画
H24(2012) H25(2013) H24(2012) H25(2013)
(各年度目標) ①波動数値解析に対
する指向性の導入手 法 確立
②フィールドにおける 検証実験
③音源指向性デ タ 法の確立 ③音源指向性データ の取得と騒音伝搬特 性に及ぼすケースス タデ
(課題と計画)
①指向性を考慮した波動数値解析手法
タディ
①指向性を考慮した波動数値解析手法 の確立と整備
②フィ ルドにおける予備実験
②フィールドにおける予備実験
②フィールドにおける検証実験
(装置等 実験準備)
③音源指向性データの取得と騒音伝搬 特性に及ぼす影響のケーススタディ
音響伝搬実験
実施日時:平成24年11月27日(火) 実施場所:大樹航空宇宙実験場
実施場所 大樹航空宇宙実験場 目的:
上空(100 500 )からの音の伝搬を実験
• 上空(100~500m)からの音の伝搬を実験 的に詳細に計測した参照データを取得すること。
[DREAMS projectの実験 (11/18-12/2) に スポット参加]
スポット参加]
• 指向性をもつ音源のデータを取得すること。
実験概要(1)
指向性音源の作成
指向性音源からの伝搬を計測するため、係留気球 に吊るす音源を作成。 GPS (位置、方位を確認)
に吊るす音源を作成。
音源 ホ ンスピ カTH 750(TOA) 音源:ホーンスピーカTH‐750(TOA)
音源の指向性-事前に無響室で計測
-5 0 0.0 0
15 30 330
250Hz 345
-5 0 0.0 0
15 30 330
500Hz 345
-20.0 -15.0 -10.0 -5.0
45
60 300
315
-20.0 -15.0 -10.0 -5.0
45
60 300
315
-40.0 -35.0 -30.0 -25.0
75
90 270
285
-40.0 -35.0 -30.0 -25.0
75
90 270
285
105
255 255 105
10dB 10dB
120
135 150 210
225
240 120
135 150 210
225 240
150 180 165
195
210 150
180 165 195 210
音源の指向性-事前に無響室で計測
-5 0 0.0 0
15 30 330
1kHz 345
-5 0 0.0 0
15 30 330
2kHz 345
-20.0 -15.0 -10.0 -5.0
45
60 300
315
-20.0 -15.0 -10.0 -5.0
45
60 300
315
-40.0 -35.0 -30.0 -25.0
75
90 270
285
-40.0 -35.0 -30.0 -25.0
75
90 270
285
105
255 255 105
10dB 10dB
120
240 240 120
実験時間:半日 音源高さ:100 音源高さ:100 m
指向性 :3ケース(直下 真横 斜め下)
指向性 :3ケ ス(直下,真横,斜め下)
実際の音源
所望の向きに設定することができなかった。
◎ロープによる係留 ◎風条件
次年度の課題
◎ロープによる係留,◎風条件,
◎装置の重心バランス
実験結果
100 150 200
高度(路面基準)
O.A. TSP-50 s Band TSP-16 s
‐50 0 50
13:51:4 13:52:0 13:52:2 13:52:4 13:53:0 13:53:2 13:53:4 13:54:0 13:54:2 13:54:4 13:55:0 13:55:2 13:55:4 13:56:0 13:56:2 13:56:4 13:57:0 13:57:2 13:57:4 13:58:0 13:58:2 13:58:4 13:59:0 13:59:2 13:59:4 14:00:0 14:00:2 14:00:4 14:01:0 14:01:2 14:01:4 14:02:0 14:02:2 14:02:4 14:03:0 14:03:2 14:03:4 14:04:0 14:04:2 14:04:4 14:05:0 14:05:2 14:05:4 14:06:0 14:06:2 14:06:4 14:07:0 14:07:2 14:07:4 14:08:0 14:08:2 14:08:4 14:09:0 14:09:214:09:314:09:514:10:0 14:10:2 14:10:4 14:11:0 14:11:2 14:11:4 14:12:0 14:12:2 14:12:4 14:13:0 14:13:2 14:13:4 14:14:0 14:14:2 14:15:0 14:15:2 14:15:4 14:16:0 14:16:2 14:16:4
200
‐150
‐100
40 00 20 40 00 20 40 00 20 40 00 20 40 00 20 40 00 20 40 00 20 40 00 20 40 00 20 40 00 20 40 00 20 40 00 20 40 00 20 40 00 20 40 00 20 40 00 20 40 00 20 40 00 20 30 50 00 20 40 00 20 40 00 20 40 00 20 40 00 20 00 20 40 00 20 40
‐200
250 Hz 500 Hz 1 kHz 2 kHz 250 Hz 500 Hz
伝搬測定 TSP(Time Stretched Pulse, Swept sine) 法
伝搬測定 ( , p ) 法
例:4s 50s, 4Oct.
実測時間:13:50-14:15(天候の変化により終了)
実測時間:13:50 14:15(天候の変化により終了)
音源高度:50m (GPSデータによる)
実験結果 気象条件・音源条件
420.0
0
270.0 300.0 330.0 360.0 390.0
風向@h=1.2m
n [deg.]
N 0
120.0 150.0 180.0 210.0
240.0 @
風向@h=10m 風向@h=5m
nd directio
270 90
15.0
風速@h 1 2 30.0 0.0
60.0
90.0 OA‐No.1
Win
180
10.0 9 0 11.0 12.0 13.0
14.0 風速@h=1.2m
風速@h=5m
OA‐No.1
[m/s]
7.0 8.0 9.0
speed
70 80 90 100
70 80 90 100
70 80 90 100
70 80 90 100
係留高度
20 30 40 50 60
20 30 40 50 60
20 30 40 50 60
20 30 40 50 60
高度はほぼ 定(50 )
360 360360 360 0 10 20 0 10 20 10 20 0 10
20
高度はほぼ一定(50m)
180 225 270 315 360
180 225 270 315 360
180 225 270 315 360
180 225 270 315 360
OA‐001 OA‐003
OA‐004
SP方位
滑走路西方向
滑走路東方向 450
90 135 180
0 45 90 135 180
0 45 90 135 180
0 45 90 135
180 OA 004
滑走路西方向
滑 路東
‐90
‐45
‐90
‐45
‐90
‐45
‐90
‐45 OA‐002
角度はバラバラ(一測定の間にも大きく変化)
角度はバラバラ( 測定の間にも大きく変化)
風速大→スピーカが大きく振られ、動揺した。
実験結果 定量化 定量化
TSP I.R. Filtering Integration
t1
ITSP ∫1 p ( ) t
2d t
t0
Unit impulse 500 Hz Unit impulse 1000 Hz Unit impulse 1000 Hz Unit impulse 2000 Hz
33 2 5 4
6
7 9 8
10 7
12 11
実験結果
1st 2nd 3rd 4th 5th
6th 7th 8th 9th 10th
130
6th 7th 8th 9th 10th
11th 12th 13th ave 10th
11th 12th 13th ave
115 120
125
250Hz
105 110 115
L[dB]
SPL 100
1st 2nd 3rd 4th 5th
6th 7th 8th 9th 10th
130
6th 7th 8th 9th 10th
11th 12th 13th ave 10th
11th 12th 13th ave
115 120
125
500Hz
105 110 115
L[dB]
90 95
SPL 100
80 85 90
500 300 100 100 300 500
‐500 ‐300 ‐100 100 300 500
Position [m]
実験結果
1st 2nd 3rd 4th 5th
6th 7th 8th 9th 10th
130
6th 7th 8th 9th 10th
11th 12th 13th ave 10th
11th 12th 13th ave
115 120
125
1kHz
105 110 115
L[dB]
90 95
SPL 100
80 85 90
500 300 100 100 300 500
‐500 ‐300 ‐100 100 300 500
Position [m]