MATLAB实现PCM13折线编码，有过经验的请进。。
%PCM编码程序function code=pcm(S) %S为输入信号
z=sign(S); %判断S的正负
MaxS=max(abs(S)); %求S的最大值
S=abs(S/MaxS); %归一化
Q=2048*S; %量化
code=zeros(length(S),8); %代码存储矩阵（全零）
%段落码判断程序
for i=1:length(S)
if (Q(i)>=128)&(Q(i)<=2048)
code(i,2)=1; %在第五段与第八段之间，段位码第一位都为“1”
end
if (Q(i)>32)&(Q(i)<128)||(Q(i)>=512)&(Q(i)<=2048)
code(i,3)=1; %在第三四七八段内，段位码第二位为“1”
end
if (Q(i)>=16)&(Q(i)<32)||(Q(i)>=64)&(Q(i)<128)||(Q(i)>=256)&(Q(i)<512)||(Q(i)>=1024)&(Q(i)<=2048)
code(i,4)=1; %在二四六八段内，段位码第三位为“1”
end
end
%段内码判断程序
N=zeros(length(S));
for i=1:length(S)
N(i)=bin2dec(num2str(code(i,2:4)))+1; %找到code位于第几段
a=; %量化间隔
b=; %除以16，得到每段的最小量化间隔
for i=1:length(S)
q=ceil((Q(i)-a(N(i)))/b(N(i))); %求出在段内的位置
if q==0
code(i,(5:8))=; %如果输入为零则输出“0”
else k=num2str(dec2bin(q-1,4)); %编码段内码为二进制
code(i,5)=str2num(k(1));
code(i,6)=str2num(k(2));
code(i,7)=str2num(k(3));
code(i,8)=str2num(k(4));
end
if z(i)>0
code(i,1)=1;
elseif z(i)<0
code(i,1)=0;
end %符号位的判断
end
end
我运行时输入:
k=90
a=pcm(k)
可以出现结果，但一旦输入多值，例如：
k=[-89 78 90 -78 24 76 -45 43]
a=pcm(k)
就提示错误：
??? Attempted to access a(0); index must be a positive integer or logical.
Error in ==> pcm at 27
q=ceil((Q(i)-a(N(i)))/b(N(i))); %求出在段内的位置
请教问题出现在哪？？ 没人看吗？自己顶下。。。 Ref: 2.常见的程序出错问题整理 http://forum.vibunion.com/thread-46001-1-1.html
2F->Subscript indices must either be real positive integers or logicals
From http://forum.vibunion.com/home-space-uid-63979-do-blog-id-18250.html
页:
[1]