function[y,dist]=mse_dist(funfcn,a,tol,p1,p2,p3)
%MSE_DIST 返回均方量化誤差.
%  [Y,DIST]=mse_dist(funfcn,a,tol,p1,p2,p3)
%  funcfcn=在一個m_file中給定的分布函數(shù).
%它最多由三個參數(shù)決定:p1,p2,p3.
%a=定義了量化區(qū)域 的邊界的向量.
%(注意:[a(1),a(length(a))]是funfcn的定義域)
%p1,p2,p3,是funfcn的參數(shù).
%tol=相對誤差.
 
args=[];
for n=1:nargin-3
     args=[args,',p',int2str(n)];
 end
 args=[args,')'];%args=[ ,p1,p2)]
 for i=1:length(a)-1%對給定的邊界求量化級數(shù)(質(zhì)心)
     y(i)=eval(['centroid(funfcn,a(i),a(i+1),tol',args]);
     %centroid(funfcn,a(i),a(i+1),tol ,p,funfcn,a,tol,)%此行錯誤
 end
 dist=0; 
 for i=1:length(a)-1%求量化信噪比
 newfun='varnormal';                                                                          
name=['varnormal(x,',num2str(p1),',',num2str(p2),',',num2str(y(i)),')'];
dist=dist+quad(name,a(i),a(i+1),tol);                                                                           
 end