第三次作

 

由上可知，

     p(a₁)=0.2 ,p(a₂)=0.3  ,p(a₃)=0.5

 因为X(a_i)=i,       X(a₁)=1，X(a₂)=2，X(a₃)=3

  F_X(0)=0，F_X(1)=0.2 ,F_X(2)=0.5  ,F_X(3)=1.0, U⁽⁰⁾=1 ,L⁽⁰⁾=0

  有公式，L⁽ⁿ⁾=L^(n-1)+(U^(n-1)-L^(n-1))F_x(x_n-1)

             u⁽ⁿ⁾=L^(n-1)+(U^(n-1)-L^(n-1))F_x(x_n)

 第一次 出现a1，时

          L⁽¹⁾=L⁽⁰⁾⁺(U⁽⁰⁾-L⁽⁰⁾)F_x(0)=0

           U⁽¹⁾=L⁽⁰⁾⁺(U⁽⁰⁾-L⁽⁰⁾)F_x(1)=0.2

第二次出现a1，时

          L⁽²⁾=L⁽¹⁾⁺(U⁽¹⁾-L⁽¹⁾)F_x(0)=0

           U⁽²⁾=L⁽¹⁾⁺(U⁽¹⁾-L⁽¹⁾)F_x(1)=0.04

第三次 出现a3，时

          L⁽³⁾=L⁽²⁾⁺(U⁽²⁾-L⁽²⁾)F_x(2)=0.02

           U⁽³⁾=L⁽²⁾⁺(U⁽²⁾-L⁽²⁾)F_x(3)=0.04

第四次 出现a2，时

          L⁽⁴⁾=L⁽³⁾⁺(U⁽³⁾-L⁽³⁾)F_x(1)=0.024

           U⁽⁴⁾=L⁽³⁾⁺(U⁽³⁾-L⁽³⁾)F_x(2)=0.03

第五次 出现a3，时

          L⁽⁵⁾=L⁽⁴⁾⁺(U⁽⁴⁾-L⁽⁴⁾)F_x(2)=0.027

           U⁽⁵⁾=L⁽⁴⁾⁺(U⁽⁴⁾-L⁽⁴⁾)F_x(3)=0.03

第六次 出现a1，时

          L⁽⁶⁾=L⁽⁵⁾⁺(U⁽⁵⁾-L⁽⁵⁾)F_x(0)=0.027

           U⁽⁶⁾=L⁽⁵⁾⁺(U⁽⁵⁾-L⁽⁵⁾)F_x(1)=0.0276

所以，序列a1a1a3a2a3a1的实值标签为：T(113231)=(L⁽⁶⁾⁺ U⁽⁶⁾)/2=0.0273;

 

解.源程序为：

#include<stdio.h>


int main()
{
    int m[10];
    double A0=0,A1=0.2,A2=0.5,A3=1.0;
    double l[10]={0.0},u[10]={1.0};//l和u分别表示下界和上界
    double tag=0.6321599,t;
    int n=10;
    for(int i=1;i<=10;i++)
    {
        t=(tag-l[i-1])/(u[i-1]-l[i-1]);
        if(t>A0&&t<A1)
        {
            l[i]=l[i-1]+(u[i-1]-l[i-1])*A0;
            u[i]=l[i-1]+(u[i-1]-l[i-1])*A1;
            m[i]=1;
        }
        else if(t>A1&&t<A2)
        {
            l[i]=l[i-1]+(u[i-1]-l[i-1])*A1;
            u[i]=l[i-1]+(u[i-1]-l[i-1])*A2;
            m[i]=2;
        }
         else
        {
            l[i]=l[i-1]+(u[i-1]-l[i-1])*A2;
            u[i]=l[i-1]+(u[i-1]-l[i-1])*A3;
            m[i]=3;
        }
         printf("%d\t",m[i]);
    }

        return 0;        
}