這是從書上copy過來的
描寫有關信噪比(SNR)的程式,程式如下
function [Y]=quan(X,N)
sgn=sign(X);
if (sgn==-1)
X=-X;
else
X=X;
end
SP=0;
SX=X-fix(X);
for i=1:1:N
SX=SX*2;
if (SX>=1)
SP=SP+2^(-i);
SX=SX-1;
C=1;
else
SP=SP;
C=0;
end
end
X=fix(X)+SP+C*2^(-N); %do round
Y=sgn*X;