pf_stats.m

来自「％ ％ set some variables in the workspace」· M 代码 · 共 17 行

function [Mean, Var] = pf_stats(pf)

	z = [pf.x.s]';
	prob = [pf.x.pi]';

	mn = sum( z .* [prob prob] );
	s2 = sum( (z - ones(pf.n,1)*mn)'.^2)';
	var = sum( s2 .* prob );

	if nargout == 0,
		fprintf('mean = %f %f, var = %f\n', mn, var);
	else
		Mean = mn;
		Var = var;
	end