哈尔滨商业大学
数学实验报告
实验题目:_ 因子分析____ ______
姓 名: 郝宇 学号: 201111530011
专 业: 数学与应用数学1班
日 期: 20##-11-22
一、实验目的
利用SPSS软件对数据进行因子分析
二、实验内容
在此实验中,通过使用SPSS对牙膏的偏好调查结果进行因子分析。
三、实验步骤及结论
(一)1、首先使用SPSS打开“牙膏偏好调查”。
2、将变量视图中名称分别为V1-V6的“标签”一栏中写入防蛀、亮白、保护牙龈、口气清新、治疗坏牙、提升魅力;
-(二)1、选择菜单栏中“分析-降维-因子分析”。
将分组V1-V6移入分组变量列表中
单击确定,执行因子分析过程。
1.该组数据是否适合做因子分析
Sig=0.000<0.05,所以适合做因子分析
2.提取多少公共因子,原因是什么
3.写出因子模型
X1=0.928X1+0.253X2+e1;
X2=-0.301X1+0.795X2+e2;
X3=0.936X1+0.131X2+e3;
X4=-0.342X1+0.789X2+e4;
X5=-0.869X1-0.351X2+e5;
X6=-0.177X1+0.871X2+e6;
4.对因子命名
1、3、5护牙因子
2、4、6美牙因子
5每个指标的共同度方差贡献
方差贡献为82.49%
四、心得体会
通过本次试验,我了解了SPSS的一些基本操作,明白了怎么使用SPSS对数据做因子分析。对于SPSS,我们有了更进一步的了解。
第二篇:数学实验报告——利用MALTAB进行数据的统计与分析
实验八 数据的统计与分析
一、滚珠
㈠问题描述
某厂从一台机床生产的滚珠中随机抽取20个,侧直径如下:
[14.6 14.7 15.1 14.9 14.8 15 15.1 15.2 14.8 14.3
15.1 14.2 14.4 14 14.6 15.1 14.9 14.7 14.5 14.7];
计算这些数据的均值、标准差、方差、极差,并画出直方图。
㈡简要分析
本题是一道简单的统计计算题,直接调用MATLAB函数进行计算、画图即可。
㈢结果与分析
平均值 dbar = 14.735;
标准差stdd = 0.332890056135181
方差 vard = 0.110815789473684
极差 ranged = 1.2
直方图如下:
㈣程序清单
clear;clc;
d = [14.6 14.7 15.1 14.9 14.8 15 15.1 15.2 14.8 14.3,...
15.1 14.2 14.4 14 14.6 15.1 14.9 14.7 14.5 14.7];
sumd = sum(d);
dbar = sumd/20;
vard = var(d);
stdd = std(d);
ranged = range(d);
hist(d)
二、炮弹射击
㈠问题描述
炮弹射击目标为一半径100m的圆形区域,弹着点亿元圆心为中心成二维正态分布,给定系数,求炮弹命中圆形区域的概率。
㈡方法与公式
使用书中266页所给的积分公式
参照例六方法,取b=d=1,f取二维正态分布:
㈢结果与分析
1、计算结果
计算求得求得积分值为p = 0.795626646433925。
2、分析
计算题目所给问题的同时,我也计算了以随机分布发射炮弹后,炮弹落入目标的概率:
p1 = 0.78449。
p1*4 = 3.13796,对应pi = 3.1416,可见虽有误差但是也还算接近。
注意,p的值应当大于p1,这是因为题目中炮弹的命中点符合而为正态分布,而二维正态分布在中心区域的概率肯定要比均匀分布在中心区域的概率大。因此p应当大于p1。
为了取得更准确的结果,这里将n扩大后再计算多次。得到:
从实验结果可以看出,命中目标的概率约为0.8。
㈣程序清单
clear;clc;
sx = 0.8;sy=0.5;p = 0.4;
n=100000;m=0;z=0;
x = unifrnd(0,1,1,n);
y = unifrnd(0,1,1,n);
for i = 1:n
if x(i)^2+y(i)^2 <= 1
u = exp(-0.5/(1-p^2)*(x(i)^2/sx^2-2*p*x(i)*y(i)/(sx*sy)+y(i)^2/sy^2));
z = z+u;
m = m+1;
end
end
p = 4*0.5*z/pi/sx/sy/sqrt(1-p^2)/n;
p1 = m/n;
三、报童问题
㈠问题描述
对报童问题,如果报纸的需求量服从正态分布,且给定批发价,建立数学模型。给定参数,求解报童每天购进的报纸数量n。
㈡简要分析
首先,对所给问题进行数学建模,考虑到本问题与原问题差异不大,可以直接在原表达式
上进行修改。
将所给a的表达式代入后,对n求导,得到:
令导数为0有:
由于,可以将上面从0开始的积分改为从负无穷开始
解此方程即可得到问题的解。
本题中,首先画出函数及导数的大致图像,选择较接近的初值,再使用fsolve求解。最后,拓展计算报童在购买特定数量的报纸后所得的收入。计算数值积分采用自适应辛普森公式完成。
㈢结果与分析
1、函数及图象
函数图像:
化简后的导数(已经约去系数):
y=b-A+2*A*n/k-(b-c)*normcdf([n],u,o);
图像:
可以看出,函数的0点在2000左右,因此初值取2000代入fzero中。
2、计算结果
(1)购买数量
n = 1.968206023024722e+03
从而,报童应当购买1968份报纸。
(2)预计收入
使用数值积分方法(自适应辛普森公式)求解积分,从而计算出报童购买1968份报纸后预计获得的收入。
计算结果:
S=37.545487871203400.
3、分析
从函数的图像可以看出收益最大值约为37.54。从导数图像可以看出该点对应的的确为最大值。
4、小结
本道题应当说是这次作业中较复杂的一道,需要先建立模型,数学推导,最终进行数学实验。当然本题带来的收获也是最大的,要求结合以前学过的非线性方程求解方法。此外,我还计算了最终的函数值,这就还需要第一次实验时数值积分的内容,很考查灵活运用知识的能力。
㈣程序清单
1、化简后的导数
function y = fun(n)
u=2000;
o=50;
A=0.5;
k=50000;
b=0.5;
c=0.35;
y=b-A+2*A*n/k-(b-c)*normcdf([n],u,o);
2、画导数图像
x = 0:5000;
y=fun(x);
plot(x,y);
3、求极值点
clear;clc;
u=2000;
o=50;
A=0.5;
k=50000;
b=0.5;
c=0.35;
[x1,fv1,ef1,out1]=fzero(@fun,2000);
4、函数计算(1)
function y = value1(x,n)
u=2000;
o=50;
A=0.5;
k=50000;
b=0.5;
c=0.35;
y=((b-A*(1-n/k))*x-(A*(1-n/k)-c)*(n-x)).*normpdf(x,u,o);
5、函数计算(2)
function y = value2(x,n)
u=2000;
o=50;
A=0.5;
k=50000;
b=0.5;
c=0.35;
y=(b-A*(1-n/k))*n*normpdf(x,u,o);
6、计算原函数
clear;clc;
u=2000;
o=50;
A=0.5;
k=50000;
b=0.5;
c=0.35;
n=1968;
x = 1:3000;
s1 = zeros(1,3000);
s2 = zeros(1,3000);
for i = 1:3000
s1(i) = quad(@(x)value1(x,i),0,i);
s2(i) = quad(@(x)value2(x,i),i,7500);
end
s = s1+s2;
plot(x,s);
四、体验与收获
这是本学期数学实验统计部分的第一次作业,总体来说比较顺利。前两道题分别考察了统计基本计算以及蒙特卡洛应用,第三道题则是应用统计理论解决实际问题,对于综合运用知识有很好的帮组。收获简要总结如下:
1、学习了数据的统计与分析等内容;
2、复习了蒙特卡洛法、大数定理等知识;
3、学习使用MATLAB解决统计相关问题;