混沌系工学特論 #5

Similar documents
Microsoft PowerPoint - 7.pptx

スライド 1

<4D F736F F F696E74202D2091E6824F82518FCD E838B C68CEB82E894AD90B B2E >

<4D F736F F F696E74202D2091E6824F82538FCD8CEB82E88C9F8F6F814592F990B382CC8CB4979D82BB82CC82505F D E95848D8682CC90B69

Information Theory

<4D F736F F F696E74202D2091E6824F82568FCD8CEB82E892F990B382CC8CF889CA82BB82CC82515F B834E838A B9797A3959C8D F A282E982C682AB82CC8CEB82E897A62E >

aaa

Microsoft PowerPoint - 14MDL.pptx

ii 3.,. 4. F. (), ,,. 8.,. 1. (75% ) (25% ) =9 7, =9 8 (. ). 1.,, (). 3.,. 1. ( ).,.,.,.,.,. ( ) (1 2 )., ( ), 0. 2., 1., 0,.

Microsoft PowerPoint - mp11-06.pptx

Probit , Mixed logit

Information Theory

: (a) ( ) A (b) B ( ) A B 11.: (a) x,y (b) r,θ (c) A (x) V A B (x + dx) ( ) ( 11.(a)) dv dt = 0 (11.6) r= θ =

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

untitled

生命情報学

Microsoft Word - 微分入門.doc

資料2-3 要求条件案.doc

Chap2

2-1 / 語問題 項書換え系 4.0. 準備 (3.1. 項 代入 等価性 ) 定義 3.1.1: - シグネチャ (signature): 関数記号の集合 (Σ と書く ) - それぞれの関数記号は アリティ (arity) と呼ばれる自然数が定められている - Σ (n) : アリ

数値計算法

Microsoft PowerPoint - SPECTPETの原理2012.ppt [互換モード]

統計的データ解析

Microsoft Word - 補論3.2

() ): (1) f(x) g(x) x = x 0 f(x) + g(x) x = x 0 lim f(x) = f(x 0 ), lim g(x) = g(x 0 ) x x 0 x x0 lim {f(x) + g(x)} = f(x 0 ) + g(x 0 ) x x0 lim x x 0

2016.

第121回関東連合産科婦人科学会総会・学術集会 プログラム・抄録

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

ii 3.,. 4. F. (), ,,. 8.,. 1. (75%) (25%) =7 20, =7 21 (. ). 1.,, (). 3.,. 1. ().,.,.,.,.,. () (12 )., (), 0. 2., 1., 0,.

Microsoft PowerPoint - 第3回2.ppt

3. 測定方法 測定系統図 測定風景写真

「産業上利用することができる発明」の審査の運用指針(案)

カイ二乗フィット検定、パラメータの誤差

講義「○○○○」

情報量と符号化

PowerPoint プレゼンテーション

Microsoft PowerPoint - mp13-07.pptx

³ÎΨÏÀ

e a b a b b a a a 1 a a 1 = a 1 a = e G G G : x ( x =, 8, 1 ) x 1,, 60 θ, ϕ ψ θ G G H H G x. n n 1 n 1 n σ = (σ 1, σ,..., σ N ) i σ i i n S n n = 1,,

春期講座 ~ 極限 1 1, 1 2, 1 3, 1 4,, 1 n, n n {a n } n a n α {a n } α {a n } α lim n an = α n a n α α {a n } {a n } {a n } 1. a n = 2 n {a n } 2, 4, 8, 16,

PowerPoint プレゼンテーション

Hara-statistics

2012-第14回.pptx

2014年度 九州大・理系数学

基礎統計

画像解析論(2) 講義内容

() x + y + y + x dy dx = 0 () dy + xy = x dx y + x y ( 5) ( s55906) 0.7. (). 5 (). ( 6) ( s6590) 0.8 m n. 0.9 n n A. ( 6) ( s6590) f A (λ) = det(a λi)

DVIOUT

オートマトン 形式言語及び演習 3. 正規表現 酒井正彦 正規表現とは 正規表現 ( 正則表現, Regular Expression) オートマトン : 言語を定義する機械正規表現 : 言語

オートマトン 形式言語及び演習 1. 有限オートマトンとは 酒井正彦 形式言語 言語とは : 文字列の集合例 : 偶数個の 1 の後に 0 を持つ列からなる集合 {0, 110, 11110,

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


Microsoft PowerPoint _SpinGlass.html

アルゴリズムとデータ構造

Microsoft PowerPoint - 13approx.pptx

そもそも情報とはなに? ある事柄に関して知識を得たり判断のより所としたりするために不可欠な 何らかの手段で伝達 ( 入手 ) された種々の事項 ( の内容 ) ( 新明解国語辞典第 6 版 三省堂より一部抜粋 ) コンピュータで取り扱う情報の定義 Claud Elwood Shannon( 情報理論

ii 3.,. 4. F. ( ), ,,. 8.,. 1. (75% ) (25% ) =7 24, =7 25, =7 26 (. ). 1.,, ( ). 3.,...,.,.,.,.,. ( ) (1 2 )., ( ), 0., 1., 0,.

4 段階推定法とは 予測に使うモデルの紹介 4 段階推定法の課題 2

1/30 平成 29 年 3 月 24 日 ( 金 ) 午前 11 時 25 分第三章フェルミ量子場 : スピノール場 ( 次元あり ) 第三章フェルミ量子場 : スピノール場 フェルミ型 ボーズ量子場のエネルギーは 第二章ボーズ量子場 : スカラー場 の (2.18) より ˆ dp 1 1 =

パソコンシミュレータの現状

Microsoft PowerPoint - H21生物計算化学2.ppt

Microsoft PowerPoint - 測量学.ppt [互換モード]

PowerPoint プレゼンテーション

Microsoft PowerPoint - 10.pptx

データ解析

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

2011年度 東京大・文系数学

τ-→K-π-π+ν τ崩壊における CP対称性の破れの探索

If(A) Vx(V) 1 最小 2 乗法で実験式のパラメータが導出できる測定で得られたデータをよく近似する式を実験式という. その利点は (M1) 多量のデータの特徴を一つの式で簡潔に表現できること. また (M2) y = f ( x ) の関係から, 任意の x のときの y が求まるので,

例 e 指数関数的に減衰する信号を h( a < + a a すると, それらのラプラス変換は, H ( ) { e } e インパルス応答が h( a < ( ただし a >, U( ) { } となるシステムにステップ信号 ( y( のラプラス変換 Y () は, Y ( ) H ( ) X (

ii 2. F. ( ), ,,. 5. G., L., D. ( ) ( ), 2005.,. 6.,,. 7.,. 8. ( ), , (20 ). 1. (75% ) (25% ). 60.,. 2. =8 5, =8 4 (. 1.) 1.,,

2015-2018年度 2次数学セレクション(整数と数列)解答解説

Microsoft Word - 02__⁄T_ŒÚ”�.doc

技術協会STD紹介

Microsoft PowerPoint - kyoto

確率分布 - 確率と計算 1 6 回に 1 回の割合で 1 の目が出るさいころがある. このさいころを 6 回投げたとき,1 度も 1 の目が出ない確率を求めよ. 5 6 /6 6 =15625/46656= (5/6) 6 = ある市の気象観測所での記録では, 毎年雨の降る

3.5 トランジスタ基本増幅回路 ベース接地基本増幅回路 C 1 C n n 2 R E p v V 2 v R E p 1 v EE 0 VCC 結合コンデンサ ベース接地基本増幅回路 V EE =0, V CC =0として交流分の回路 (C 1, C 2 により短絡 ) トランジスタ

工業数学F2-04(ウェブ用).pptx

Microsoft Word - 8章(CI).doc

日心TWS

森林水文 水資源学 2 2. 水文統計 豪雨があった時, 新聞やテレビのニュースで 50 年に一度の大雨だった などと報告されることがある. 今争点となっている川辺川ダムは,80 年に 1 回の洪水を想定して治水計画が立てられている. 畑地かんがいでは,10 年に 1 回の渇水を対象として計画が立て

(1) (2) (1) (2) 2 3 {a n } a 2 + a 4 + a a n S n S n = n = S n

Microsoft PowerPoint - DA2_2019.pptx

2010 II / y = e x y = log x = log e x 2. ( e x ) = e x 3. ( ) log x = 1 x 1.2 Warming Up 1 u = log a M a u = M a 0

newmain.dvi

第6章 実験モード解析

スライド 1

SAP11_03

arduino プログラミング課題集 ( Ver /06/01 ) arduino と各種ボードを組み合わせ 制御するためのプログラミングを学 ぼう! 1 入出力ポートの設定と利用方法 (1) 制御( コントロール ) する とは 外部装置( ペリフェラル ) が必要とする信号をマイ

( 最初の等号は,N =0, 番目は,j= のとき j =0 による ) j>r のときは p =0 から和の上限は r で十分 定義 命題 3 ⑵ 実数 ( 0) に対して, ⑴ =[] []=( 0 または ) =[6]+[] [4] [3] [] =( 0 または ) 実数 に対して, π()

布に従う しかし サイコロが均質でなく偏っていて の出る確率がひとつひとつ異なっているならば 二項分布でなくなる そこで このような場合に の出る確率が同じであるサイコロをもっている対象者をひとつのグループにまとめてしまえば このグループの中では回数分布は二項分布になる 全グループの合計の分布を求め

コロイド化学と界面化学

(1) 3 A B E e AE = e AB OE = OA + e AB = (1 35 e ) e OE z 1 1 e E xy e = 0 e = 5 OE = ( 2 0 0) E ( 2 0 0) (2) 3 E P Q k EQ = k EP E y 0

1999年度 センター試験・数学ⅡB

医系の統計入門第 2 版 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. このサンプルページの内容は, 第 2 版 1 刷発行時のものです.

Microsoft PowerPoint - 小路田俊子 [互換モード]


2001 年度 『数学基礎 IV』 講義録

Microsoft Word ã‡»ã…«ã‡ªã…¼ã…‹ã…žã…‹ã…³ã†¨åłºæœ›å•¤(佒芤喋çfl�)

18 ( ) I II III A B C(100 ) 1, 2, 3, 5 I II A B (100 ) 1, 2, 3 I II A B (80 ) 6 8 I II III A B C(80 ) 1 n (1 + x) n (1) n C 1 + n C

Transcription:

混沌系工学特論 #5 情報科学研究科井上純一 URL : htt://chaosweb.comlex.eng.hokudai.ac.j/~j_inoue/ Mirror : htt://www5.u.so-net.ne.j/j_inoue/index.html 平成 17 年 11 月 14 日第 5 回講義

デジタルデータの転送と復元再考 P ({ σ} ) = ex σ ( σσ ) < ij> i j ( σ ) iσ j ex < ij> 画像の復元では画像の特質を事前分布で表現できる 劣化されたデジタル画像 デジタル画像 確率モデルで表現 伝送路 { τ} 符号 ビット列ノイズの乗ったビット列復元されたデジタル画像 1111111 復号 111111111 { ξ} ビット列の場合は? { σ} 復号 11111111 復元されたビット列

雑音のある状況でのデータ送信誤り訂正符号 誤り率 の 元対称通信路 ビット反転率ただ1 度だけでなく 何回か同じ記号を繰り返し送信する 受信者側は多数決で送信された記号を確定する 1 1-1- 1 元対称通信路 5 回送信した場合 : 11 1111 1 等 送信回数 n = 5 のときの誤り確率は f ( ) = C (1 ) + C (1 ) + (5) 3 4 5 e 5 3 5 4 = 1 (1 ) + 5 (1 ) + 3 4 5 3 ビットの誤りが 5 ビットの何番目にくるかという場合の数 繰り返し回数 n を大きくするとどうなるか?

伝送速度と誤り確率 R log M 1 = = n n 伝送速度 ( レート ) 送信記号の種類の数 1, ならば ビット誤り率 繰り返し回数を限りなく多くとると 誤り確率は限りなくゼロに近づくが 同時に伝送速度もゼロになってしまう ( n) lim fe = n 実用的ではない しかし R < C 通信路容量 ならば 限りなく小さな誤り確率での情報伝送が可能 通信路での反転確率

通信路容量 C = max I( X; Y) P X I( X; Y) = H( Y) H( Y X ) X : 入力 Y : 出力 通信路容量の定義 相互情報量 : 入力 X を知ったときに 出力 Y に対して得られる知識の増加分 入力の確率分布に関して 相互情報量を最大化したものが通信路容量 通信路容量は 通信路が伝送できる最大の情報量という意味を持つ

通信路容量の計算例 誤りの無い 元対称通信路 1- P ( ) = P (1 1) = 1 YX YX P ( 1) = P (1 ) = YX YX 1 1-1 条件付き確率での表現 グラフ表現 HY ( X) = log (1 ) log(1 ) { } C = I( X; Y) = max H( Y) H( Y X) = 1 h( ) P X = 1+ log + (1 )log(1 ) 入力分布によらない BSC の通信路容量 同様にガウス通信路に対して 通信路容量は C 1 = max I( X; Y) = log 1+ a a

通信路符号化定理 通信路符号化定理 i) R< Cなる任意の速度 R= log M n に対し ii) 任意に小さな復号誤り率 R> CなるRに対し 任意に小さな E の符号が存在する E の符号は存在しない ゼロでなくてもよい しかし 符号長 n を無限大ととることは必要 情報源の記号数 M = nr この定理を ランダム符号 と呼ばれる符号に対して証明していく

ランダム符号 情報源の記号 :,,..., ( 等確率で生成される ) 符号化 : 1 S1 S S M S, S,..., S の各々にランダムに,1を n個並べた系列を割りあてる S S S S 1 3 M M 1 1 1 1 1 1 1 符号長 n 各ビットに ½ の確率で,1 を割り振る M nr = < n R 1 として議論を進める これらを 元対称通信路を介して転送する状況下で定理を証明していく

通信路符号化定理の証明 #1 一つの符号語が転送により拡大する大きさの評価 入力出力 1- 転送によりオリジナルな符号語が取りうる場合の数 1 1 1-1 ( ) nh 典型的な系列の個数 z = = 符号語は 元対称通信路を介して転送される ある符号語 S 1 = は通信路の雑音により1 11として受信される 全 nビットのなかで誤りの生じるビット数はおおよそnと見積もられる (1 11), (111 1),, (11 111) n ビットの誤り n ビットの誤り n ビットの誤り 通信路出力の典型的な系列 典型的な系列のなかの一つが現れる確率 n n(1 ) nlog n+ n(1 )log(1 ) nh( ) = (1 ) = = 値エントロピー関数

通信路符号化定理の証明 # 復号誤り率の評価 S3 復号誤りが生じるのは S 受信される可能性のある 3 nh( ) S, S,..., S の ( M 1) 個のいずれかに間違って復号されるとき M 個の各々が 実際に送信された符号語 S 以外の ( M 1) M = = = nh( ) nh( ) n( R + 1 h( )) n( R C) E n n 受信される可能性の (M-1) 個のいずれかある系列の個数が選ばれる確率 i) の証明終わ 1 ( C= 1 h( ) > R)

通信路符号化定理の証明 #3 ii) の証明 受信系列の全ての可能性 1つの符号語に対する通信路出力の典型列を収める 箱 の数 = z = n M n 箱 箱 元対称通信路により実際に得られる典型列の数 は全てこの 箱 に収まらなければならない n nh( ) n( C R) > > 1 M これはR> Cでは満たされない nh( ) ii) の証明終わり

スピンモデルを用いた符号 / 復号化 ビット列を送るのではなく 結合 = ξ ξ ij i j を送信 エネルギー関数 H = ij ξ ξσσ i j i j エネルギー関数の最小状態を選ぶことで復号する = ξ ξ ij i j ノイズが無い場合 = ξ ξ ij i j ノイズのある場合 フラストレーション無し エネルギー最小を与える配列が非自明 ( スピングラス ) フラストレーションフラストレーション ( 笑ってはいるが顔色が悪い ) ( 笑ってはいるが顔色が悪い ) ノイズにより反転したボンド

パリティと画素の同時送信 = ξ ξ ij i j 画素だけではなく 隣接画素対も同時に送信する ( 正方格子で N 個ある ) 尤度 P ({ },{ τ} { σ} ) = パリティ送信部分 β ijσiσ j+ h ij τ i iσi e [cosh ] [cosh ] NB ( β ) ( h) 画素送信部分 N 事後分布とエネルギー関数 P ({ σ} { τ} ) = e H eff Heff e { σ} 事前分布の部分 H = β σσ σσ h τσ eff ij ij i j ij i j i i j MAP 復元 MPM 復元を行う

復元結果 ビット誤り率 パリティ有り パリティ無し 原画像 % 劣化画像 パリティ項の割合 パリティ無し MPM 復号 パリティ有り MPM 復号 冗長性を付加することにより より効果的な復元が実現される

Sourlas 符号 N 個のビットの中から 個をピックアップして送信する i i = ξi ξi 1 1 オリジナル ビットの 個の積を送る 速度 R N = C N! 1 N ガウス通信路で送信 P i 通信路 ({ } { }) ( ) ( ) 1 1 N N! ξ = ex i i ξ 1 i ξi C π!! N i < < i ( a ) a = log 1+ 1 1! N log 1 1 R C シャノンの限界式 = 1 log

P 事後確率 { σ} { } ( ) 平均場理論による性能評価 = e β σ σ { σ} i1 < < i i1 i i1 i e β σ σ i1 < < i i1 i i1 i 磁化 1 H ( ) m= σ i = Dxtanh G β q 1 ( ) q = σ i = Dxtanh G G = x + β m エネルギー関数 レプリカ法から得られる状態方程式 スピングラス秩序変数 = β σ σ eff i1 < < i i 1 i i 1 i セルフ コンシステントに解く ( ) ( ( ) ) 1 ξ1 sgn σ 1 sgn ( ) i 1 1 Pb = = Dx G ビット誤り率

解析結果 #1 磁化 スピングラス秩序変数 強磁性 常磁性相転移 = 3 シグナルノイズ比 =.8 ビット誤り率 復号可能相 復号不可能相転移

解析結果 # 温度に関する最良性 ビット誤り率 = 3 =.8 西森温度 = 3 =.8 T =

解析結果 #3 磁化 低温ではスピングラス スピングラス秩序変数 高温では常磁性 = 3 =.6 シグナルノイズ比 ビット誤り率 誤り訂正不可能

解析結果 #4 の増加につれて減少 : log, = 強磁性相の存在限界 R C Sourlas 符号は漸近的にシャノン限界を達成する 次回は CDMA に移ります