matlab程序求反应谱,MATLAB程序精确法求解反应谱

MATLAB程序精确法求解反应谱

”

ResponSespectrumProgram(精确法求解)

%%%%%%%%%%%%%%%%%%%%% 反应谱 精确法 程序 Begin With

matlab6.5%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% clear

% ***********读入地震记录***********

fid = fopen('CHI010.txt');

[Accelerate,count] = fscanf(fid,'%g'); %count 读入的记录的量

Accelerate=9.8*Accelerate'; %单位统一为 m和s

time=0:0.005:(count-1)*0.005; %单位 s

% ***********精确法计算各反应***********

%初始化各储存向量

Displace=zeros(1,count); %相对位移

Velocity=zeros(1,count); %相对速度

AbsAcce=zeros(1,count); %绝对加速度

% ***********A,B矩阵***********

DampA=[0,0.05,0.1]; %三个阻尼比

TA=0.0:0.05:6; %TA=0.000001:0.02:6; %结构周期

Dt=0.005; %地震记录的步长

%记录计算得到的反应，MDis为某阻尼时最大相对位移，MVel为某阻尼