实验一_离散信源的编程实现.docVIP

实验报告 一、实验名称 离散信源的编程实现 二、实验目的 了解不同离散信源产生消息的原理，利用程序进行实现。 三、实验原理 在matlab中，要随机产生给定概率的A、B、C、D、E五个离散消息，可以利用随机数函数unidrnd（），它每次随机产生一个特定范围内的数，例如：r unidrnd（100），则每执行一次，r的值随机为1到100之间的某一个整数。 如果要以给定的概率产生消息，我们可以通过限制r的范围，来对应输出消息即可。例如：A、B、C、D、E产生的概率相等，均为0.2，可以通过编程实现当0 r 20时，输出A，当20 r 40时，输出B, 当40 r 60时，输出C, 当60 r 80时，输出D, 当80 r 100时，输出E。 当要产生其他给定概率的消息时，对应改变r的范围即可。 四、实验内容 1．实现产生A、B、C、D、E五个消息的信源，各消息的产生概率相等，均为0.2，试连续产生该信源的100个消息，根据产生的消息序列统计每个消息的相应后验概率。 流程图： 源程序： clear clc A 0;B 0;C 0;D 0;E 0; for i 1:100 r unidrnd 100 ; %产生一个1到100之间的随机数 if r 80 message E;E 0.01+E; %80 r 100时，message为字符E,E增加0.01 elseif r 60 message D; D 0.01+D; %60 r 80时，message为字符D,D增加0.01 elseif r 40 message C;C 0.01+C; %40 r 60时，message为字符C,C增加0.01 elseif r 20 message B;B 0.01+B; %20 r 40时，message为字符B,B增加0.01 elseif r 0 message A;A 0.01+A; %0 r 20时，message为字符A,A增加0.01 end output i message ; %将message的内容赋给output i end fprintf 产生的100个离散消息为： ; %输出 fprintf \n ; %换行 for i 1:100 %循环 fprintf %c,output i ; %输出output i 的内容（即输出消息） end %循环结束 fprintf \n ; %换行 fprintf \n ; fprintf A后验概率为: %f\n,A ; %输出A的后验概率 fprintf B后验概率为: %f\n,B ; %输出B的后验概率 fprintf C后验概率为: %f\n,C ; %输出C的后验概率 fprintf D后验概率为: %f\n,D ; %输出D的后验概率 fprintf E后验概率为: %f\n,E ; %输出E的后验概率 运行结果： 2．实现产生A、B、C、D、E五个消息的信源，各消息的产生概率分别为0.4,0.1,0.2,0.2,0.1，试连续产生该信源的100个消息，根据产生的消息序列统计每个消息的相应后验概率。 流程图： 源程序： clear clc A 0;B 0;C 0;D 0;E 0; %定义变量 for i 1:100 %循环开始 r unidrnd 100 ; %产生一个1到100之间的随机数 if r 90 %判断语句 message E;E 0.01+E; %80 r 100时，message为字符E,E增加0.01 elseif r 70 message D; D 0.01+D; %60 r 80时，message为字符D,D增加0.01 elseif r 50 message C;C 0.01+C; %40 r 60时，message为字符C,C增加0.01 elseif r 40 message B;B 0.01+B; %20 r 40时，message为字符B,B增加0.01 elseif r 0 message A;A 0.01+A; %0 r 20时，message为字符A,A增加0.01 end %if语句结束 output i message; %将message的内容赋给output i end %循环结束 fprintf 产生的100个离散消息为： ; %输出 fprintf \n ; %换行 for i 1:100 %循环 fprintf %c,output i ; %输出output i 的内容（即输出消息） end %循环结束 fprintf \n ; %换行 fprintf \n ; fprintf A后验概率为: %f\n,A ; %输出A的后验概率 fprint