%最大積雪深度與灌溉面積之間的關系
%繪制散點圖并添加趨勢線
x=[15.2	10.4?21.2?18.6?26.4	23.4?13.5?16.7?24?19.1];
y=[28.6	19.3?40.5?35.6?48.9	45?29.2	34.1?46.7?37.4];
plot（xy‘Or‘）;????????????????????????%繪制散點圖
xlabel（‘最大積雪深度x‘）;????????????????%添加橫軸標簽
ylabel（‘灌溉面積y‘）;????????????????????%添加縱軸標簽
hold?on????????????????????????????????%保持圖形
lsline?????????????????????????????????%添加趨勢線

%基于矩陣運算計算回歸系數并添加趨勢線
[mn]=size（x）;?????????????????????????%計算自變量矩陣行列數
X=[ones（n1）x‘];??????????????????????%在自變量矩陣中添加常數向量并轉置
Y=y‘;??????????????????????????????????%因變量向量轉置
B=inv（X‘*X）*X‘*Y;??????????????????????%計算回歸系數
f=B（1）+B（2）*x;?????????????????????????%建設模型
plot（xf）;?????????????????????????????%基于模型預測值添加趨勢線
hold?off???????????????????????????????%第一次繪圖結束添加趨勢
%計算用于模型檢驗的主要統計量
R2=（abs（B‘*X‘*Y）-n*mean（y）^2）/（abs（Y‘*Y）-n*mean（y）^2）;?%計算擬合優度
s=sqrt（（Y‘*Y-B‘*X‘*Y）/（n-m-1））;????????%計算標準誤差
v=s/mean（y）;???????????????????????????%計算變異系數
F=（abs（B‘*X‘*Y）-n*mean（y）^2）/（m*s^2）;??%計算F統計量
p=inv（X‘*X）;???????????????????????????%計算一個方陣的逆矩陣
t1=B（1）/（sqrt（p（11））*s）;??????????????%計算截距的t統計量
t2=B（2）/（sqrt（p（22））*s）;??????????????%計算斜率的t統計量
e=y-f;?????????????????????????????????%計算殘差
i=1:n-1;???????????????????????????????%殘差編號
DW=sumsqr（e（i+1）-e（i））/sumsqr（e）;??????%計算Durbin-Watson統計量
BR2svFt1t2DW????????????????????%給出參數和統計量的計算值

%將F值和t值轉換為概率值
Sig=1-fcdf（Fmn-m-1）
P1=2*（1-tcdf（abs（t1）n-m-1））
P2=2*（1-tcdf（abs（t2）n-m-1））

%F檢驗、t檢驗和相關系數檢驗的等價性
F1=t2^2
F2=（（n-2）*R2）/（1-R2）
R=sqrt（F/（F+n-2））
%計算相關系數的臨界值
a=0.05;
Rc=sqrt（finv（1-a1n-2）/（finv（1-a1n-2）+n-2））
%預測
x=27.5;
y=B（1）+B（2）*x

%基于多項式擬合計算回歸系數并添加趨勢線
clear
x=[15.2	10.4?21.2?18.6?26.4	23.4?13.5?16.7?24?19.1];
y=[28.6	19.3?40.5?35.6?48.9	45?29.2	34.1?46.7?37.4];
plot（xy‘+r‘）;????????????????????????%繪制散點圖
xlabel（‘最大積雪深度x‘）;????????????????%添加橫軸標簽
ylabel（‘灌溉面積y‘）;????????????????????%添加縱軸標簽
title（‘最大積雪深度與灌溉面積的關系‘）;???%添加圖名
hold?on????????????????????????????????%保持圖形
B=polyfit（xy1）;??????????????????????%多項式擬合
lx=linspace（1028）;????????????????????%限定趨勢線的長度范圍
z=polyval（Blx）;???????????????????????%計算分割點上多項式的函數值
plot（lxz‘b-‘）;???????????????????????%將趨勢線添加到散點圖中
hold?off???????????????????????????????%繪圖結束
B??????????????????????????????????????%給出參數估計值

%計算用于檢驗的主要統計量
%計算相關系數和相關系數平方
C=corrcoef（xy）;???????????????????????%計算相關系數矩陣
R=C（12）;??????????????????????????????%提取相關系數
R2=R^2;????????????????????????????????%計算擬合優度
%計算標準誤差和變異系數
[mn]=size（x）;?????????????????????????%計算自變量矩陣行列數
f=B（1）*x+B（2）;?????????????????????????%建設模型
e=y-f;?????????????????????????????????%計算殘差
s=sqrt（（sumsqr（e）/（n-m-1）））;???????????%計算標準差
v=s/mean（y）;???????????????????????????%計算變異系數
%計算Durbin-Watson統計量
i=1:n-1;???????????????????????????????%殘差編號
DW=sumsqr（e（i+1）-e（i））/sumsqr（e）;??????%計算DW值
RsDW?????????????????????????????????%給出主要的統計量

%基于多項式擬合計算回歸系數并添加趨勢線
clear
x=[15.2	10.4?21.2?18.6?26.4	23.4?13.5?16.7?2