% 目标函数:高斯混合模N[x(i),mu,sm]=0.8N[x(i),mu1,sm1]+0.2N[x(i),mu2,sm2]

%                                           mu1=2,sm1=1;    mu2=1,sm2=1

clc

clear all

N=30;                               %变量x数目为30个

for i=1:N                            %构造目标函数

if rand(1)>0.2    x(i)=randn(1)+2;

else              x(i)=randn(1)+1;

end

end

a(1,1)=0.8;mu(1,1)=3;sm(1,1)=2;      %给定初始值

a(2,1)=0.2;mu(2,1)=-2;sm(2,1)=7;

for n=1:30                          %EM迭代次数为30次

for j=1:2                        %高斯分量2个

for i=1:N

p1(j)=a(j,n)*normpdf(x(i),mu(j,n),sm(j,n));

P=0;  %清零

for s=1:2

P=a(s,n)*normpdf(x(i),mu(s,n),sm(s,n))+P;  %累加

end

p(j,i)=p1(j)/P;

end

temp1=0; temp2=0;temp3=0;

for i=1:N

temp1=p(j,i)+temp1;  %累加

temp2=p(j,i)*x(i)+temp2;  %累加

temp3=p(j,i)*((x(i)-(temp2/temp1)).^2)+temp3;  %累加

end

a(j,n+1)=temp1/N;  %得到后面的项

mu(j,n+1)=temp2/temp1;

sm(j,n+1)=sqrt(temp3/temp1);

end

end

matlab编写的程序输入参数怎么写,用matlab编写一个用EM算法估计参数的程序相关推荐

  1. EM算法 估计混合高斯模型参数 Python实现

    EM算法 估计混合高斯模型参数 Python实现 EM算法是一种用来解决含有隐变量问题的算法,混合高斯模型中对于某个数据我们并不知道是来自于哪个模型,因此可以视为隐 变量,故可以采用隐含高斯模型来求解 ...

  2. EM算法估计GMM (二维高斯混合模型)参数 python 2.7(西瓜书版)

    #!/usr/bin/python # coding:utf-8 # 19-5-23 上午10:07 # @File : EM.py import numpy as np import math im ...

  3. 使用EM算法估计GMM参数的原理及matlab实现

    相关数学概念 协方差矩阵 多维高斯分布 其中k=n,即x的维度. GMM的原理 GMM,高斯混合模型,是一种聚类算法. 1.GMM概念: -将k个高斯模型混合在一起,每个点出现的概率是几个高斯混合的结 ...

  4. 拼团小程序源码_纯小白如何做一个摄影电商拼团小程序?

    过去,摄影行业都是线下实体店为主,宣传手段单一.推广效率差,客户复购率更是不怎样.随着微信小程序的出现,各个摄影行业商家开始通过小程序来进行推广,连接线上线下场景,这样便能触达更多消费者. 小程序拥有 ...

  5. 可以进行单元测试么_为什么程序员都讨厌写单元测试?有一个词叫“相爱相杀”!...

    面对现实吧!没有人真的喜欢做单元测试.有很多人向我讲述他们超级讨厌单元测试.尽管有些人擅长于此,但对于我们大多数人而言,不管有多少抱怨.多少反感,单元测试都是一件必不可少的事情.今天,我将探讨为什么我 ...

  6. c语言程序输入自己班级学号,C语言中设计一个程序就是通过输入学号就能知道班级的程序...

    最后,说说电池,d80的电池据说跟d200是通用,跟d3000的电池相比也更重一点,其贮电量达到了1500毫安,要比d3000稍微耐用一些,而且d80还内置了对电池的检测系统,能够测试出电池的总体使用 ...

  7. 小程序能用vue写么_仿网易云音乐APP的微信小程序【小程序和Vue版本】

    小程序版本: 首先是网易云的音乐接口: npm 源码获取见文章底部. 目前实现功能 用户 歌单 FM 播放 评论 MV 专辑 歌手 登录 歌曲红心,FM trash,收藏单曲至歌单 收听记录 歌单歌曲 ...

  8. java语言编程三角形图形_编程题:编写程序输入三角形的3条边长,计算并输出... 求助一道JAVA编程题:编写一个类似记事本的图形用......

    导航:网站首页 > 编程题:编写程序输入三角形的3条边长,计算并输出... 求助一道JAVA编程题:编写一个类似记事本的图形用... 编程题:编写程序输入三角形的3条边长,计算并输出... 求助 ...

  9. em算法matlab图像应用,em算法matlab程序

    EM 算法作业 EM 算法简单 介绍及应用 EM 算法是当存在数据缺失问题时,极... Matlab 实现根据以上推导,可以很容易实现 EM 算法估计 GMM 参数.现... 题目:matlab 实现 ...

最新文章

  1. DM 源码阅读系列文章(二)整体架构介绍
  2. 整合Swagger2文档api
  3. 「权威发布」2019年电赛最全各类题目细节问题解答汇总
  4. 关于document.cookie的使用
  5. QT每日一练day27:绘制不规则窗体
  6. mui mui.plusReady() 事件中的变量问题;
  7. lixnu 启动php,好奇怪的代码
  8. unity 观察者模式
  9. 聚苯硫醚的全球与中国市场2022-2028年:技术、参与者、趋势、市场规模及占有率研究报告
  10. 用LoopBack接口配置EBGP邻居
  11. Spark执行流程详解
  12. 回收站里的文件删除了怎么恢复,其实很简单
  13. 车载诊断数据库ODX——ODX参数解析类型(上)
  14. 将Google Analytics for Firebase数据导入BigQuery
  15. jquery/javascript画心爱心方法
  16. 近来学习的一些东西和思考(2)
  17. 【Scala学习】之函数式风格编程
  18. C++实现判断一个数是否为素数
  19. 熊猫TV构建差异化竞争壁垒 游戏直播现新突破
  20. 谈谈软件项目合作开发

热门文章

  1. [Math]理解卡尔曼滤波器 (Understanding Kalman Filter)
  2. Delphi中的消息截获
  3. 欧盟科学家突破纳米尺度高清晰红外显微成像技术
  4. MySQL主从同步问题集
  5. CentOS6.6+Puppet3.7.3分布式部署Zabbix2.4监控系统
  6. Java魔法堂:深入正则表达式API
  7. 升级SharePoint场的时候, 运行Configuration Wizard需要有什么顺序么?
  8. mysql 唯一索引_MySQL学会用索引,让你数据库的查询速度起飞
  9. 【python教程入门学习】Pandas库下载和安装
  10. cmd find 命令 多个 或者 关系 +windows_Chocolatey -Windows系统的yum||apt 软件安装工具...