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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 求助 大虾 看下这个程序
31/1返回列表
查看: 439  |  回复: 2
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

vikinglife

金虫 (正式写手)

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

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


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);
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-06-10 11:12:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gbl404173223

银虫 (小有名气)

【答案】应助回帖

★
感谢参与，应助指数 +1
臭水沟: 金币+1, 谢谢交流~~ 2012-06-11 22:09:11
这个程序描述的是主成分分析法，输出结果是可选分布的累积分布函数和逆累积分布函数，具体理论见教材《主成分分析》
赞一下 回复此楼
没什么大不了的！
2楼2012-06-11 17:17:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vikinglife

金虫 (正式写手)
引用回帖:
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 的主题更新
31/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭