日志
上传一个自己做的气候场代码
|
fid=fopen('C:\Users\Desktop\H500.txt','r');
fp=fopen('C:\Users\Desktop\qihou1.dat','w');
fp=fopen('C:\Users\Desktop\qihou5.dat','w');
data=zeros(17,37,48);m=zeros(37,17);
geo=zeros(17*37,2);
qihou1=zeros(17,37);qihou5=zeros(17,37);
qihou11=zeros(37*17,1);qihou55=zeros(37*17,1);
juping1=zeros(17,37);juping5=zeros(17,37);
juping11=zeros(37*17,1);juping55=zeros(37*17,1);
std1=zeros(17,37);std5=zeros(17,37);
std11=zeros(37*17,1);std55=zeros(37*17,1);
for k=1:48
m1=fscanf(fid,’%f %f’);
m=fscanf(fid,'%f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f',[37,17]);
data(:,:,k)=m';
end
%生成坐标
for i=1:17
for j=1:37
geo(37*(i-1)+j,1)=0+2.5*(i-1);geo(37*(i-1)+j,2)=60+2.5*(j-1);
end
end
%求一月气候场
for k=1:12:37
for i=1:17
for j=1:37
qihou1(i,j)=qihou1(i,j)+data(i,j,k)/4;qihou11(37*(i-1)+j,1)=qihou1(i,j);
end
end
end
fprintf(fp,'%f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f',qihou1);
%求五月气候场
for k=5:12:41
for i=1:17
for j=1:37
qihou5(i,j)=qihou5(i,j)+data(i,j,k)/4;qihou55(37*(i-1)+j)=qihou5(i,j);
end
end
end
fprintf(fp,'%f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f',qihou5);
%求1982年1月距平场
for i=1:17
for j=1:37
juping1(i,j)=data(i,j,1)-qihou1(i,j);juping11(37*(i-1)+j,1)=juping1(i,j);
end
end
%求1982年5月距平场
for i=1:17
for j=1:37
juping5(i,j)=data(i,j,5)-qihou5(i,j);juping55(37*(i-1)+j,1)=juping5(i,j);
end
end
%1月均方差场
for k=1:12:37
for i=1:17
for j=1:37
std1(i,j)=std1(i,j)+(data(i,j,k)-qihou1(i,j)).^2/4;std11(37*(i-1)+j,1)=std1(i,j);
end
end
end
%5月均方差场
for k=5:12:41
for i=1:17
for j=1:37
std5(i,j)=std5(i,j)+(data(i,j,k)-qihou5(i,j)).^2/4;std55(37*(i-1)+j,1)=std5(i,j);
end
end
end