置信区间与假设检验matlab程序

版权所有严禁侵犯不得用于商业用途统计学专用程序---基于MATLAB7.0开发---置信区间与假设检验2013年8月1日版权所有严禁侵犯不得用于商业用途第1页共60页置信区间与假设检验程序【开发目的】众所周知，统计工作面对的数据量繁琐而且庞大，在统计的过程中和计算中容易出错，并统计决定着国民经济的命脉，开发此软件就是为了进行验证统计的准确性以及理论可行性，减少统计工作中的错误，使统计工作者更好地进行工作与学习；所以特意开发此程序来检验统计中的参数估计和假设检验。【开发特色】本软件基于matlab7.0进行运算，对于样本的输入采用行矩阵的形式，并且开发了样本频数输入，对于多样本的输入可以减缓工作量，对于显著性水平本程序自带正态分布函数，t分布函数，F分布函数，2分布函数的计算公式，不用再为查表和计算而苦恼，只需输入显著性水平即可，大大的简化了计算量。【关键技术】矩阵输入进行频数判断条件循环语句的使用等【程序界面】版权所有严禁侵犯不得用于商业用途第2页共60页版权所有严禁侵犯不得用于商业用途第3页共60页【程序代码】此程序采用多文件结构，在建立文件时不能改变文件名；以下是各个文件的代码：（Zhucaidan.m)：clc;disp('统计学专用');disp('1.假设检验');disp('2.置信区间');disp('3.使用说明');disp('4.打开代码');disp('0.退出程序');disp('请进行选择：');a=input('');ifa==0exit;elseifa==1jiashejianyan;elseifa==2zhixinqujian;elseifa==3help1;elseifa==4open('zhucaidan');disp('菜单选项');disp('1.返回主菜单');disp('2.退出程序！');p=input('');ifp==1zhucaidan;elseifp==2disp('正在退出，请稍候。。。');(exit);endendendendendendend(Zhixinqujian.m)：clc;disp('置信区间');disp('菜单选项');disp('0.退出程序!!!');disp('1.返回主菜单');disp('2.方差已知，待估参数为u');disp('3.方差未知，待估参数为u');disp('4.均值已知，待估参数为方差');disp('5.均值未知，待估参数为方差');disp('请进行选择：');版权所有严禁侵犯不得用于商业用途第4页共60页c=input('');ifc==0exit;elseifc==1zhucaidan;elseifc==2uguji;elseifc==3uguji1;elseifc==4fguji;elseifc==5fguji1;endendendendendend(Jiashejianyan.m)：clc;disp('*********假设检验**********');disp('*********菜单选项**********');disp('********0.退出程序!!!******');disp('********1.返回主菜单*******');disp('***2.σ已知，对u进行检验***');disp('***3.σ未知，对u进行检验***');disp('***4.u已知，对σ进行检验***');disp('***5.u未知，对σ进行检验***');disp('******请进行选择：');c=input('');ifc==0exit;elseifc==1zhucaidan;elseifc==2ujianyan;elseifc==3ujianyan1;elseifc==4fangcha;elseifc==5fangcha1;endendendendendendHelp.m：版权所有严禁侵犯不得用于商业用途第5页共60页clc;disp('--假设检验与置信区间--');disp('是否使用帮助（1/0)');c=input('');ifc==1clc;disp('--假设检验与置信区间--');disp('********************帮助与操作说明*************');disp('1.请严格按照说明与提示进行操作，否则可能会出现错误！！');disp('2.对于样本的输入，如果频数都为1时直接回车就可，不必输入！！');disp('3.版权所有，仿冒必究！！');disp('0.退出说明进入主菜单！！');t=input('');ift==0zhucaidan;elseift~=0exit;endendelseifc==0zhucaidan;elseifc~=1|c~=0exit;endendend(Help1.m)：clc;disp('--假设检验与置信区间--');disp('********************帮助与操作说明*************');disp('1.请严格按照说明与提示进行操作，否则可能会出现错误！！');disp('2.对于样本的输入，如果频数都为1时直接回车就可，不必输入！！');disp('3.版权所有，仿冒必究！！');disp('0.退出说明进入主菜单！！');t=input('');ift==0zhucaidan;elseift~=0exit;endend(Ujianyan.m):clc;disp('菜单选项');disp('0.返回主菜单');disp('1.返回上一层');disp('2.退出程序！！！');disp('3.方差已知，单个母体，u检验');版权所有严禁侵犯不得用于商业用途第6页共60页disp('4.方差已知，两个母体，t检验');disp('请进行选择：');c=input('');ifc==0zhucaidan;elseifc==1jiashejianyan;elseifc==2exit;elseifc==3clc;disp('u检验，方差已知，服从N（0,1）');disp('检验假设选择');disp('0.返回主菜单');disp('1.返回上一层');disp('2.退出程序！！！');disp('3.H0:μ=μ0;H1:μμ0');disp('4.H0:μ=μ0;H1:μμ0');disp('5.H0:μ=μ0;H1:μ≠μ0');disp('请进行选择：');c=input('');ifc==0zhucaidan;elseifc==1jiashejianyan;elseifc==2exit;elseifc==3clc;disp('H0:μ=μ0;H1:μμ0');disp('请输入子样样本，以行矩阵的形式输入');A=input('');disp('请输入各个子样样本的频数，以行矩阵的形式输入');B=input('');disp('请输入均值：');a=input('');disp('请输入标准差：');b=input('');disp('请输入显著性水平α：');d=input('');n=numel(A);sum=0;sun=0;o=numel(B);ifn==ofori=1:nsum=sum+A(1,i)*B(1,i);endforj=1:nsun=sun+B(1,j);endu=(sum/sun-a)*sqrt(sun)/belseifonfori=1:n版权所有严禁侵犯不得用于商业用途第7页共60页sum=sum+A(1,i);endu=(sum/n-a)*sqrt(n)/bendendu1=norminv(1-d,0,1)ifu=u1disp('****该假设检验拒绝H0，接受H1****');elseifuu1disp('****该假设检验拒绝H1，接受H0****');endendelseifc==4clc;disp('H0:μ=μ0;H1:μμ0');disp('请输入子样样本，以行矩阵的形式输入');A=input('');disp('请输入各个子样样本的频数，以行矩阵的形式输入');B=input('');disp('请输入均值：');a=input('');disp('请输入标准差：');b=input('');disp('请输入显著性水平α：');d=input('');n=numel(A);sum=0;sun=0;o=numel(B);ifn==ofori=1:nsum=sum+A(1,i)*B(1,i);endforj=1:nsun=sun+B(1,j);endu=(sum/sun-a)*sqrt(sun)/belseifonfori=1:nsum=sum+A(1,i);endu=(sum/n-a)*sqrt(n)/bendendu1=norminv(1-d,0,1)ifu=u1disp('****该假设检验拒绝H0，接受H1****');elseifuu1disp('****该假设检验拒绝H1，接受H0****');endendelseifc==5clc;disp('H0:μ=μ0;H1:μ≠μ0');disp('请输入子样样本，以行矩阵的形式输入');版权所有严禁侵犯不得用于商业用途第8页共60页A=input('');disp('请输入各个子样样本的频数，以行矩阵的形式输入');B=input('');disp('请输入均值：');a=input('');disp('请输入标准差：');b=input('');disp('请输入显著性水平α：');d=input('');n=numel(A);sum=0;sun=0;o=numel(B);ifn==ofori=1:nsum=sum+A(1,i)*B(1,i);endforj=1:nsun=sun+B(1,j);endu=(sum/sun-a)*sqrt(sun)/belseifonfori=1:nsum=sum+A(1,i);endu=(sum/n-a)*sqrt(n)/bendendu1=norminv(1-0.5*d,0,1);u1ifabs(u)=u1disp('****该假设检验拒绝H0，接受H1****');elseifabs(u)u1disp('****该假设检验拒绝H1，接受H0****');endendendendendendendendelseifc==4clc;disp('u检验，方差已知，服从N（0,1）');disp('检验假设选择');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:μ1≠μ2');disp('请进行选择：');c=input('');版权所有严禁侵犯不得用于商业用途第9页共60页ifc==0zhucaidan;elseifc==1jiashejianyan;elseifc==2exit;elseifc==3clc;disp('H0:μ1=μ2;H1:μ1μ2');disp('请输入子样样本1');A=input('');disp('请输入样本1的频数，以行矩阵的形式输入');C=input('');disp('请输入子样样本2');B=input('');disp('请输入样本2的频数，以行矩阵的形式输入');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);ifm==o&&n==zfori=1:msum=sum+A(1,i)*C(1,i