协方差驱动SSI模态参数识别教程：适合初学者的文本详解

本人也在学习的过程，如果有不对的地方请谅解！希望大家共同进步
第一步将时程数据加入白噪声后保存，第二步通过奇异值分解图判断状态矩阵的阶数，和系统阶数。第三步将加入白噪声的数据选中打开，并输入i和N，即可得到模态参数。

%协方差驱动SSI模态参数识别程序;  TimeVy需为偶数  此时选择用加入白噪声的TimeVY2 ，第一、二部分是给这一部分做铺垫。
clc;clear
[FileName,PathName]=uigetfile('*.txt','选择文件');  %显示选择时程文件的对话框
File=strcat(PathName,FileName);  %把内容横向连接成单个字符串  将括号内路径和文件名连成一条字符串/多个字符串连接成单个字符串
fip=fopen(File,'r');       %打开文件，r:读出   fip=+N表示文件打开成功，文件代号是N
TimeVy=load('-ascii',File);  %  load函数默认是读取.mat格式（二进制文件）用的，如果你是文本文件，请加上'-ascii'参数，用load的时候第一行注释要删除，直接从数据开始。 本句是载入时程信息
TimeVy=TimeVy(1:6000,1:6);  %得到添加白噪声的时间序列的行、列值
[TimePointNum,MesNodeTotalNum]=size(TimeVy); %[r,c]=size（A）将矩阵A的行数返回到第一个输出变量r，将矩阵的列数返回到第二个输出变量c  
MesNodeTotalNum=MesNodeTotalNum-1;   %将第一列设置为时间序列    （ 数据中第一列是时间序列 ）  测点数据就是除去第一列其他列
TimePointNum=TimePointNum;    %重新赋值一下时间点，   如100HZ采样60S，一共有6000个时间点，所以有6000行
prompt={'TOEPLITZ矩阵块数M(>=采样频率/(2*结构基频)):','系统阶数N(<i*实际测点数):'}; %所谓结构基频就是结构本身最小的那个固有频率 建立一个提示字符串  系统阶数可以通过第二部跌出的图来判断
dlg_title='Input Parameter of Calculate'; %赋值 输入计算参数  一会用于对话框的名称
num_lines=1; %赋值  一会用于对话框显示行数
def={'50','30'};       %赋值 用于对话框默认值    %系统阶数N=2*n
answer=inputdlg(prompt,dlg_title,num_lines,def); %输入对话框  inputdlg
M=str2num(answer{1,1});%将字符串类型转换为数值类型  例如50
N=str2num(answer{2,1});%将字符串类型转换为数值类型  例如30
j=TimePointNum-2*M;  % Hankel矩阵最后一个矩阵块对应坐标为（2*i，j）,即2*i+j需<=TimePointNum      相当于输入M=i   i需要人为确定
Hankel=zeros((2*M+1)*MesNodeTotalNum,j); %建立一个零矩阵分配内存，后期填数，因为矩阵每一列比前一列递增一个值，所以到最后一个值的时候等于2*i+j  自己可以举个简单例子，在纸上列举一下就明白了
for ti=1:j
      m=0;
      for tii=1:(2*M+1)     %这个循环就是构造一个矩阵，矩阵的每个元素都是一个列向量，列向量个数等于你监测点数，比如你安装了5个传感器，那么每个元素就是同一时间点的5个值。
           km=ti+tii-1;
            m=m+1;
            beginy=(m-1)*MesNodeTotalNum+1;%计算每次赋值在所在行的起始位置
              endy=m*MesNodeTotalNum;%计算每次赋值在所在行的结束位置
          Hankel(beginy:endy,ti)=(TimeVy(km,2:(MesNodeTotalNum+1)))'; %提取同一时刻的不同测点的监测值 即每一行块由L行组成，L为监测通道。