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

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

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

3、 else if a=3 help1; else if a=4 open(zhucaidan);disp( 菜單選項(xiàng));disp(1.返回主菜單);disp(2.退出程序！);p=input( ); if p=1 zhucaidan; else if p=2 disp(正在退出，請(qǐng)稍候。); (exit);end end end end end endend(Zhixinqujian.m) ：clc;disp( 置信區(qū)間);disp( 菜單選項(xiàng));disp(0.退出程序!);disp(1.返回主菜單);disp(2.方差已知，待估參數(shù)為u);disp(3.方差未知，待估參數(shù)為u);disp(4

4、.均值已知，待估參數(shù)為方差);disp(5.均值未知，待估參數(shù)為方差);disp(請(qǐng)進(jìn)行選擇：);c=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è)檢驗(yàn)*);disp(*菜單選項(xiàng)*);disp(*0.退出程序!*);disp(*1.返回主菜單*);disp(*2.已知，對(duì)u進(jìn)行檢驗(yàn)*);disp(*

5、3.未知，對(duì)u進(jìn)行檢驗(yàn)*);disp(*4.u已知，對(duì)進(jìn)行檢驗(yàn)*);disp(*5.u未知，對(duì)進(jìn)行檢驗(yàn)*);disp(*請(qǐng)進(jìn)行選擇：);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(-假設(shè)檢驗(yàn)與置信區(qū)間-);disp(是否使用幫助（1/0);c=input();if c=1 clc; disp( -假設(shè)

6、檢驗(yàn)與置信區(qū)間-); disp(*幫助與操作說明*); disp(1.請(qǐng)嚴(yán)格按照說明與提示進(jìn)行操作，否則可能會(huì)出現(xiàn)錯(cuò)誤！); disp(2.對(duì)于樣本的輸入，如果頻數(shù)都為1時(shí)直接回車就可，不必輸入！); disp(3.版權(quán)所有，仿冒必究！); disp(0.退出說明進(jìn)入主菜單！); t=input(); 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è)檢驗(yàn)與置信區(qū)間-); disp(*幫助與操作說明*);

7、 disp(1.請(qǐng)嚴(yán)格按照說明與提示進(jìn)行操作，否則可能會(huì)出現(xiàn)錯(cuò)誤！); disp(2.對(duì)于樣本的輸入，如果頻數(shù)都為1時(shí)直接回車就可，不必輸入！); disp(3.版權(quán)所有，仿冒必究！); disp(0.退出說明進(jìn)入主菜單！); t=input(); if t=0 zhucaidan; else if t=0 exit; end end(Ujianyan.m):clc;disp( 菜單選項(xiàng));disp(0.返回主菜單);disp(1.返回上一層);disp(2.退出程序！);disp(3.方差已知，單個(gè)母體，u檢驗(yàn));disp(4.方差已知，兩個(gè)母體，t檢驗(yàn));disp(請(qǐng)進(jìn)行選擇： );c=i

8、nput();if c=0 zhucaidan;else if c=1 jiashejianyan;else if c=2 exit; else if c=3clc;disp(u檢驗(yàn)，方差已知，服從N（0,1）);disp( 檢驗(yàn)假設(shè)選擇);disp(0.返回主菜單);disp(1.返回上一層);disp(2.退出程序！);disp(3.H0:=0;H1:0);disp(4.H0:=0;H1:0); disp(請(qǐng)輸入子樣樣本，以行矩陣的形式輸入); A=input(); disp(請(qǐng)輸入各個(gè)子樣樣本的頻數(shù)，以行矩陣的形式輸入); B=input(); disp(請(qǐng)輸入均值：); a=input

9、(); disp(請(qǐng)輸入標(biāo)準(zhǔn)差：); b=input(); disp(請(qǐng)輸入顯著性水平：); 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=u1 disp(*該假設(shè)檢驗(yàn)拒絕H0，接受H1*); else if uu1 disp(*該假設(shè)檢驗(yàn)拒絕H1，接受H0*); end end else if c=4 clc; disp(H0:=

10、0;H1:0); disp(請(qǐng)輸入子樣樣本，以行矩陣的形式輸入); A=input(); disp(請(qǐng)輸入各個(gè)子樣樣本的頻數(shù)，以行矩陣的形式輸入); B=input(); disp(請(qǐng)輸入均值：); a=input(); disp(請(qǐng)輸入標(biāo)準(zhǔn)差：); b=input(); disp(請(qǐng)輸入顯著性水平：); 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

11、else if on 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 uu1 disp(*該假設(shè)檢驗(yàn)拒絕H1，接受H0*); end endelse if c=5 clc; disp(H0:=0;H1:0); disp(請(qǐng)輸入子樣樣本，以行矩陣的形式輸入); A=input(); disp(請(qǐng)輸入各個(gè)子樣樣本的頻數(shù)，以行矩陣的形式輸入); B=input(); disp(請(qǐng)輸入均值：); a=input(); disp(請(qǐng)輸入標(biāo)準(zhǔn)差：); b=input(); disp(請(qǐng)輸入

12、顯著性水平：); 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=u1 disp(*該假設(shè)檢驗(yàn)拒絕H0，接受H1*); else if abs(u)2);disp(4.H0:1=2;H1:12); disp(請(qǐng)輸入子樣樣本1); A=input(); disp(請(qǐng)輸入樣本1的頻數(shù)，以行矩陣的形式輸入); C=input(); disp

13、(請(qǐng)輸入子樣樣本2); B=input(); disp(請(qǐng)輸入樣本2的頻數(shù)，以行矩陣的形式輸入); D=input(); disp(請(qǐng)輸入標(biāo)準(zhǔn)差1：); a=input(); disp(請(qǐng)輸入標(biāo)準(zhǔn)差2：); b=input(); disp(請(qǐng)輸入顯著性水平：); 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 su

14、n=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 om&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&zn for i=1:m sum=sum+A(1,i)*C

15、(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 om&z=u1 disp(*該假設(shè)檢驗(yàn)拒絕H0，接受H1*); else if uu1 disp(*該假設(shè)檢驗(yàn)拒絕H1，接受H0*); end end else if c=4 clc; disp(H0:1=2;H1:12); disp(請(qǐng)輸入子樣樣本1); A=input(); disp(請(qǐng)輸入樣本1的頻數(shù)，以行矩陣的形式輸入); C=input(); di

16、sp(請(qǐng)輸入子樣樣本2); B=input(); disp(請(qǐng)輸入樣本2的頻數(shù)，以行矩陣的形式輸入); D=input(); disp(請(qǐng)輸入標(biāo)準(zhǔn)差1：); a=input(); disp(請(qǐng)輸入標(biāo)準(zhǔn)差2：); b=input(); disp(請(qǐng)輸入顯著性水平：); 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

17、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 om&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&zn for i=1:m sum=sum+A(1,i)

18、*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 om&zn 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 uu1 disp(*該假設(shè)檢驗(yàn)拒絕H1，接受H0*); end en

19、delse if c=5 clc; disp(H0:1=2;H1:12); disp(請(qǐng)輸入子樣樣本1); A=input(); disp(請(qǐng)輸入樣本1的頻數(shù)，以行矩陣的形式輸入); C=input(); disp(請(qǐng)輸入子樣樣本2); B=input(); disp(請(qǐng)輸入樣本2的頻數(shù)，以行矩陣的形式輸入); D=input(); disp(請(qǐng)輸入標(biāo)準(zhǔn)差1：); a=input(); disp(請(qǐng)輸入標(biāo)準(zhǔn)差2：); b=input(); disp(請(qǐng)輸入顯著性水平：); d=input(); m=numel(A);n=numel(B);sum=0;sun=0;jiegua=0;jieguo=

20、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 om&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=ji

21、egua+D(1,j); end u=(sum/m)-(sun/jiegua)/(sqrt(a2)/m)+(b2)/jiegua) else if o=m&zn 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 om&z=u1 disp(*該假設(shè)檢驗(yàn)拒絕H0，接受H1*); else if abs(u)0,);disp(4.H0:=0;H1:0); dis

22、p(請(qǐng)輸入子樣樣本,以行向量的形式輸入); A=input(); disp(請(qǐng)輸入樣本的頻數(shù)，以行矩陣的形式輸入); C=input(); disp(請(qǐng)輸入均值：); a=input(); disp(請(qǐng)輸入顯著性水平：); 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

23、-a)*sqrt(rong2-rong)/jieguo) t1=tinv(1-b,rong-1) else if o=t1 disp(該假設(shè)檢驗(yàn)拒絕H0，接受H1); else if tt1 disp(該假設(shè)檢驗(yàn)拒絕H1，接受H0); end end else if c=4 clc; disp(方差未知，H0:=0;H1:0); disp(請(qǐng)輸入子樣樣本,以行向量的形式輸入); A=input(); disp(請(qǐng)輸入樣本的頻數(shù)，以行矩陣的形式輸入); C=input(); disp(請(qǐng)輸入均值：); a=input(); disp(請(qǐng)輸入顯著性水平：); b=input(); n=numel(

24、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)/jieguo) t1=tinv(1-b,rong-1) else if on for i=1:n sum=sum+A(1,i); end sun=sum/n; for j=1:n jieguo=jieguo+(A(1,j

25、)-sun)2; end t=(sun-a)*sqrt(n2-n)/jieguo) t1=tinv(1-b,n-1) end end if tt1 disp(該假設(shè)檢驗(yàn)拒絕H1，接受H0); end endelse if c=5 clc; disp(方差未知，H0:=0;H1:0); disp(請(qǐng)輸入子樣樣本,以行向量的形式輸入); A=input(); disp(請(qǐng)輸入樣本的頻數(shù)，以行矩陣的形式輸入); C=input(); disp(請(qǐng)輸入均值：); a=input(); disp(請(qǐng)輸入顯著性水平：); b=input(); n=numel(A);sum=0;jieguo=0;rong=

26、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)/jieguo) t1=tinv(1-0.5*b,rong-1) else if o=t1 disp(該假設(shè)檢驗(yàn)拒絕H0，接受H1); else if abs(t)2);disp(4.H0:1=2;H1:12); disp(請(qǐng)輸入子樣樣本1); A=input();

27、disp(請(qǐng)輸入樣本1的頻數(shù)，以行矩陣的形式輸入); C=input(); disp(請(qǐng)輸入子樣樣本2); B=input(); disp(請(qǐng)輸入樣本2的頻數(shù)，以行矩陣的形式輸入); D=input(); disp(請(qǐng)輸入顯著性水平：); b=input(); m=numel(A);n=numel(B);sum=0;jieguo=0;sun=0;jiegua=0;rong=0;rong1=0;o=numel(C);z=numel(D); if o=m&z=n for i=1:m sum=sum+A(1,i)*C(1,i); rong=rong+C(1,i); end sump=sum/rong; for j=1:n jieguo=jieguo+B(1,j)*D(1,j); rong1=rong1+D(1,j); end jieguop=jieguo/rong1; for k=1:m sun=s