统计计算R软件课程设计

发布 2022-10-03 19:28:28 阅读 5846

2016 年 7 月 8 日。

一. 课程设计目的。

统计计算课程设计是在学习了《统计计算》课程之后,进行此课程设计,是对这门课程的全面复习,也是这门课程理论知识的实践,是整个教学工作的重要环节。

通过统计计算课程设计教学所要达到的目的是:以统计计算课程和理论知识为基础,通过课程设计的实践,加强学生对所学相关课程的理解、掌握,训练并提高r软件的使用、统计蒙特卡洛方法、独立解决问题的能力。

二. 设计内容。

1. 运用r软件完成《习题一》1.5题。

2. 编制梯形求积公式和抛物线求积公式计算的程序。

3. 编制二分法和牛顿法的求解程序。

4. 用直接抽样法产生[1,10]上的均匀随机数和指数分布的随机数。

5. 用二维变换抽样法产生正态分布的随机数。

6. 用gfi法产生gamma分布的随机数。

7. 用舍选抽样法产生beta分布的随机数。

8. 用复合抽样法产生密度函数为的随机数。

9. 用随机投针试验方法求π的近似值。

三. 程序及结果。

1.题目:运用r软件完成《习题一》1.5题。

程序**。x1<-c(23,20,18,29,43,35,32,40,29,26,24,26);/x2<-c(1000,1000,500,500,500,100,100,100,100,100,100,100,100,100,100);/x3<-c(70,70,70,70,10);

mean(x1);#均值。

var(x1);#方差。

sd(x1);#标准差。

median(x1);#中位数

which(table(x1)==max(table(x1)))#众数。

cv=sd(x1)/mean(x1);#变异系数。

print(cv);

n=length(x1);

g1=sqrt(1/6*n)*sum(((x1-mean(x1))/sd(x1))^3);#标准偏度系数。

print(g1);

g2=sqrt(n/24)*(1/n)*sum(((x1-mean(x1))/sd(x1))^4)-3);#标准峰度系数,sqrt为平方根函数。

print(g2) ;

结果。1)x1<-c(23,20,18,29,43,35,32,40,29,26,24,26)

mean(x1);#均值。

var(x1);#方差。

sd(x1);#标准差。

median(x1);#中位数。

which(table(x1)==max(table(x1)))#众数。

cv=sd(x1)/mean(x1);#变异系数。

print(cv);

n=length(x1);

g1=sqrt(1/6*n)*sum(((x1-mean(x1))/sd(x1))^3);#标准偏度系数。

print(g1);

g2=sqrt(n/24)*(1/n)*sum(((x1-mean(x1))/sd(x1))^4)-3);#标准峰度系数。

print(g2) ;

2)x2<-c(1000,1000,500,500,500,100,100,100,100,100,100,100,100,100,100)

mean(x1);#均值。

var(x1);#方差。

sd(x1);#标准差。

median(x1);#中位数。

which(table(x1)==max(table(x1)))#众数。

cv=sd(x1)/mean(x1);#变异系数。

print(cv);

n=length(x1);

g1=sqrt(1/6*n)*sum(((x1-mean(x1))/sd(x1))^3);#标准偏度系数。

print(g1);

g2=sqrt(n/24)*(1/n)*sum(((x1-mean(x1))/sd(x1))^4)-3);#标准峰度系数。

print(g2) ;

3)x3<-c(70,70,70,70,10)

mean(x1);#均值。

var(x1);#方差。

sd(x1);#标准差。

median(x1);#中位数。

which(table(x1)==max(table(x1)))#众数。

cv=sd(x1)/mean(x1);#变异系数。

print(cv);

n=length(x1);

g1=sqrt(1/6*n)*sum(((x1-mean(x1))/sd(x1))^3);#标准偏度系数。

print(g1);

g2=sqrt(n/24)*(1/n)*sum(((x1-mean(x1))/sd(x1))^4)-3);#标准峰度系数。

print(g2) ;

2.题目:编制梯形求积公式计算的程序。

程序**。integral=function(x) #abs 返回整形数据的绝对值。

print(n[1:i])

结果: print(a[1:i])

5. 题目:编制二分法的求解程序。

程序**:ef=function(x)

if(fb*fc<0)

if(abs(a[i]-b[i])<1*exp(-5))

print(c(i,a[i+1],b[i+1],c[i]))

c[i+1]=(a[i+1]+b[i+1])/2;

print(c[i+1])结果:

c[i+1]=(a[i+1]+b[i+1])/2;

print(c[i+1])

6. 题目:用直接抽样法产生[1,10]上的均匀随机数的随机数。

程序**:junyun=function(n,a,b){

统计计算课程设计作业

班级 数学121班。姓名 林孝涛。学号 1920125079 2015年1月13日星期二。作业一。必做题一。步骤 1.打开word文件和 准备填入数据。2.首先填序号列,这时候要注意单元格的格式,所以设置单元格,要设置成文本格式如图。3.填入001,然后下拉自动填充序列 4.紧接着填第二列 书名。4...

统计计算课程设计作业

报告书。实验一 统计学中有关概念。1 基本要求 了解统计学中常用的一些基本概念,重点掌握几个集中指标和变异指标。2 完成时间 2012 06 11 3 实验内容 2 根据下表求平均数 区域 四分位数 平均绝对离差 标准差 方差。操作过程 1 建立excel,把问题中的数据输入,下面计算所求值。2 在...

软件学院课程设计

2011年软件学院c 课程设计。课程设计目的 1 熟悉利用面向对象的方法以及c 的编程思想来完成系统的设计 2 锻炼学生在设计的过程中,建立清晰的类层次,应用继承和多态等面向对象的编程思想 3 通过本课程设计,加深对面向对象程序设计课程所学知识的理解,熟练掌握和巩固c 语言的基本知识和语法规范,深刻...