pc_analysis.asv

统计模式识别工具包(统计模式识别工具包

PCA_data=cell(1214,2);
for i=1:1:1214
    X=data{i,1}';
    model=pca(X,1);
    Z = linproj(X,model);
    XR = pcarec(X,model);
    [dummy,mn] = min(Z);
    [dummy,mx] = max(Z);
    PCA_data{i,1}(1,1)=XR(1,mn);
    PCA_data{i,1}(2,1)=XR(1,mx);
    PCA_data{i,1}(1,2)=XR(2,mn);
    PCA_data{i,1}(2,2)=XR(2,mx);
end

for i=1:1:1214
    if PCA_data{i,1}(1,1)==PCA_data{i,1}(2,1)
        k=inf;
    else k=abs((PCA_data{i,1}(2,2)-PCA_data{i,1}(1,2))./(PCA_data{i,1}(2,1)-PCA_data{i,1}(1,1)));
    end
    SIZE_PCA=size(data{i,1},1);
    if (data{i,1}(1,2)-data{i,1}(SIZE_PCA,2))./(data{i,1}(1,1)-data{i,1}(SIZE_PCA,1))>=0
        PCA_data{i,2}(1,1)=atand(k);
end