本科学生综合性实验
报告
软件系统测试报告下载sgs报告如何下载关于路面塌陷情况报告535n,sgs报告怎么下载竣工报告下载
项目组长：杨 维
学 号：0113239
成 员： 杨维
专 业：电子信息工程
班 级：电信111班
实验项目名称：统计方法回归
分析
定性数据统计分析pdf销售业绩分析模板建筑结构震害分析销售进度分析表京东商城竞争战略分析
建模实验
指导教师及职称：党建武（教授）
开课学期 2013至 2014 学年 第一 学期
上课时间 2013 年 11 月 26 日
一、实验
设计
领导形象设计圆作业设计ao工艺污水处理厂设计附属工程施工组织设计清扫机器人结构设计
方案
气瓶 现场处置方案 .pdf气瓶 现场处置方案 .doc见习基地管理方案.doc关于群访事件的化解方案建筑工地扬尘治理专项方案下载
实验名称：统计方法回归分析建模实验
实验时间：2013.11.26
小组合作：是○ 否●
小组成员：无
实验目的及要求：
学习统计方法回归分析的思想和基本原理；掌握建立回归模型的基本步骤，明确回归分析的主要任务；熟悉MATLAB软件进行回归模型的各种统计分析；通过范例学习，熟悉统计分析思想和建立回归模型的基本要素。
通过该实验的学习，使学生掌握回归分析的统计思想，认识面对什么样的实际问题可以建立回归模型，并且对回归模型作统计分析，同时使学生学会使用MATLAB软件进行回归分析和计算的基本命令，了解统计软件的功能和作用；熟悉处理大量数据的要领和方法。
2.实验场地及仪器、设备和材料：
实验场地：H123
仪器、设备和材料：计算机，MATLAB软件
3.实验学时数与实验类型
3个学时，综合性实验
4. 实验思路（实验
内容
财务内部控制制度的内容财务内部控制制度的内容人员招聘与配置的内容项目成本控制的内容消防安全演练内容
、数据处理方法及实验步骤等）：
a.实验内容：
1．线性回归模型的建立与分析步骤（问题假设→模型→参数估计→模型检验→确定最优回归方程→预测）；
2．非线性回归模型的建立与分析步骤；
3．使用MATLAB命令对回归模型进行计算与分析（包括模型检验与预测）；
4．利用某些数值与图形对统计特征作定性分析。
b.数据处理方法:
通过matlab软件的使用，实现一些线性回归方程问题求解的功能，得出数据及分析结果
c.实验步骤：
1．开启软件平台——MATLAB，开启MATLAB编辑窗口；
2．根据微分方程求解步骤编写M文件
3．保存文件并运行；
4．观察运行结果(数值或图形)；
5．根据观察到的结果和体会写出实验报告。
d实验要求与任务：
1.某校60名学生的一次考试成绩如下：
93 75 83 93 91 85 84 82 77 76
77 95 94 89 91 88 86 83 96 81
79 97 78 75 67 69 68 84 83 81
75 66 85 70 94 84 83 82 80 78
74 73 76 70 86 76 90 89 71 66
86 73 80 94 79 78 77 63 53 55
（1） 计算均值、标准差、极差、偏度、峰度，画出直方图；
（2） 检验分布的正态性；
（3） 若检验符合正态分布，估计正态分布的参数并检验参数。
2．混凝土的抗压强度随养护时间的延长而增加，现将一批混凝土作成12个试块，记录了养护日期x（日）及抗压强度y（kg/cm2）的数据：
养护时间x
2
3
4
5
7
9
12
14
17
21
28
56
抗压强度y
35
42
47
53
59
65
68
73
76
82
86
99
试求
型回归方程。
3.电影院调查电视广告费用和报纸广告费用对每周收入的影响，得到下面的数据，试建立回归模型，并进行检验（写出模型检验的依据），并预测电视广告费用为1，报纸广告费用为6时的周收入（写出预测的程序指令）。
每周收入
96
90
95
92
95
95
94
94
电视广告费
1.5
2.0
1.5
2.5
3.3
2.3
4.2
2.5
报纸广告费
5.0
2.0
4.0
2.5
3.0
3.5
2.5
3.0
二、实验结果与分析
1、实验目的、场地及仪器、设备和材料、实验思路等见实验设计方案
2、实验现象、数据及结果
第一题：某校60名学生的一次考试成绩如下：
93 75 83 93 91 85 84 82 77 76
77 95 94 89 91 88 86 83 96 81
79 97 78 75 67 69 68 84 83 81
75 66 85 70 94 84 83 82 80 78
74 73 76 70 86 76 90 89 71 66
86 73 80 94 79 78 77 63 53 55
1) 计算均值、标准差、极差、偏度、峰度，画出直方图；
2) 检验分布的正态性；
3) 若检验符合正态分布，估计正态分布的参数并检验参数。
解：输入数据
s1=[93 75 83 93 91 85 84 82 77 76];
s2=[77 95 94 89 91 88 86 83 96 81];
s3=[79 97 78 75 67 69 68 84 83 81];
s4=[75 66 85 70 94 84 83 82 80 78];
s5=[74 73 76 70 86 76 90 89 71 66];
s6=[86 73 80 94 79 78 77 63 53 55];
s=[s1 s2 s3 s4 s5 s6];
save data s
保存在当前工作目录下
调用load data
1）程序
mean(s) %均值
std(s) %标准差
d=max(s)-min(s)%极差
skewness(s)%偏度
kurtosis(s)%峰度
hist(s,10)%直方图
结果：
ans =
80.1000
ans =
9.7106
d =
44
ans =
-0.4682
ans =
3.1529
2）normplot(s) %检验分布的正态性
初步可以断定为正态分布
3）[muhat,sigmahat,muci,sigmaci]=normfit(s)%估计正态分布的参数
muhat =
80.1000
sigmahat =
9.7106
muci =
77.5915
82.6085
sigmaci =
8.2310
11.8436
[h,sig,ci]=ttest(s,80.1)%检验参数
h =
0
sig =
1
ci =
77.5915 82.6085
可以认为学生的考试成绩服从正态分布，平均成绩为80.1
2．混凝土的抗压强度随养护时间的延长而增加，现将一批混凝土作成12个试块，记录了养护日期x（日）及抗压强度y（kg/cm2）的数据：
养护时间x
2
3
4
5
7
9
12
14
17
21
28
56
抗压强度y
35
42
47
53
59
65
68
73
76
82
86
99
试求
型回归方程。
解一：对将要拟合的非线性模型
，建立M文件volum.m如下
function yhat=volum(beta,x)
yhat=beta(1)+beta(2)*log(x);
输入数据
x=[2 3 4 5 7 9 12 14 17 21 28 56];
y=[35 42 47 53 59 65 68 73 76 82 86 99];
beta0=[5 1]';
求回归系数
[beta,r,J]=nlinfit(x',y','volum',beta0);
Beta
结果：
beta =
21.0053
19.5287
即得回规模型为
预测和作图
解二：令ln x=u，则x=e
.原回归方程为
输入数据
x=[2 3 4 5 7 9 12 14 17 21 28 56]';
u=log(x);
u=[ones(12,1) u];
y=[35 42 47 53 59 65 68 73 76 82 86 99]';
回归分析及检测
[b,bint,r,rint,stats]=regress(y,u);
b,bint,stats
结果：
b =
21.0058
19.5285
bint =
19.4463 22.5653
18.8943 20.1627
stats =
1.0e+003 *0.0010 4.7069 0.0000
作残差图 rcoplot(r,rint)%作残差图
预测及作图
z=b(1)+b(2)*log(x)%预测及作图
plot(x,y,'k+',x,z,'r')
z =
34.5419
42.4600
48.0780
52.4357
59.0065
63.9143
69.5323
72.5426
76.3342
80.4608
86.0788
99.6149
3.电影院调查电视广告费用和报纸广告费用对每周收入的影响，得到下面的数据，试建立回归模型，并进行检验（写出模型检验的依据），并预测电视广告费用为1，报纸广告费用为6时的周收入（写出预测的程序指令）。
每周收入
96
90
95
92
95
95
94
94
电视广告费
1.5
2.0
1.5
2.5
3.3
2.3
4.2
2.5
报纸广告费
5.0
2.0
4.0
2.5
3.0
3.5
2.5
3.0
解：选择纯二次模型，即
方法一用多元二项式回归命令
数据输入
x1=[1.5 2.0 1.5 2.5 3.3 2.3 4.2 2.5];
x2=[5.0 2.0 4.0 2.5 3.0 3.5 2.5 3.0];
y=[96 90 95 92 95 95 94 94]';
x=[x1' x2'];
回归、检测及预测
rstool(x,y,'purequadratic')
beta,rmse
结果：
beta =
76.3019
1.5280
6.6454
-0.0779
-0.6252
rmse =
0.2496
故回归模型为
剩余标准差为0.2496，显著性很好。
预测电视广告费用为1，报纸广告费用为6时的周收入为95.1186
修改rstool函数：maxx = max(x)+1 minx = min(x)-1
在左边图形下方的方框中输入1，在右边图形下方的方框中输入6，则画面左边的“predicted Y1”下方的数据即是周收入。
方法二化为多元线性回归
程序
x1=[1.5 2.0 1.5 2.5 3.3 2.3 4.2 2.5];
x2=[5.0 2.0 4.0 2.5 3.0 3.5 2.5 3.0];
y=[96 90 95 92 95 95 94 94]';
x=[ones(8,1) x1' x2' (x1.^2)' (x2.^2)'];
[b,bint,r,rint,stats]=regress(y,x);
b,stats
结果
b =
76.3019
1.5280
6.6454
-0.0779
-0.6252
stats =
0.9930 107.0540 0.0014 0.0623
2．混凝土的抗压强度随养护时间的延长而增加，现将一批混凝土作成12个试块，记录了养护日期x（日）及抗压强度y（kg/cm2）的数据：
养护时间x
2
3
4
5
7
9
12
14
17
21
28
56
抗压强度y
35
42
47
53
59
65
68
73
76
82
86
99
试求
型回归方程。
解一：对将要拟合的非线性模型
，建立M文件volum.m如下
function yhat=volum(beta,x)
yhat=beta(1)+beta(2)*log(x);
输入数据
x=[2 3 4 5 7 9 12 14 17 21 28 56];
y=[35 42 47 53 59 65 68 73 76 82 86 99];
beta0=[5 1]';
求回归系数
[beta,r,J]=nlinfit(x',y','volum',beta0);
Beta
结果：
beta =
21.0053
19.5287
即得回规模型为
预测和作图
解二：令ln x=u，则x=e
.原回归方程为
输入数据
x=[2 3 4 5 7 9 12 14 17 21 28 56]';
u=log(x);
u=[ones(12,1) u];
y=[35 42 47 53 59 65 68 73 76 82 86 99]';
回归分析及检测
[b,bint,r,rint,stats]=regress(y,u);
b,bint,stats
结果：
b =
21.0058
19.5285
bint =
19.4463 22.5653
18.8943 20.1627
stats =
1.0e+003 *0.0010 4.7069 0.0000
作残差图 rcoplot(r,rint)%作残差图
预测及作图
z=b(1)+b(2)*log(x)%预测及作图
plot(x,y,'k+',x,z,'r')
z =
34.5419
42.4600
48.0780
52.4357
59.0065
63.9143
69.5323
72.5426
76.3342
80.4608
86.0788
99.6149
3.电影院调查电视广告费用和报纸广告费用对每周收入的影响，得到下面的数据，试建立回归模型，并进行检验（写出模型检验的依据），并预测电视广告费用为1，报纸广告费用为6时的周收入（写出预测的程序指令）。
每周收入
96
90
95
92
95
95
94
94
电视广告费
1.5
2.0
1.5
2.5
3.3
2.3
4.2
2.5
报纸广告费
5.0
2.0
4.0
2.5
3.0
3.5
2.5
3.0
解：选择纯二次模型，即
方法一用多元二项式回归命令
数据输入
x1=[1.5 2.0 1.5 2.5 3.3 2.3 4.2 2.5];
x2=[5.0 2.0 4.0 2.5 3.0 3.5 2.5 3.0];
y=[96 90 95 92 95 95 94 94]';
x=[x1' x2'];
回归、检测及预测
rstool(x,y,'purequadratic')
beta,rmse
结果：
beta =
76.3019
1.5280
6.6454
-0.0779
-0.6252
rmse =
0.2496
故回归模型为
剩余标准差为0.2496，显著性很好。
预测电视广告费用为1，报纸广告费用为6时的周收入为95.1186
修改rstool函数：maxx = max(x)+1 minx = min(x)-1
在左边图形下方的方框中输入1，在右边图形下方的方框中输入6，则画面左边的“predicted Y1”下方的数据即是周收入。
方法二化为多元线性回归
程序
x1=[1.5 2.0 1.5 2.5 3.3 2.3 4.2 2.5];
x2=[5.0 2.0 4.0 2.5 3.0 3.5 2.5 3.0];
y=[96 90 95 92 95 95 94 94]';
x=[ones(8,1) x1' x2' (x1.^2)' (x2.^2)'];
[b,bint,r,rint,stats]=regress(y,x);
b,stats
结果
b =
76.3019
1.5280
6.6454
-0.0779
-0.6252
stats =
0.9930 107.0540 0.0014 0.0623
2．混凝土的抗压强度随养护时间的延长而增加，现将一批混凝土作成12个试块，记录了养护日期x（日）及抗压强度y（kg/cm2）的数据：
养护时间x
2
3
4
5
7
9
12
14
17
21
28
56
抗压强度y
35
42
47
53
59
65
68
73
76
82
86
99
试求
型回归方程。
解一：对将要拟合的非线性模型
，建立M文件volum.m如下
function yhat=volum(beta,x)
yhat=beta(1)+beta(2)*log(x);
输入数据
x=[2 3 4 5 7 9 12 14 17 21 28 56];
y=[35 42 47 53 59 65 68 73 76 82 86 99];
beta0=[5 1]';
求回归系数
[beta,r,J]=nlinfit(x',y','volum',beta0);
Beta
结果：
beta =
21.0053
19.5287
即得回规模型为
预测和作图
解二：令ln x=u，则x=e
.原回归方程为
输入数据
x=[2 3 4 5 7 9 12 14 17 21 28 56]';
u=log(x);
u=[ones(12,1) u];
y=[35 42 47 53 59 65 68 73 76 82 86 99]';
回归分析及检测
[b,bint,r,rint,stats]=regress(y,u);
b,bint,stats
结果：
b =
21.0058
19.5285
bint =
19.4463 22.5653
18.8943 20.1627
stats =
1.0e+003 *0.0010 4.7069 0.0000
作残差图 rcoplot(r,rint)%作残差图
预测及作图
z=b(1)+b(2)*log(x)%预测及作图
plot(x,y,'k+',x,z,'r')
z =
34.5419
42.4600
48.0780
52.4357
59.0065
63.9143
69.5323
72.5426
76.3342
80.4608
86.0788
99.6149
3.电影院调查电视广告费用和报纸广告费用对每周收入的影响，得到下面的数据，试建立回归模型，并进行检验（写出模型检验的依据），并预测电视广告费用为1，报纸广告费用为6时的周收入（写出预测的程序指令）。
每周收入
96
90
95
92
95
95
94
94
电视广告费
1.5
2.0
1.5
2.5
3.3
2.3
4.2
2.5
报纸广告费
5.0
2.0
4.0
2.5
3.0
3.5
2.5
3.0
解：选择纯二次模型，即
方法一用多元二项式回归命令
数据输入
x1=[1.5 2.0 1.5 2.5 3.3 2.3 4.2 2.5];
x2=[5.0 2.0 4.0 2.5 3.0 3.5 2.5 3.0];
y=[96 90 95 92 95 95 94 94]';
x=[x1' x2'];
回归、检测及预测
rstool(x,y,'purequadratic')
beta,rmse
结果：
beta =
76.3019
1.5280
6.6454
-0.0779
-0.6252
rmse =
0.2496
故回归模型为
剩余标准差为0.2496，显著性很好。
预测电视广告费用为1，报纸广告费用为6时的周收入为95.1186
修改rstool函数：maxx = max(x)+1 minx = min(x)-1
在左边图形下方的方框中输入1，在右边图形下方的方框中输入6，则画面左边的“predicted Y1”下方的数据即是周收入。
方法二化为多元线性回归
程序
x1=[1.5 2.0 1.5 2.5 3.3 2.3 4.2 2.5];
x2=[5.0 2.0 4.0 2.5 3.0 3.5 2.5 3.0];
y=[96 90 95 92 95 95 94 94]';
x=[ones(8,1) x1' x2' (x1.^2)' (x2.^2)'];
[b,bint,r,rint,stats]=regress(y,x);
b,stats
结果
b =
76.3019
1.5280
6.6454
-0.0779
-0.6252
stats =
0.9930 107.0540 0.0014 0.0623
3、对实验现象、数据及观察结果的分析与讨论：
对于第1题的实验结果的分析，第1题是这4道题中最简单的，我是根据书上的例题及自己的理解将图形及运行代码写出来，直接输入命令语句以及画图语句就可将图形输出，操作较简单。
对于第2道题的实验结果及数据的分析，虽然第2题的程序我自己按照书上的要求写的，但是系统总是显示有错误，原来我将主函数及输入命令放在了同一个窗口，所以系统总是显示错误，我只好将程序重新书写并做了修改，将他们分别放入不同的窗口，最后得出了实验结果，并正确得出了图形。
对于第3题是这几道题目中算较难的，因为它所给的微分方程中含有3个参数，而它只给了我们其中的两个，第三个参数需要我们自己根据它的范围来了解方程解得变化，我在书本上都没有找到这样的例子，不知道该如何做，后来问了旁边的同学，才了解了这道题的意义，写出了代码，画出了图形，由于是需要需找规律，所以需要我们输入多个a的取值，得到多个图形，最后也得出了实验结果：当a增大时，x，y，z慢慢不收敛而发散，当到一定值后而趋于稳定。
对于第4道题
4、结论：通过本次实验的学习，让我掌握微分方程(组)求解方法（解析法、欧拉法、梯度法、改进欧拉法等），对常微分方程的数值解法有一个初步了解，同时学会使用MATLAB软件求解微分方程的基本命令，学会建立微分方程方面的数学模型。这对于我想深入理解微分、积分的数学概念，掌握数学的分析思维方法，熟悉处理大量的工程计算问题的方法是很有帮助的。
5、实验总结
（1）本次实验成败之处及其原因分析： 本实验的成功之处在于对微分方程(组)模型及求解微分方程的基本命令等的知识的运用；失败之处在于对matlab软件掌握的还不是特别的熟练以及对函数调用的方法还不是特别熟悉，尤其是函数调用这方面的知识还不是特别的熟悉，因此要加强这方面的练习。
（2）本实验的关键环节及改进措施：
做好本实验需要把握的关键环节：
将求解微分方程的主函数和输入命令语句分开编写，否则容易出现错误，得不出实验结果；了解各个库函数，如ode45,ode23,ode113,ode15s,ode23s等求解微分数值函数的使用方法，另外掌握一些基本函数简化程序，以便在使用过程中方便使用；能够根据图形的变换的出一定的数学规律，并掌握绘制图形的一般方法及微分方程数值解得方法在实际问题中的应用，解决实际生活中的问题。
若重做本实验，为实现预期效果，仪器操作和实验步骤应如何改善：
我认为要实现预期的效果，关键是掌握matlab软件的使用方法，对于实验步骤我觉得应该进一步完善，比如使用一些常用的函数，并掌握它们的用法。首先应该把课本的例题在机子上操作一遍，增加自己对基础知识的理解，这样也许就可以做的快点，做好实验之后理解得也要更深。实在做不出来时，可以请教老师或者同学，相互交流，以便得出正确的实验结果。另外还应该尽可能多的使用matlab软件，做到熟能生巧，灵活应用。并且要识记一些常用的函数命令，掌握方法，加强运用。
（3）对实验的自我评价：
通过做本次实验我认识到了matlab软件在实际运用中对于我们的帮助很大，不仅能够解决一些数学方面的问题，对于我们实际生活中的问题也能够解决。对于本实验我基本上算是完成了，但有一些代码及图形出来的效果不是特别的理想，我认为我还需要进一步的对matlab软件熟悉，不仅是一些基本的操作，还应该熟记一些函数的用法，这样就不需要每次写代码的时候都翻课本了，可以节约很多的时间了。本次实验比较难，较上几次的实验难度加大了很多，但我基本上完成了实验的要求，不懂的就翻翻书，重温下知识，或是问下周围的同学和老师，总的来说，这次实验让我对matlab软件有了新的看法，它不仅能够帮助我们解决一些简单的问题，也能够帮助我们解决一些实际生活中的问题，对我们有很大的帮助，对于实验的结果我还是比较满意的。
指导老师评语及得分：
成绩： 指导教师签名：
批阅日期：
浙公网安备 33021202002483