%?第二步，根據(jù)雷達(dá)數(shù)據(jù)時(shí)間找到分別在該時(shí)間前后的兩個(gè)ECMWF數(shù)據(jù)文件，如14:05的雷達(dá)數(shù)據(jù)，則找到12:00和18：00的文件
%?第三步，根據(jù)雷達(dá)數(shù)據(jù)的經(jīng)緯度，分別在兩個(gè)文件中找到距離該位置最近的四個(gè)角經(jīng)緯度數(shù)據(jù)，做反距離加權(quán)，這樣可以得到兩個(gè)值
%?第四步，將這兩個(gè)值根據(jù)時(shí)間做反距離加權(quán)，就可以得到與雷達(dá)數(shù)據(jù)匹配的ECMWF數(shù)據(jù)（并沒有完全按照這個(gè)要求來啊，好像也沒事吧，反正結(jié)果就是好了）
%?close?all
%?clear%還有專門的netcdf工具包。。。2021b之后的matlab版本可以不用安裝了吧，只要用ncread就能讀取了
%?clc%nc文件轉(zhuǎn)換成mat文件——save出去就是mat文件
%?file2??=dir（‘E:\m\nc\‘）;%第一個(gè)第二個(gè)出來是.和..
%?buoy_4=load（‘buoy_012.dat‘）;%包含了時(shí)分秒，波高m，公尺波浪周期s，最大浪高的浪向or風(fēng)向or波向，用360°表示？
%?波高、波長、波周期、波向

fid=fopen（‘match_data2010_1.dat‘‘w‘）;
%?folder=‘E:\m‘;
%?files=dir（[folder?‘\*.hdf‘]）;%也可以改寫成folder=‘E:\m\‘;files=dir（[folder?‘*.hdf‘]）;的吧
%?h=length（files）;
filename=‘2A25.20100102.69110.7.hdf‘;%降雨雷達(dá)文件，先選用的是第2天的文件
hour=hdfread（filename‘Hour‘）;
minute=hdfread（filename‘Minute‘）;
lati=hdfread（filename‘Latitude‘）;
longi=hdfread（filename‘Longitude‘）;
[row?~]?=?size（longi）;%此處賦值給‘column’的值似乎未使用。請(qǐng)考慮將其替換為～，因?yàn)楹竺鎲为?dú)定義了一個(gè)column
%?for?i=1:length（longi（:1））??
day=hdfread（filename‘DayOfMonth‘）;

filename1=‘E:\m\nc\2010-1-00.nc‘;
u?=?ncread（filename1‘u10‘）;
v?=?ncread（filename1‘v10‘）;
ec00=sqrt（u.^2+v.^2）;
EC00=ec00（::1）;

filename2=‘E:\m\nc\2010-1-06.nc‘;
u6?=?ncread（filename2‘u10‘）;
v6?=?ncread（filename2‘v10‘）;
ec06=sqrt（u6.^2+v6.^2）;
EC06=ec06（::1）;

filename3=‘E:\m\nc\2010-1-12.nc‘;
u12?=?ncread（filename3‘u10‘）;
v12?=?ncread（filename3‘v10‘）;
ec12=sqrt（u12.^2+v12.^2）;
EC12=ec12（::1）;

filename4=‘E:\m\nc\2010-1-18.nc‘;
u18?=?ncread（filename4‘u10‘）;
v18?=?ncread（filename4‘v10‘）;
ec18=sqrt（u18.^2+v18.^2）;
EC18=ec18（::1）;

time=hour+minute/60;%這個(gè)時(shí)間是個(gè)數(shù)組，在執(zhí)行到69行那個(gè)與或的語句時(shí)，就出錯(cuò)了，必須為一個(gè)整的數(shù)值，加個(gè)循環(huán)吧？？？
%?if?0%-------------------------------------------------------------------------------------------------------------??
%%這部分是循環(huán)的，可以先不看吧。。。
%?column?=?[17181921222425262829313233];%
a=hdfread（filename‘rainAve‘）;
column=a（:[17181921222425262829313233]2）;%寫成這種格式嗎？
%?rain=a（:[17181921222425262829313233]2）;%經(jīng)緯度類似吧，下面原先寫的是錯(cuò)的
%?la=lati（:[17181921222425262829313233]）;
????for?i=1:row%把經(jīng)緯度和時(shí)間分別取出來進(jìn)行時(shí)空插值
????????for?j?=?1:13%對(duì)應(yīng)上面的13列之前是1：12
%?????????????lat=
????????????lat?=?lati（icolumn（j））;%索引超出矩陣維度，就是前面入射角對(duì)應(yīng)的行列數(shù)沒有寫完整導(dǎo)致的
%下標(biāo)索引必須為正整數(shù)類型或邏輯類型，問題就出在column（j）這里，看一看變量空間就發(fā)現(xiàn)要不是0（不是正數(shù)），要不是0.222（不是整數(shù)）
????????????lon?=?longi（icolumn（j））;%？？？這里應(yīng)該是longi吧，之前寫的是lon?=?lati（icolumn（j））;
%?????????????警告：賦值給變量‘lat’‘lon’的值可能未使用，這是因?yàn)楹竺嬗侄x了一個(gè)具體的經(jīng)緯度值，這個(gè)循環(huán)是之后各個(gè)點(diǎn)的，可暫時(shí)不看
????????????time1?=?time（i）;%先不加分號(hào)吧，這樣就能使時(shí)間變成了一個(gè)單獨(dú)的數(shù)的形式了，而不是一整個(gè)數(shù)組了
%?????????????if?time1>=12?&&?time1<18
%?????????????end

if?time1<6%這里的time表示的是pr的啊，跟ecmwf的不一樣，所以這部分還是不能省略掉
????????????dt1=time1;
????????????dt2=6-time1;
????????????ecdata1=EC00;
????????????ecdata2=EC06;
???????????elseif?time1>=6?&&?time1<12%||?和?&&?運(yùn)算符的操作數(shù)必須能夠轉(zhuǎn)換為邏輯標(biāo)量值time必須是一個(gè)數(shù)，你這里都成數(shù)組了
???????????????dt1=time1-6;
???????????????dt2=12-time1;
???????????????ecdata1=EC06;???
???????????????ecdata2=EC12;
???????????elseif?time1>=12?&&?time1<18
???????????????dt1=time1-12;

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件????????4894??2018-07-25?11:33??shiyishi.m