版块导航: 正在加载中...

客户端APP下载

登录注册

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

24小时热门版块排行榜

vikinglife

金虫 (正式写手)

应助: 0 (幼儿园)
金币: 1262.8
散金: 970
红花: 1
帖子: 391
在线: 88.8小时
虫号: 1420964
注册: 2011-09-28
专业: 电力电子学

[求助] 求助大虾看下这个程序

麻烦大家帮我看下这个程序实现什么功能输出值各代表什么含义谢谢

function [contribution,T2,SPE,t2cl,s_cl] = PCA_model(Xtrain,Xtest)
X_mean = mean(Xtrain);
X_std = std(Xtrain);
[X_row ,X_col]= size(Xtrain);
for i = 1:X_col
Xtrain(:,i) = (Xtrain(:,i)-X_mean(i))./X_std(i);
Xtest(:,i) = (Xtest(:,i)-X_mean(i))./X_std(i);
end
[U,S,V]=svd(Xtrain./sqrt(size(Xtrain,1)-1),0);
D= S^2;
lamda=diag(D);
num_pc=1;
while sum(lamda(1:num_pc))/sum(lamda)<0.9
num_pc=num_pc+1;
end
D=diag(lamda);
P=V(:,1:num_pc);
[a,b]=size(Xtest);
[r,y]=size(P*P');
I=eye(r,y);
e=Xtest*(I-P*P');
for i=1:a
T2(i)=Xtest(i,

*P*inv(D(1:num_pc,1:num_pc))*P'*Xtest(i,

';
end
for l=1:a
SPE(l)=e(l,

*e(l,

';
end

for j=1:b
contribution(j)=(norm(e(:,j)))^2;
end

t2cl=num_pc*(X_row-1)*(X_row+1)*icdf('f',0.99,num_pc,X_row-num_pc)/(X_row*(X_row-num_pc));
for i=1:3
   theta(i)=trace((D(num_pc+1:X_col,num_pc+1:X_col))^i);
end
% 另一种SPE控制线算法
% h=(theta(1)^2)/theta(2);
% g=theta(2)/theta(1);
% conf=0.95;
% df=round(h);
% delta2a1=g*pinv(df,2);

h0=1-2*theta(1)*theta(3)/(3*theta(2)^2);
ca=icdf('norm',0.99,0,1);
s_cl=theta(1)*(ca*sqrt(2*theta(2)*h0^2)/theta(1)+1+theta(2)*h0*(h0-1)/theta(1)^2)^(1/h0);

» 猜你喜欢

拟解决的关键科学问题还要不要写已经有11人回复
救命帖已经有5人回复
限项规定已经有5人回复
为什么nbs上溴没有产物点出现呢已经有9人回复
招博士已经有3人回复
存款400万可以在学校里躺平吗已经有35人回复
最失望的一年已经有18人回复
求推荐博导已经有4人回复
求推荐英文EI期刊已经有5人回复
疑惑？已经有5人回复

» 本主题相关价值贴推荐，对您同样有帮助:

求助各位大虾，孪生金相中这个形貌是什么组织，要如何解释？已经有6人回复
求助SAS 多元非线性回归拟合编程已经有3人回复
【求助】请哪位大虾帮忙查下这个单晶结构有没有报道过已经有6人回复
【求助】关于未知树状结构存储的问题？请大虾指教！已经有11人回复
【求助/交流】纳豆激酶制备过程中蛋白含量测定结果令我很费解，希望各位大虾多多指教� 已经有14人回复
【求助/交流】那位大虾帮忙好好解释一下CDS和cdna的区别和联系？？已经有17人回复
【求助】看文献中遇到一个问题，请教各位大虾! 已经有10人回复
【求助】请各位大虾帮我看一下，我的图谱和物质对的上吗已经有5人回复
【讨论】求助：大虾们觉得这个怎么回比较好已经有4人回复

1楼 2012-06-10 11:12:17

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

gbl404173223

银虫 (小有名气)

应助: 3 (幼儿园)
金币: 262.7
帖子: 117
在线: 44.6小时
虫号: 1811707
注册: 2012-05-12
性别: GG
专业: 控制理论与方法

【答案】应助回帖

★
感谢参与，应助指数 +1
臭水沟: 金币+1, 谢谢交流~~ 2012-06-11 22:09:11

这个程序描述的是主成分分析法，输出结果是可选分布的累积分布函数和逆累积分布函数，具体理论见教材《主成分分析》

没什么大不了的！

2楼2012-06-11 17:17:23

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

vikinglife

金虫 (正式写手)

应助: 0 (幼儿园)
金币: 1262.8
散金: 970
红花: 1
帖子: 391
在线: 88.8小时
虫号: 1420964
注册: 2011-09-28
专业: 电力电子学

引用回帖:

2楼: Originally posted by gbl404173223 at 2012-06-11 17:17:23
这个程序描述的是主成分分析法，输出结果是可选分布的累积分布函数和逆累积分布函数，具体理论见教材《主成分分析》

嗯是的，是主成分分析。。。但我一直不知道它做异常数据检测时
[contribution,T2,SPE,t2cl,s_cl] = PCA_model(Xtrain,Xtest) 输出参数各有什么意义

3楼2012-06-13 14:17:54

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 vikinglife 的主题更新