matlab 离散无记忆信源建模

单符号离散无记忆二元信源

%离散无记忆信源的MATLAB建模
%假设有X个球，红球占30%，白球为70%，抓取后再放回，为每次抓球的结果做建模
P0=0.3;%每次抓取红球的概率
P1=0.7;%每次抓取白球的概率
N=10;%抓球的次数
x=randsrc(1, N, [0  1; P0  P1]);%生成一个 1行N列 的矩阵， 其中的元素只有0（抓取出红球）和1（抓取出白球），元素0出现的概率为P0，元素1出现的概率为P1
disp(x);
N0=length(find(x==0));%find(x==0)的含义是找到所有 x==0 的位置，将这些位置生成一个数组
disp(find(x==0));
P0x=N0/N;%N次抓取的结果中，红球出现的概率
P1x=1-P0x;%N次抓取的结果中，白球出现的概率
disp(P0x);
disp(P1x);

将 N 改为 1000 后的结果

单符号离散无记忆四元信源

输出的符号数为4， p(0)=0.4, p(1)=0.3, p(2)=0.2, p(3)=0.1，长度为N=10000

%单符号离散无记忆四元信源的MATLAB建模
%输出符号数为4， p(0)=0.4, p(1)=0.3, p(2)=0.2, p(3)=0.1，长度为10000
P0=0.4;
P1=0.3;
P2=0.2;
P3=0.1;

N=10000;
x=randsrc(1, N, [0  1  2  3; P0  P1  P2  P3]);%生成一个 1行N列 的矩阵， 其中的元素只有0,1,2,3，出现的概率分别是P0  P1  P2  P3
disp(x);

N0=length(find(x==0));%find(x==0)的含义是找到所有 x==0 的位置，将这些位置生成一个数组
disp(find(x==0));
N1=length(find(x==1));%find(x==1)的含义是找到所有 x==1 的位置，将这些位置生成一个数组
disp(find(x==1));
N2=length(find(x==2));%find(x==2)的含义是找到所有 x==2 的位置，将这些位置生成一个数组
disp(find(x==2));
N3=length(find(x==3));%find(x==3)的含义是找到所有 x==3 的位置，将这些位置生成一个数组
disp(find(x==3));

P0x=N0/N;
P1x=N1/N;
P2x=N2/N;
P3x=N3/N;
disp(P0x);
disp(P1x);
disp(P2x);
disp(P3x);