clear clc close all hidden fs=10000; Data601=load('C:\Users\蓝紫\Desktop\01.txt'); Data602=load('C:\Users\蓝紫\Desktop\02.txt'); Data603=load('C:\Users\蓝紫\Desktop\03.txt'); Data604=load('C:\Users\蓝紫\Desktop\11.txt'); Data605=load('C:\Users\蓝紫\Desktop\12.txt'); Data606=load('C:\Users\蓝紫\Desktop\13.txt'); D601=Data601(:,3); D602=Data602(:,3); D603=Data603(:,3); detx=0.1; %归零 for i=1:1 %1/3倍频db值 [san601(:,i),f(:,i)]=sfzybp_db0(D601(:,i),fs); %计算测点1数据垂向加速度三分之一倍频程值,分别计算测点2,测点3的 [san602(:,i),f(:,i)]=sfzybp_db0(D602(:,i),fs); [san603(:,i),f(:,i)]=sfzybp_db0(D603(:,i),fs); A=[san601(:,i),san602(:,i),san603(:,i)]; %编制一个计算振动衰减率子程序 [DR(:,i),f(:,i)]=DRR(A,detx); end %画图 plot(DR,f);
时间: 2024-04-01 21:37:35 浏览: 14
这段代码的作用是读取六个文本文件中的数据,并进行数据处理和绘图。具体来说,它的主要步骤包括:
1. 使用 `load` 函数读取六个文本文件中的数据,分别存储在 `Data601` 到 `Data606` 中。
2. 从读取的数据中提取出第三列数据,分别存储在 `D601` 到 `D603` 中。
3. 使用函数 `sfzybp_db0` 计算三个测点的垂向加速度三分之一倍频程值,分别存储在 `san601` 到 `san603` 中。
4. 将三个测点的三分之一倍频程值组成一个矩阵 `A`。
5. 使用自编的 `DRR` 函数计算振动衰减率,存储在 `DR` 中。
6. 绘制振动衰减率与频率的变化关系图,使用 `plot` 函数完成。
需要注意的是,这段代码只对第一个测点的数据进行了处理,因为 `for` 循环只执行了一次。如果需要对多个测点进行处理,需要修改循环的终止条件和处理的数据。