功能需求：
     1、可自动读取多条txt格式的地震波的加速度（地震波的文件名可以用1-n的数字命名）。

     2、可自动按比例对各条地震波进行放缩和调整，但持续时间和频谱特性不作修改，仅对峰值加速度PGA进行调幅。

     3、三线型滞回曲线模型图中，初始刚度、第一阶段强度（如混凝土开裂点）、第二阶段强度（如钢筋屈服点）、第三阶段强度（如混凝土压碎）均已知。第二阶段、第三阶段的刚度退化系数为手动输入已知值。常用模型中，卸载段、反向卸载段终点连线的斜率均为0，此次模型计算中，卸载段终点连线的斜率为≥0（为手动输入已知值），反向卸载段终点连线的斜率为≤0（为手动输入已知值）。

     4、分别采用Wilson- θ法和Newmark-β法（两种方法可以分别编制MATLAB程序）按照三线型滞回曲线模型（模型如下图）在目标加速度a^\' (t)下对单自由度体系进行动力时程分析，并计算在目标加速度a^\' (t)下的最大层间位移角θ_max。a^\' (t)的起始值为5cm/s2，步长可5cm/s2。对于每一条地震波，当结构分析异常或k_i/k_1=10%时计算停止。
   
     5、将每条地震波各目标加速度a^\' (t)下的θ_max，即（a^\' (t)，θ_max）保存于统一的指定Excel文件中。

     6、求得每条地震波在指定〖θ_1〗_max、〖θ_2〗_max、〖θ_3〗_max、〖θ_4〗_max、〖θ_5〗_max下对应的目标加速度a^\' (t)，并保存于统一的指定Excel文件中。

     7、程序应具有尽量详细的注释。

功能需求：
     1、可自动读取多条txt格式的地震波的加速度（地震波的文件名可以用1-n的数字命名）。

     2、可自动按比例对各条地震波进行放缩和调整，但持续时间和频谱特性不作修改，仅对峰值加速度PGA进行调幅。

     3、三线型滞回曲线模型图中，初始刚度、第一阶段强度（如混凝土开裂点）、第二阶段强度（如钢筋屈服点）、第三阶段强度（如混凝土压碎）均已知。第二阶段、第三阶段的刚度退化系数为手动输入已知值。常用模型中，卸载段、反向卸载段终点连线的斜率均为0，此次模型计算中，卸载段终点连线的斜率为≥0（为手动输入已知值），反向卸载段终点连线的斜率为≤0（为手动输入已知值）。

     4、分别采用Wilson- θ法和Newmark-β法（两种方法可以分别编制MATLAB程序）按照三线型滞回曲线模型（模型如下图）在目标加速度a^\' (t)下对单自由度体系进行动力时程分析，并计算在目标加速度a^\' (t)下的最大层间位移角θ_max。a^\' (t)的起始值为5cm/s2，步长可5cm/s2。对于每一条地震波，当结构分析异常或k_i/k_1=10%时计算停止。
   
     5、将每条地震波各目标加速度a^\' (t)下的θ_max，即（a^\' (t)，θ_max）保存于统一的指定Excel文件中。

     6、求得每条地震波在指定〖θ_1〗_max、〖θ_2〗_max、〖θ_3〗_max、〖θ_4〗_max、〖θ_5〗_max下对应的目标加速度a^\' (t)，并保存于统一的指定Excel文件中。

     7、程序应具有尽量详细的注释。