置信區(qū)間與假設(shè)檢驗matlab程序

1、 版權(quán)所有 嚴禁侵犯 不得用于商業(yè)用途統(tǒng)計學專用程序 -基于MATLAB 7.0開發(fā) -置信區(qū)間與假設(shè)檢驗 2013年8月1日 置信區(qū)間與假設(shè)檢驗程序【開發(fā)目的】眾所周知，統(tǒng)計工作面對的數(shù)據(jù)量繁瑣而且龐大，在統(tǒng)計的過程中和計算中容易出錯，并統(tǒng)計決定著國民經(jīng)濟的命脈，開發(fā)此軟件就是為了進行驗證統(tǒng)計的準確性以及理論可行性，減少統(tǒng)計工作中的錯誤，使統(tǒng)計工作者更好地進行工作與學習；所以特意開發(fā)此程序來檢驗統(tǒng)計中的參數(shù)估計和假設(shè)檢驗?！鹃_發(fā)特色】本軟件基于matlab7.0進行運算，對于樣本的輸入采用行矩陣的形式，并且開發(fā)了樣本頻數(shù)輸入，對于多樣本的輸入可以減緩工作量，對于顯著性水平本程序自帶正態(tài)分布函

2、數(shù)，t分布函數(shù)，F(xiàn)分布函數(shù)，分布函數(shù)的計算公式，不用再為查表和計算而苦惱，只需輸入顯著性水平即可，大大的簡化了計算量?！娟P(guān)鍵技術(shù)】矩陣輸入 進行頻數(shù)判斷 條件循環(huán)語句的使用等【程序界面】【程序代碼】此程序采用多文件結(jié)構(gòu)，在建立文件時不能改變文件名；以下是各個文件的代碼：（Zhucaidan.m)：clc;disp('統(tǒng)計學專用');disp('1.假設(shè)檢驗');disp('2.置信區(qū)間');disp('3.使用說明');disp('4.打開代碼');disp('0.退出程序');disp('

3、請進行選擇：');a=input('');if a=0 exit;else if a=1 jiashejianyan ; else if a=2 zhixinqujian ; else if a=3 help1; else if a=4 open('zhucaidan');disp(' 菜單選項');disp('1.返回主菜單');disp('2.退出程序！');p=input(' '); if p=1 zhucaidan; else if p=2 disp('正在退出，請稍候。&#

4、39;); (exit);end end end end end endend(Zhixinqujian.m) ：clc;disp(' 置信區(qū)間');disp(' 菜單選項');disp('0.退出程序!');disp('1.返回主菜單');disp('2.方差已知，待估參數(shù)為u');disp('3.方差未知，待估參數(shù)為u');disp('4.均值已知，待估參數(shù)為方差');disp('5.均值未知，待估參數(shù)為方差');disp('請進行選擇：');c=

5、input('');if c=0 exit;else if c=1 zhucaidan; else if c=2 uguji; else if c=3 uguji1; else if c=4 fguji; else if c=5 fguji1; end end end end endend(Jiashejianyan.m)：clc;disp('*假設(shè)檢驗*');disp('*菜單選項*');disp('*0.退出程序!*');disp('*1.返回主菜單*');disp('*2.已知，對u進行檢驗*'

6、;);disp('*3.未知，對u進行檢驗*');disp('*4.u已知，對進行檢驗*');disp('*5.u未知，對進行檢驗*');disp('*請進行選擇：');c=input('');if c=0 exit;else if c=1 zhucaidan; else if c=2 ujianyan; else if c=3 ujianyan1; else if c=4 fangcha; else if c=5 fangcha1; end end end end endendHelp.m：clc;disp(&#

7、39;-假設(shè)檢驗與置信區(qū)間-');disp('是否使用幫助（1/0)');c=input('');if c=1 clc; disp(' -假設(shè)檢驗與置信區(qū)間-'); disp('*幫助與操作說明*'); disp('1.請嚴格按照說明與提示進行操作，否則可能會出現(xiàn)錯誤！'); disp('2.對于樣本的輸入，如果頻數(shù)都為1時直接回車就可，不必輸入！'); disp('3.版權(quán)所有，仿冒必究！'); disp('0.退出說明進入主菜單！'); t=input(&

8、#39;'); if t=0 zhucaidan; else if t=0 exit; endendelse if c=0 zhucaidan; else if c=1|c=0 exit; end endend(Help1.m)：clc; disp(' -假設(shè)檢驗與置信區(qū)間-'); disp('*幫助與操作說明*'); disp('1.請嚴格按照說明與提示進行操作，否則可能會出現(xiàn)錯誤！'); disp('2.對于樣本的輸入，如果頻數(shù)都為1時直接回車就可，不必輸入！'); disp('3.版權(quán)所有，仿冒必究！'

9、;); disp('0.退出說明進入主菜單！'); t=input(''); if t=0 zhucaidan; else if t=0 exit; end end(Ujianyan.m):clc;disp(' 菜單選項');disp('0.返回主菜單');disp('1.返回上一層');disp('2.退出程序！');disp('3.方差已知，單個母體，u檢驗');disp('4.方差已知，兩個母體，t檢驗');disp('請進行選擇： ');c=i

10、nput('');if c=0 zhucaidan;else if c=1 jiashejianyan;else if c=2 exit; else if c=3clc;disp('u檢驗，方差已知，服從N（0,1）');disp(' 檢驗假設(shè)選擇');disp('0.返回主菜單');disp('1.返回上一層');disp('2.退出程序！');disp('3.H0:=0;H1:>0');disp('4.H0:=0;H1:<0');disp('5

11、.H0:=0;H1:0');disp('請進行選擇： ');c=input('');if c=0 zhucaidan;else if c=1 jiashejianyan;else if c=2 exit; else if c=3 clc; disp('H0:=0;H1:>0'); disp('請輸入子樣樣本，以行矩陣的形式輸入'); A=input(''); disp('請輸入各個子樣樣本的頻數(shù)，以行矩陣的形式輸入'); B=input(''); disp('請

12、輸入均值：'); a=input(''); disp('請輸入標準差：'); b=input(''); disp('請輸入顯著性水平：'); d=input(''); n=numel(A);sum=0;sun=0;o=numel(B); if n=o for i=1:n sum=sum+A(1,i)*B(1,i); end for j=1:n sun=sun+B(1,j); end u=(sum/sun-a)*sqrt(sun)/b else if o<n for i=1:n sum=sum+A(1,

13、i); end u=(sum/n-a)*sqrt(n)/b end end u1=norminv(1-d,0,1) if u>=u1 disp('*該假設(shè)檢驗拒絕H0，接受H1*'); else if u<u1 disp('*該假設(shè)檢驗拒絕H1，接受H0*'); end end else if c=4 clc; disp('H0:=0;H1:<0'); disp('請輸入子樣樣本，以行矩陣的形式輸入'); A=input(''); disp('請輸入各個子樣樣本的頻數(shù)，以行矩陣的形式輸入&

14、#39;); B=input(''); disp('請輸入均值：'); a=input(''); disp('請輸入標準差：'); b=input(''); disp('請輸入顯著性水平：'); d=input(''); n=numel(A);sum=0;sun=0;o=numel(B); if n=o for i=1:n sum=sum+A(1,i)*B(1,i); end for j=1:n sun=sun+B(1,j); end u=(sum/sun-a)*sqrt(sun)/

15、b else if o<n for i=1:n sum=sum+A(1,i); end u=(sum/n-a)*sqrt(n)/b end end u1=norminv(1-d,0,1) if u<=u1 disp('*該假設(shè)檢驗拒絕H0，接受H1*'); else if u>u1 disp('*該假設(shè)檢驗拒絕H1，接受H0*'); end endelse if c=5 clc; disp('H0:=0;H1:0'); disp('請輸入子樣樣本，以行矩陣的形式輸入'); A=input('');

16、 disp('請輸入各個子樣樣本的頻數(shù)，以行矩陣的形式輸入'); B=input(''); disp('請輸入均值：'); a=input(''); disp('請輸入標準差：'); b=input(''); disp('請輸入顯著性水平：'); d=input(''); n=numel(A);sum=0;sun=0;o=numel(B); if n=o for i=1:n sum=sum+A(1,i)*B(1,i); end for j=1:n sun=sun+B(

17、1,j); end u=(sum/sun-a)*sqrt(sun)/b else if o<n for i=1:n sum=sum+A(1,i); end u=(sum/n-a)*sqrt(n)/b end end u1=norminv(1-0.5*d,0,1); u1 if abs(u)>=u1 disp('*該假設(shè)檢驗拒絕H0，接受H1*'); else if abs(u)<u1 disp('*該假設(shè)檢驗拒絕H1，接受H0*'); end end end endendendendend else if c=4 clc;disp('u

18、檢驗，方差已知，服從N（0,1）');disp(' 檢驗假設(shè)選擇');disp('0.返回主菜單');disp('1.返回上一層');disp('2.退出程序！');disp('3.H0:1=2;H1:1>2');disp('4.H0:1=2;H1:1<2');disp('5.H0:1=2;H1:12');disp('請進行選擇： ');c=input('');if c=0 zhucaidan;else if c=1 jiashej

19、ianyan;else if c=2 exit; else if c=3 clc; disp('H0:1=2;H1:1>2'); disp('請輸入子樣樣本1'); A=input(''); disp('請輸入樣本1的頻數(shù)，以行矩陣的形式輸入'); C=input(''); disp('請輸入子樣樣本2'); B=input(''); disp('請輸入樣本2的頻數(shù)，以行矩陣的形式輸入'); D=input(''); disp('請輸入標

20、準差1：'); a=input(''); disp('請輸入標準差2：'); b=input(''); disp('請輸入顯著性水平：'); d=input(''); m=numel(A);n=numel(B);sum=0;sun=0;jiegua=0;jieguo=0;o=numel(C);z=numel(D); if m=o&&n=z for i=1:m sum=sum+A(1,i)*C(1,i); jieguo=jieguo+C(1,i); end for j=1:n sun=sun+

21、B(1,j)*D(1,j); jiegua=jiegua+D(1,j); end u=(sum/jieguo)-(sun/jiegua)/(sqrt(a2)/jieguo)+(b2)/jiegua) else if o<m&&n=z for i=1:m sum=sum+A(1,i); end for j=1:n sun=sun+B(1,j)*D(1,j); jiegua=jiegua+D(1,j); end u=(sum/m)-(sun/jiegua)/(sqrt(a2)/m)+(b2)/jiegua) else if o=m&&z<n for i=

22、1:m sum=sum+A(1,i)*C(1,i); jieguo=jieguo+C(1,i); end for j=1:n sun=sun+B(1,j); end u=(sum/jieguo)-(sun/n)/(sqrt(a2)/jieguo)+(b2)/n) else if o<m&&z<n for i=1:m sum=sum+A(1,i); end for j=1:n sun=sun+B(1,j); end u=(sum/m)-(sun/n)/(sqrt(a2)/m)+(b2)/n) end end end end u1=norminv(1-d,0,1) if

23、 u>=u1 disp('*該假設(shè)檢驗拒絕H0，接受H1*'); else if u<u1 disp('*該假設(shè)檢驗拒絕H1，接受H0*'); end end else if c=4 clc; disp('H0:1=2;H1:1<2'); disp('請輸入子樣樣本1'); A=input(''); disp('請輸入樣本1的頻數(shù)，以行矩陣的形式輸入'); C=input(''); disp('請輸入子樣樣本2'); B=input('

24、9;); disp('請輸入樣本2的頻數(shù)，以行矩陣的形式輸入'); D=input(''); disp('請輸入標準差1：'); a=input(''); disp('請輸入標準差2：'); b=input(''); disp('請輸入顯著性水平：'); d=input(''); m=numel(A);n=numel(B);sum=0;sun=0;jiegua=0;jieguo=0;o=numel(C);z=numel(D); if m=o&&n=z

25、for i=1:m sum=sum+A(1,i)*C(1,i); jieguo=jieguo+C(1,i); end for j=1:n sun=sun+B(1,j)*D(1,j); jiegua=jiegua+D(1,j); end u=(sum/jieguo)-(sun/jiegua)/(sqrt(a2)/jieguo)+(b2)/jiegua) else if o<m&&n=z for i=1:m sum=sum+A(1,i); end for j=1:n sun=sun+B(1,j)*D(1,j); jiegua=jiegua+D(1,j); end u=(sum

26、/m)-(sun/jiegua)/(sqrt(a2)/m)+(b2)/jiegua) else if o=m&&z<n for i=1:m sum=sum+A(1,i)*C(1,i); jieguo=jieguo+C(1,i); end for j=1:n sun=sun+B(1,j); end u=(sum/jieguo)-(sun/n)/(sqrt(a2)/jieguo)+(b2)/n) else if o<m&&z<n for i=1:m sum=sum+A(1,i); end for j=1:n sun=sun+B(1,j); end

27、u=(sum/m)-(sun/n)/(sqrt(a2)/m)+(b2)/n) end end end end u1=norminv(1-d,0,1) if u<=u1 disp('*該假設(shè)檢驗拒絕H0，接受H1*'); else if u>u1 disp('*該假設(shè)檢驗拒絕H1，接受H0*'); end endelse if c=5 clc; disp('H0:1=2;H1:12'); disp('請輸入子樣樣本1'); A=input(''); disp('請輸入樣本1的頻數(shù)，以行矩陣的形式輸

28、入'); C=input(''); disp('請輸入子樣樣本2'); B=input(''); disp('請輸入樣本2的頻數(shù)，以行矩陣的形式輸入'); D=input(''); disp('請輸入標準差1：'); a=input(''); disp('請輸入標準差2：'); b=input(''); disp('請輸入顯著性水平：'); d=input(''); m=numel(A);n=numel(B);s

29、um=0;sun=0;jiegua=0;jieguo=0;o=numel(C);z=numel(D); if m=o&&n=z for i=1:m sum=sum+A(1,i)*C(1,i); jieguo=jieguo+C(1,i); end for j=1:n sun=sun+B(1,j)*D(1,j); jiegua=jiegua+D(1,j); end u=(sum/jieguo)-(sun/jiegua)/(sqrt(a2)/jieguo)+(b2)/jiegua) else if o<m&&n=z for i=1:m sum=sum+A(1,i

30、); end for j=1:n sun=sun+B(1,j)*D(1,j); jiegua=jiegua+D(1,j); end u=(sum/m)-(sun/jiegua)/(sqrt(a2)/m)+(b2)/jiegua) else if o=m&&z<n for i=1:m sum=sum+A(1,i)*C(1,i); jieguo=jieguo+C(1,i); end for j=1:n sun=sun+B(1,j); end u=(sum/jieguo)-(sun/n)/(sqrt(a2)/jieguo)+(b2)/n) else if o<m&

31、&z<n for i=1:m sum=sum+A(1,i); end for j=1:n sun=sun+B(1,j); end u=(sum/m)-(sun/n)/(sqrt(a2)/m)+(b2)/n) end end end end u1=norminv(1-0.5*d,0,1) if abs(u)>=u1 disp('*該假設(shè)檢驗拒絕H0，接受H1*'); else if abs(u)<u1 disp('*該假設(shè)檢驗拒絕H1，接受H0*'); end end end endendendendend end end end end

32、enddisp('菜單選項');disp('1.繼續(xù)運算');disp('2.返回主菜單');disp('3.打開代碼');disp('4.退出程序！');p=input(' ');if p=1 (ujianyan); else if p=2 zhucaidan; else if p=3 open('ujianyan'); disp(' 菜單選項');disp('1.返回主菜單');disp('2.退出程序！');p=input(

33、9; '); if p=1 zhucaidan; else if p=2 disp('正在退出，請稍候。'); (exit);end endelse if p=4 disp('正在退出，請稍候。'); (exit);endendendend Ujianyan1.m ：clc;disp(' 菜單選項');disp('0.返回主菜單');disp('1.返回上一層');disp('2.退出程序！');disp('3.方差未知，單個母體，t檢驗');disp('4.方差未知

34、，兩個母體，t檢驗');disp('請進行選擇： ');c=input('');if c=0 zhucaidan;else if c=1 jiashejianyan;else if c=2 exit; else if c=3clc;disp('t檢驗，方差未知，服從t（n-1）');disp(' 檢驗假設(shè)選擇');disp('0.返回主菜單');disp('1.返回上一層');disp('2.退出程序！');disp('3.H0:=0;H1:>0,');

35、disp('4.H0:=0;H1:<0,');disp('5.H0:=0;H1:0,');disp('請進行選擇： ');c=input('');if c=0 zhucaidan;else if c=1 jiashejianyan;else if c=2 exit; else if c=3 clc; disp('方差未知，H0:=0;H1:>0'); disp('請輸入子樣樣本,以行向量的形式輸入'); A=input(''); disp('請輸入樣本的頻數(shù)，以行

36、矩陣的形式輸入'); C=input(''); disp('請輸入均值：'); a=input(''); disp('請輸入顯著性水平：'); b=input(''); n=numel(A);sum=0;jieguo=0;rong=0;o=numel(C); if n=0 for i=1:n sum=sum+A(1,i)*C(1,i); rong=rong+C(1,i); end sun=sum/rong; for j=1:n jieguo=jieguo+C(1,j)*(A(1,j)-sun)2; end

37、t=(sun-a)*sqrt(rong2-rong)/jieguo) t1=tinv(1-b,rong-1) else if o<n for i=1:n sum=sum+A(1,i); end sun=sum/n; for j=1:n jieguo=jieguo+(A(1,j)-sun)2; end t=(sun-a)*sqrt(n2-n)/jieguo) t1=tinv(1-b,n-1) end end if t>=t1 disp('該假設(shè)檢驗拒絕H0，接受H1'); else if t<t1 disp('該假設(shè)檢驗拒絕H1，接受H0');

38、end end else if c=4 clc; disp('方差未知，H0:=0;H1:<0'); disp('請輸入子樣樣本,以行向量的形式輸入'); A=input(''); disp('請輸入樣本的頻數(shù)，以行矩陣的形式輸入'); C=input(''); disp('請輸入均值：'); a=input(''); disp('請輸入顯著性水平：'); b=input(''); n=numel(A);sum=0;jieguo=0;rong=0

39、;o=numel(C); if n=0 for i=1:n sum=sum+A(1,i)*C(1,i); rong=rong+C(1,i); end sun=sum/rong; for j=1:n jieguo=jieguo+C(1,j)*(A(1,j)-sun)2; end t=(sun-a)*sqrt(rong2-rong)/jieguo) t1=tinv(1-b,rong-1) else if o<n for i=1:n sum=sum+A(1,i); end sun=sum/n; for j=1:n jieguo=jieguo+(A(1,j)-sun)2; end t=(sun-

40、a)*sqrt(n2-n)/jieguo) t1=tinv(1-b,n-1) end end if t<=t1 disp('該假設(shè)檢驗拒絕H0，接受H1'); else if t>t1 disp('該假設(shè)檢驗拒絕H1，接受H0'); end endelse if c=5 clc; disp('方差未知，H0:=0;H1:0'); disp('請輸入子樣樣本,以行向量的形式輸入'); A=input(''); disp('請輸入樣本的頻數(shù)，以行矩陣的形式輸入'); C=input('

41、;'); disp('請輸入均值：'); a=input(''); disp('請輸入顯著性水平：'); b=input(''); n=numel(A);sum=0;jieguo=0;rong=0;o=numel(C); if n=0 for i=1:n sum=sum+A(1,i)*C(1,i); rong=rong+C(1,i); end sun=sum/rong; for j=1:n jieguo=jieguo+C(1,j)*(A(1,j)-sun)2; end t=(sun-a)*sqrt(rong2-rong)/

42、jieguo) t1=tinv(1-0.5*b,rong-1) else if o<n for i=1:n sum=sum+A(1,i); end sun=sum/n; for j=1:n jieguo=jieguo+(A(1,j)-sun)2; end t=(sun-a)*sqrt(n2-n)/jieguo) t1=tinv(1-0.5*b,n-1) end end if abs(t)>=t1 disp('該假設(shè)檢驗拒絕H0，接受H1'); else if abs(t)<t1 disp('該假設(shè)檢驗拒絕H1，接受H0'); end end e

43、nd endendendendend else if c=4 clc;disp('t檢驗，方差未知');disp(' 檢驗假設(shè)選擇');disp('0.返回主菜單');disp('1.返回上一層');disp('2.退出程序！');disp('3.H0:1=2;H1:1>2');disp('4.H0:1=2;H1:1<2');disp('5.H0:1=2;H1:12');disp('請進行選擇： ');c=input('')

44、;if c=0 zhucaidan;else if c=1 jiashejianyan;else if c=2 exit; else if c=3 clc; disp('方差未知，H0:1=2;H1:1>2'); disp('請輸入子樣樣本1'); A=input(''); disp('請輸入樣本1的頻數(shù)，以行矩陣的形式輸入'); C=input(''); disp('請輸入子樣樣本2'); B=input(''); disp('請輸入樣本2的頻數(shù)，以行矩陣的形式輸入'); D=input(''); disp('請輸入顯著性水平：'); b=input(''); m=numel(A);n=numel(B);sum=0;j