联系QQ:

2181264433

新闻资讯
联系我们
联系:张女士
微信:扫一扫右侧二维码
QQ :2181264433
邮箱:2181264433@qq.com
地址:浙江省-嘉兴市-南湖区
网址:www.mhslogic.com
您当前位置:首页 > 国内资讯 > 正文国内资讯
MATLAB代做|FPGA代做|simulink代做——高斯混合模型参数估计
添加时间:2020-7-5 来源:本站整理
% 高斯混合模型参数估计:EM算法
 
function [omega mu nu  t  delta]=sunMOG(x, J, maxit, pre)
% function [omega mu nu  t  delta]=sunMOG(x, J, maxit, pre)
% x为一维信号, J为高斯函数个数,maxit为最大迭代次数,pre为精度
% p(x)=g(x,mu(1),nu(1))*w(1) + ... + g(x,mu(J),nu(J))*w(J)
% g为高斯函数,[w(1) ... w(J)]为线性组合系数
% [mu(1)... mu(J)]为各个高斯函数的均值,[nu(1)..nu(J)]为各个高斯函数的均值sigma^2
if nargin<4
    pre=1e-3;
end
if nargin<3
    maxit=388;
end
if nargin<2
    J=3;
end

N=size(x,2);
omega=rand(1,J);
omega=omega/sum(omega);
mu=rand(1,J);
nu=rand(1,J);
delta=1;
for t=1:maxit
    oldomega=omega;
    oldmu=mu;
    oldnu=nu;
    % calculate p_x
    p_x=zeros(1,N);
    for i=1:N
        p_x(i)=px(x(i),oldomega,oldmu,oldnu);
    end
    % calculate omegam, mu, nu
    omega=zeros(1,J);
    mu=zeros(1,J);
    nu=zeros(1,J);
    for j=1:J
        for n=1:N
            omega(j)=omega(j)+oldomega(j)*g(x(n),oldmu(j),oldnu(j))/p_x(n);
            mu(j)=mu(j)+oldomega(j)*g(x(n),oldmu(j),oldnu(j))*x(n)/p_x(n);
        end
    end
    mu=mu./omega;
    for j=1:J
        for n=1:N
            nu(j)=nu(j)+oldomega(j)*g(x(n),oldmu(j),oldnu(j))*(x(n)-mu(j))^2/p_x(n);
        end
    end
    nu=nu./omega;
    omega=omega/N;
    delta=norm(mu-oldmu)/norm(oldmu)+norm(nu-oldnu)/norm(oldnu);
    if delta<pre        
        break;        
    end    
end
[t delta]
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 高斯函数
function p=g(x,mu,nu)
    p=exp(-(x-mu)^2/nu/2) / sqrt(2*pi*nu);
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%一维信号x的pdf
function p=px(x,omega,mu,nu)
    J=size(omega,2);
    p=0;
    for j=1:J
        p=p + g(x,mu(j),nu(j))*omega(j);
    end
end

联系:张女士

QQ :2181264433

微信:lovemike121

邮箱:2181264433@qq.com

网站:http://www.mhslogic.com/

------------------------------------------------------------------------------------------------

交易可直接通过淘宝店铺


(可直接通过旺旺与技术员联系)




https://shop185632181.taobao.com/?spm=2013.1.w5001-22653403434.6.27f92e0a6xR1sB&scene=taobao_shop

扫一扫,关注我们