2009年6月22日 星期一

三、使用Matlab進行分析

(2)使用ttest(t-檢定)

使用第一段和第二段的礫石成份進行比較,首先將兩者所佔的成份畫成直方圖
十號橋剖面礫石成份依橫軸依序從1-10為變質砂岩、板岩、脈石英、泥岩、蛇紋岩、輝長岩、玄武岩、基性岩、安山岩、安山岩質凝灰岩,縱軸為百分比。
最後利用Ttest(t檢定)選擇信心程度為95%,後得到兩者為同一岩層成立, h為0,然而其significance接近於1,因此相似度極大,信心範圍為[-0.0203 0.0203]
程式碼
bb=load('sediments2.txt')%在入資料%
first=bb(:,1);%第一段成份百分比%
second=bb(:,2);%第二段成份百分比%
h1= bar(bb(:,1));%畫第一段直方圖%
hold on%留住第一張圖,並把下面的也畫到第一張圖上
h2 = bar(bb(:,2));%畫第二段直方圖%
set(h1, 'FaceColor', 'none','EdgeColor', 'r')%讓直方圖內部中空,邊框為紅色%
set(h2, 'FaceColor', 'none','EdgeColor', 'b')%讓直方圖內部中空,邊框為藍色%
[h,significance,ci]=ttest(first, second, 0.05)%算ttest

沒有留言:

張貼留言