24小时热门版块排行榜    

查看: 2842  |  回复: 12
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

张小小旭

铜虫 (小有名气)


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

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

对于光谱数据来说,一般样本数几百个,每个样本光谱数据超过一千个,我用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函数计算报错如上。
谢谢大家共同讨论!
回复此楼

» 猜你喜欢

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

查看全部散金贴

已阅   回复此楼   关注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的回帖
查看全部 13 个回答

张小小旭

铜虫 (小有名气)


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

wlpl

木虫 (小有名气)



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

发自小木虫Android客户端
3楼2019-12-27 17:38:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
2019-12-27 23:53   回复  
张小小旭(金币+1): 谢谢参与
发自小木虫Android客户端
LWJL20198楼
2019-12-28 10:58   回复  
张小小旭(金币+1): 谢谢参与
发自小木虫Android客户端
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 收调剂 +6 调剂的考研学生 2026-03-10 6/300 2026-03-10 15:30 by 2204230311
[考研] 0856材料与化工309分求调剂 +4 ZyZy…… 2026-03-10 4/200 2026-03-10 15:20 by houyaoxu
[考研] 2026考研求调剂-材料类-本科211一志愿985-初试301分 +9 虫友233 2026-03-07 9/450 2026-03-10 15:06 by renwoxing13
[考研] 材料工程专硕调剂自荐信,初试323 +5 李白26 2026-03-07 5/250 2026-03-10 14:57 by ruiyingmiao
[考研] 材料与化工求调剂 +11 与冬清宁 2026-03-07 14/700 2026-03-10 14:23 by Linda Hu
[考研] 接受26届调剂生 +21 猪猪猪毛 2026-03-06 21/1050 2026-03-10 13:39 by Jack_333
[考研] 一志愿郑州大学,学硕,物理化学, 333求调剂 +8 李魔女斗篷 2026-03-04 9/450 2026-03-10 10:33 by allen-yin
[硕博家园] 2026级硕士研究生招生/调剂 +4 lbj6746988 2026-03-03 7/350 2026-03-10 08:15 by lbj6746988
[考研] 320求调剂 +4 魏zy 2026-03-08 4/200 2026-03-09 16:14 by ruiyingmiao
[基金申请] 第二个面上项目的获批几率如何? +3 wawm9296 2026-03-04 5/250 2026-03-09 13:00 by xujie4132
[考研] 求0856调剂 +5 小力气珂珂 2026-03-08 5/250 2026-03-09 06:42 by houyaoxu
[考研] 求调剂,一志愿江南大学,食品科学与工程,总分,320 +3 yyyyyukino 2026-03-07 3/150 2026-03-08 23:07 by 清风月
[考研] 材料工程307,求调剂 +4 我要燃烧你的梦 2026-03-08 4/200 2026-03-08 21:07 by 535743368
[基金申请] 进入个人成果库好难,一下午都没进去 +6 mi_dilee 2026-03-05 6/300 2026-03-08 11:43 by lan111
[考研] 化学工程与技术316求调剂 +7 Sixuan wang 2026-03-06 7/350 2026-03-08 11:23 by 我的船我的海
[考研] 材料化工求调剂 +7 博斯特525 2026-03-06 8/400 2026-03-08 03:51 by lfhuang
[考研] 材料科学(0805)338 求调剂 +7 xiaokang3286 2026-03-07 7/350 2026-03-08 03:42 by lfhuang
[考研] 282求调剂 +7 夕~日 2026-03-05 8/400 2026-03-05 21:31 by zzpnuli111
[考研] 322分 085600求调剂,有互联网+国金及主持省级大创经历 +6 熊境喆 2026-03-04 6/300 2026-03-04 20:32 by kakakapanpan
[考研] 085601 材料工程 320 +6 和乐瑶 2026-03-03 6/300 2026-03-04 16:01 by chixmc
信息提示
请填处理意见