Python实现的计算马氏距离算法示例


Posted in Python onApril 03, 2018

本文实例讲述了Python实现的计算马氏距离算法。分享给大家供大家参考,具体如下:

Python实现的计算马氏距离算法示例

我给写成函数调用了

python实现马氏距离源代码:

# encoding: utf-8
from __future__ import division
import sys
reload(sys)
sys.setdefaultencoding('utf-8')
import numpy as np
def mashi_distance(x,y):
  print x
  print y
  #马氏距离要求样本数要大于维数,否则无法求协方差矩阵
  #此处进行转置,表示10个样本,每个样本2维
  X=np.vstack([x,y])
  print X
  XT=X.T
  print XT
  #方法一:根据公式求解
  S=np.cov(X)  #两个维度之间协方差矩阵
  SI = np.linalg.inv(S) #协方差矩阵的逆矩阵
  #马氏距离计算两个样本之间的距离,此处共有4个样本,两两组合,共有6个距离。
  n=XT.shape[0]
  d1=[]
  for i in range(0,n):
    for j in range(i+1,n):
      delta=XT[i]-XT[j]
      d=np.sqrt(np.dot(np.dot(delta,SI),delta.T))
      print d
      d1.append(d)
if __name__ == '__main__':
  # 第一列
  x = [3, 5, 2, 8]
  # 第二列
  y = [4, 6, 2, 4]
  mashi_distance(x,y)

运行结果:

Python实现的计算马氏距离算法示例

希望本文所述对大家Python程序设计有所帮助。

Python 相关文章推荐
Python遍历zip文件输出名称时出现乱码问题的解决方法
Apr 08 Python
Python3里的super()和__class__使用介绍
Apr 23 Python
Python字符串转换成浮点数函数分享
Jul 24 Python
Python 数据结构之旋转链表
Feb 25 Python
Pytorch实现神经网络的分类方式
Jan 08 Python
python绘制动态曲线教程
Feb 24 Python
keras中的History对象用法
Jun 19 Python
Python 制作查询商品历史价格的小工具
Oct 20 Python
关于pycharm 切换 python3.9 报错 ‘HTMLParser‘ object has no attribute ‘unescape‘ 的问题
Nov 24 Python
python 逆向爬虫正确调用 JAR 加密逻辑
Jan 12 Python
Python中json.load()和json.loads()有哪些区别
Jun 07 Python
Pandas数据结构之Series的使用
Mar 31 Python
python2.7无法使用pip的解决方法(安装easy_install)
Apr 03 #Python
快速解决安装python没有scripts文件夹的问题
Apr 03 #Python
Python中将变量按行写入txt文本中的方法
Apr 03 #Python
python使用tensorflow深度学习识别验证码
Apr 03 #Python
使用python获取csv文本的某行或某列数据的实例
Apr 03 #Python
python中实现数组和列表读取一列的方法
Apr 03 #Python
numpy.delete删除一列或多列的方法
Apr 03 #Python
You might like
编写PHP的安全策略
2006/10/09 PHP
超级简单的发送邮件程序
2006/10/09 PHP
php输出表格的实现代码(修正版)
2010/12/29 PHP
对淘宝URL中ID提取的PHP代码
2013/09/01 PHP
php中的路径问题与set_include_path使用介绍
2014/02/11 PHP
不常用但很实用的PHP预定义变量分析
2019/06/25 PHP
JS 密码强度验证(兼容IE,火狐,谷歌)
2010/03/15 Javascript
JQuery开发的数独游戏代码
2010/10/29 Javascript
JavaScript实现经典排序算法之插入排序
2016/12/28 Javascript
微信小程序实战之上拉(分页加载)效果(2)
2017/04/17 Javascript
bootstrap-table实现服务器分页的示例 (spring 后台)
2017/09/01 Javascript
vue中遇到的坑之变化检测问题(数组相关)
2017/10/13 Javascript
Node.js Stream ondata触发时机与顺序的探索
2019/03/08 Javascript
浅谈javascript中的prototype和__proto__的理解
2019/04/07 Javascript
Fetch超时设置与终止请求详解
2019/05/18 Javascript
微信小程序 scroll-view的使用案例代码详解
2020/06/11 Javascript
vue实现简单全选和反选功能
2020/09/15 Javascript
[00:58]他们到底在电话里听到了什么?
2017/11/21 DOTA
Python微信操控itchat的方法
2019/05/31 Python
django中使用事务及接入支付宝支付功能
2019/09/15 Python
Weekendesk意大利:探索多种引人入胜的周末主题
2016/10/14 全球购物
Under Armour澳大利亚官网:美国知名的高端功能性运动品牌
2018/02/22 全球购物
意大利制造的西装、衬衫和针对男士量身定制的服装:Lanieri
2018/04/08 全球购物
美国踏板车和轻便摩托车销售网站:Mega Motor Madness
2020/02/26 全球购物
县优秀教师事迹材料
2014/01/31 职场文书
买房协议书
2014/04/11 职场文书
石油工程专业毕业生求职信
2014/04/13 职场文书
农行心得体会
2014/09/02 职场文书
就业协议书
2014/09/12 职场文书
学生逃课检讨书1000字
2014/10/20 职场文书
党校党性分析材料
2014/12/19 职场文书
2015年教研工作总结
2015/05/23 职场文书
黄河绝恋观后感
2015/06/08 职场文书
新入职员工工作总结
2015/10/15 职场文书
MySQL令人咋舌的隐式转换
2021/04/05 MySQL
用python开发一款操作MySQL的小工具
2021/05/12 Python