以Gaussian p.d.f.求方差的一元一次方程为例,已知均值(MU) / 样本点(x,y),求方差(sigma2)。
实参赋值
MU= 4.5e-8; % 均值 x_thres= 5.1e-4; % 找一个好的门限->对于solve(sigma2)->很重要。 y_thres= mvnpdf(x_thres, noise_mu, noise_sigma);
建立方程
syms x y sigma2 MU eqn= ( 1/sqrt(2*pi*sigma2)* exp( -(x-MU)^2/(2*sigma2) )==y ); % 建立方程:Gaussian p.d.f. sol_sigma2 = solve(eqn, sigma2); % 得到2个解:solx(1), solx(2)
将 MU/ (x,y) 的实参带入方程,得到sigma2的数值解
Sigma_1= subs( sol_sigma2(1), {MU, x, y}, {mu_val, x_thres, y_thres}); % 带入实参 sigma_1= eval(Sigma_1); % 得到方差的数值解
感谢 https://www.cnblogs.com/sunny99/ sumoier对本文的帮助!
标签:方差,MU,solve,matlab,解方程,实参,thres,sigma2 From: https://www.cnblogs.com/Yi-920259020/p/16869735.html