function [kd,s,b] = drillselection(a)
%統計地層可鑽性應用；
%根據抽樣測定的乙個地層整體的岩石可鑽性級值資料；
% 計算出1、地層整體的統計地層可鑽性kd;
% 2、地層級值波動範圍s;
% 3、驗證地層可鑽性級別統計分布圖和理論分布圖等；
% 例：輸入測得實驗資料
% a=load('result.txt');
% 輸出岩石可鑽性等級統計表
% b=[l;l+1;k;ck;fk;fn;pk;f;d]』；
% 其中l和l+1兩列表示等級區間，k表示頻數，ck表示累計頻數；
% fk表示頻率，fn表示累計頻率，pk表示概率，f表示累計概率，d表示abs(fn-f);
a=load('result.txt');
a=a(:)';
%% 求均值和方差
kd=mean(a); %均值
s=std(a); %方差
%% 求各級值區間頻數，累計頻數，頻率和累計頻率
for i=1:10
n=find(a>=(i-1)&ak(i)=size(n,2); %頻數
endk;
ck=cumsum(k); %累計頻數
sk=sum(k);
fk=k./sk; %頻率
fn=cumsum(fk); %累計頻率
%% 求各區間概率和累計概率以及d
for j=1:10
pk(j) = normpdf(j,kd,s); %概率
endpk;
f=cumsum(pk); %累計概率
d=abs(fn-f);
%% 輸出岩石可鑽性等級統計表及分布圖
l=[0,1,2,3,4,5,6,7,8,9];
b=[l;l+1;k;ck;fk;fn;pk;f;d]';
figure,bar(fk,1)
set(gca,'xticklabel',);
xlabel('地層可鑽性級值');
ylabel('頻率概率');
hold on;
x=1:10;
values = spcrv([[x(1) x x(end)];[pk(1) pk pk(end)]],3);
plot(values(1,:),values(2,:))
hold off;
end
MATLAB概率統計
一 產生隨機變數 二項分布隨機資料產生 n1 10 10 60 a1 binornd n1,1.n1 b1 binornd n1,1.n1,1,6 一行六列 c1 binornd n1 n1 1.n1 1.n1 2,6 兩行六列 正態分佈隨機資料產生 a2 normrnd 0,1,1,5 標準正態分...
Matlab之統計函式
年紀越大發現記憶越差了，看來好記憶確實不如爛筆頭。用matlab已經有些年了，但是一些常用函式時間長了不用還是記不住。於是乎決定花點心思整理下 基本統計函式 均值 mean x 注意預設是對列求均值 mean x,2 對行求均值 方差 tabulate x 統計頻數 由於本人職業所需，經常需要統計考...
matlab雜記2概率統計
1，函式binornd可以產生二項分布隨機資料 r binornd n,p n,p為二項分布的兩個引數，n,p 大小相同。r binornd n,p,m m指的是隨機數的個數，與r同維數，r bonornd n,p,m,n m,n分別表示r的行數和列數 2，normrnd可以產生正態分佈的隨機資料。...