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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 马氏距离剔除光谱数据异常值
131/1返回列表
查看: 3123  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 时间戳今天，20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 2026国自然函评费到账 +9 羊腰板 2026-08-21 9/450 2026-08-22 00:03 by anji999
[基金申请] 建议基金发布提前给出明确的时间点 +10 kulium 2026-08-21 13/650 2026-08-21 21:58 by alongwaytogo
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 让我中一个面上吧！ +11 大萍1987 2026-08-20 13/650 2026-08-21 14:48 by 20120902066
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 只有每年这种时候来逛逛小木虫 +23 yaoyewhu2008 2026-08-20 25/1250 2026-08-21 12:19 by Poppy1104
[基金申请] 什么时候开奖？ +7 CrisMessi 2026-08-18 7/350 2026-08-21 12:15 by weiyin
[基金申请] 今天基金会出结果吗？20260819 +15 kkkl_v 2026-08-19 16/800 2026-08-21 11:22 by 365372687
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 今天系统多次维护，明天很可能放榜！ +10 zju2000 2026-08-16 11/550 2026-08-20 20:02 by cl479861084
[基金申请] 时间戳变了，能看出什么问题？ +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 重要消息，中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[教师之家] 为什么余额宝的年化利率越来越低？主要原因有哪些？ +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
关闭
请填处理意见
关闭