为了正常的体验网站,请在浏览器设置里面开启Javascript功能!

统计计算源程序代码

2011-05-14 13页 doc 65KB 36阅读

用户头像

is_942467

暂无简介

举报
统计计算源程序代码////////////////////////////////////////////////////////////////////////////////////////////// //样本均值和样本方差的计算公式-------P10-------3.2-------- #include #include void main() { int n; cout>n; double *a=new double [n+1]; //记录样本值 double *x=new double [n+1]; //记录均值 double *...
统计计算源程序代码
////////////////////////////////////////////////////////////////////////////////////////////// //样本均值和样本方差的-------P10-------3.2-------- #include #include void main() { int n; cout<<"请输入样本值的个数n:"<>n; double *a=new double [n+1]; //记录样本值 double *x=new double [n+1]; //记录均值 double *s=new double [n+1]; //记录方差 cout<<"请输入样本值:"<>a[i]; } //初始化 x[1]=a[1]; s[1]=0; //开始计算其余的均值和方差 for (int k=1;k #include const int n=10,p=4; //定义样本值的个数为n=10,矩阵为4阶矩阵 void main() { double m[n+1][p]; //记录样本值 double x[n+1][p]; //记录均值向量 double A[n+1][p][p]; //记录样本离差阵 int i,j; cout<<"请输入"<>m[i][j]; } } //初始化 for(j=0;j
数----------P57--4.1--- #include #include double Calculate(int i,double x) //定义一个函数用来计算后面的级数展开式 { double px=x,pa=1,pb=1; for (int j=1;j>x>> precision; while(IsGoOn) //开始求级数展开式的值 { sumall+=Calculate(i,x); //判断精度是否达到要求 if(Calculate(i+1,x)-Calculate(i,x)< precision) { IsGoOn=false; } i++; } erf=2*exp(-x*x)*sumall/sqrt(3.1415926); //最后求函数的值 cout< #include void main() { //已知系数数组 double a[]={0.0705230784,0.0422820123,0.0092705272, 0.0001520143,0.0002765672,0.0000430638}; double x; //result用来记录中间变量,erf用来记录最后函数的值 double result=0,erf; cin>>x; //输入初始值 for(int i=0;i<6;i++) { result+=a[i]*pow(x,i+1); //调用数学函数求n次方 } result+=1; erf=1-pow(result,-16); cout< #include void main() { double x,y,u,result=0; double a[11]={...}; //此处没有把全部的值给写出来 cin>>x; y=-log(4*x*(1-x)); for(int i=0;i<11;i++) { result+=a[i]*pow(y,i); //调用数学函数 } result*=y; u=sqrt(result); cout< #include double I(double a,double b,double x); //声明两个叠代函数 double U(double a,double b,double x); //这个叠代函数只能进行初始值为(0.5,0.5),(0.5,1),(1,0.5),(1,1)的计算,下同 double I(double a,double b,double x) //函数的定义部分 { double pa=a,pb=b,px=x; if(pa==0.5 && pb==0.5) //根据公式写出各种情况的返回值 { return (1-2/(3.1415926*tan(sqrt((1-px)/px)))); } else if(pa==0.5 && pb==1) { return sqrt(px); } else if(pa==1 && pb==0.5) { return (1-sqrt(1-px)); } else if(pa==1 && pb==1) { return px; } else { if(pa>1) return I(pa-1,pb,px)-U(pa-1,pb,px)/(pa-1); //不断分解,实现递归调用 if(pb>1) return I(pa,pb-1,px)+U(pa,pb-1,px)/(pb-1); } } double U(double a,double b,double x) { double pa=a,pb=b,px=x; if(pa==0.5 && pb==0.5) { return (sqrt(px*(1-px))/3.1415926); } else if(pa==0.5 && pb==1) { return sqrt(px)*(1-px)/2; } else if(pa==1 && pb==0.5) { return px*sqrt(1-px)/2; } else if(pa==1 && pb==1) { return px*(1-px); } else { if(pa>1) return (pa-1+pb)*px*U(pa-1,pb,px)/(pa-1); if(pb>1) return (pa+pb-1)*(1-px)*U(pa,pb-1,px)/(pb-1); } } void main() { double a,b,x; cout<<"请输入所求I(a,b)中的a,b的值(请注意:a,b的值只能是1/2的整数倍!):"<>a>>b; if(fmod(a*2,1)!=0 || fmod(b*2,1)!=0) //出错控制 { cout<<"a,b的值不是1/2的整数倍!"<>x; if(x>1 || x<0) { cout<<"X不在0到1之间!"<>a0>>b0>>c0>>d0; cout<<"所求方程的解的区间(a,b)中a,b各自的值:"<>a>>b; cout<<"请输入计算结果要达到的精度percision:"<>percision; if(f(a)*f(b)>0) { cout<<"区间不符合要求,方程在此区间内无解!请重新确定区间!"<>n; double *X=new double[n]; cout<<"请输入初始值:"<>X[0]; for(int i=1;i #include void main() { double M=pow(2,31)-1; double a=pow(7,5); int b=M/a; int c=fmod(M,a); int n; cout<<"请输入要求的次数n"<>n; double *X=new double[n]; //动态定义一个一维数组 cout<<"请输入初始值:"<>X[0]; int k0,k1; for(int i=1;i #include const int n=3; //定义全局变量,表示数组的维数 double a[n][n]; //定义一个二维数组来表示矩阵 int L[n]; //记录主元的列标 //此函数用来对矩阵A施行以(i,j)为主元的消去变换,即定义6.1的实现 void sweep(int i,int j) { double b[n][n]; for(int t=0;t>a[i][j]; } } for(i=0;i
/
本文档为【统计计算源程序代码】,请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑, 图片更改请在作品中右键图片并更换,文字修改请直接点击文字进行修改,也可以新增和删除文档中的内容。
[版权声明] 本站所有资料为用户分享产生,若发现您的权利被侵害,请联系客服邮件isharekefu@iask.cn,我们尽快处理。 本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用。 网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。

历史搜索

    清空历史搜索