随机子空间法的修正之后代码,,,,,%1、噪声添加程序 y向加速度
%读取数据文件
clc;clear
[FileName,PathName]=uigetfile('*.txt','Select the TimeVy File'); %显示选择时程文件的对话框
File=strcat(PathName,FileName); %把数横向连接成单个字符串
fip=fopen(File,'r');
TimeVy=load('-ascii',File); %载入时程信息
[TimePointNum,MesNodeTotalNum]=size(TimeVy); %这个矩阵与hankel块方向不同,下面为转置
MesNodeTotalNum=MesNodeTotalNum-1; %第一列为时间序列
TimePointNum=TimePointNum;
Baizaosheng=wgn(TimePointNum,MesNodeTotalNum,1);
2022-08-25 18:51:12
342KB
ssi
1