第三次作业

 参考书《数据压缩导论(第4版)》  Page 100   5, 6

 解:

从概率模型可知:

Fx(k)=0, k≤0, Fx(1)=0.2, Fx(2)=0.5, Fx(3)=1, k>3.

我们可以利用公式确定标签所在的上下限。将u(0)初始化为1,将l(0)初始化为0。

该序列的第1个元素为a1,得到以下更新:

l(1) =0+(1-0)Fx(0)=0

u(1) =0+(1-0)Fx(1)=0.2 

也就是说,该标签包含在区间[0,0.2)中,该序列的第2个元素为a1,得到以下更新:

l(2) =0+(0.2-0)Fx(0)=0

u(2) =0+(0.2-0)Fx(1)=0.04

因此,序列a1a1的标签所在的区间为[0,0.04)。该序列的第3个元素为a3,得到以下更新:

l(3) =0+(0.04-0)Fx(2)=0.02

u(3) =0+(0.04-0)Fx(3)=0.04

因此,序列a1a1a3的标签所在的区间为[0.02,0.04)。该序列的第4个元素为a2,得到以下更新:

l(4) =0.02+(0.04-0.02)Fx(1)=0.024

u(4) =0.02+(0.04-0.02)Fx(2)=0.03

因此,该标签所在的区间为[0.024,0.03)。该序列的第5个元素为a3,得到以下更新:

l(5) =0.024+(0.03-0.024)Fx(2)=0.027

u(5) =0.024+(0.03-0.024)Fx(3)=0.03

因此,该标签所在的区间为[0.027,0.03)。该序列的第6个元素为a1,得到以下更新:

l(6) =0.027+(0.03-0.027)Fx(0)=0.027

u(6) =0.027+(0.03-0.027)Fx(1)=0.0276

可以生成序列a1a1a3a2a3a1的标签如下:

Tx(a1a1a3a2a3a1)=(0.027+0.0276)/2=0.0273

    

解:

由题可知:Tx=0.63215699,Fx(k)=0, k≤0, Fx(1)=0.2, Fx(2)=0.5, Fx(3)=1, k>3.

                    将u(0)=1,l(0)=0。

且l(k) =l(k-1) +(u(k-1) -l(k-1) )Fx(xk-1)

   u(k) =l(k-1) +(u(k-1) -l(k-1) )Fx(xk)

   t*=(tag-l(k-1))/(u(k-1) -l(k-1))

所以:

t*=(0.63215699-0)/(1 -0)=0.63215699

Fx(2)=0.5≤t*≤1=Fx(3)

l(1) =l(0) +(u(0) -l(0) )Fx(2)=0+(1-0)*0.5=0.5

u(1) =l(0) +(u(0) -l(0) )Fx(3)=0+(1-0)*1=1

因此,得到的第1个序列为:a3

t*=(0.63215699-0.5)/(1 -0.5)=0.26431398

Fx(1)=0.2≤t*≤0.5=Fx(2)

l(2) =l(1) +(u(1) -l(1) )Fx(1)=0.5+(1-0.5)*0.2=0.6

u(2) =l(1) +(u(1) -l(1) )Fx(2)=0.5+(1-0.5)*0.5=0.75

因此,得到的第2个序列为:a2

t*=(0.63215699-0.5)/(1 -0.5)=0.26431398

Fx(1)=0.2≤t*≤0.5=Fx(2)

l(3) =l(2) +(u(2) -l(2) )Fx(1)=0.6+(0.75-0.6)*0.2=0.63

u(3) =l(2) +(u(2) -l(2) )Fx(2)=0.6+(0.75-0.6)*0.5=0.675

因此,得到的第3个序列为:a2

t*=(0.63215699-0.63)/(0.675 -0.63)=0.04793311

Fx(0)=0≤t*≤0.2=Fx(1)

l(4) =l(3) +(u(3) -l(3) )Fx(0)=0.63+(0.675-0.63)*0=0.63

u(4) =l(3) +(u(3) -l(3) )Fx(1)=0.63+(0.675-0.63)*0.2=0.639

因此,得到的第4个序列为:a1

t*=(0.63215699-0.63)/(0.639 -0.63)=0.23966556

Fx(1)=0.2≤t*≤0.5=Fx(2)

l(5) =l(4) +(u(4) -l(4) )Fx(1)=0.63+(0.639-0.63)*0.2=0.6318

u(5) =l(4) +(u(4) -l(4) )Fx(2)=0.63+(0.639-0.63)*0.5=0.6345

因此,得到的第5个序列为:a2

t*=(0.63215699-0.6318)/(0.6345 -0.6318)=0.13221852

Fx(0)=0≤t*≤0.2=Fx(1)

l(6) =l(5) +(u(5) -l(5) )Fx(0)=0.6318+(0.6345-0.6318)*0=0.6318

u(6) =l(5) +(u(5) -l(5) )Fx(1)=0.6318+(0.6345-0.6318)*0.2=0.63234

因此,得到的第6个序列为:a1

t*=(0.63215699-0.6318)/(0.63234 -0.6318)=0.66109259

Fx(2)=0.5≤t*≤1=Fx(3)

l(7) =l(6) +(u(6) -l(6) )Fx(2)=0.6318+(0.63234-0.6318)*0.5=0.63207

u(7) =l(6) +(u(6) -l(6) )Fx(3)=0.6318+(0.63234-0.6318)*1=0.63234

因此,得到的第7个序列为:a3

t*=(0.63215699-0.63207)/(0.63234 -0.63207)=0.66109259

Fx(1)=0.2≤t*≤0.5=Fx(2)

l(8) =l(7) +(u(7) -l(7) )Fx(1)=0.63207+(0.63234-0.63207)*0.2=0.632124

u(8) =l(7) +(u(7) -l(7) )Fx(2)=0.63207+(0.63234-0.63207)*0.5=0.632205

因此,得到的第8个序列为:a2

t*=(0.63215699-0.63207)/(0.63234 -0.63207)=0.66109259

Fx(1)=0.2≤t*≤0.5=Fx(2)

l(9) =l(8) +(u(8) -l(8) )Fx(1)=0.632124+(0.632205-0.632124)*0.2=0.6321402

u(9) =l(8) +(u(8) -l(8) )Fx(2)=0.632124+(0.632205-0.632124)*0.5=0.6321645

因此,得到的第9个序列为:a2

t*=(0.63215699-0.6321402)/(0.6321645 -0.6321402)=0.69835391

Fx(2)=0.5≤t*≤1=Fx(3)

因此,得到的第10个序列为:a3

所以,序列为a3a2a2a1a2a1a3a2a2a3

原文地址:https://www.cnblogs.com/wushuangfei/p/4812491.html