mula_sqnr.m

通原实验课的作业

close all
clear all
echo on
mu=255;
m=randn(1,500);
n=10;%
[f,maximum]=mulaw(m,mu);%
[sqnr,y_q,code]=u_pcm(f,n);%
m_quan=invmulaw(y_q,mu);%
m_quan=maximum*m_quan;%
sqnr=20*log10(norm(m)/norm(m-m_quan))%