当前位置:首页 >> 信息与通信 >>

MATLAB基础性实验报告4


一、实验目的:
学习统计方法回归分析的思想和基本原理;掌握建立回归模型的基本步骤,明 确回归分析的主要任务;熟悉 MATLAB 软件进行回归模型的各种统计分析;通过范 例学习,熟悉统计分析思想和建立回归模型的基本要素。 通过该实验的学习,使学生掌握回归分析的统计思想,认识面对什么样的实际 问题可以建立回归模型,并且对回归模型作统计分析,同时使学生学会使用 MATLAB

软件进行回归分析和计算的基本命令,了解统计软件的功能和作用;熟悉处理大量 数据的要领和方法。

二、实验仪器、设备或软件:
电脑,MATLAB 软件

三、实验内容:
1.线性回归模型的建立与分析步骤(问题假设→模型→参数估计→模型检验 →确定最优回归方程→预测) ; 2.非线性回归模型的建立与分析步骤; 3.使用 MATLAB 命令对回归模型进行计算与分析(包括模型检验与预测) ; 4.利用某些数值与图形对统计特征作定性分析。

四、实验过程(实验步骤、记录、数据、分析)
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) 若检验符合正态分布,估计正态分布的参数并检验参数。 解:输入数据 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

12

10

8

6

4

2

0 50

55

60

65

70

75

80

85

90

95

100

2)normplot(s)

%检验分布的正态性
Normal Probability Plot

0.99 0.98 0.95 0.90 0.75

Probability

0.50 0.25 0.10 0.05 0.02 0.01 55 60 65 70 75 Data 80 85 90 95

初步可以断定为正态分布

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 抗压强度 y 2 35 3 42 4 47 5 53 7 59 9 65 12 68 14 73 17 76 21 82 28 86 56 99

? ? a ? b ln x 型回归方程。 试求 y ? ? a ? b ln x ,建立 M 文件 volum.m 如下 解一:对将要拟合的非线性模型 y
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

? ? 21.0053? 19.5287ln x 即得回规模型为 y
预测和作图

100

90

80

70

60

50

40

30

0

10

20

30

40

50

60

解二:令 ln x=u,则 输入数据

x=e

u

? ? a ? bu .原回归方程为 y

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 18.8943

22.5653 20.1627

stats =

1.0e+003 *0.0010 作残差图

4.7069

0.0000

rcoplot(r,rint)%作残差图
Residual Case Order Plot

3

2

1

Residuals

0

-1

-2

-3 2 4 6 8 Case Number 10 12

预测及作图

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

100

90

80

70

60

50

40

30

0

10

20

30

40

50

60


相关文章:
MATLAB基础性实验报告4
一、实验目的:学习统计方法回归分析的思想和基本原理;掌握建立回归模型的基本步骤,明 确回归分析的主要任务;熟悉 MATLAB 软件进行回归模型的各种统计分析;通过范 例...
MATLAB基础性实验报告4
MATLAB基础性实验报告4_信息与通信_工程科技_专业资料。学生实验报告 一、实验目的 1. 归纳和学习求解常微分方程(组)的基本原理和方法; 2. 掌握解析、数值解法,...
MATLAB基础性实验报告5
MATLAB基础性实验报告5_调查/报告_表格/模板_应用文书。一、实验要求:了解插值与...四、实验步骤: 1.开启软件平台——MATLAB,开启 MATLAB 编辑窗口; 2.根据各种...
MATLAB基础性实验报告3
MATLAB,开启 MATLAB 编辑窗口; 2.根据微分方程求解步骤编写 M 文件 3.保存文件并运行; 4.观察运行结果(数值或图形); 5.根据观察到的结果和体会写出实验报告...
2014《MATLAB及应用》实验报告4 -
实验报告4_MATLAB基础 5页 免费 matlab实验报告4 2页 1下载券 matlab高级编程...MATLAB 及应用 上机实验 2014-04 -16 盛义发 班学姓成 级号名绩 2013 级...
上机实验报告-操作性实验4
上机实验报告-操作性实验4_学习总结_总结/汇报_实用文档。昆明理工大学理学院 ...实验环境:Matlab 运行环境 实验要求: 1.掌握正态分布随机数的产生与操作命令 2...
Matlab实验报告4
Matlab实验报告4_计算机软件及应用_IT/计算机_专业资料。Matlab实验报告4,西北工业大学一、实验内容 1. 输入20个数,求其中最大数和最小数。要求分别用循环结构和调...
南华大学MATLAB实验报告4
南华大学MATLAB实验报告4_调查/报告_表格/模板_实用文档。2015-2016南华大学...plot(x, y):基本格式,x 和 y 可为向量或矩阵。 plot(x1, y1, x2, y...
机械CAD技术基础实验报告4
机械CAD技术基础实验报告4_机械/仪表_工程科技_专业资料。实验报告 课程名称: 机械...实验报告 课程名称: 机械 CAD 技术基础 实验项目: 用 MATLAB 进行插值和拟合 ...
信号与系统实验报告4
实验报告课程名称:信号与系统 实验项目名称:实验 4 信号抽样及抽样定理 实验目的: 学会运用 MATLAB 完成信号抽样及对抽样信号的频谱进行分析;学会运用 MATLAB 改变...
更多相关标签:
计算机基础实验报告4 | matlab实验报告 | matlab实验报告总结 | matlab上机实验报告 | matlab实验报告模板 | matlab仿真实验报告 | matlab数学实验报告 | matlab实验报告答案 |