MATLAB中SSI参数识别源程序

%数据驱动驱动随机子空间方法
clear all;
 clc;
 fid=fopen（‘Abc1#1.txt‘）;        
 x=fscanf（fid‘%f‘[1inf]）;    %读入速度数据，原始数据文件名在程序中输入
sf=79.9;  %采样频率

 for lll=1:1           %%%l个测点不包括参考点 
     imfh（lll:）=x（1:）;
 end
 
 n=length（imfh（1:））;
 c=50;l=1;             %%%%%c=i;定义的块行数＝2*c，l为测点数
 for i=0:2*c-1         %%%%%块行＝1:2i
 for j=1:n-2*c+1       %%%%%列＝1:n-2i+1  
 for k=1:l             %%%%%块行里的行＝测点数
 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加权法

 [USV]=svd（PP）;         %%%%%奇异值分解
 n=rank（S）;               %%%%理论上是按照此条命令，根据奇异值矩阵S的对角线非零元素个