variance.m

u率pcm编码,还有误差分析,通信原理实验的一个例子

function v=variance(funfcn,a,b,tol,p1,p2,p3)
%EXPECTED求出概率密度函数为''funfcn,定义域为[a,b]的随机变量的方差.
%Y=VARIANCE(F',A,B,TOL.P1,P2,P3),求出在m_file中给出的密度函数F的方差
%该函数最多可包含三个参数:p1,p2,p3
%tol=相对误差

args=[];
for n=1:nargin-4
    args=[args,',p',int2str(n)];
end
args=[args,')'];
funfcn1=['x',funfcn];
e=eval(['quad(funfcn1,a,b,tol,[]',args]);
% newfun=['(x-(',num2str(e),')).^2.*',funfcn];%此行程序错误
newfun='varnormal';
v=eval(['quad(newfun,a,b,tol,[]',args]);%对varnomal求积分,也即得到给定区域的方差