1、 大氣綜合觀測實習報告1.實驗目的：掌握超聲風速儀數(shù)據(jù)后處理過程；掌握渦度相關法計算湍流通量；掌握基本湍流統(tǒng)計特征量的日、月變化過程。掌握基本氣象要素、輻射分量、能量通量的日、月變化過程。2.試驗原理：渦度相關法原理超聲風速儀數(shù)據(jù)質量控制渦度相關法計算湍流通量第一題：能量通量：clear;clf;year0=2007;year1=2007;month1=6;day1=1;year2=2008;month2=6;day2=30;stepxlabel=1;meandaystep=7; load datasetflux75861.dat;load datasetft75861.dat;load da

2、tasethour.txt;load datasetdoy.txt;load datasetrain7586.dat;load datasett7120.txt;load datasetp7120.txt;load datasetu7120.txt;load datasetprec7120.txt;load datasetdoy_int.txt;%history_DOYload datasetrain.txt;load datasetsun.txt;load datasetg0.txt;data=flux(:, 6);a= flux75861;b=size(a);c=a(9793:9840,6

3、:7);figureholdonplot(c(:,1)plot(c(:,2)d=a(17857:17904,12:13);figureholdonplot(d(:,1)plot(d(:,2)（參見圖能量通量-1）Co2通量：clear;clf;year0=2007;year1=2007;month1=6;day1=1;year2=2008;month2=6;day2=30;stepxlabel=1;meandaystep=7; load datasetflux75861.dat;load datasetft75861.dat;load datasethour.txt;load datasetd

4、oy.txt;load datasetrain7586.dat;load datasett7120.txt;load datasetp7120.txt;load datasetu7120.txt;load datasetprec7120.txt;load datasetdoy_int.txt;%history_DOYload datasetrain.txt;load datasetsun.txt;load datasetg0.txt;data=flux(:, 6);a= flux75861;b=size(a);c=a(9793:9840,12:13);figureholdonplot(c(:,

5、1)plot(c(:,2)d=a(17857:17904,12:13);figureholdonplot(d(:,1)plot(d(:,2)（參見圖co2-1,co2-2）短波輻射輻照度：clear;clf;year0=2007;year1=2007;month1=6;day1=1;year2=2008;month2=6;day2=30;stepxlabel=1;meandaystep=7; load datasetflux75861.dat;load datasetft75861.dat;load datasethour.txt;load datasetdoy.txt;load datase

6、train7586.dat;load datasett7120.txt;load datasetp7120.txt;load datasetu7120.txt;load datasetprec7120.txt;load datasetdoy_int.txt;%history_DOYload datasetrain.txt;load datasetsun.txt;load datasetg0.txt;data=flux(:, 6);a= ft75861b=size(a);c=a(9793:9840,3:4:7);figureholdonplot(c(:,1)plot(c(:,2)(參見圖短波輻射

7、輻照度-1)長波輻射輻照度：clear;clf;year0=2007;year1=2007;month1=6;day1=1;year2=2008;month2=6;day2=30;stepxlabel=1;meandaystep=7; load datasetflux75861.dat;load datasetft75861.dat;load datasethour.txt;load datasetdoy.txt;load datasetrain7586.dat;load datasett7120.txt;load datasetp7120.txt;load datasetu7120.txt;

8、load datasetprec7120.txt;load datasetdoy_int.txt;%history_DOYload datasetrain.txt;load datasetsun.txt;load datasetg0.txt;data=flux(:, 6);a= ft75861b=size(a);c=a(9793:9840, 11:4:15);figureholdonplot(c(:,1)plot(c(:,2)(參見圖長波輻射輻照度-1)總輻射與凈輻射：clear;clf;year0=2007;year1=2007;month1=6;day1=1;year2=2008;mont

9、h2=6;day2=30;stepxlabel=1;meandaystep=7;load datasetflux75861.dat;load datasetft75861.dat;load datasethour.txt;load datasetdoy.txt;load datasetrain7586.dat;load datasett7120.txt;load datasetp7120.txt;load datasetu7120.txt;load datasetprec7120.txt;load datasetdoy_int.txt;%history_DOYload datasetrain.

10、txt;load datasetsun.txt;load datasetg0.txt;data=flux(:, 6);a= ft75861;b=size(a);c=a(9793:9840, 11:4:15);figureholdonplot(c(:,1)plot(c(:,2)（參見圖總輻射與凈輻射-1）地面有效輻射：clear;clf;year0=2007;year1=2007;month1=6;day1=1;year2=2008;month2=6;day2=30;stepxlabel=1;meandaystep=7;load datasetflux75861.dat;load dataset

11、ft75861.dat;load datasethour.txt;load datasetdoy.txt;load datasetrain7586.dat;load datasett7120.txt;load datasetp7120.txt;load datasetu7120.txt;load datasetprec7120.txt;load datasetdoy_int.txt;%history_DOYload datasetrain.txt;load datasetsun.txt;load datasetg0.txt;data=flux(:, 6);a= ft75861;b=size(a

12、);c=a(9793:9840, 11:4:15);d=a(17857:17904, 11:4:15);figureholdonplot(c(:,1)plot(c(:,2)(參見圖地面有效輻射-1)反照率：clear;clf;year0=2007;year1=2007;month1=6;day1=1;year2=2008;month2=6;day2=30;stepxlabel=1;meandaystep=7;load datasetflux75861.dat;load datasetft75861.dat;load datasethour.txt;load datasetdoy.txt;loa

13、d datasetrain7586.dat;load datasett7120.txt;load datasetp7120.txt;load datasetu7120.txt;load datasetprec7120.txt;load datasetdoy_int.txt;%history_DOYload datasetrain.txt;load datasetsun.txt;load datasetg0.txt;data=flux(:, 6);bb= ft75861;up=bb(9793:9840,7);down=bb(9793:9840,3);for i=1:48albedo(i)=up(

14、i)/down(i)endzz=albedo;figureplot(zz)（參見圖反照率-1）光合有效輻射：clear;clf;year0=2007;year1=2007;month1=6;day1=1;year2=2008;month2=6;day2=30;stepxlabel=1;meandaystep=7; load datasetflux75861.dat;load datasetft75861.dat;load datasethour.txt;load datasetdoy.txt;load datasetrain7586.dat;load datasett7120.txt;load

15、 datasetp7120.txt;load datasetu7120.txt;load datasetprec7120.txt;load datasetdoy_int.txt;%history_DOYload datasetrain.txt;load datasetsun.txt;load datasetg0.txt;data=flux(:, 6);a= ft75861;b=size(a);c=a(9793:9840, 25);d=a(17857:17904, 25);figureholdonplot(c(:,1)plot(c(:,2)（參見光合有效輻射-1）地表反照率：clear;clf;

16、year0=2007;year1=2007;month1=6;day1=1;year2=2008;month2=6;day2=30;stepxlabel=1;meandaystep=7; load datasetflux75861.dat;load datasetft75861.dat;load datasethour.txt;load datasetdoy.txt;load datasetrain7586.dat;load datasett7120.txt;load datasetp7120.txt;load datasetu7120.txt;load datasetprec7120.txt

17、;load datasetdoy_int.txt;%history_DOYload datasetrain.txt;load datasetsun.txt;load datasetg0.txt;data=flux(:, 6);a= ft75861;b=size(a);c=a(9793:9840, 3:4:7);d=a(17857:17904, 3:4:7);figureholdonplot(c(:,1)plot(c(:,2)（參見地表反照率-1）第二題：clear;clf;year0=2007;year1=2007;month1=6;day1=1;year2=2008;month2=6;day

18、2=30;stepxlabel=1;meandaystep=7; load datasetflux75861.dat;load datasetft75861.dat;load datasethour.txt;load datasetdoy.txt;load datasetrain7586.dat;load datasett7120.txt;load datasetp7120.txt;load datasetu7120.txt;load datasetprec7120.txt;load datasetdoy_int.txt;%history_DOYload datasetrain.txt;loa

19、d datasetsun.txt;load datasetg0.txt;data=flux(:, 6);a= flux75861;b=ft75861;x=aa(:,34) aa(:,35) aa(:,36) aa(:,43) aa(:,44) aa(:,45); m6=x(1:30*48,:);m6=m6(abs(m6(:,1)0 & abs(m6(:,2)0 & abs(m6(:,3)0 & abs(m6(:,4)0 & abs(m6(:,5)0 & abs(m6(:,6)0,:);m7=x(1+30*48:61*48,:); m7=m7(abs(m7(:,1)0 & abs(m7(:,2)

20、0 & abs(m7(:,3)0 & abs(m7(:,4)0 & abs(m7(:,5)0 & abs(m7(:,6)0,:);m8=x(1+61*48:92*48,:);m8=m8(abs(m8(:,1)0 & abs(m8(:,2)0 & abs(m8(:,3)0 & abs(m8(:,4)0 & abs(m8(:,5)0 & abs(m8(:,6)0,:);m9=x(1+92*48:122*48,:);m9=m9(abs(m9(:,1)0 & abs(m9(:,2)0 & abs(m9(:,3)0 & abs(m9(:,4)0 & abs(m9(:,5)0 & abs(m9(:,6)0

21、,:);m10=x(1+122*48:153*48,:);m10=m10(abs(m10(:,1)0 & abs(m10(:,2)0 & abs(m10(:,3)0 & abs(m10(:,4)0 & abs(m10(:,5)0 & abs(m10(:,6)0,:);m11=x(1+122*48:153*48,:);m11=m11(abs(m11(:,1)0 & abs(m11(:,2)0 & abs(m11(:,3)0 & abs(m11(:,4)0 & abs(m11(:,5)0 & abs(m11(:,6)0,:);m12=x(1+183*48:214*48,:);m12=m12(abs

22、(m12(:,1)0 & abs(m12(:,2)0 & abs(m12(:,3)0 & abs(m12(:,4)0 & abs(m12(:,5)0 & abs(m12(:,6)0,:);m01=x(1+214*48:245*48,:);m01=m01(abs(m01(:,1)0 & abs(m01(:,2)0 & abs(m01(:,3)0 & abs(m01(:,4)0 & abs(m01(:,5)0 & abs(m01(:,6)0,:);m02=x(1+245*48:274*48,:);m02=m02(abs(m02(:,1)0 & abs(m02(:,2)0 & abs(m02(:,3

23、)0 & abs(m02(:,4)0 & abs(m02(:,5)0 & abs(m02(:,6)0,:);m03=x(1+274*48:305*48,:);m03=m03(abs(m03(:,1)0 & abs(m03(:,2)0 & abs(m03(:,3)0 & abs(m03(:,4)0 & abs(m03(:,5)0 & abs(m03(:,6)0,:);m04=x(1+305*48:335*48,:);m04=m04(abs(m04(:,1)0 & abs(m04(:,2)0 & abs(m04(:,3)0 & abs(m04(:,4)0 & abs(m04(:,5)0 & abs

24、(m04(:,6)0,:);m05=x(1+335*48:366*48,:);m05=m05(abs(m05(:,1)0 & abs(m05(:,2)0 & abs(m05(:,3)0 & abs(m05(:,4)0 & abs(m05(:,5)0 & abs(m05(:,6)0,:);m06=x(1+366*48:19001,:);m06=m06(abs(m06(:,1)0 & abs(m06(:,2)0 & abs(m06(:,3)0 & abs(m06(:,4)0 & abs(m06(:,5)0 & abs(m06(:,6)0,:);ww7m6=mean(m6,1);ww7m7=mean

25、(m7,1);ww7m8=mean(m8,1);ww7m9=mean(m9,1);ww7m10=mean(m10,1);ww7m11=mean(m11,1);ww7m12=mean(m12,1);ww8m01=mean(m01,1);ww8m02=mean(m02,1);ww8m03=mean(m03,1);ww8m04=mean(m04,1);ww8m05=mean(m05,1);ww8m06=mean(m06,1);data=ww7m6;ww7m7;ww7m8;ww7m9;ww7m10;ww7m11;ww7m12;ww8m01;ww8m02;ww8m03;ww8m04;ww8m05;ww8

26、m06;figureplot(data(:,1)figureplot(data(:,2)figureplot(data(:,3)figureplot(data(:,4)figureplot(data(:,5)figureplot(data(:,6)第三題：能量閉合度線型擬合：clear;clf;year0=2007;year1=2007;month1=6;day1=1;year2=2008;month2=6;day2=30;stepxlabel=1;meandaystep=7;load datasetflux75861.dat;load datasetft75861.dat;load data

27、sethour.txt;load datasetdoy.txt;load datasetrain7586.dat;load datasett7120.txt;load datasetp7120.txt;load datasetu7120.txt;load datasetprec7120.txt;load datasetdoy_int.txt;%history_DOYload datasetrain.txt;load datasetsun.txt;load datasetg0.txt;aa=flux75861;bb=ft75861;LE=aa(:,6);H=aa(:,7);G=bb(:,44);

28、Rn=bb(:,19);x=Rn-G;y=H+LE;cc=x y;cc=cc(abs(cc(:,2)0&abs(cc(:,1)0,:)p=polyfit(cc(:,1),cc(:,2),1);y1=polyval(p,cc(:,1)scatter(cc(:,1),cc(:,2)hold onplot(cc(:,1),cc(:,1)hold onplot(cc(:,1),y1(:,1)（參見圖3）第四題：aa= flux75861;bb=ft75861;xx=aa(:,6) aa(:,7) aa(:,8) bb(:,3) bb(:,7) bb(:,11) bb(:,15) bb(:,19) bb

29、(:,25);m6=xx(1:30*48,:);m6=m6(abs(m6(:,1)0 & abs(m6(:,2)0 & abs(m6(:,3)0 & abs(m6(:,4)0 & abs(m6(:,5)0 & abs(m6(:,6)0&abs(m6(:,7)0& abs(m6(:,8)0& abs(m6(:,9)0,:);m7=xx(1+30*48:61*48,:);m7=m7(abs(m7(:,1)0 & abs(m7(:,2)0 & abs(m7(:,3)0 & abs(m7(:,4)0 & abs(m7(:,5)0 & abs(m7(:,6)0 & abs(m7(:,7)0&abs(m7

30、(:,8)0 & abs(m7(:,9)0,:);m8=xx(1+61*48:92*48,:);m8=m8(abs(m8(:,1)0 & abs(m8(:,2)0 & abs(m8(:,3)0 & abs(m8(:,4)0 & abs(m8(:,5)0 & abs(m8(:,6)0& abs(m8(:,7)0& abs(m8(:,8)0& abs(m8(:,9)0,:);m9=xx(1+92*48:122*48,:);m9=m9(abs(m9(:,1)0 & abs(m9(:,2)0 & abs(m9(:,3)0 & abs(m9(:,4)0 & abs(m9(:,5)0 & abs(m9(:

31、,6)0&abs(m9(:,7)0&abs(m9(:,8)0&abs(m9(:,9)0,:);m10=xx(1+122*48:153*48,:);m10=m10(abs(m10(:,1)0 & abs(m10(:,2)0 & abs(m10(:,3)0 & abs(m10(:,4)0 & abs(m10(:,5)0 & abs(m10(:,6)0& abs(m10(:,7)0& abs(m10(:,8)0& abs(m10(:,9)0,:);m11=xx(1+122*48:153*48,:);m11=m11(abs(m11(:,1)0 & abs(m11(:,2)0 & abs(m11(:,3

32、)0 & abs(m11(:,4)0 & abs(m11(:,5)0 & abs(m11(:,6)0& abs(m11(:,7)0& abs(m11(:,8)0& abs(m11(:,9)0,:);m12=xx(1+183*48:214*48,:);m12=m12(abs(m12(:,1)0 & abs(m12(:,2)0 & abs(m12(:,3)0 & abs(m12(:,4)0 & abs(m12(:,5)0 & abs(m12(:,6)0& abs(m12(:,7)0& abs(m12(:,8)0& abs(m12(:,9)0,:);m1=xx(1+214*48:245*48,:);

33、m1=m1(abs(m1(:,1)0 & abs(m1(:,2)0 & abs(m1(:,3)0 & abs(m1(:,4)0 & abs(m1(:,5)0 & abs(m1(:,6)0& abs(m1(:,7)0& abs(m1(:,8)0& abs(m1(:,9)0,:);m2=xx(1+245*48:274*48,:);m2=m2(abs(m2(:,1)0 & abs(m2(:,2)0 & abs(m2(:,3)0 & abs(m2(:,4)0 & abs(m2(:,5)0 & abs(m2(:,6)0& abs(m2(:,7)0& abs(m2(:,8)0& abs(m2(:,9)0,

34、:);m3=xx(1+274*48:305*48,:);m3=m3(abs(m3(:,1)0 & abs(m3(:,2)0 & abs(m3(:,3)0 & abs(m3(:,4)0 & abs(m3(:,5)0 & abs(m3(:,6)0& abs(m3(:,7)0& abs(m3(:,8)0& abs(m3(:,9)0,:);m4=xx(1+305*48:335*48,:);m4=m4(abs(m4(:,1)0 & abs(m4(:,2)0 & abs(m4(:,3)0 & abs(m4(:,4)0 & abs(m4(:,5)0 & abs(m4(:,6)0& abs(m4(:,7)0& abs(m4(:,8)0& abs(m4(:,9)0,:);m