1、clearclcfile1=TEST20180626(955-1005).xlsx; %打开截取好时间范围的 excel 文件sheet=1; %选取 excel 文件中的工作表项T=2;M=8;N=9; %数据中初始行数 record1,1=时间; %创建一个广义矩阵用于存放数据，并设置标题栏分别为时间，流速，流向record1,2=流速m/s;record1,3=流向;for i=1:1:172; %构造循环选取每个采样里的 Xm 及邻近对应的数据a=num2str(T); %将选取区域所在的行数变量从数字变为字符串 b=num2str(M); %将选取区域左上角的行数变量从数字变为字符串
2、c=num2str(N); %将选取区域右下角的行数变量从数字变为字符串xlRange=strcat(D,a,:,F,a); %将选取区域代表的各字符串组成一个整体t_n=xlsread(file1,sheet,xlRange); %随 i 变化，选取区域中的时间数据hh=num2str(t_n(1)+8); %将选取的 小时 从格林尼治时间变为北京时间，并变为字符串mm=num2str(t_n(2); %将选取的 分 变为字符串ss=num2str(t_n(3); %将选取的 秒 变为字符串combt=strcat(hh,:,mm,:,ss); %将小时，分，秒字符串合并recordi+1,
3、1=combt; %将时间字符串储存到广义矩阵中xlRange=strcat(B,b,:,D,c); %将选取区域代表的各字符串组成一个整体f_n=xlsread(file1,sheet,xlRange); %随 i 变化，选取不同位置中的区域的数据name_string = f num2str(i) =f_n; %创建随提取数据位置改变而改变的字符串变量eval(name_string); %将 name_string 的字符串视为语句并运行x_mod=(f_n(2)*cos(f_n(3)+f_n(5)*cos(f_n(6)/2 %将当前选取区域中速度分解到横坐标y_mod=(f_n(2)*sin(f_n(3)+f_n(5)*sin(f_n(6)/2 %将当前选取区域中速度分解到纵坐标vector=complex(x_mod,y_mod); %通过 x_mod 和 y_mod 构建复数来简化构造矢量velocity=abs(vector)/100; %求 SI 下流速大小的模recordi+1,2=velocity; %记录流速大小recordi+1,3=angle(vector); %记录流向M=M+41;N=N+41;T=T+41; %改变行数end
