三、水文資料分析(一)

1.月均流量Histogarm

程式碼
clear all; close all
bio=dlmread('flow.txt');%讀取資料


flow=bio(1:end,2);figure(1);%選擇資料區間
[n1,x1]=hist(flow);h1=bar(x1,n1);%繪製長條圖
title('Histogram of 1969~2008 立霧溪月均流量');%加標題
xlabel('水流量 cms')%x軸標示
ylabel('數量 month')%y軸標示
mean(flow)%平均
std(flow)%標準差

資料區間 1959~2008，每月一筆。
mean = 3.2319e+004
std = 3.0804e+004

月均流量資料的初步呈現，
約70%的月均流量小於16570cms。

2.含沙量Histogram

程式碼
clear all; close all
bio=dlmread('sand.txt');%讀取資料


sand=bio(1:end,2);figure(1);%選擇資料區間
[n1,x1]=hist(sand);h1=bar(x1,n1);%繪製長條圖
title('Histogram of 1986~2008 立霧溪含沙量');%顯示標題
xlabel('含沙量 ppm')%x軸標示
ylabel('次數')%y軸標示
mean(sand)%平均
std(sand)%標準差

mean = 1.4952e+003
std = 4.0156e+003
資料區間1986~2008，記錄無固定間隔，

超過90%的含沙量小於5000ppm，
含沙量資料初步呈現。

3.1959~2008立霧溪月均流量分佈圖

程式碼
clear all; close all
f=load('flow.txt')
x=f(1:end, 1)
y=f(1:end, 2)
z=log(f(:,2))
figure(1); plot(x,y);
xlabel('時間序列1959~2008'); % X軸註解,一月一筆資料,以1959年1月為0、2月為1/12，
1960年1月為1，以此類推
ylabel('流量 cms');% Y軸註解
title('1959～2008年立霧溪月平均流量圖');
figure(2);plot(x,z);
xlabel('時間序列1959~2008');
ylabel('Log 流量 cms');% Y軸註解
title('Log 1959～2008年立霧溪月平均流量圖');

第一張圖為月均流量分佈圖，第二張圖為y軸取log。

一個月一筆資料以1959年1月為0、2月為1/12，1960年則為1，以此類推，
水流量的最大值出現在1999年，50年來並沒有特別乾旱的現象，
僅2003年的數值較往年稍小，表示立霧溪水量較為穩定，
水流量的分佈似遵循某個規律的循環，因為台灣乾濕季分明，
降雨明顯集中於7~10月的颱風季節，可能與此種季節降雨相關。