氣象統計 第三次實習

一、實習題目

二、實習程序

import pandas as pd import numpy as np##讀文件 data1=np.array(pd.read_table('mon1.txt',header=None)) data2=np.array(pd.read_table('mon2.txt',header=None)) data12=np.array(pd.read_table('mon12.txt',header=None)) ##簡單相關系數,注意x,y長度相等 def xgxs(x,y):n1=len(x) ave1=np.mean(x)ave2=np.mean(y)s1= np.std(x)s2= np.std(y)sum=0for i in range(0,n1):sum=(x[i]-ave1)*(y[i]-ave2)+sumr=sum/s1/s2/n1return r ###落后交叉相關系數函數,注意x,y長度相等 def lhjcxgxs(x,y):n1=len(x)N=int(n1/2)ave1=np.mean(x)ave2=np.mean(y)s1= np.std(x)s2= np.std(y)r=np.zeros(N)for i in range(0,N):sum=0for j in range(0,n1-i-1):sum=((x[j]-ave1)/s1)*((y[i+1+j]-ave2)/s2)+sumr[i]=1/(n1-i-1)*sumreturn r##12月與1月的落后交叉相關系數 print('12月與1月的落后交叉相關系數') r_xy=lhjcxgxs(data12,data1) for i in range(0,10):print('t=%d r_xy=%lf'%((i+1),r_xy[i])) ##12月與1月的偏相關系數 r=(xgxs(data12,data1)-xgxs(data1,data2)*xgxs(data12,data2))/((1-xgxs(data1,data2)**2)*(1-xgxs(data12,data2)**2))**0.5 print('12月與1月的偏相關系數:%lf'%(r)) ##12月與2月的落后交叉相關系數 print('12月與2月的落后交叉相關系數') r_xy=lhjcxgxs(data12,data2) for i in range(0,10):print('t=%d r_xy=%lf'%((i+1),r_xy[i])) ##12月與2月的偏相關系數 r=(xgxs(data12,data2)-xgxs(data1,data2)*xgxs(data12,data1))/((1-xgxs(data12,data1)**2)*(1-xgxs(data1,data2)**2))**0.5 print('12月與2月的偏相關系數:%lf'%(r))

三、例圖

四、代碼及源文件

實習三

程序中如有疏漏歡迎指正

總結

以上是生活随笔為你收集整理的气象统计 第三次实习 落后交叉相关系数和偏相关系数的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。