24小时热门版块排行榜    

查看: 3214  |  回复: 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的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

wlpl

木虫 (小有名气)


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

发自小木虫Android客户端
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): 谢谢参与
2019-12-27 23:53   回复  
张小小旭(金币+1): 谢谢参与
。 发自小木虫Android客户端
LWJL20198楼
2019-12-28 10:58   回复  
张小小旭(金币+1): 谢谢参与
。 发自小木虫Android客户端
2019-12-30 08:10   回复  
张小小旭(金币+1): 谢谢参与
nono200910楼
2019-12-30 18:41   回复  
张小小旭(金币+1): 谢谢参与
·
2019-12-31 12:00   回复  
张小小旭(金币+1): 谢谢参与
赞
youngen12楼
2019-12-31 12:07   回复  
张小小旭(金币+1): 谢谢参与
2020-04-08 11:34   回复  
张小小旭(金币+1): 谢谢参与
相关版块跳转 我要订阅楼主 张小小旭 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] 同一个分子,一条来自文献,一条来自AI——不告诉你答案,你会选哪条? +5 tianxiaoxian 2026-09-28 9/450 2026-10-02 19:54 by yunshui01
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-10-01 21:59 by Equinoxhua
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 tqUTOxClQUMF 2026-09-28 4/200 2026-10-01 21:07 by Equinoxhua
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +7 IDs3scOF0tjC 2026-09-28 8/400 2026-10-01 19:08 by Equinoxhua
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 CfXuS1rDhLYN 2026-09-28 4/200 2026-10-01 09:23 by Equinoxhua
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 16:52 by Z9YWQ5EAO3qp
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 16:45 by Z9YWQ5EAO3qp
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 6/300 2026-09-30 16:26 by Z9YWQ5EAO3qp
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 13:35 by C3vABopYbZJS
[教师之家] 某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意! +10 zju2000 2026-09-26 13/650 2026-09-30 10:26 by yexuqing
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:37 by etmYJ6d2rquH
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:30 by etmYJ6d2rquH
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:46 by etmYJ6d2rquH
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 tqUTOxClQUMF 2026-09-28 4/200 2026-09-29 15:15 by etmYJ6d2rquH
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 14:58 by etmYJ6d2rquH
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 5/250 2026-09-29 13:50 by etmYJ6d2rquH
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:42 by ZvyPK8n6Nfki
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 6/300 2026-09-28 22:15 by ez6fGg9abYaj
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-28 19:08 by 9lS3ad5oOymn
信息提示
请填处理意见