信号处理多媒体通信 matlab实验算术编码/译码

disp('程序限定字符为：a b c d e'); str=input('请输入编码的字符串：'); [j,k]=size(str); l=0; r=1; d=1; p=[0.2 0.3 0.1 0.15 0.25]; pa=0.2; pb=0.3; pc=0.1; pd=0.15; pe=0.25; n=length(str); disp('a b c d e') disp(num2str(p)) for i=1:k if i==1 switch str(i) case 'a' m=1;a1=0; a2=pa; case 'b' m=2; a1=pa; a2=pa+pb; case 'c' m=3; a1=pa+pb; a2=pa+pb+pc; case 'd' m=4; a1=pa+pb+pc; a2=pa+pb+pc+pd; case 'e' m=5; a1=pa+pb+pc+pd; a2=pa+pb+pc+pd+pe; end l=a2-a1; end if (i>=2)&(i<=k) switch str(i) case 'a' aa=a1; ab=a1+l*pa; a1=aa;a2=ab; case 'b' aa=a1+l*pa; ab=a1+l*(pa+pb); a1=aa; a2=ab; case 'c' aa=a1+l*(pa+pb); ab=a1+l*(pa+pb+pc); a1=aa; a2=ab; case 'd' aa=a1+l*(pa+pb+pc); ab=a1+l*(pa+pb+pc+pd); a1=aa; a2=ab; case 'e' aa=a1+l*(pa+pb+pc+pd); ab=a1+l*(pa+pb+pc+pc+pd+pe); a1=aa; a2=ab; end l=a2-a1; end strl=strcat('input',int2str(i),'编码区间起始值为：'); disp(strl); disp(a1) strl=strcat('input',int2str(i),'编码区间的终止值为：'); disp(strl); disp(a2) strl=strcat('input',int2str(i),'本程序选择区间中点作为编码是：'); disp(strl); disp((a1+a2)/2) end disp('是否要译码，请选择') disp('1、译码') disp('2、不译码？') disp('说明：输入相应的数字进行选择') YN=input('您选择的是：'); if YN==1 bm=(a1+a2)/2; SuansuJiema(bm,k) end if YN==2 disp('你选择了不译码。') end