第二.三章习题：
2.1为评价家电行业售后服务的质量，随机抽取了由100家庭构成的一个样本。服务质量的
等级分别表示为：A.好；B.较好；C.一般；D.差；E.较差。调查结果如下：
B E
C C A
D C B A E
D A C B C D
E C E E
A D
B
C C A E
D C B
B A
C
D
E A B D D C
C B C E
D B C C B C
D A C B C D
E C E B
B E
C C A
D C B A E
B A
C
D
E A B D D C
A D
B
C C A E
D C B
C B C E
D B C C B C
(1)指出上面的数据属于什么类型；
属于顺序数据。
(2)用Excel制作一张频数分布表；
（3）绘制一张条形图，反映评价等级的分布。
汇总
5
10
15
20
25
30
35
A
B
C D
E
2.2 某行业管理局所属40个企业2002年的产品销售收入数据如下（单位：万元）：
152 124 129 116 100 103 92 95 127 104 105 119 114 115 87 103 118 142 135 125 117 108 105 110 107 137 120 136 117 108 97
88
123
115
119
138
112
146
113
126
(2) 如果按规定：销售收入在125万元以上为先进企业，115万～125万元为良好企业，
105万～115万元为一般企业，105万元以下为落后企业，按先进企业、良好企业、一般企业、落后企业进行分组。 按销售收入分组（万元） 企业数（个）
频率（%） 先进企业 良好企业 一般企业 落后企业 11 11 9 9 27.5 27.5 22.5 22.5 合计
40
100.0
2.3 某百货公司连续40天的商品销售额如下（单位：万元）：
41 25 29 47 38 34
30 38 43 40
46 36 45 37 37 36 45 43 33 44 35 28 46 34 30 37 44 26 38 44 42 36 37 37 49 39 42 32 36 35
根据上面的数据进行适当的分组，编制频数分布表，并绘制直方图。按销售额分组（万元）频数（天）频率（%）
25～30 30～35 35～40 40～45 45～50
4
6
15
9
6
10.0
15.0
37.5
22.5
15.0
合计40 100.0
在图中显示的统计结果中，可以看见输出的内容分为两部分，一部分是数据表示形式，一部分是直方图形式。在数据表部分，显示每个区间中的人数及累计百分率数值。通过该统计结果，我们可以知道，在40天的商品销售额中有27天销售额在35^50万元之间，11人在25~30万元之间，1人在25万及之下。2.4为了确定灯泡的使用寿命（小时），在一批灯泡中随机抽取100只进行测试，所得结果
如下：
700 716 728 719 685 709 691 684 705 718
706 715 712 722 691 708 690 692 707 701
708 729 694 681 695 685 706 661 735 665
668 710 693 697 674 658 698 666 696 698
706 692 691 747 699 682 698 700 710 722
694 690 736 689 696 651 673 749 708 727
688 689 683 685 702 741 698 713 676 702
701 671 718 707 683 717 733 712 683 692
693 697 664 681 721 720 677 679 695 691
713 699 725 726 704 729 703 696 717 688
(1)利用计算机对上面的数据进行排序；
(2)以组距为10进行等距分组，整理成频数分布表，并绘制直方图；
按使用寿命分组（小时）灯泡个数（只）频率（%）
650~660 2 2
660~670 5 5
670~680 6 6
680~690 14 14
690~700 26 26
700~710 18 18
710~720 13 13
720~730 10 10
730~740 3 3
740~750 3 3
合计100 100
2.51997年我国几个主要城市各月份的平均相对湿度数据如下表，试绘制箱线图，并分析
资料来源：《中国统计年鉴1998》，中国统计出版社1998，第10页。
2.6 某百货公司6月份各天的销售额数据如下（单位：万元）：
257 276 297 252 238 310 240 236 265 278 271 292 261 281 301 274 267 280 291 258 272
284
268
303
273
263
322
249
269
295
（1） 计算该百货公司日销售额的均值、中位数和四分位数； x =274.1（万元）；Me =272.5 ；Q L =260.25；Q U =291.25。 （2）计算日销售额的标准差。
17.21 s （万元）
2.7
比较哪个企业的总平均成本高？并分析其原因。
。。甲企业平均成本＝19.41（元），乙企业平均成本＝18.29（元）；原因：尽管两个企业的单位成本相同，但单位成本较低的产品在乙企业的产量中所占比重较大，因此拉低了总平均成本。
2.8 在某地区抽取的120家企业按利润额进行分组，结果如下：
按利润额分组（万元）
企业数（个）
200～300
19
300～400 30 400～500 42 500～600 18 600以上 11 合计
120
计算120家企业利润额的均值和标准差。
。、x =426.67（万元）；48.116=s （万元）
第四章练习题
4.1 一个具有100=n 个观察值的随机样本选自于30=μ、16=σ的总体。试求下列概率的近似值：
⑴ 0.8944； ⑵ 0.0228； ⑶ 0.1292； ⑷ 0.9699。
4.2 美国汽车联合会（AAA ）是一个拥有90个俱乐部的非营利联盟，它对其成员提供旅行、
金融、保险以及与汽车相关的各项服务。1999年5月，AAA 通过对会员调查得知一个4口之家出游中平均每日餐饮和住宿费用大约是213美元（《旅行新闻》Travel News ，1999年5月11日）。假设这个花费的标准差是15美元，并且AAA 所报道的平均每日消费是总体均值。又假设选取49个4口之家，并对其在1999年6月期间的旅行费用进行记录。
⑴ 描述x （样本家庭平均每日餐饮和住宿的消费）的抽样分布。特别说明x 服从怎样
的分布以及x 的均值和方差是什么？ 正态分布, x =213, s=4.5918； ⑵ 对于样本家庭来说平均每日消费大于213美元的概率是什么？大于217美元的概率
呢？在209美元和217美元之间的概率呢？ P= 0.5, p= 0.031, p=0.938。
4.3 技术人员对奶粉装袋过程进行了质量检验。每袋的平均重量标准为406=μ克、标准差
为1.10=σ克。监控这一过程的技术人者每天随机地抽取36袋，并对每袋重量进行测量。现考虑这36袋奶粉所组成样本的平均重量x 。
（1）描述x 的抽样分布，并给出x μ和x σ的值，以及概率分布的形状；
x μ=406, x σ=1.68, 属于正态分布
=0.001
（2） 假设某一天技术人员观察到8.400=x ，这是否意味着装袋过程出现问题了
呢，为什么？
是，因为小概率发生了
4.4 某制造商为击剑运动员生产安全夹克，这些夹克是以剑锋刺入其中时所需的最小力量
（以牛顿为单位）来定级的。如果生产工艺操作正确，则他生产的夹克级别应平均840牛顿，标准差15牛顿。国际击剑管理组织（FIE ）希望这些夹克的最低级别不小于800牛顿。为了检查其生产过程是否正常，某检验人员从生产过程中抽取了50个夹克作为
一个随机样本进行定级，并计算x，即该样本中夹克级别的均值。她假设这个过程的标准差是固定的，但是担心级别均值可能已经发生变化。
⑴如果该生产过程仍旧正常，则x的样本分布为何？
正态分布
⑵假设这个检验人员所抽取样本的级别均值为830牛顿，则如果生产过程正常的话，
样本均值x≤830牛顿的概率是多少？
P(x≤830)约等于0
⑶在检验人员假定生产过程的标准差固定不变时，你对b部分有关当前生产过程的现
状有何看法（即夹克级别均值是否仍为840牛顿）？
不正常
⑷现在假设该生产过程的均值没有变化，但是过程的标准差从15牛顿增加到了45牛
顿。在这种情况下x的抽样分布是什么？当x具有这种分布时，则x≤830牛顿的概率是多少？
正态分布, p(x≤830)=0.06
4.5 某大学为了解学生每天上网的时间，在全校7500名学生中采取不重复抽样方法随机抽取36人，调查他们每天上网的时间，得到下面的数据（单位：小时）：
3.3 3.1 6.2 5.8 2.3
4.1
5.4 4.5 3.2
4.4 2.0
5.4 2.6
6.4 1.8 3.5 5.7 2.3
2.1 1.9 1.2 5.1 4.3 4.2
3.6 0.8 1.5
4.7 1.4 1.2 2.9 3.5 2.4 0.5 3.6 2.5
求该校大学生平均上网时间的置信区间，置信水平分别为90%、95%和99%。
90%时置信区间为；
95%时置信区间；
99%时置信区间。
4.6从两个正态总体中分别抽取两个独立的随机样本，它们的均值和标准差如下表：
来自总体1的样本
来自总体2的样本
141=n 72=n 2.531=x
4.432=x
8.9621=s
0.10222=s
（1） 求21μμ-90%的置信区间；=（1.86,17.74）
（2） 求21μμ-95%的置信区间。=（0.19,19.41）
4.7从两个正态总体中分别抽取两个独立的随机样本，它们的均值和标准差如下表：
来自总体1的样本
来自总体2的样本
251=x
232=x
1621=s
2022=s
（3） 设10021==n n ，求21μμ-95%的置信区间；=2±1.176
（4） 设1021==n n ，2
221σσ=，求21μμ-95%的置信区间；=2±3.986 （5） 设1021==n n ，2
221σσ≠，求21μμ-95%的置信区间；=2±3.986
（6） 设20,1021==n n ，2
221σσ=，求21μμ-95%的置信区间；=2±3.587
（7） 设20,1021==n n ，2
221σσ≠，求21μμ-95%的置信区间。=2±3.364 4.8从两个总体中各抽取一个25021==n n 的独立随机样本，来自总体1的样本比率为
%401=p ，来自总体2的样本比率为%302=p 。
（8） 构造21ππ-90%的置信区间；=10%±6.98% （9） 构造21ππ-95%的置信区间。=10%±8.32%
4.9生产工序的方差是共需质量的一个重要度量。当方差较大时，需要对共需进行改进以减
构造两个总体方差比2221σσ95%的置信区间
置信区间为（4.06,14.35）
4.10假定两个总体的标准差分别为：121=σ，152=σ，若要求误差范围不超过5，相应的置信水平为95%，假定21n n =，估计两个总体均值之差21μμ-时所需的样本容量为多
大？ 57
生物统计学课后习题解答-李春喜汇总
第一章概论 解释以下概念：总体、个体、样本、样本容量、变量、参数、统计数、效应、互作、随机误差、系统误差、准确性、精确性。 第二章试验资料的整理与特征数的计算习题 2.1 某地 100 例 30 ～ 40 岁健康男子血清总胆固醇(mol · L -1 ) 测定结果如下： 4.77 3.37 6.14 3.95 3.56 4.23 4.31 4.71 5.69 4.12 4.56 4.37 5.39 6.30 5.21 7.22 5.54 3.93 5.21 6.51 5.18 5.77 4.79 5.12 5.20 5.10 4.70 4.74 3.50 4.69 4.38 4.89 6.25 5.32 4.50 4.63 3.61 4.44 4.43 4.25 4.03 5.85 4.09 3.35 4.08 4.79 5.30 4.97 3.18 3.97 5.16 5.10 5.85 4.79 5.34 4.24 4.32 4.77 6.36 6.38 4.88 5.55 3.04 4.55 3.35 4.87 4.17 5.85 5.16 5.09 4.52 4.38 4.31 4.58 5.72 6.55 4.76 4.61 4.17 4.03 4.47 3.40 3.91 2.70 4.60 4.09 5.96 5.48 4.40 4.55 5.38 3.89 4.60 4.47 3.64 4.34 5.18 6.14 3.24 4.90 计算平均数、标准差和变异系数。 【答案】=4.7398, s=0.866, CV =18.27 ％ 2.2 试计算下列两个玉米品种 10 个果穗长度 (cm) 的标准差和变异系数，并解释所得结果。 24 号： 19 ， 21 ， 20 ， 20 ， 18 ， 19 ， 22 ， 21 ， 21 ， 19 ； 金皇后： 16 ， 21 ， 24 ， 15 ， 26 ， 18 ， 20 ， 19 ， 22 ， 19 。 【答案】 1 =20, s 1 =1.247, CV 1 =6.235% ； 2 =20, s 2 =3.400, CV 2 =17.0% 。 2.3 某海水养殖场进行贻贝单养和贻贝与海带混养的对比试验，收获时各随机抽取 50 绳测其毛重(kg) ，结果分别如下：
数值分析上机实验报告 选题：曲线拟合的最小二乘法 指导老师： 专业： 学号： 姓名：
课题八曲线拟合的最小二乘法 一、问题提出 从随机的数据中找出其规律性，给出其近似表达式的问题，在生产实践和科学实验中大量存在，通常利用数据的最小二乘法求得拟合曲线。 在某冶炼过程中，根据统计数据的含碳量与时间关系，试求含碳量y 与时间t 的拟合曲线。 二、要求 1、用最小二乘法进行曲线拟合； 2、近似解析表达式为()33221t a t a t a t ++=?； 3、打印出拟合函数()t ?，并打印出()j t ?与()j t y 的误差，12,,2,1 =j ； 4、另外选取一个近似表达式，尝试拟合效果的比较； 5、*绘制出曲线拟合图*。 三、目的和意义 1、掌握曲线拟合的最小二乘法； 2、最小二乘法亦可用于解超定线代数方程组； 3、探索拟合函数的选择与拟合精度间的关系。 四、计算公式 对于给定的测量数据(x i ,f i )(i=1,2,…，n )，设函数分布为 ∑==m j j j x a x y 0)()(? 特别的，取)(x j ?为多项式 j j x x =)(? (j=0, 1,…，m )
则根据最小二乘法原理，可以构造泛函 ∑∑==-=n i m j i j j i m x a f a a a H 1 10))((),,,(? 令 0=??k a H (k=0, 1,…，m ) 则可以得到法方程 ???? ??????? ?=????????????????????????),(),(),(),(),(),(),(),(),(),(),(),(1010101111000100m m m m m m m m f f f a a a ????????????????????? 求该解方程组，则可以得到解m a a a ,,,10 ，因此可得到数据的最小二乘解 ∑=≈m j j j x a x f 0)()(? 曲线拟合：实际工作中，变量间未必都有线性关系，如服药后血药浓度与时间的关系；疾病疗效与疗程长短的关系；毒物剂量与致死率的关系等常呈曲线关系。曲线拟合是指选择适当的曲线类型来拟合观测数据，并用拟合的曲线方程分析两变量间的关系。 五、结构程序设计 在程序结构方面主要是按照顺序结构进行设计，在进行曲线的拟合时，为了进行比较，在程序设计中，直接调用了最小二乘法的拟合函数polyfit ，并且依次调用了plot 、figure 、hold on 函数进行图象的绘制，最后调用了一个绝对值函数abs 用于计算拟合函数与原有数据的误差，进行拟合效果的比较。
计算方法上机报告 姓名： 学号： 班级： 上课班级：
说明： 本次上机实验使用的编程语言是Matlab 语言，编译环境为MATLAB 7.11.0，运行平台为Windows 7。 1. 对以下和式计算： ∑ ∞ ? ?? ??+-+-+-+=0681581482184161n n n n S n ，要求： ① 若只需保留11个有效数字，该如何进行计算； ② 若要保留30个有效数字，则又将如何进行计算； （1） 算法思想 1、根据精度要求估计所加的项数，可以使用后验误差估计，通项为： 1421114 16818485861681 n n n a n n n n n ε??= ---<< ?+++++??； 2、为了保证计算结果的准确性，写程序时，从后向前计算； 3、使用Matlab 时，可以使用以下函数控制位数： digits(位数)或vpa(变量，精度为数) （2）算法结构 1. ;0=s ?? ? ??+-+-+-+= 681581482184161n n n n t n ； 2. for 0,1,2,,n i =??? if 10m t -≤ end; 3. for ,1,2,,0n i i i =--??? ;s s t =+
（3）Matlab源程序 clear; %清除工作空间变量 clc; %清除命令窗口命令 m=input('请输入有效数字的位数m='); %输入有效数字的位数 s=0; for n=0:50 t=(1/16^n)*(4/(8*n+1)-2/(8*n+4)-1/(8*n+5)-1/(8*n+6)); if t<=10^(-m) %判断通项与精度的关系break; end end; fprintf('需要将n值加到n=%d\n',n-1); %需要将n值加到的数值 for i=n-1:-1:0 t=(1/16^i)*(4/(8*i+1)-2/(8*i+4)-1/(8*i+5)-1/(8*i+6)); s=s+t; %求和运算 end s=vpa(s,m) %控制s的精度 （4）结果与分析 当保留11位有效数字时，需要将n值加到n=7， s =3.1415926536； 当保留30位有效数字时，需要将n值加到n=22, s =3.14159265358979323846264338328。 通过上面的实验结果可以看出，通过从后往前计算，这种算法很好的保证了计算结果要求保留的准确数字位数的要求。
统计学作业2 单项选择题 第1题某地区有10万人口，共有80个医院。平均每个医院要服务1250人，这个指标是（）。 A、平均指标 B、强度相对指标 C、总量指标 D、发展水平指标 答案：B 第2题某企业2002年工业总产值比1992年增长了3倍，则该公司1992-2002年间工业总产值平均增长速度为（） A、11.61% B、14.87% C、13.43% D、16.65% 答案：A 第3题某工业企业的某种产品成本，第一季度是连续下降的。1月份产量750件，单位成本20元；2月份产量1000件，单位成本18元；3月份产量1500件，单位成本15元。则第一季度的平均成本为（）。 A、17.67 B、17.54 C、17.08 D、16.83 答案：C 第4题已知4个水果商店苹果的单价和销售额，要求计算4个商店苹果的平均单价，应该采用（）。 A、简单算术平均数 B、加权算术平均数 C、加权调和平均数 D、几何平均数
答案：C 第5题如果分配数列把频数换成频率，那么方差（）。 A、不变 B、增大 C、减小 D、无法预期变化 答案：A 第6题某厂5年的销售收入如下：200万、220万、250万、300万、320万，则平均增长量为（）。 A、120/5 B、120/4 C、320/200的开5次方 D、320/200的开4次方 答案：B 第7题直接反映总体规模大小的指标是（）。 A、平均指标 B、相对指标 C、总量指标 D、变异指标 答案：C 第8题计算结构相对指标时，总体各部分数值与总体数值对比求得的比重之和（）。 A、小于100% B、大于100% C、等于100% D、小于或大于100% 答案：C 多项选择题
2.用下列方法求方程e^x+10x-2=0的近似根，要求误差不超过5*10的负4次方，并比较计算量 （1）二分法 （局部，大图不太看得清，故后面两小题都用局部截图） （2）迭代法
（3）牛顿法 顺序消元法 #include for(k=p+1;k 昆明理工大学工科研究生《数值分析》上机实验 学院：材料科学与工程学院 专业：材料物理与化学 学号：2011230024 姓名: 郑录 任课教师：胡杰 P277-E1 1.已知矩阵A= 10787 7565 86109 75910 ?? ?? ?? ?? ?? ??,B= 23456 44567 03678 00289 00010 ?? ?? ?? ?? ?? ?? ?? ?? ,错误！未找到引用源。 = 11/21/31/41/51/6 1/21/31/41/51/61/7 1/31/41/51/61/71/8 1/41/51/61/71/81/9 1/51/61/71/81/91/10 1/61/71/81/91/101/11?????????????????? （1）用MA TLAB函数“eig”求矩阵全部特征值。 （2）用基本QR算法求全部特征值（可用MA TLAB函数“qr”实现矩阵的QR分解）。解：MA TLAB程序如下： 求矩阵A的特征值： clear; A=[10 7 8 7;7 5 6 5;8 6 10 9;7 5 9 10]; E=eig(A) 输出结果： 求矩阵B的特征值： clear; B=[2 3 4 5 6;4 4 5 6 7;0 3 6 7 8;0 0 2 8 9;0 0 0 1 0]; E=eig(B) 输出结果： 求矩阵错误！未找到引用源。的特征值： clear; 错误！未找到引用源。=[1 1/2 1/3 1/4 1/5 1/6; 1/2 1/3 1/4 1/5 1/6 1/7; 1/3 1/4 1/5 1/6 1/7 1/8; 1/4 1/5 1/6 1/7 1/8 1/9;1/5 1/6 1/7 1/8 1/9 1/10; 1/6 1/7 1/8 1/9 1/10 1/11]; E=eig(错误！未找到引用源。) 输出结果： （2）A= 10 7877565861097 5 9 10 第一步：A0=hess(A);[Q0,R0]=qr(A0);A1=R0*Q0 返回得到： 第二部：[Q1,R1]=qr(A1);A2=R1*Q1 《数值计算方法》上机实验报告华北电力大学 实验名称数值il?算方法》上机实验课程名称数值计算方法专业班级：电力实08学生姓名：李超然学号:200801001008 成绩: 指导教师:郝育黔老师实验日期：2010年04月华北电力大学实验报告数值计算方法上机实验报吿一. 各算法的算法原理及计算机程序框图1、牛顿法求解非线性方程 *对于非线性方程，若已知根的一个近似值，将在处展开成一阶 xxfx ()0, fx ()xkk 泰勒公式 "f 0 / 2 八八,fxfxfxxxxx 0 0 0 0 0 kkkk2! 忽略高次项，有 ，fxfxfxxx 0 ()()()，，, kkk 右端是直线方程，用这个直线方程来近似非线性方程。将非线性方程的 **根代入，即fx ()0, X ，* fxfxxx 0 0 0 0, ,, kkk fx 0 fx 0 0, 解出 fX 0 *k XX,, k' fx 0 k 水将右端取为，则是比更接近于的近似值，即xxxxk, Ik, Ik fx ()k 八XX, Ikk* fx()k 这就是牛顿迭代公式。 ,2,计算机程序框图:，见, ,3,输入变量、输出变量说明: X输入变量:迭代初值，迭代精度，迭代最大次数，\0 输出变量:当前迭代次数，当前迭代值xkl ,4,具体算例及求解结果: 2/16 华北电力大学实验报吿 开始 读入 l>k /fx()0?,0 fx 0 Oxx，，01* fx ()0 XX,,，?10 kk, ,1，kN, ?xx, 10 输出迭代输出X输出奇异标志1失败标志 ,3,输入变量、输出变量说明: 结束 例：导出计算的牛顿迭代公式，并il ?算。（课本P39例2-16） 115cc （0）, 求解结果: 10. 750000 10.723837 10. 723805 10. 723805 2、列主元素消去法求解线性方程组，1,算法原理: 高斯消去法是利用现行方程组初等变换中的一种变换，即用一个不为零的数乘 -个 方程后加只另一个方程，使方程组变成同解的上三角方程组，然后再自下而上 对上三角 3/16 华北电力大学实验报告方程组求解。 列选主元是当高斯消元到第步时，从列的以下（包括）的各元素中选出绝 aakkkkkk 对值最大的，然后通过行交换将其交换到的位置上。交换系数矩阵中的 两行（包括常ekk 数项），只相当于两个方程的位置交换了，因此，列选主元不影响求解的结 ,2,计算机程序框图:，见下页, 输入变量:系数矩阵元素，常向量元素baiji 输出变量:解向量元素bbb,,12n 计算方法第四次上机报告 2.用欧拉方法解初值 y’=10x(1-y) 0<=x<=1 Y(0)=0 取步长h=0.1,保留5位有效数字，并与准确解相比较 分析：该题目考察欧拉方法解初值问题 程序如下： function Heun(a,b,y0,n) h=(b-a)/n;x=a:h:b; y=y0*ones(1,n+1); for j=2:n+1 yp=y(j-1)+h*f(x(j-1),y(j-1)); yc=y(j-1)+h*f(x(j),yp); y(j)=1/2*(yp+yc); end for k=1:n+1 fprintf('x[%d]=%f\ty[%d]=%f\n',k-1,x(k),k-1,y(k)); end function z=f(xx,yy) z=10*xx*(1-yy); 运行结果： >> Heun(0,1,0,10) x[0]=0.000000 y[0]=0.000000 x[1]=0.100000 y[1]=0.050000 x[2]=0.200000 y[2]=0.183000 x[3]=0.300000 y[3]=0.362740 x[4]=0.400000 y[4]=0.547545 x[5]=0.500000 y[5]=0.705905 x[6]=0.600000 y[6]=0.823543 x[7]=0.700000 y[7]=0.901184 x[8]=0.800000 y[8]=0.947627 x[9]=0.900000 y[9]=0.973290 x[10]=1.000000 y[10]=0.986645 >> 分析： 该结果与准确结果相比比较接近，但是有一定的误差。 6.用四阶龙格—库塔公式解第三题中的初值问题，取步长h=0.2，保留五位有效数字。 题目目的分析： 该题考查四阶龙格-库塔方法和改进欧拉方法求解精确度问题。 程序： 改进欧拉法： function Heun(a,b,y0,n) h=(b-a)/n;x=a:h:b; y=y0*ones(1,n+1); for j=2:n+1 yp=y(j-1)+h*f(x(j-1),y(j-1)); yc=y(j-1)+h*f(x(j),yp); y(j)=1/2*(yp+yc); end for k=1:n+1 fprintf('x[%d]=%f\ty[%d]=%f\n',k-1,x(k),k-1,y(k)); end 统计学论文错误 一、科研设计 1.论文题目：ＨＯ-１在阿尔兹海默病额叶及海马组织中的表达及对神经的保护作用 摘要：方法选择2012年7月至2013年12月该院收集的30例神经活检脑标本作为研究对象。结合病理诊断选择其中肯定的散发型AD病例记为AD组，共16 例，其余14例有中枢神经系统等病症且无明显脑病理变化者作为对照组。以免疫组化法以及免疫荧光染色法分别检测HO-1在额叶及海马组织中的表达，对比两组HO-1IHC阳性情况以及两组额叶及海马神经元中HO-1表达IA值水平。结论HO-1在AD疾病患者额叶及海马组织中均有表达，且可能是通过启动机体内源性的神经保护机制从而对大脑发挥一定的保护作用。 错误之处： 在科研设计中，病例与对照选择的基本原则是所调查的病例足以代表总体中该病的病例，对照足以代表产生病例的总体。 此调查选取了来源于某一所医院住院部的病例在一定时期内得神经活检标本作为病例。较易进行，省经费；但是带有选择性，容易产生选择偏倚，仅反映该机构的病人特点，而不是全人群该病的特点。 对照的选择关系到病例对照研究的成败。该实验从本医院的其他病人中选对照，即在选择病例的医院内选择有中枢神经系统等病症且无明显脑病理变化者作为对照组，但并未提及是何病种，也并知道该病种是否对研究阿尔兹海默病是否产生影响，病种以愈复杂愈好。这样比较方便，且这种对照的应答率和信息的质量均较高。 建议： 1.得出结论是可以给出范围，如某地区某医院HO-1在AD疾病患者额叶及海马组织中均有表达，且可能是通过启动机体内源性的神经保护机制从而对大脑发挥一定的保护作用。 2.可在方法中给出对照组疾病对研究并未产生统计学作用。 3.可以同时选取两种对照，即从一般人口中选择对照，又可以从住院病人中选择对照。研究结果一致，则能增强评价的依据。如结果不一致，则需分析其原因，可能有偏倚。 4.所选病例必须是患同一种疾病的病人，诊断标准、病例的年龄、性别、种族、职业等，选择时要有一个明确的规定。如本例可选取某排除标准将具有其他痴呆相关的神经系统疾病、抑郁症等排除入组等。 2.论文题目：补肾益智颗粒联合盐酸多奈哌齐片治疗阿尔兹海默病的临床疗效研究 摘要：目的:观察补肾益智颗粒联合盐酸多奈哌齐片治疗阿尔兹海默病的临床疗效。方法:选取2012年11月至2013年11月期间,在广西中医药大学第一附属医院脑病科及干部科门诊或住院治疗符合诊断标准的136例阿尔兹海默病患者,按就诊先后随机分为治疗组和对照组,每组各68例。对照组单用盐酸多奈哌齐片治疗,治疗组在对照组的基础上加用补肾益智颗粒治疗,4周为一疗程,共治疗24周。治疗前及治疗后采用临床总体印象量表(CGI)评价患者的整体情况,采用简易智能精神状态量表(MMSE)评价患者的认知功能及痴呆的严重程度,采用AD评定量表的认知次级量表(ADAS-Cog)评价患者认知功能及精神行为,采用日常生活能力量表(ADL)评价患者的日常生活自理能力,采用神经精神量表(NPI)评价患者的精神状态,采用中医证候评分评价患者症状、体征的变化。所有数据均采用SPSS17.0 《生物统计学》教学大纲 课程名称：生物统计学 课程类型：范围选修课-基础课 学时：56学时，3.5学分 适用对象：农学、植物保护、生物技术、生物科学、草业科学等本科专业先修课程：高等数学、线性代数、概率论与数理统计、植物学、植物生理学、遗传学等课程 一、课程性质、目的与任务以及对先开课程要求 统计学是论述收集、分析并解释数字信息的科学，生物统计学则是一门运用统计学的原理和方法，研究生物学数据资料的一般统计学。统计方法是现代生物学研究不可缺少的工具。正确的统计分析能够帮助我们正确认识事物客观存在的规律性。 概率论与数理统计等先开课程的重点是讲述没有量纲或单位抽象的数量规律，为生物学科应用这些规律打基础。 二、教学重点及难点 本课程教学的全过程可以看成是一个生物信息搜集、处理、分析，从而提炼新的生物信息的过程。教学重点是通过生物现象的数量观察、对比、归纳和分析，揭示那些困惑费解的生物学问题，从偶然性的剖析中，发现事物的必然性，指导生物科学的理论和实践。本课程的难点是概念较多、理论抽象、系统严密、实践性强、公式复杂、符号繁多、计算量大，因此，教学安排上除精讲48学时外，有针对性的安排上机操作8学时。 三、与其他课程关系 生物统计学与数学有密切关系，现代统计学用到了较多的数学知识，研究理论生物统计学的人需要有较深的数学功底，应用统计方法的人也应具备良好的数学基础。统计学又是一门应用性很强的学科，几乎生物学科所有的门类都要研究和分析数据，掌握生物学类学科专业基础课和专业课程知识有利于对统计分析的结果做出合理的解释和分析。 四、教学内容、学时分配及基本要求 绪论（1学时） 基本要求：理解什么是统计？什么是统计学；统计数据与统计学的关系，描述统计与推断统计内涵；统计方法能解决生物学科中哪些问题，了解生物 统计学的产生与发展。 计算方法 一、填空题 1.假定x ≤1，用泰勒多项式?+??+++=! !212n x x x e n x ,计算e x 的值，若要求截断误差不超过0.005，则n=_5___ 2. 解 方 程 03432 3=-+x － x x 的牛顿迭代公式 )463/()343(121121311+--+--=------k k k k k k k x x x x x x x 3.一阶常微分方程初值问题 ?????= ='y x y y x f y 0 0)() ,(,其改进的欧拉方法格式为)],(),([21 1 1 y x y x y y i i i i i i f f h +++++= 4.解三对角线方程组的计算方法称为追赶法或回代法 5. 数值求解初值问题的四阶龙格——库塔公式的局部截断误差为o(h 5 ) 6.在ALGOL 中，简单算术表达式y x 3 + 的写法为x+y ↑3 7.循环语句分为离散型循环,步长型循环,当型循环. 8.函数)(x f 在[a,b]上的一次（线性）插值函数= )(x l )()(b f a b a x a f b a b x --+-- 9.在实际进行插值时插值时，将插值范围分为若干段，然后在每个分段上使用低阶插值————如线性插值和抛物插值，这就是所谓分段插值法 10、数值计算中，误差主要来源于模型误差、观测误差、截断误差和舍入误差。 11、电子计算机的结构大体上可分为输入设备 、 存储器、运算器、控制器、 输出设备 五个主要部分。 12、算式2 cos sin 2x x x +在ALGOL 中写为))2cos()(sin(2↑+↑x x x 。 13、ALGOL 算法语言的基本符号分为 字母 、 数字 、 逻辑值、 定义符四大 计算方法实验报告 班级： 学号： 姓名： 成绩： 1 舍入误差及稳定性 一、实验目的 （1）通过上机编程，复习巩固以前所学程序设计语言及上机操作指令； （2）通过上机计算，了解舍入误差所引起的数值不稳定性 二、实验内容 1、用两种不同的顺序计算10000 21n n -=∑，分析其误差的变化 2、已知连分数() 1 01223//(.../)n n a f b b a b a a b =+ +++，利用下面的算法计算f ： 1 1 ,i n n i i i a d b d b d ++==+ (1,2,...,0 i n n =-- 0f d = 写一程序，读入011,,,...,,,...,,n n n b b b a a 计算并打印f 3、给出一个有效的算法和一个无效的算法计算积分 1 041 n n x y dx x =+? (0,1,...,1 n = 4、设2 2 11N N j S j == -∑ ，已知其精确值为1311221N N ?? -- ?+?? （1）编制按从大到小的顺序计算N S 的程序 （2）编制按从小到大的顺序计算N S 的程序 （3）按两种顺序分别计算10001000030000,,,S S S 并指出有效位数 三、实验步骤、程序设计、实验结果及分析 1、用两种不同的顺序计算10000 2 1n n -=∑，分析其误差的变化 （1）实验步骤： 分别从1~10000和从10000~1两种顺序进行计算，应包含的头文件有stdio.h 和math.h （2）程序设计： a.顺序计算 #include 《统计学原理》作业（二） （第四章） 一、判断题 1、总体单位总量和总体标志总量是固定不变的，不能互相变换。（×） 2、相对指标都是用无名数形式表现出来的。（×） 3、能计算总量指标的总体必须是有限总体。（×） 4、按人口平均的粮食产量是一个平均数。（×） 5、在特定条件下，加权算术平均数等于简单算术平均数。（√） 6、用总体部分数值与总体全部数值对比求得的相对指标。说明总体内部的组成状况，这个相对指标是比例相对指标。（×） 7、国民收入中积累额与消费额之比为1：3，这是一个比较相对指标。（×） 8、总量指标和平均指标反映了现象总体的规模和一般水平。但掩盖了总体各单位的差异情况，因此通过这两个指标不能全面认识总体的特征。（√） 9、用相对指标分子资料作权数计算平均数应采用加权算术平均法。（×） 10、标志变异指标数值越大，说明总体中各单位标志值的变异程度就越大，则平均指标的代表性就越小。（√） 二、单项选择 1、总量指标数值大小（A） A、随总体范围扩大而增大 B、随总体范围扩大而减小 C、随总体范围缩小而增大 D、与总体范围大小无关 2、直接反映总体规模大小的指标是（C） A、平均指标 B、相对指标 C、总量指标 D、变异指标 3、总量指标按其反映的时间状况不同可以分为（D） A、数量指标和质量指标 B、实物指标和价值指标 C、总体单位总量和总体标志总量 D、时期指标和时点指标 4、不同时点的指标数值（B） A、具有可加性 B、不具有可加性 C、可加或可减 D、都不对 5、由反映总体各单位数量特征的标志值汇总得出的指标是（B） A、总体单位总量 B、总体标志总量 C、质量指标 D、相对指标 1. 资料可以分为哪几类？它们有何区别？ 答：（1）资料一般可以分为数量性状资料、质量性状资料、半定量资料三大类，其中 数量性状资料又包括计量资料和计数资料。 （2）区别：数量性状资料是能够以量测或计数的方式获得的资料，质量性状资料是只能观察而不能直接测量的资料，半定量资料既有计数资料的特点又有程度或量的不同。联系 三种不同类型的资料有时可根据研究目的和统计方法的要求将一种类型资料转化成另一种类型的资料。 2. 为什么要对资料进行整理？对于计量资料，整理的基本步骤是什么？ 答：（1）由调查或试验收集来的原始资料往往是零乱的，无规律可循。只有通过统计整理，才能发现其内部的联系和规律性，从而揭示事物的本质。资料整理是进行统计分析的基础。（2）计量资料整理的基本步骤包括：①求全距，全距即为资料中最大值与最小值之差。②确定组数，一般根据样本含量及资料的变动范围大小确定组数。③确定组距，通常根据等距离分组的原则，组距等于全距除以组数。④确定组限和组中值，各组的最大值为组上限，最小值为组下限 每一组的中点值称为组中值。⑤归组划线计数，作次数分布表。 3. 在对计量资料进行整理时，为什么第一组的组中值以接近或等于资料中的最小值为好？答：在对计量资料进行整理时，第一组的组中值以接近或等于资料中的最小值可以避免第一组中观察值过多的情况，同时也确保资料中最小值不会遗漏。 4. 统计表与统计图有何用途？常用统计表、统计图有哪些？列统计表、绘统计图时，应注意什么？ 答：（1）统计表用表格形式来表示数量关系 统计图用几何图形来表示数量关 系。用统计表和统计图可以把研究对象的特征、内部构成、相互关系等简明、形 象地表达出来，便于比较分析。 （2）常用的统计图有长条图、圆图、线图、直方图和折线图等。 （3）常用的统计表有简单表和复合表两大类。 （4）列统计表的注意事项： ①标题要简明扼要、准确地说明表的内容，有时须注明时间、地点。 ②标目分横标目和纵标目两项，横标目列在表的左侧用以表示被说明事物的主要标志，纵标目列在表的上端说明横标目各统计指标内容，并注明计算单位。 ③数字一律用阿拉伯数字，数字小数点对齐，小数位数一致，无数字的用“—”表示，数字是“0”的须写“0”。 ④表的上下两条边线略粗、纵、横标目间及合计用细线分开 表的左右边线可以省去表的左上角一般不用斜线。 （5）绘统计图的注意事项： ①标题简明扼要并列于图的下方。 ②纵、横两轴应有刻度，注明单位。 ③横轴由左至右，纵轴由上而下，数值由小到大，图形长宽比例约为5：4或6：5。 ④图中需用不同颜色或线条表示不同事物时应有图例说明。 5. 生物统计中常用的平均数有几种？各在什么情况下应用？ 答：生物统计中常用的平均数有算术平均数、几何平均数、调和平均数、中位数和众数。算术平均数较常用，简称平均数，当资料呈正态分布时可用算术平均数描述其中心位置。几何均数主要应用于畜牧、水产业的动态分析，畜禽疾病及药物效价的统计分析，如畜禽、水产养殖的增长率，抗体的滴度，药物的效价，畜禽疾病的潜伏期等。调和均数主要用于反映畜群不同阶段的平均增长率或畜群不同规模的平均规模。当所获得的数据资料呈偏态分布时中位数的代表性优于算术平均数。众数也适用于资料呈偏态分布的情况。 1． 设?+=1 05dx x x I n n ， （1） 由递推公式n I I n n 1 51+-=-，从0I 的几个近似值出发，计算20I ； 解：易得：0I =ln6-ln5=0.1823, 程序为： I=0.182; for n=1:20 I=(-5)*I+1/n; end I 输出结果为：20I = -3.0666e+010 （2） 粗糙估计20I ，用n I I n n 51 5111+- =--，计算0I ； 因为 0095.05 6 0079.01020 201 020 ≈<<≈??dx x I dx x 所以取0087.0)0095.00079.0(2 1 20=+= I 程序为：I=0.0087; for n=1:20 I=(-1/5)*I+1/(5*n); end I 0I = 0.0083 （3） 分析结果的可靠性及产生此现象的原因(重点分析原因)。 首先分析两种递推式的误差；设第一递推式中开始时的误差为000I I E '-=，递推过程的舍入误差不计。并记n n n I I E '-=，则有01)5(5E E E n n n -==-=-Λ。因为=20E 20020)5(I E >>-，所此递推式不可靠。而在第二种递推式中n n E E E )5 1(5110-==-=Λ，误差在缩小， 所以此递推式是可靠的。出现以上运行结果的主要原因是在构造递推式过程中，考虑误差是否得到控制， 即算法是否数值稳定。 2． 求方程0210=-+x e x 的近似根，要求4 1105-+?<-k k x x ，并比较计算量。 （1） 在[0，1]上用二分法； 程序：a=0;b=1.0; while abs(b-a)>5*1e-4 c=(b+a)/2; 统计学作业 Document number：WTWYT-WYWY-BTGTT-YTTYU-2018GT 《统计学》作业 适用层次：专科培养类型：经管类 专业班级：姓名：学号： 作业要求：作业可打印，答案要求手写；考试时交作业。 第三章综合指标 1、某厂按计划规定，第一季度的单位产品成本比去年同期降低10%，实际执行结果， 单位产品成本较去年同期降低8％。问该厂第一季度产品单位成本计划的完成程度如 何 2、某乡甲、乙两个村的粮食生产情况如下： 试分别计算甲、乙两个村的平均亩产。根据表列资料及计算结果，比较分析哪一个村 的生产经营管理工作做得好，并简述作出这一结论的理由。 3、甲、乙两单位工人的生产资料如下： 试分析：（1）哪个单位工人的生产水平高 （2）哪个单位工人的生产水平整齐 第四章时间数列 1、某企业2010年第一季度职工人数及产值资料如下： 要求：（1）计算第一季度的月平均劳动生产率。 （2）计算第一季度的劳动生产率。 2、某炼钢厂连续5年钢产量资料如下： 要求：（1）平均发展水平（2）平均增减量（3）平均发展速度（4）平均增减速度3、某企业2002～2008年的产品销售额资料如表所示。 要求：用最小平方法配合直线趋势方程，并预测2009年的趋势值。 第五章统计指数 1、试根据某企业三种产品总产值和产量动态的资料，计算三种产品产量总指数，以及由于产量增加使企业所增加的总产值。 2、某企业三种产品总产值、出厂价格增长速度资料如下表： 要求：计算出厂价格总指数和由于价格变化而增加的总产值； 3、某商场销售三种商品，销售量及销售价格如下： 要求用总量指标变动中的因素分析法，对销售额的变动作出分析。 生物统计学学习心得 一、《生物统计学》这一门课。你学到什么？谈谈你学习这一门课的心得体会。 （一）、《生物统计学》这门课，首先，我不仅学到了很多生物统计方面的基础知识、基本概念和相关的应用，还学习了如何设计试验。 在第一章，我学了统计数据的收集与整理。首先学习的是总体与样本的概念，统计学研究的核心问题是如何通过样本推断总体，因此，总体与样本是生物统计学中的两个最基本概念。总体是我们研究的全部对象。构成总体的一个研究单位称为个体。样本是总体的一部分，样本内包含的个体数目称为样本含量。接着学习了数据类型及频数分布。生物统计学中经常遇到的数据有两种类型，一种是连续型数据，指与某种标准做比较所得到的数据，采用变量的方法进行分析。另一种是离散型数据，指由记录不同类别的个体的数目所得到的数据，采用属性的方法进行分析。最后学习了样本的几个特征数,平均数、标准差、方差。 在第二章，我学了概率和概率分布。概率是事件所固有的，且不随人的主观意识而改变。总体分布是建立在概率这一概念基础之上的，因此在研究总体分布之前首先应对概率的基本知识有所了解。试验的每一最基本的结果称为基本事件，指不能再分的事件。复合事件指由若干个基本事件组合而成的事件。概率的基本运算法则包括概率加法法则、条件概率、概率乘法法则、独立事件。概率分布包括离散型概率分布和连续型概率分布。 在第三章，我学了几种常见的概率分布律。首先学了二项分布，二项分布的基本情况是：设有一随机试验，每次试验都有两种不同的结果，如成功的（事件A）和失败的（事件A’）；生男孩（事件A）和生女孩（事件A’）。显然这两种可能的结果是互不相容的，独立地将此试验重复做n次，求在n次试验中，一种结果出现y次的概率。接着学了泊松分布、超几何分布、负二项分布、正态分布、指数分布等。 在第四章，我学了抽样分布。首先学了从一个正态总体中抽取的样本统计量的分布，学了一些基本概念，如标准误差、样本标准误差、自由度、查表。然后学了从两个正态总体中抽取的样本统计量的分布，包括标准差已知时两个平均数的和与差的分布、标准未知但相等时两个平均数的和与差的分布、两个样本方差比的分布----F分布。 在第五章，我学了统计推断。对总体做统计推断可以通过两条途径进行，一是首先对所估计的总体提出一个假设，称为统计假设检验，二是通过样本统计量估计总体参数，称为总体参数估计。首先学习单个样本的统计假设检验，检验的基本步骤：1.提出假设。2.构造并计算检验统计量：利用原假设所提供的信息，而且抽样分布已知。3.确定否定域（临界值）：根据小概率事件原理，比较检验统计量和临界值的关系，确定其落在否定域还是接受域。主要学了t检验，u检验、x2检验。接着学了两个样本的差异显著性检验，包括两个方差的检验----F检验，标准差已知时两个平均数间差异显著性的检验，标准差未知但相等时，两平均数之间差异显著性的检验，标准差未知且可能不等时两平均数之间差异显著性的检验，配对数据的显著性检验-----配对数据的t检验，二项分布数据的显著性检验。 在第六章，我学了参数估计，即由样本统计量估计总体参数。估计量是估计总体参数的统计量，一个好的估计量应该满足三个条件：无偏性、有效性、相容性。对总体参数的估计，可分为点估计和区间估计。区间估计是指在一定概率保证下指出总体参数的可能范围，所给出的可能范围叫置信区间，本章我学习了μ的置信区间、σ的置信区间、平均数差的置信区间、配对数据的置信区间、标准差比的置信区间二项分布总体的置信区间。 在第七章，我学了拟合优度检验，拟合优度检验是用来检验实际观测数与依照某种假设或模型计算出来的理论数之间的一致性，以便判断该假设或模型是否与观测数相配合。做拟合优度检验一般需一下各步：1.对数据进行分组。2.计算理论数Ti。3分别合并两个尾区的理论数。4.零假设。5.计算出x2与x2临界值（查附表6）做比较。 第一次&第二次上机作业 上机作业： 1.在Matlab上执行：>> 5.1-5-0.1和>> 1.5-1-0.5 给出执行结果，并简要分析一下产生现象的原因。 解：执行结果如下： 在Matlab中，小数值很难用二进制进行描述。由于计算精度的影响，相近两数相减会出现误差。 2.（课本181页第一题） 解：（1）n=0时，积分得I0=ln6-ln5，编写如下图代码 从以上代码显示的结果可以看出，I 20的近似值为0.7465 （2）I I =∫I I 5+I 10dx,可得∫I I 610dx ≤∫I I 5+I 10dx ≤∫I I 510dx,得 16(I +1)≤I I ≤15(I +1),则有1126≤I 20≤1105, 取I 20=1 105 ,以此逆序估算I 0。代码段及结果如下图所示 （3）从I20估计的过程更为可靠。首先根据积分得表达式是可知，被积函数随着n的增大，其所围面积应当是逐步减小的，即积分值应是随着n的递增二单调减小的，（1）中输出的值不满足这一条件，（2）满足。设I I表示I I的近似值，I I-I I=(?5)I(I0?I0)(根据递推公式可以导出此式),可以看出，随着n的增大，误差也在增大，所以顺序估计时，算法不稳定性逐渐增大，逆序估计情况则刚好相反，误差不断减小，算法逐渐趋于稳定。 2.（课本181页第二题） （1）上机代码如图所示 求得近似根为0.09058 （2）上机代码如图所示 得近似根为0.09064； （3）牛顿法上机代码如下 计算所得近似解为0.09091 第三次上机作业上机作业181页第四题 线性方程组为 [1.13483.8326 0.53011.7875 1.16513.4017 2.53301.5435 3.4129 4.9317 1.23714.9998 8.76431.3142 10.67210.0147 ][ I1 I2 I3 I4 ]=[ 9.5342 6.3941 18.4231 16.9237 ] （1）顺序消元法 A=[1.1348,3.8326,1.1651,3.4017;0.5301,1.7875,2.5330,1.5435; 3.4129, 4.9317,8.7643,1.3142;1.2371,4.9998,10.6721,0.0147]; b=[9.5342;6.3941;18.4231;16.9237]; 上机代码（函数部分）如下 function [b] = gaus( A,b )%用b返回方程组的解 B=[A,b]; n=length(b); RA=rank(A); RB=rank(B);数值分析上机作业
《数值计算方法》上机实验报告
计算方法上机作业
统计学作业完整版本
《生物统计学》教学大纲
计算方法试题库讲解
计算方法上机实习题大作业(实验报告).
统计学原理作业2答案(新)
生物统计学课堂作业及答案
(完整版)数值计算方法上机实习题答案
统计学作业
生物统计学学习心得
计算方法上机作业集合
相关主题
文本预览