1、读入资料read_origin_data

关键命令：textscan? fopen fclose

===================================================================
% function [time_array,on,temp]=read_origin_data(filename)
% % read_origin_data
% % Function name:[time_axis,on,temp]=read_origin_data(filename)
% % 读取原始asc文件，包括特殊资料 brk015，brk31，brk40
% % 得到时间轴time_array，向岸速度分量on，温度temp
% % 空缺值设为nan
% % input=
% %?? filename
% % output=
% %?? time_array(N*6), onshore speed arry(N*1), temperature array(N*1)
% % Last version: 2009-10-09
% % For test:
% % clear
% % clc
% % filename='1016brk011.asc';?????????? %普通格式
% % filename='1016brk015.asc';??????? %特殊格式
%
%
% directory='d:\study\mooring\origin_data\';
% fid=fopen([directory,filename],'r'); %打开mooring文件
%
% esp_filename={'1016brk015.asc','1016brk031.asc','1016brk040.asc'};
% format=strcmp(filename,esp_filename);
%
% % 读取mooring资料组成的元胞数组
% % 数据格式1：普通格式
% % HHMM:SS DDMMYY??? DIR'N??? SPEED?? E/COMP?? N/COMP???? TEMP
% %????????????????? (deg)?? (cm/s)?? (cm/s)?? (cm/s)? (deg C)
% % 1258: 0? 9 292?? 331.48??? 28.27?? -13.50??? 24.84??? 22.09
% % C1C2 C3C4C5C6???? C7?????? C8?????? C9?????? C10????? C11
%
% % 数据格式2：特殊格式
% %? Time?? Date??? Directn?? Speed??? East???? North???? Temp?
% %????????????????????????????? ?? ???? Comp???? Comp??????????
% % (hhmm) (d m y)?? (deg)?? (cm/s)?? (cm/s)?? (cm/s)?? (deg C)
% %? 1246? 9? 2 92?? 143.63??? 30.05??? 17.82?? -24.20??? 26.52
% %?? C1? C2 C3 C4??? C5??????? C6?????? C7????? C8??????? C9
%
% if? format(1)||format(2)||format(3); % 特殊格式
%???? C=textscan(fid,'%n%n%n%n%f%f%f%f%f','headerlines',10);
%???? %定义各量
%???? year=C{4}; ? ?? %年
%???? month=C{3};?? %月
%???? day=C{2}; ? ?? %日
%???? time=C{1};???? %时刻
%???? dir=C{5}; ? ? ? %流速方向,由北向顺时针旋转的角度
%???? speed=C{6}; %流速总量
%???? temp=C{9}; ? %温度
% else??????????????????????????????? % 普通格式
%????
%???? C=textscan(fid,'%n%c%s%n%1n%2n%f%f%f%f%f','headerlines',10);
%???? %定义各量
%???? year=C{6}; ? ? ? %年
%???? month=C{5}; ? %月
%???? day=C{4};? ? ??? %日
%???? time=C{1};???? ? %时刻
%???? dir=C{7};? ? ? ?? %流速方向,由北向顺时针旋转的角度
%???? speed=C{8}; ? %流速总量
%???? temp=C{11};?? %温度???????
% end
%
% % 时间序列，例如 （92+1900，2，3，2，23,0）年月日时分秒
% min=mod(time,100);
% hh=(time-min)/100;
% ss=zeros(length(min),1);
% time_array=[year+1900,month,day,hh,min,ss];
%
% % 筛选NaN值(-8888.0)
% speed(speed==-8888.0)=nan;
% temp(temp==-8888.0)=nan;
%
% % 坐标轴顺时针转66度后，流速方向dir2(弧度)=dir-66
% % 向岸为on 沿岸为along
% dir2=(dir-66)*pi/180;
% on=speed.*sin(dir2);
% % al=speed.*cos(dir2);
%
% fclose(fid);
% return