实验七美国七种犯罪的主成分分析
一、实验目的
1.掌握如何使用SAS软件来进行主成分分析和聚类分析；
2.看懂和理解SAS输出的结果，并学会以此来作出分析；
3.掌握对实际数据如何来进行主成分分析；
4.对同一组数据使用五种系统聚类方法，学会对各种聚类效果的比较，获取重要经验；
5.掌握使用主成分进行聚类
二、实验内容
下表中给出的是美国50个州每天100000个人中七种犯罪的比率数据。
这七种犯罪是：
X1:杀人犯 x2：强奸罪 x3：抢劫罪
X4：斗殴罪 x5：夜盗罪 x6：偷盗罪 x7：汽车犯罪
试对表中犯罪数据进行主成分分析。
state x1x2x3x4x5x6x7 Alabama14.225.296.8278.31135.51881.9280.7 Alaska10.851.696.82841331.73369.8753.3 Arizona9.534.2138.2312.32346.14467.4439.5 Arkansas8.827.683.2203.4972.61862.1183.4 California11.549.42873582139.43499.8663.5 Colorado 6.342170.7292.91935.23903.2477.1 Connecticut 4.216.8129.5131.813462620.7593.2 Delaware624.9157194.21682.63678.4467 Florida10.239.6187.9449.11859.93840.5351.4 Georgia11.731.1140.5256.51351.12170.2297.9 Hawaii7.225.512864.11911.53920.4489.4 Idaho 5.519.439.6172.51050.82599.6237.6 Illinois9.921.8211.320910852828.5528.6 Indiana7.426.5123.2153.51086.22498.7377.4 Iowa 2.310.641.289.8812.52685.1219.9 Kansas 6.622100.7180.51270.42739.3244.3 Kentucky10.119.181.1123.3872.21662.1245.4 Louisiana15.530.9142.9335.51165.52469.9337.7 Maine 2.413.538.71701253.12350.7246.9 Maryland834.8292.1358.914003177.7428.5 Massachuset 3.120.8169.1231.61532.22311.31140.1 Michigan9.338.9261.9274.61522.73159545.5 Minnesota 2.719.585.985.81134.72559.3343.1 Mississippi14.319.665.7189.1915.61239.9144.4 Missouri9.628.3189233.51318.32424.2378.4 Montana 5.416.739.2156.8804.92773.2309.2 Nebraska 3.918.164.7112.77602316.1249.1 Nevada15.849.1323.13552453.14212.6559.2 New Hampshi 3.210.723.2761041.72343.9293.4 New Jersey 5.621180.4185.11435.82774.5511.5 New Mexico8.839.1109.6343.41418.73008.6259.5 New York10.729.4472.6319.117282782745.8 North Carol10.61761.3318.31154.12037.8192.1 North Dakot0.9913.343.8446.11843144.7 Ohio7.827.3190.5181.112162696.8400.4 Oklahoma8.629.273.82051288.22228.1326.8 Oregon 4.939.9124.1286.91636.43506.1388.9 Pennsylvani 5.619130.3128877.51624.1333.2 Rhode Islan 3.610.586.52011489.52844.1791.4 South Carol11.933105.9485.31613.62342.4245.1 South Dakot213.517.9155.7570.51704.4147.5 Tennessee10.129.7145.8203.91259.71776.5314 Texas13.333.8152.4208.21603.12988.7397.6 Utah 3.520.368.8147.31171.63004.6334.5 Vermont 1.415.930.8101.21348.22201265.2 Virginia923.392.1165.7986.22521.2226.7 Washington 4.339.6106.2224.81605.63386.9360.3 West Virgin613.242.290.9597.41341.7163.3 Wisconsin 2.812.952.263.7846.92614.2220.7 Wyoming 5.421.939.7173.9811.62772.2282
四、实验指导
1.直接编程进行主成分分析
proc princomp data=work.aa out=prin;
var x1-x7;
proc sort;
by prin1;
proc print;
id state;
var prin1 prin2;
proc sort;
by prin2;
proc print;
id state;
var prin1 prin2;
proc plot;
plot prin2*prin1=alphabet/haxis=-4 to6 by 1vaxis=-3.5 to 2.5 by 0.5;
run;
2.分析家进行主成分分析
在t中打开数据集aa，见图1。选菜单过程如下：
在图1中选统计?多元分析?在变量框中选x1，x2，x3，x4，x5，x6，x7，x8 (见
图2) ?variables?ok
图1
图2 得到如图6所示的结果：
图6
从图5可以看出，前两个和前三个主成分的累计贡献率分别达到78.5％和86.8％，第一
主成分1?y 在所有变量(除在*
2x 上的载荷稍偏小外)上都有近似相等的正载荷，反映了综合消费
性支出的水平，因此第一主成分可称为综合性犯罪成分。第二主成分2?y
在变量x3，x7上有很高的正载荷，在变量x1上有中等的正载荷，而在其余变量上有负载荷或很小的正载荷。可以认为这个主成分度量了暴力犯罪成分(主要是抢劫罪x3，其次是汽车犯罪x7)在所有犯罪成分中占的比重(也可理解为一种犯罪倾向)，第二主成分可称为暴力倾向成分。前三个主成分的总贡献率能达到86%，但第三主成分的特征值小于1，因此我们只取前面两个主成分。
Pcr1=0.3*x1+0.43*x2+0.397*x3+0.397*x4+0.44*x5+0.36*x6+0.3*x7 Pcr2=-0.63*x1-0.17*x2+0.04*x3-0.34*x4+0.2*x5+0.4*x6+0.5*x7 散点图：
在图1中选分析?散点图（Y X ）?选 PCR1→X ；PCR2→Y ；region →标签变量 (见图7)?确定
随即出现如图8所示的散点图，接下来我们对该图进行设置、调整，使之符合我们的要求。
图7
浙江万里学院 期末考试1 课程名称:统计学原理(专升本) 1、(单选题)某地区有10万人口,共有80个医院。平均每个医院要服务1250人,这个指标就是( )。(本题 2、0分) A、平均指标 B、强度相对指标 C、总量指标 D、发展水平指标 答案:B、 解析:无、 2、(单选题)总量指标数值大小( )(本题2、0分) A、随总体范围扩大而增大 B、随总体范围扩大而减小 C、随总体范围缩小而增大 D、与总体范围大小无关 答案:A、 解析:无、 3、(单选题)数据筛选的主要目的就是( )。(本题2、0分) A、发现数据的错误 B、对数据进行排序 C、找出所需要的某类数据 D、纠正数据中的错误 答案:C、 解析:无、 4、(单选题)抽样调查所特有的误差就是( )。(本题2、0分) A、A 由于样本的随机性而产生的误差 B、B 登记误差 C、C 系统性误差 D、D ABC都错 答案:A、 解析:无、 5、(单选题)某工厂有100名职工,把她们的工资加总除以100,这就是对100个( )求平均数(本题2、0分) A、A 变量
B、B 标志 C、C 变量值 D、D 指标 答案:C、 解析:无、 6、(单选题)“增长1%的绝对值”反映的就是同样的增长速度在不同条件下所包含的绝对水平。(本题2、0分) A、计量单位 B、数据类型 C、时间 D、调查方法 答案:C、 解析:无、 7、(单选题)在其她条件不变的情况下,当总体数据的离散程度较大时,总体均值的置信区间( )。(本题2、0分) A、可能变宽也可能变窄 B、变窄 C、变宽 D、保持不变 答案:C、 解析:无、 8、(单选题)将某企业职工的月收入依次分为2000元以下、2000元～3000元、3000元～4000元、4000元～5000元、5000元以上几个组。第一组的组中值近似为( )(本题2、0分) A、2000 B、1000 C、1500 D、2500 答案:C、 解析:无、 9、(单选题)统计指标有如下特点( )。(本题2、0分) A、A 数量性、综合性、具体性 B、B 准确性、及时性、完整性 C、C 大量性、同质性、差异性 D、D 科学性、客观性、社会性 答案:A、 解析:无、 10、(单选题)设 zc 为检验统计量的计算值,检验的假设为H 0 : μ =μ0 , H 1 : μ ≠μ0 , zc = - 1、96 时,计算出的 P 值为( )。(本题2、0
Matlab作业3（数值分析） 机电工程学院（院、系）专业班组 学号姓名实验日期教师评定 1.计算多项式乘法(x2+2x+2)(x2+5x+4)。 答： 2. （1）将(x-6)(x-3)(x-8)展开为系数多项式的形式。（2）求解在x=8时多项 式(x-1)(x-2) (x-3)(x-4)的值。 答：（1） （2）
3. y=sin(x)，x从0到2π，?x=0.02π，求y的最大值、最小值、均值和标准差。 4.设ｘ＝[0.00.30.8 1.1 1.6 2.3]'，y=[0.500.82 1.14 1.25 1.35 1.40]'，试求二次多项式拟合系数，并据此计算x1=[0.9 1.2]时对应的y1。解：x=[0.0 0.3 0.8 1.1 1.6 2.3]'; %输入变量数据x y=[0.50 0.82 1.14 1.25 1.35 1.40]'; %输入变量数据y p=polyfit(x,y,2) %对x,y用二次多项式拟合,得到系数p x1=[0.9 1.2]; %输入点x1 y1=polyval(p,x1) %估计x1处对应的y1 p = -0.2387 0.9191 0.5318 y1 = a) 1.2909
5.实验数据处理：已知某压力传感器的测试数据如下表 p为压力值，u为电压值，试用多项式 d cp bp ap p u+ + + =2 3 ) ( 来拟 合其特性函数，求出a,b,c,d，并把拟合曲线和各个测试数据点画在同一幅图上。解： >> p=[0.0,1.1,2.1,2.8,4.2,5.0,6.1,6.9,8.1,9.0,9.9]; u=[10,11,13,14,17,18,22,24,29,34,39]; x=polyfit(p,u,3) %得多项式系数 t=linspace(0,10,100); y=polyval(x,t); %求多项式得值 plot(p,u,'*',t,y,'r') %画拟和曲线 x = 0.0195 -0.0412 1.4469 9.8267
数据分析 实验指导书 理学院实验中心数学专业实验室编写
实验一SAS系统的使用 【实验类型】（验证性） 【实验学时】2学时 【实验目的】使学生了解SAS系统，熟练掌握SAS数据集的建立及一些必要的SAS语句。 【实验内容】 1. 启动SAS系统，熟悉各个菜单的内容；在编辑窗口、日志窗口、输出窗口之间切换。 2. 建立数据集 表1 Name Sex Math Chinese English Alice f908591 Tom m958784 Jenny f939083 Mike m808580 Fred m848589 Kate f978382 Alex m929091 Cook m757876 Bennie f827984 Hellen f857484 Wincelet f908287 Butt m778179 Geoge m868582 Tod m898484 Chris f898487 Janet f866587 1）通过编辑程序将表1读入数据集sasuser.score; 2）将下面记事本中的数据读入SAS数据集，变量名为code name scale share price: 000096 广聚能源8500 0.059 1000 13.27 000099 中信海直6000 0.028 2000 14.2 000150 ST麦科特12600 -0.003 1500 7.12 000151 中成股份10500 0.026 1300 10.08 000153 新力药业2500 0.056 2000 22.75
3)将下面Excel表格中的数据导入SAS数据集work.gnp； name x1 x2 x3 x4 x5 x6 北京190.33 43.77 7.93 60.54 49.01 90.4 天津135.2 36.4 10.47 44.16 36.49 3.94 河北95.21 22.83 9.3 22.44 22.81 2.8 山西104.78 25.11 6.46 9.89 18.17 3.25 内蒙古128.41 27.63 8.94 12.58 23.99 3.27 辽宁145.68 32.83 17.79 27.29 39.09 3.47 吉林159.37 33.38 18.37 11.81 25.29 5.22 黑龙江116.22 29.57 13.24 13.76 21.75 6.04 上海221.11 38.64 12.53 115.65 50.82 5.89 江苏144.98 29.12 11.67 42.6 27.3 5.74 浙江169.92 32.75 21.72 47.12 34.35 5 安徽153.11 23.09 15.62 23.54 18.18 6.39 福建144.92 21.26 16.96 19.52 21.75 6.73 江西140.54 21.59 17.64 19.19 15.97 4.94 山东115.84 30.76 12.2 33.1 33.77 3.85 河南101.18 23.26 8.46 20.2 20.5 4.3 湖北140.64 28.26 12.35 18.53 20.95 6.23 湖南164.02 24.74 13.63 22.2 18.06 6.04 广东182.55 20.52 18.32 42.4 36.97 11.68 广西139.08 18.47 14.68 13.41 20.66 3.85 四川137.8 20.74 11.07 17.74 16.49 4.39 贵州121.67 21.53 12.58 14.49 12.18 4.57 云南124.27 19.81 8.89 14.22 15.53 3.03 陕西106.02 20.56 10.94 10.11 18 3.29 甘肃95.65 16.82 5.7 6.03 12.36 4.49 青海107.12 16.45 8.98 5.4 8.78 5.93 宁夏113.74 24.11 6.46 9.61 22.92 2.53 新疆123.24 38 13.72 4.64 17.77 5.75 4)使用VIEWTABLE格式新建数据集earn,输入如表所示数据Year earn 1981 125000 1982 136000 1983 122350 1984 65200 1985 844600 1986 255000 1987 265000 1988 280000 1989 136000
浙江万里学院实验报告 课程名称：数据通信与计算机网络及实践 实验名称：OSPF路由协议配置 专业班级：姓名：小组学号：2012014048实验日期：6.6
再测试。要求写出两台路由器上的ospf路由配置命令。
[RTC-rip-1]import ospf [RTC-rip-1]quit [RTC]ospf [RTC-ospf-1]import rip [RTC-ospf-1]quit
结合第五步得到的路由表分析出现表中结果的原因： RouteB 通过RIP学习到C和D 的路由情况，通过OSPF学习到A 的路由信息 实验个人总结 班级通信123班本人学号后三位__048__ 本人姓名_ 徐波_ 日期2014.6.06 本次实验是我们的最后一次实验，再次之前我们已经做了很多的有关于华为的实验，从一开始的一头雾水到现在的有一些思路，不管碰到什么问题，都能够利用自己所学的知识去解决或者有一些办法。这些华为实验都让我受益匪浅。 实验个人总结 班级通信123班本人学号后三位__046__ 本人姓名_ 金振宁_ 日期2014.6.06 这两次实验都可以利用软件在寝室或者去其他的地方去做，并不拘泥于实验室，好好的利用华为的模拟机软件对我们来说都是非常有用的。 实验个人总结 班级通信123班本人学号后三位__044_ 本人姓名_ 陈哲日期2014.6.06
理解OSPF路由协议，OSPF协议具有如下特点： 适应范围：OSPF 支持各种规模的网络，最多可支持几百台路由器。 快速收敛：如果网络的拓扑结构发生变化，OSPF 立即发送更新报文，使这一变化在自治系统中同步。 无自环：由于OSPF 通过收集到的链路状态用最短路径树算法计算路由，故从算法本身保证了不会生成自环路由。 实验个人总结 班级通信123班本人学号后三位__050 本人姓名_ 赵权日期2014.6.06 通过本次实验学会了基本的在路由器上配置OSPF路由协议，组建一个简单的路由网络。想必以后的生活中有可能会用到。
实验一 误差分析 实验（病态问题） 实验目的：算法有“优”与“劣”之分，问题也有“好”与“坏”之别。对数值方法的研究而言，所谓坏问题就是问题本身对扰动敏感者，反之属于好问题。通过本实验可获得一个初步体会。 数值分析的大部分研究课题中，如线性代数方程组、矩阵特征值问题、非线性方程及方程组等都存在病态的问题。病态问题要通过研究和构造特殊的算法来解决，当然一般要付出一些代价（如耗用更多的机器时间、占用更多的存储空间等）。 问题提出：考虑一个高次的代数多项式 )1.1() ()20()2)(1()(20 1∏=-=---=k k x x x x x p 显然该多项式的全部根为1,2,…,20共计20个，且每个根都是单重的。现考虑该多项式的一个扰动 )2.1(0 )(19=+x x p ε 其中ε是一个非常小的数。这相当于是对（）中19x 的系数作一个小的扰动。我们希望比较（）和（）根的差别，从而分析方程（）的解对扰动的敏感性。 实验内容：为了实现方便，我们先介绍两个Matlab 函数：“roots ”和“poly ”。 roots(a)u = 其中若变量a 存储n+1维的向量，则该函数的输出u 为一个n 维的向量。设a 的元素依次为121,,,+n a a a ，则输出u 的各分量是多项式方程 01121=+++++-n n n n a x a x a x a 的全部根；而函数 poly(v)b =
的输出b 是一个n+1维变量，它是以n 维变量v 的各分量为根的多项式的系数。可见“roots ”和“poly ”是两个互逆的运算函数。 ;000000001.0=ess );21,1(zeros ve = ;)2(ess ve = ))20:1((ve poly roots + 上述简单的Matlab 程序便得到（）的全部根，程序中的“ess ”即是（）中的ε。 实验要求： （1）选择充分小的ess ，反复进行上述实验，记录结果的变化并分析它们。 如果扰动项的系数ε很小，我们自然感觉（）和（）的解应当相差很小。计算中你有什么出乎意料的发现表明有些解关于如此的扰动敏感性如何 （2）将方程（）中的扰动项改成18x ε或其它形式，实验中又有怎样的现象 出现 （3）（选作部分）请从理论上分析产生这一问题的根源。注意我们可以将 方程（）写成展开的形式， ) 3.1(0 ),(1920=+-= x x x p αα 同时将方程的解x 看成是系数α的函数，考察方程的某个解关于α的扰动是否敏感，与研究它关于α的导数的大小有何关系为什么你发现了什么现象，哪些根关于α的变化更敏感 思考题一：（上述实验的改进） 在上述实验中我们会发现用roots 函数求解多项式方程的精度不高，为此你可以考虑用符号函数solve 来提高解的精确度，这需要用到将多项式转换为符号多项式的函数poly2sym,函数的具体使用方法可参考Matlab 的帮助。
实验课程：数据分析 信息与计算科学 业： 专 级： 班 号：学 姓名： 中北大学理学院．
实验一 SAS系统的使用 【实验目的】 了解SAS系统，熟练掌握SAS数据集的建立及一些必要的SAS语句。 【实验内容】 1. 将SCORE数据集的内容复制到一个临时数据集test。 SCORE数据集 English Math Sex Chinese Name 91 90 f 85 Alice 95 Tom m 87 84 93 90 Jenny f 83 80 85 80 Mike m 84 85 89 m Fred 97 83 f 82 Kate 92 Alex 90 m 91 75 Cook m 78 76 82 f Bennie 79 84 85 Hellen f 74 84 90 82 Wincelet f 87 77 Butt m 81 79 86 85 Geoge m 82 89 Tod m 84 84 89 Chris f 84 87 86 65 f 87 Janet math的高低拆分到3个不同的数据集：SCORE2．将数据集中的记录按照math大于等于90的到good数据集，math在80到89之间的到normal数据集，math 在80以下的到bad数据集。 3．将3题中得到的good，normal，bad数据集合并。 【实验所使用的仪器设备与软件平台】SAS 【实验方法与步骤】 1： DATA SCORE; INPUT NAME $ Sex $ Math Chinese English; CARDS; 2
91 85 Alice f 90 84 Tom m 95 87 83 f 93 90 Jenny 80 80 85 Mike m 89 85 m Fred 84 82 83 Kate f 97 91 Alex m 92 90 76 Cook m 78 75 84 82 79 f Bennie 84 74 Hellen f 85 87 82 Wincelet f 90 79 Butt m 77 81 82 m 86 85 Geoge 84 89 84 Tod m 87 84 f Chris 89 87 Janet f 86 65 ; ; Run PROC PRINT DATA=SCORE; DATA test; SET SCORE; ：2 good normal bad; DATA SCORE; SET; SELECT) output good; 90when(math>=) output normal; 80when(math>=&math<90) output bad; when(math<80; end; Run=good; DATA PRINT PROC=normal; DATA PRINT PROC=bad; DATA PRINT PROC ：3 All; DATA good normal bad; SET=All; DATA PROC PRINT;Run 3 【实验结果】 结果一：
ospf协议,实验报告 篇一：实验7 OSPF路由协议配置实验报告 浙江万里学院实验报告 课程名称：数据通信与计算机网络及实践 实验名称： OSPF路由协议配置专业班级：姓名：小组学号：XX014048 实验日期： 再测试。要求写出两台路由器上的ospf路由配置命令。 第页共页 [RTC-rip-1]import ospf [RTC-rip-1]quit [RTC]ospf [RTC-ospf-1]import rip [RTC-ospf-1]quit 结合第五步得到的路由表分析出现表中结果的原因： RouteB 通过RIP学习到C和D 的路由情况，通过OSPF 学习到A 的路由信息 实验个人总结 班级通信123班本人学号后三位__048__ 本人姓名_徐波_ 日期 本次实验是我们的最后一次实验，再次之前我们已经做了很多的有关于华为的实验，从一开始的一头雾水到现在的有一些思路，不管碰到什么问题，都能够利用自己所学的知识去解决或者有一些办法。这些华为实验都让我受益匪浅。
实验个人总结 班级通信123班本人学号后三位__046__ 本人姓名_金振宁_ 日期 这两次实验都可以利用软件在寝室或者去其他的地方去做，并不拘泥于实验室，好好的利用华为的模拟机软件对我们来说都是非常有用的。 实验个人总结 班级通信123班本人学号后三位本人姓名_陈哲日期 第页共页 篇二：单区域的OSPF协议配置实验报告 学生实验报告 *********学院 篇三：OSPF实验报告 计算机学院 实验报告 （ XX 年春季学期） 课程名称：局域网设计与管理 主讲教师：李辉 指导教师：学生姓名： 学 年郑思楠号： XX012019 级： XX级
一、填空 1、江河水系或某一河段，监测时要求设置背景断面、对照断面、控制断面、监测断面等 断面，断面选定后，应根据水面的宽度确定断面上的采样垂线，再根据采样垂线的深度确定采样点的位置。监测断面和采样点确定后，应设置固定天然或人工的标志物以保证样品的代表性和可比性。 2、地表水的采样时间和采样频次的确定原则是：对较大的水系干流和中小河流，全年采样时间为丰 水期、枯水期、平水期，每期采样两次，流经城市工业区、污染较重的河流、游览水域，饮用水源地全年采样不少于 12 次，采样时间为每月一次或视具体情况而定。 3、为了克服污染物随时间的变化，通常采用优化频次的办法来解决，克服在空间上的不同， 通常采用优化布点的办法来解决。 4、测定COD的水样，采样时用硬质玻璃瓶容器贮存，保存时加硫酸。 5、硝酸-高氯酸消解法，一般先加入硝酸氧化水样中的羟基化合物，稍冷后再加另一种酸 进行消解。 6、水的颜色分为真色和表色，真色是指除去悬浮物后水的颜色，反之则为表色；采用铂钴标准比 色法是先配制标准色列再与水样的颜色进行目视比色确定，稀释倍数法是 7、取某水样50ml于坩埚中蒸干后移入103-105℃烘箱内烘至恒重，冷却后称重为32.9682g，然后再 移入马弗炉中高温燃烧，冷却后称重为31.5349g，已知坩埚的重量为31.3259g，则该水样的TSS(总 悬浮物)= 1.6423 ，TVSS（总可挥发物）= 1.4333 ，总固定性固体浓度= 32846mg/L 。 8、大气监测必测项目有 PM10 、 PM2.5 、 SO2/NO2/CO/O3 、硫酸盐化速率、灰尘 自然降尘量。 9、采集降水的采样器要高于基础面 1.2cm 以上，采集雪水时的采样器要求是塑料容器且上口比 采集降雨时的采样器大10cm以上。 10、固定污染源采样位置要求采样断面设置在阻力构件下游方向大于6倍管道直径处或上游方 向大于3倍管道直径处，难于满足时，应≮1.5d，并适当增加采样点数目和采样频 率，最好在5m/S的断面上设置，优先考虑在垂直管道上设置。 11、测定咽气烟尘浓度必须采用等速采样采样，当采样速度大于采样点的咽气速度 时，导致测量结果偏低，反之结果则偏高。
实验2.1 多项式插值的振荡现象 问题提出：考虑在一个固定的区间上用插值逼近一个函数。显然Lagrange 插值中的使用节点越多，插值多项式的次数就越高.我们自然就关心插值多项式的次数增加时，()n L x 是否也更加靠近被逼近的函数.Runge 给出的一个例子就是极著名并富有启发性的.设区间[-1,1]上的函数 ()2 1 125f x x = +. 实验内容：考虑区间[-1,1]的一个等距划分，分点为 21,0,1,2,,,i i x i n n =-+ = 则拉格朗日插值多项式为 ()()201 125n n i i i L x l x x ==+∑ . 其中，(),0,1,2,,i l x i n =是n 次Lagrange 插值基函数. 实验要求： （1） 选择不断增大的分点数目2,3, n =，画出原函数()f x 及插值多项式()n L x 在 [-1,1]上的图像，并比较分析实验结果. （2） 选择其他的函数，例如定义在区间[-5,5]上的函数 ()()4 ,arctan ,1x h x g x x x = =+ 重复上面的实验看结果如何. 解：matlab 程序代码
实验结果： f（x）结果如下
h（x）结果如下 g（x）结果如下
结果分析：适当提高插值多项式的次数，可以提高逼近的精度，但次数太高反而会产生不良效果。主要是次数越高，计算工作量大，积累的误差也大；在整个区间上做高次多项式，但局部插值节点处的值有微笑偏差时，可能会影响整个区间上函数值的很大变化，使计算很不稳定。从上图可以看出，高次插值不准确。 实验3.1 编制以函数{} n k k x =为基的多项式最小拟合程序，并对表3.11中的数据作3次多项式最小二 乘拟合. 表3.11 取权数1i w ≡，求拟合曲线0 n k k k x ?α * *== ∑中的参数{}k α、平方误差2δ，并作离散数据 {},i i x y 的拟合函数()y x ?*=的图形. 解：matlab 程序代码 实验结果：
实验一分析太阳黑子数序列 一、实验目的：了解时间序列分析的基本步骤，熟悉SAS/ETS软件使用方法。 二、实验内容：分析太阳黑子数序列。 三、实验要求：了解时间序列分析的基本步骤，注意各种语句的输出结果。 四、实验时间：2小时。 五、实验软件：SAS系统。 六、实验步骤 1、开机进入SAS系统。 2、创建名为exp1的SAS数据集，即在窗中输入下列语句： 3、保存此步骤中的程序，供以后分析使用（只需按工具条上的保存按钮然后填写完提问 后就可以把这段程序保存下来即可）。 4、绘数据与时间的关系图，初步识别序列，输入下列程序： ods html; ods listing close; 5、run;提交程序，在graph窗口中观察序列，可以看出此序列是均值平稳序列。
6、识别模型，输入如下程序。 7、提交程序，观察输出结果。初步识别序列为AR(2)模型。 8、估计和诊断。输入如下程序： 9、提交程序，观察输出结果。假设通过了白噪声检验，且模型合理，则进行预测。 10、进行预测，输入如下程序： 11、提交程序，观察输出结果。
12、退出SAS系统，关闭计算机。总程序： data exp1; infile "D:\"; input a1 @@;
year=intnx('year','1jan1742'd,_n_-1); format year year4.; ; proc print;run; ods html; ods listing close; proc gplot data=exp1 ; symbol i=spline v=dot h=1 cv=red ci=green w=1; plot a1*year/autovref lvref=2 cframe=yellow cvref=black ; title "太阳黑子数序列"; run; proc arima data=exp1; identify var=a1 nlag=24 minic p=(0:5) q=(0:5); estimate p=3; forecast lead=6 interval=year id=year out=out; run; proc print data=out; run; 选取拟合模型的规则: 1.模型显著有效(残差检验为白噪声)
实验报告一 题目：非线性方程求解 摘要：非线性方程的解析解通常很难给出，因此线性方程的数值解法就尤为重要。本实验采用两种常见的求解方法二分法和Newton法及改进的Newton法。 前言：（目的和意义） 掌握二分法与Newton法的基本原理和应用。 数学原理： 对于一个非线性方程的数值解法很多。在此介绍两种最常见的方法：二分法和Newton法。 对于二分法，其数学实质就是说对于给定的待求解的方程f(x)，其在[a,b]上连续，f(a)f(b)<0，且f(x)在[a,b]内仅有一个实根x*，取区间中点c，若，则c恰为其根，否则根据f(a)f(c)<0是否成立判断根在区间[a,c]和[c,b]中的哪一个，从而得出新区间，仍称为[a,b]。重复运行计算，直至满足精度为止。这就是二分法的计算思想。
Newton法通常预先要给出一个猜测初值x0，然后根据其迭代公式 产生逼近解x*的迭代数列{x k}，这就是Newton法的思想。当x0接近x*时收敛很快，但是当x0选择不好时，可能会发散，因此初值的选取很重要。另外，若将该迭代公式改进为 其中r为要求的方程的根的重数，这就是改进的Newton法，当求解已知重数的方程的根时，在同种条件下其收敛速度要比Newton法快的多。 程序设计： 本实验采用Matlab的M文件编写。其中待求解的方程写成function的方式，如下 function y=f(x); y=-x*x-sin(x); 写成如上形式即可，下面给出主程序。 二分法源程序： clear %%%给定求解区间 b=1.5; a=0;
%%%误差 R=1; k=0;%迭代次数初值 while (R>5e-6) ; c=(a+b)/2; if f12(a)*f12(c)>0; a=c; else b=c; end R=b-a;%求出误差 k=k+1; end x=c%给出解 Newton法及改进的Newton法源程序：clear %%%% 输入函数 f=input('请输入需要求解函数>>','s') %%%求解f(x)的导数 df=diff(f);
90-08年人民消费能力分析 一、问题提出 改革开放以来中国经济飞速发展，GDP连续超过德国、日本，现以成为世界上第二大经济体，人民生活水平不断提高，但受金融危机的影响，近几年来物价持续上涨，本月CPI创历史新高，人民的消费能力是否随着GDP的增加而增加呢？本文以中国经济年鉴中的“人民消费支出构成”的数据为依据利用统计软件SAS 进行了相关分析。数据如下 食品衣着居住家庭设备用品及服务交通通讯文教娱乐用品及服务医疗保健其他商品及服务 1990 58.8000 7.7700 17.3400 5.2900 1.4400 5.3700 3.2500 0.7400 1995 58.6200 6.8500 13.9100 5.2300 2.5800 7.8100 3.2400 1.7600 2000 49.1300 5.7500 15.4700 4.5200 5.5800 11.1800 5.2400 3.1400 2005 45.4800 5.8100 14.4900 4.3600 9.5900 11.5600 6.5800 2.1300 2007 43.0800 6.0000 17.8000 4.6300 10.1900 9.4800 6.5200 2.3000 2008 43.6700 5.7900 18.5400 4.7500 9.8400 8.5900 6.7200 2.0900 二、问题分析 1、通过对消费种类进行主成分分析判断人民的消费情况。 2、对主成分标准化后在分析各年的消费能力排名。 三、解决问题 3.1 SAS程序： data examp4_4; input id x1-x8; cards; 1990 58.8000 7.7700 17.3400 5.2900 1.4400 5.3700 3.2500 0.7400 1995 58.6200 6.8500 13.9100 5.2300 2.5800 7.8100 3.2400 1.7600 2000 49.1300 5.7500 15.4700 4.5200 5.5800 11.1800 5.2400 3.1400 2005 45.4800 5.8100 14.4900 4.3600 9.5900 11.5600 6.5800 2.1300 2007 43.0800 6.0000 17.8000 4.6300 10.1900 9.4800 6.5200 2.3000 2008 43.6700 5.7900 18.5400 4.7500 9.8400 8.5900 6.7200 2.0900 ; run; proc corr cov nosimple data=examp4_4; var x1-x8; run; proc princomp data=examp4_4 out=bb; var x1-x8; run; data score1; /*以下程序是对各年按第一主成分得分进行排名并打印结果*/ set bb; keep id prin1;
浙江万里学院实验报告 专业班级：信管111 姓名：学号：实验日期：2014.05.06
五、实验中遇到的问题及相应的解决方案 如果需要删除已经生成的单据或发票，必须先删除凭证，然后在“应收单审核”窗口中取消审核操作，通过执行“应收单审核\应收单列表”命令，在“应收单列表”窗口中删除。 存货核算系统必须执行正常单据记账后，才能确认销售成本，并生成结转销售成本凭证。 有时候自己填写的发货单不能保存，是由于之前做采购管理的时候采购的数量不够，这时候要么在做销售的时候少卖一些，要么再去采购一些。 业务单据中录入项目不能选择业务单据（比如采购入库单、其他出库单等）在表体中已经有项目，但是在添置业务单据时，不能对表体中的项目进行选择录入。只能对项目编号进行选择。在基础设置/单据设计中，打开业务单据，然后增加“项目编码”，保存即可。 六、实验心得 亲身实践后，我的总体感觉是，很多实验前面进行了很多单据的相关操作后，为的就是最后生成一张凭证，实验操作的过程进展的也并不总是那么顺利，只要稍有失误，如漏了其中某一步骤，就导致最后凭证不能生成。所以感觉业务员操作还是很不好做的，更不用说是开发这个系统的人了，进行每一步操作都要谨慎小心。在上课的过程中，不敢有一丝丝的大意，只要在刚开始的过程中有一个小地方没有做完全，就会让稍后的业务采购中有一部分实现不了。每一堂实验课都动手按书本上的步骤实际操作。在实验过程中还是会犯一些小错误，总是需要老师的帮助。 在这段时间中，我对这个系统有了一定的了解。在实验中学习理论知识，这使我对理论中的ERP有了更直观的认识。然而，我总是会存在一些小问题。总是莫名其妙的会出现一些错误，总会让我很无可奈何。还好有同学和老师愿意帮助我。我觉得我在这个学习中学到很多实践课上学不到的内容，感觉十分的充实。
A卷 一、单项选择题（每题1分） 1．借贷记帐法起源于（ B ）。 A、英国 B、意大利 C、美国 D、希腊 2．提取无形资产减值准备金这一做法体现的原则是（ C ）。 A、配比原则 B、重要性原则 C、谨慎原则 D、客观性原则 3．在借贷记帐法下，帐户的贷方登记（ B ）。 A、资产的增加 B、负债的增加 C、费用的增加 D、收入的减少 4．定期对外提供会计报表，这是（ D ）的要求。 A、会计主体 B、货币计量 C、持续经营 D、会计分期 5．会计科目是（ A ） A、账户的名称 B、账簿的名称 C、报表项目的名称 D、会计要素的名称 6．“四柱清册”是我国古代会计模式，四柱中的“实在”是指（ D ）。 A、期初结存 B、本期收入 C、本期付出 D、期末结存 7．（ B ）是记录经济业务，明确经济责任，作为记账依据的书面证明。 A、记账凭证 B、会计凭证 C、汇总记账凭证 D、原始凭证 8．下列不属于资产类账户的是（ D ）。 A、预付账款 B、原材料 C、应收账款 D、实收资本 9．按照权责发生制原则，下述各项目中属于本年收入的是（ D ）。 A、收到上年所销产品的货款85000元 B、预收下年度仓库租金4500元 C、预付下年度财产保险费6000元 D、本年销售产品一批，价款20000元，货款将于下年收到 10．原始凭证按取得的来源不同，可分为（ B ）。 A、一次凭证和累计凭证 B、自制原始凭证和外来原始凭证 C、原始凭证和汇总原始凭证 D、收款凭证、付款凭证和转账凭证 11．下列支出中属于资本性支出的是（ C ）。 A、材料采购支出 B、工资费用支出 C、购置设备支出 D、广告宣传费支出 12．复式记账法是指任何一笔经济业务都必须用相等的金额在两个或两个以上的有关账户中（ D ）。 A、一个记增加，另一个记减少 B、两个都记增加
专业班级：金融106姓名：周吉利1222朱宁宁1224杨程琤1212周孟杰1207实验日期：2012.3.27 浙江万里学院实验报告 课程名称：2011/2012学年第二学期统计实验 实验名称：备择实验专业班级：金融105-106姓名：叶美君1219胡志晖1206黄世杰1208崔 迦楠1175 实验日期：2012.3.29 成绩： 教师：
专业班级：金融106姓名：周吉利1222朱宁宁1224杨程琤1212周孟杰1207实验日期：2012.3.27 一、实验目的：统计分析的目的在于研究总体特征。但是，由于各种各样的原因，我们能够得到的往往只能是从总体中随机抽取的一部分观察对象，他们构成了样本，只有通过对样本的研究，我们才能对总体的实际情况作出可能的推断。因此描述性统计分析是统计分析的第一步，做好这一步是进行正确统计推断的先决条件。通过描述性统计分析可以大致了解数据的分布类型和特点、数据分布的集中趋势和离散程度，或对数据进行初步的探索性分析（包括检查数据是否有错误，对数据分布特征和规律进行初步观察）。 本试验旨在于：引到学生利用正确的统计方法对数据进行适当的整理和显示， 描述并探索出数据内在的数量规律性，掌握统计思想，培养学生学习统计学的兴趣，为继续学习推断统计方法及应用各种统计方法解决实际问题打下必要而坚实的基础。 二、实验内容: 1.表 2.7为某班级16位学生的身高数据，对其进行频数分析，并对实验报告作出说明。 表2.7 某班16位学生的身高数据 学号性别身高（cm ）学号性别身高（cm ） 1 M 170 9 M 150 2 F 17 3 10 M 157 3 F 169 11 F 177 4 M 15 5 12 M 160 5 F 174 13 F 169 6 F 178 14 M 154 7 M 156 15 F 172 8 F 171 16 F 180 三、实验过程: 1、输入某班级16位学生的身高数据。 2、然后选择分析，描述统计，频率，并选择统计量。
一、判断题 （ X ） 1. 支持价格是政府规定的某种产品的最高价格。 （√ ） 2. 经济分析中一般假定消费者的目标是效用最大化。 （√） 3.完全竞争市场的含义:指不受任何阻碍和干扰的市场结构. （√ ） 4. 在一条直线型的需求曲线上每一点的需求弹性系数都不一样。 （√ ） 5. 消费者均衡的实现条件是消费者花在每一元钱上的商品的边际效用都相等。 （ X ） 6. 在同一平面内，两条无差异曲线是可以相交的。 （√ ） 7.需求规律的含义是, 在其他条件不变的情况下，某商品的需求量与价格之间成反方向变动。 （ X ） 8. 规模经济和边际收益递减规律所研究的是同一个问题。 （×） 9. 基尼系数越大，收入分配越平均。 （√ ） 10、在一个国家或家庭中,食物支出在收入中所占比例随着收入的增加而减少。 （√ ） 11. 边际成本曲线总是交于平均成本曲线的最低点。 （√） 12. 完全竞争市场中厂商所面对的需求曲线是一条水平线。 （√） 13. 一个博弈中存在纳什均衡，但不一定存在占优均衡。 （√） 14.边际替代率是消费者在获得相同的满足程度时,每增加一种商品的数量与放弃的另一种商品的数量之比。 （√） 15. 微观经济学的中心理论是价格理论。 （ X ） 16. 经济学按研究内容的不同可分为实证经济学和规范经济学。 （√） 17. 需求弹性这一概念表示的是需求变化对影响需求的因素变化反应程度 （√） 18. 消费者均衡表示既定收入下效用最大化。 （ X ） 19. 经济学按研究内容的不同可分为实证经济学和规范经济学。 （ X ） 20. 平均产量曲线可以和边际产量曲线在任何一点上相 交。 （√） 21. 在长期分析中，固定成本与可变成本的划分是不存在的。 （√） 22. 吉芬商品的需求曲线与正常商品的需求曲线是不相同的。 （√ ） 23. 任何经济决策都必须考虑到机会成本。 （√ ） 24、指经济变量之间存在函数关系时,弹性表示因变量的变动率对自变量的变动 率的反应程度。需求价格弹性指价格变动的比率所引起的需求量变动的比率，即 求量变动对价格变动的反应程度。 （√ ） 25. 在完全竞争市场中每一厂商长期均衡时都没有超额利润。 （X ） 26. 平均产量曲线可以和边际产量曲线在任何一点上相交。 （√ ） 27. 寡头垄断市场的最大特征就是寡头厂商之间的相互依赖性。（√） 28. 微观经济学的研究对象是包括单个消费者、单个生产者、单个市
内容包括： 实验题目1：算法的数值稳定性实验 实验题目2：LU分解实验 实验题目3：三次样条插值外推样条实验 实验题目4：第二类Fredholm 积分方程实验实验题目5：M级显式R_K法
实验题目：算法的数值稳定性实验 实验内容：计算积分()1 0()d 1515n x I n x a x ==+? (n=1,2,…,20) 易得到下面递推公式 ()()1 1I n aI n n =--+ 并有估计式 ()() ()() 1 1 111I n a n a n << +++ 计算方法： 算法一:采用下面递推公式计算: ()()1 1I n aI n n =--+ ()1,2,,20 n = 取初值()116 0ln ln 15a I a +== 算法二: 采用下面递推公式计算: ()()111I n I n a n ??-= -+???? ()20,19,,1 n =
结果分析：（分析哪个好哪个不好，原因是什么） 我觉得算法二比较好, 原因一:根据式 ()() ()() 1 1 111I n a n a n << +++得知,I(n)不可能小于 零,而算法一的计算结果有部分结果小于零。原因二:对算法一记初始误差 ε0=/I 0-I(0)/>0; 则εn =/I n -I(n)/=a/I n-1-I(n-1)/=a n *ε0 由此可知,当n=20时, ε20把ε0放大了a 20倍，其结果造成严重的。 而对于算法二^ ^ 11n n a εε-= ,…, ^ ^ 01 n n a εε=,尽管有初始误差^ 20ε,但随着计算的进程,这个误差的影响不断减小。 附：源程序：（把源程序附上） 算法一程序: >> format long >> a=15;I=log(16/15); for n=1:20 n I=-a*I+1/n end 算法二程序: >> format long >> a=15;I=31/10080; >> for n=20:-1:1 n I I=1/a*(-I+1/n); End
浙江万里学院生物与环境学院 化学工程实验技术实验报告 实验名称：乙酸乙酯皂化反应 姓名成绩 班级学号 同组姓名实验日期 指导教师签字批改日期年月日
一、实验预习（30分） 1．实验装置预习（10分）＿＿＿＿＿年＿＿＿＿月＿＿＿＿日 指导教师＿＿＿＿＿＿（签字）成绩 2．实验仿真预习（10分）＿＿＿＿＿年＿＿＿＿月＿＿＿＿日 指导教师＿＿＿＿＿＿（签字）成绩 3．预习报告（10分） 指导教师＿＿＿＿＿＿（签字）成绩 （1）实验目的 1．用电导率仪测定乙酸乙酯皂化反应进程中的电导率。 2．掌握用图解法求二级反应的速率常数，并计算该反应的活化能。 3．学会使用电导率仪和超级恒温水槽。 （2）实验原理 乙酸乙酯皂化反应是个二级反应，其反应方程式为 CH3COOC2H5+Na++OH-→CH3COO-+Na++C2H5OH 当乙酸乙酯与氢氧化钠溶液的起始浓度相同时，如均为a，则反应速率表示为 （1）式中，x为时间t时反应物消耗掉的浓度，k为反应速率常数。将上式积分得 （2） 起始浓度a为已知，因此只要由实验测得不同时间t时
的x值，以对t作图，应得一直线，从直线的斜率便可求出k值。 乙酸乙酯皂化反应中，参加导电的离子有OH-、Na+和CH3COO-，由于反应体系是很稀的水溶液，可认为CH3COONa是全部电离的，因此，反应前后Na+的浓度不变，随着反应的进行，仅仅是导电能力很强的OH-离子逐渐被导电能力弱的CH3COO-离子所取代，致使溶液的电导逐渐减小，因此可用电导率仪测量皂化反应进程中电导率随时间的变化，从而达到跟踪反应物浓度随时间变化的目的。 令G0为t=0时溶液的电导，G t为时间t时混合溶液的电导，G∞为t=∞（反应完毕）时溶液的电导。则稀溶液中，电导值的减少量与CH3COO-浓度成正比，设K为比例常数，则 由此可得 所以（2）式中的a-x和x可以用溶液相应的电导表示，将其代入（2）式得： 重新排列得: （3） 因此，只要测不同时间溶液的电导值G t和起始溶液的电导值G0，然后 以G t对作图应得一直线，直线的斜率为，由此便求出某温 度下的反应速率常数k值。由电导与电导率κ的关系式:G=κ代入（3）式得: （4） 通过实验测定不同时间溶液的电导率κt和起始溶液 的电导率κ0，以κt，对作图，也得一直线，从直线的斜率也可求出反应速率数k值。如果知道不同温度下的反应速率常数k(T2)和k(T1)，根据Arrhenius公式，可计算出该反应的活化能E和反应半衰期。 （5）