%數(shù)據(jù)驅(qū)動驅(qū)動隨機子空間方法
clear?all;
?clc;
?fid=fopen（‘Abc1#1.txt‘）;????????
?x=fscanf（fid‘%f‘[1inf]）;????%讀入速度數(shù)據(jù)，原始數(shù)據(jù)文件名在程序中輸入
sf=79.9;??%采樣頻率

?for?lll=1:1???????????%%%l個測點不包括參考點?
?????imfh（lll:）=x（1:）;
?end
?
?n=length（imfh（1:））;
?c=50;l=1;?????????????%%%%%c=i;定義的塊行數(shù)＝2*c，l為測點數(shù)
?for?i=0:2*c-1?????????%%%%%塊行＝1:2i
?for?j=1:n-2*c+1???????%%%%%列＝1:n-2i+1??
?for?k=1:l?????????????%%%%%塊行里的行＝測點數(shù)
?hank（i*l+kj）=imfh（ki+j）/（（n-2*c+1）^0.5）;?
?end
?end
?end
?
?%Yf=hank（c*l+1:2*c*l:）;
?%W1=（Yf*Yf‘）^（-0.5）;
?
?[QR]=qr（hank‘）;
?Q=Q‘;?R=R‘;?
?R21=R（l*c+1:end1:l*c）;
?Q1=Q（1:l*c:）;
?PP=R21*Q1;
?%PP=W1*PP;????????????????%%%%%CAV加權(quán)法

?[USV]=svd（PP）;?????????%%%%%奇異值分解
?n=rank（S）;???????????????%%%%理論上是按照此條命令，根據(jù)奇異值矩陣S的對角線非零元素個