应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 马氏距离剔除光谱数据异常值
131/1返回列表
查看: 2855  |  回复: 12
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

张小小旭

铜虫 (小有名气)

[交流] 马氏距离剔除光谱数据异常值

请问大家，我在应用马氏距离剔除光谱数据异常值，我看到关于马氏距离有这样的描述“求在计算马氏距离过程中，要求总体样本数大于样本的维数，否则得到的总体样本协方差矩阵逆矩阵不存在，这种情况下，用欧式距离计算即可”

对于光谱数据来说，一般样本数几百个，每个样本光谱数据超过一千个，我用pdist函数报错显示
错误使用 pdist (line 190)
the covariance matrix for the mahalanobis metric must be symmetric and positive definite.
出错 shishi (line 132)
dist= pdist(x,'mahalanobis')

我看到很多论文是应用马氏距离剔除光谱数据异常值，这种样本数小于样本的维数情况下，是否可以实现马氏距离计算呢？
我的代码如下：

meanx=mean(x); %求均值
[m,n]=size(x);
covx=cov(x);%求协方差矩阵
cov_x=inv(covx);%逆矩阵
for i=1:m
dist(i)=(x(i,: )-meanx)*cov_x*(x(i,: )-meanx)'%求出每个样本到u的马氏距离,
end
dist= pdist(x,'mahalanobis')

x数据120样本，1900吸光度，用定义计算出马氏距离平方有负数，用pdist函数计算报错如上。
谢谢大家共同讨论！
回复此楼

» 猜你喜欢

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2019-12-27 16:17:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

wlpl

木虫 (小有名气)
★
张小小旭(金币+1): 谢谢参与
要是不用pca降维的话，可以Covx加一个很小的单位阵，保证可逆。

发自小木虫Android客户端
赞一下(1人) 回复此楼
3楼2019-12-27 17:38:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

张小小旭

铜虫 (小有名气)
可能是这样的方法《PCA结合马氏距离法剔除近红外异常样品 陈 斌 , 邹贤勇 , 朱文静》
赞一下 回复此楼
2楼2019-12-27 17:35:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

张小小旭

铜虫 (小有名气)
引用回帖:
3楼: Originally posted by wlpl at 2019-12-27 17:38:14
要是不用pca降维的话，可以Covx加一个很小的单位阵，保证可逆。

谢谢回复，请问是先求协方差矩阵cov，再加上单位阵，然后求逆矩阵吗？具体是如何加单位阵呢，多谢~
赞一下 回复此楼
5楼2019-12-27 19:21:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
tzynew4楼
2019-12-27 17:39   回复  
张小小旭(金币+1): 谢谢参与
psylhh6楼
2019-12-27 21:11   回复  
张小小旭(金币+1): 谢谢参与
唉尼摸尚靓7楼
2019-12-27 23:53   回复  
张小小旭(金币+1): 谢谢参与
。 发自小木虫Android客户端
LWJL20198楼
2019-12-28 10:58   回复  
张小小旭(金币+1): 谢谢参与
。 发自小木虫Android客户端
miaojiabing9楼
2019-12-30 08:10   回复  
张小小旭(金币+1): 谢谢参与
nono200910楼
2019-12-30 18:41   回复  
张小小旭(金币+1): 谢谢参与
·
星城石墨11楼
2019-12-31 12:00   回复  
张小小旭(金币+1): 谢谢参与
赞
youngen12楼
2019-12-31 12:07   回复  
张小小旭(金币+1): 谢谢参与
reminiscence13楼
2020-04-08 11:34   回复  
张小小旭(金币+1): 谢谢参与
相关版块跳转 我要订阅楼主 张小小旭 的主题更新
131/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085701环境工程319分,本科985，求调剂 +3 liuziqi_0622 2026-03-07 3/150 2026-03-14 05:31 by JourneyLucky
[考研] 考研一志愿南京理工大学化学工程专硕311分求调剂 +3 宁倾衍 2026-03-07 3/150 2026-03-14 05:13 by JourneyLucky
[考研] 347求调剂 +3 vevee 2026-03-08 3/150 2026-03-14 03:38 by JourneyLucky
[考研] 265求调剂 +9 小木虫085600 2026-03-09 12/600 2026-03-14 01:11 by JourneyLucky
[考研] 调剂 +3 13853210211 2026-03-10 3/150 2026-03-14 00:47 by JourneyLucky
[考研] b区环境工程求调剂 +4 Maps1 2026-03-10 6/300 2026-03-14 00:23 by JourneyLucky
[考研] 285 求调剂 资源与环境 一志愿北京化工大学 +3 未名考生 2026-03-10 3/150 2026-03-13 23:04 by JourneyLucky
[考研] 材料专硕288分求调剂 一志愿211 +4 在家想你 2026-03-11 4/200 2026-03-13 22:49 by JourneyLucky
[考研] 材料工程调剂 +9 咪咪空空 2026-03-12 9/450 2026-03-13 22:05 by 星空星月
[考研] 材料工程调剂 +4 咪咪空空 2026-03-11 4/200 2026-03-13 19:57 by JourneyLucky
[考研] 328求调剂 +3 5201314Lsy！ 2026-03-13 5/250 2026-03-13 16:05 by 5201314Lsy！
[考研] 求调剂 +3 程雨杭 2026-03-12 3/150 2026-03-13 15:06 by JourneyLucky
[考研] 材料301分求调剂 +5 Liyouyumairs 2026-03-12 5/250 2026-03-13 14:42 by JourneyLucky
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
[考研] 材料调剂，307分 +13 张泳铭1 2026-03-09 17/850 2026-03-13 11:09 by 薛云鹏
[考研] 288求调剂085600材料与化工 +14 Daunrin 2026-03-07 16/800 2026-03-12 16:19 by songkelly
[考研] 420求调剂 +4 莫向外求11 2026-03-10 6/300 2026-03-12 14:41 by ruiyingmiao
[考研] 收调剂 +7 调剂的考研学生 2026-03-10 7/350 2026-03-10 17:57 by 麦茶汤圆
[考研] 考研调剂，一志愿山东大学材料与化工，328分，政治51 +5 关你西红柿929 2026-03-08 6/300 2026-03-09 13:50 by 新篇章DFSS
[考研] 293一志愿华东理工 0817化学工程与技术 调剂 +5 fjj0912 2026-03-07 5/250 2026-03-09 09:13 by 30660438
信息提示
关闭
请填处理意见
关闭