西安交通大学计算方法A实验报告
实验一 矩阵的分解
一、实验目的
掌握矩阵的分解原理和一般方法，学会利用矩阵分解直接求解线性方程组。
二、实验内容
求矩阵()
2020
=ij
A α?的T
LDL 分解与Cholesky 分解，其中
,min(,),ij i i j
i j i j α=?=?
≠?
。
三、问题分析
1． Cholesky 分解
Cholesky 分解是针对被分解矩阵为对称正定的情况给出的。
分解步骤如下：
11g =1111/y b g =，1111i i g g α=
2i n = ；
DO
2j n =
jj g = IF 0jj g < STOP ，JUMP TO (5)
DO 1i j n =+
1
1j ij ik kj k ij jj
g g g g α-=??- ?
?
?=∑
ji ij g g =
1
1j i ik k k i jj b g y y g -=??- ?
?
?=∑
END DO
END DO
2． T
LDL 分解
T LDL 分解是针对Cholesky 分解中的开平方运算进行的改进。
分解步骤如下：
11i i r α=，1111/i i r r r =，11y b = 1i n =
DO
2i n =
DO
j i n =
1
1i ij ij ik kj k r l r α-=??=- ???
∑
/ji ij ii l r r =
1
1i i i ik k k y b l b -=??=- ???
∑
END DO
END DO
四、matlab 求解
分别写出T
LDL 分解和Cholesky 分解的函数程序gaijinsqrt.m 和.cholesky m ，调用格
式如下：
1. [index,x,r]=gaijinsqrt(A,b) 参数说明：
A 和b 分别是线性代数方程组Ax =b 的系数矩阵和右端向量；输出x 为解向量。
[index,x,g]=Cholesky(A,b)
参数说明：
A 和b 分别是线性代数方程组Ax =b 的系数矩阵和右端向量；输出x 为解向量。 然后写出主程序2homework .m 如下： %生成矩阵A
A=zeros(20,20);
for i=1:20
for j=1:20
if i~=j
if i>j
A(i,j)=j; else A(i,j)=i;
end
else A(i,j)=i;
end
end
end
b=ones(20,1);
for i=1:10
b(i)=i; b(21-i)=i;
end
%LDLt 分解
[index1,x1,r]=gaijinsqrt(A,b) %Cholesky 分解
[index2,x2,g]=Cholesky(A,b)
gaijinsqrt.m 和.cholesky m 见附件
五、实验结果
选取b=[1 2 3 4 5 6 7 8 9 10 10 9 8 7 6 5 4 3 2 1]。实验结果如下：
202011=11D ?????????????????
，2020000010=10001110L ????????????????? ，2020
111=111111G ???????
????
??
????
实验二 三对角方程组Tx f =的求解
一、 实验目的
掌握三对角方程组Tx f =求解的原理和方法。
二、 实验内容
求三对角方程组Tx f =的解，其中：
1 1 0 0 01
2 1 0 00 1
3 1 00 0 1
4 10 0 0 1 5T ????????=????????
，()12345,,,,T x x x x x x =，()3,8,15,24,29T
f =。
三、 问题分析
追赶法的算法组织如下：
(1) 输入三对角矩阵T 和右端向量f ；
(2) 将T 压缩为四个一维数组{}{}{}{}i i i i a b c d 、、、，将分解矩阵压缩为三个一维数组
{}{}{}i i i l r y 、、
(3) 对T 做Crout 分解（也可以用Doolittle 分解）导出追赶法的计算步骤如下：
1
1111,d l b y l ==
DO
2i n =
()11111, ,i i i i i i i i i i i i
r c l l b a r y d a y l -----==-=-
END DO
(4) 回代求解x
n n x y =
DO 11i n =-
1i i i i x x r x +=-
END DO
(5) 停止，输出结果
四、matlab求解
zhuigan m，调用格式如下：
编写追赶法的函数程序.
[x,r1]=zhuigan(A,b)
参数说明：
A和b分别是线性代数方程组的系数矩阵和右端向量，x是求得的解向量，r1是分解矩阵。
homework3m如下：
编写主程序.
T=[1 1 0 0 0;1 2 1 0 0;0 1 3 1 0;0 0 1 4 1;0 0 0 1 5];
f=[3;8;15;24;29];
%调用函数zhuigan.m
[x,r]=zhuigan(T,f);
zhuigan m见附件。
函数文件.
五、实验结果
x 。
(1.0000 2.0000 3.0000 4.0000 5.0000)
实验三 线性代数方程组的迭代求解
一、 实验目的
掌握线性代数方程组的迭代求解的原理和基本方法Jacobi 迭代法和Gauss Seidel -法。
二、 实验内容
用Jacobi 迭代法和Gauss Seidel -迭代法求解下列方程组：
1234 0.78 -0.02 -0.12 -0.140.76-0.02 0.86 -0.04 0.06 0.08-0.12 -0.04 -0.72 -0.08 1.12-0.14 0.06 -0.08 0.74 0.68x x x x ??????????????????=????????????????
?? 并比较迭代次数，使误差小于3
10-。
三、 问题分析
将原始线性代数方程组Ax b =改写为()x x ?= 的形式，其中()x ? 为x 的矩阵函数。
于是可以得到迭代格式：(1)()()k k x x ?+= ，此即为Jacobi 迭代法的迭代格式。如果在计算(1)k x
+ 时，将已经算出的分量立即代换()
k x 对应分量，则得到Gauss Seidel -迭代法的迭代格式。
1. Jacobi 迭代法的算法组织如下：
(1) 给出迭代格式(1)()
()k k x x ?+=
(2) 给出迭代初始向量0x 、允许误差ε和最大迭代次数N
(3) 按照迭代格式(1)()()k k x x ?+= 进行迭代，直至达满足迭代停止条件
(1)()1max k k i i i n
x x ε+≤≤-<
(4) 停止，输出结果
2. Gauss Seidel -迭代法的算法组织如下：
(1) 给出迭代格式(1)()
()k k x x ?+=
(2) 给出迭代初始向量0x 、允许误差ε和最大迭代次数N
(3) 按照迭代格式(1)()()k k x x ?+= ，并且将已经算出的分量立即代换()k x 对应分量进行迭
代，直至达满足迭代停止条件(1)
()1max k k i
i i n
x x ε+≤≤-<
(4) 停止，输出结果
四、 matlab 求解
编写Jacobi 迭代法的函数.Jacobi m 和Gauss Seidel -迭代法函数.Gauss Seidel m -，他们的调用格式如下：
1、 [x,k]=jacobi(A,b,x0,N,epsilon)
参数说明：
A 和b 分别是线性代数方程组的系数矩阵和右端向量，x0是迭代初始向量，N 是最
大迭代次数，epsilon 是迭代允许误差，x 是求得的解向量，k 是程序结束时迭代次数。 2、 [x,k]=gaosisidel(A,b,x0,N,epsilon)
参数说明：
A 和b 分别是线性代数方程组的系数矩阵和右端向量，x0是迭代初始向量，N 是最
大迭代系数，epsilon 是迭代允许误差，x 是求得的解向量，k 是程序结束时迭代次数。
编写主程序homework3.m 如下：
A=[0.78 -0.02 -0.12 -0.14;-0.02 0.86 -0.04 0.06;-0.12 -0.04 -0.72 -0.08;-0.14 0.06 -0.08 0.74]; b=[0.76;0.08;1.12;0.68]; x0=[0 0 0 0]; N=1e+3;
epsilon=1e-3;
%Jacobi 迭代法
[x1,k1]=jacobi(A,b,x0,N,epsilon); %Gauss-Sidel 迭代法
[x2,k2]=gaosisidel(A,b,x0,N,epsilon);
五、 实验结果
1(0.8570 -0.0326 -1.7953 0.8896) , =6x k = 2(0.8569 -0.0326 -1.7954 0.8896) , 4x k ==
实验四 多项式插值
一、实验目的
掌握多项式插值的原理和基本方法。
二、实验内容
已知2
1()(11)125f x x x =
-≤≤+，对5,10,20n =
1、 计算函数()f x 在点2
1,(0,1,2,,)i x i i n n
=-+= 处的值()i f x 2、 求插值数据点
(){},
(0,1,2,,)i i
x y i n = 的Newton 插值多项式()n
N x 和三次
样条插值多项式()n S x 3、 对5,20n =，计算2
1,(110,9099)100
k x k k =-+
= 和相应的函数值()
,(),()k n k n k k y f x N x S x =
4、 计算()()max n n k
k k
E N y N x =-，()()
max n
n k k
k
E S y
S x =-，解释所得到结果。
三、问题分析
1. Newton 插值
Newton 插值多项式是实质是Lagrange 插值多项式的零次式和一次式的推广，一般形
式为：
()()()()()()()010201011n n n N x c c x x c x x x x c x x x x x x -=+-+--++---
2. 三次样条插值
在节点0121()i n n x x x x x x b -<<<<= 分成的每个小区间1[,]i i x x -上利用其节点处
的二次导数值（即弯矩值）进行线性插值，再在此区间上积分两次，利用节点处的函数值求
得两个积分常数，从而得到三次样条函数的M 表达式。
要确定三次样条插值函数，用到()S x 的光滑性，利用三次样条函数在节点处的一阶导
数连续得到i M 满足1n -个方程：
112,1,2,,1i i i i i i M M M d i n μλ-+++==-
其中
[]111111116
,1,()6,,.
i i i i i i i i i i i i i i i i i i i i i
h h y y y y d f x x x h h h h h h h h μλμ+--+++++--=
==-=-=+++
三弯矩方程组只有1n -个方程，不能确定1n +未知量(0,1,2,)i M i n = 。为此，可
以给出区间[,]a b 的端点处的一阶导数值''0n y y 、。
由()'
'
00S x y =和()'
'
n n S x y =可以导出
01012,2n n n M M d M M d -+=+=
其中''101001166
(),()n n n n n n
y y y y d y d y h h h h ---=
-=- 联合先前得到的三弯矩方程组和两个补充方程得到恰定的三弯矩方程组。它是严格对角
占优的三对角方程组，可以用追赶法求解。
四、matlab 求解
1. Newton 插值
编写Newton 插值的函数Ninterpolation.m ，其调用格式如下：
[y,y0,ddt]=Ninterpolation(xx,yy,x0)
参数说明：
xx,yy 为输入数据的横坐标向量和纵坐标向量；x0为待求点处的横坐标；y 为插值后
得到的多项式函数，此函数会给出插值多项式的系数；y0为对应于x0的函数值；ddt 为求得的差商表。 2. 三次样条插值
编写求解给出边界点的函数值和一阶导数值的三次样条插值的函数.myspline2m ，其
调用格式如下：
[M,y0]=myspline2(x,y,m1,mn,x0)
参数说明：
x,y 为输入数据的横坐标向量和纵坐标向量；x0为待求点处的横坐标；m1,mn 为边界
点的一阶导数值；y0为对应于x0的函数值；M 为各节点处的弯矩值组成的弯矩向量。
编写主程序homework4.m 如下： n=20;h=2/n;
x=zeros(n,1);y=zeros(n,1); for i=1:n+1
x(i)=-1+(i-1)*h;
y(i)=1/(1+25*x(i)^2); end
%输出第一问答案
fprintf('\n 函数f(x)在各点处的函数值为：\n'); y %
x0=zeros(20,1); for i=1:10
x0(i)=-1+2*i/100; x0(21-i)=-x0(i); end
y0=1./(1+25*x0.^2);
%Newton 插值多项式
[y1,y10]=Ninterpolation(x,y,x0); hold on
%三次样条插值多项式
[y2,y20]=myspline2(x,y,0,50/26^2,x0); hold on
%第四问的答案
EN=max(abs(y0-y10)); fprintf('\nE(N)=%d\n',EN); ES=max(abs(y0-y20));
fprintf('\nE(S)=%d\n',ES);
五、实验结果
1.
计算函数()f x 在点2
1,(0,1,2,,)i x i i n n
=-+
= 处的值()i f x ，()i f x 按i 的大小从左到右，再从上到下依次排列如下： 20n =时
0.03850.04710.05880.07550.10000.13790.20000.30770.50000.80001.00000.80000.50000.30770.20000.13790.10000.07550.05880.04710.0385
10n =时
0.03850.05880.10000.20000.50001.00000.50000.20000.10000.05880.0385
5n =时
0.03850.10000.50000.50000.10000.0385
2. 求得插值数据点
(){},
(0,1,2,,)i i
x y i n = 的Newton 插值多项式()n
N x 和三次样
条插值多项式()n S x 如下： (1) 20n =
201816201412108642()260178.6309717141012094.874479971639177.4108481442944.96343925+757287.0927749245249.28392450749317.54291293926119.2199494214470.84622675944524.14347962478121
N x x x x x x x x x x x =-++
--+-+-+
当10.9x -≤≤-时
3220() 5.24996414.36516712.980442-3.826777S x x x x =---
当0.90.8x -≤≤-时
3220() 1.722742+4.461140+3.963234 1.256325S x x x x =+
当0.80.7x -≤≤-时
3220()0.0750370.506647+0.7996400.412700S x x x x =++
当0.70.6x -≤≤-时
3220()0.9737802.394008+2.120793+0.720969S x x x x =+
当0.60.5x -≤≤-时
3220() 1.552441 3.435598 2.7457470.845960S x x x x =+++
当0.50.4x -≤≤-时
3220() 3.551653+6.434415+4.245155+1.095861S x x x x =
当0.40.3x -≤≤-时
3220() 5.726359+9.044062+5.289014+1.235043S x x x x =
当0.30.2x -≤≤-时
3220()12.534955+15.171799+7.127335+1.418875S x x x x =
当0.20.1x -≤≤-时
3220()32.78925612.022728+1.688430+1.056281S x x x x =--
当0.10x -≤≤时
3220()89.07023928.907022+1.000000S x x x =--
当00.1x ≤≤时
3220()89.070210-28.907022+1.000000S x x x =
当0.10.2x ≤≤时
3220()32.78939712.022728 1.688424+1.056281S x x x x =--
当0.20.3x ≤≤时
3220()12.535491+15.1721557.127411+1.418880S x x x x =--
当0.30.4x ≤≤时
3220() 5.724355+9.042132 5.288404+1.234979S x x x x =--
当0.40.5x ≤≤时
3220() 3.559131+6.443864 4.249097+1.096405S x x x x =--
当0.50.6x ≤≤时
3220() 1.524530 3.391962 2.7231460.842080S x x x x =-+-+
当0.60.7x ≤≤时
3220() 1.077945 2.588109 2.240834+0.745617S x x x x =-+-
当0.70.8x ≤≤时
3220()0.3137120.3343720.1950980.268279S x x x x =--+
当0.80.9x ≤≤时
3220() 3.173574+8.035114 6.890687 2.053769S x x x x =--+
当0.91x ≤≤时
3220()10.66454129.327795+26.7359328.034216S x x x x =--
(2) 10n =
1086
1042()220.941742081448+494.909502262443381.433823529411123.35972850678716.85520361990951
N x x x x x x =--+-+
当10.8x -≤≤-时
3210() 1.271419 3.050923 2.2875890.469624S x x x x =----
当0.80.6x -≤≤-时
3210() 1.325569 3.181847+2.6986270.860034S x x x x =++
当0.60.4x -≤≤-时
3210()0.720275 2.092319+2.0449100.729290S x x x x =++
当0.40.2x -≤≤-时
3210()13.44038917.356456+8.150564 1.543378S x x x x =++
当0.20x -≤≤时
3210()54.48183223.3968770.000102 1.000000S x x x x =---+
当00.2x ≤≤时
3210()54.48183223.3968770.000102 1.000000S x x x x =--+
当0.20.4x ≤≤时
3210()13.46593017.3748458.154447 1.543623S x x x x =-+-+
当0.40.6x ≤≤时
3210()0.623222 1.963596 1.9899470.721690S x x x x =-+-+
当0.60.8x ≤≤时
3210() 1.688240 3.880628 3.1401660.951734S x x x x =-+-+
当0.81x ≤≤时
3210() 2.625051 6.471270+5.141352 1.256671S x x x x =--
(3) 5n =
4210() 1.20192307692308 1.730769230769230.567307692307692N x x x =-+
当10.6x -≤≤-时
325() 2.274939+6.299457+5.77409 1.788040S x x x x =+
当0.60.2x -≤≤-时
325() 3.459433 4.0224120.4190250.549416S x x x x =---+
当0.20.2x -≤≤时
325()0.024330 1.9321540.0009730.577286S x x x x =--+
当0.20.6x ≤≤时
325() 3.362111 3.934823+0.3995600.550584S x x x x =-+
当0.20.6x ≤≤时
325() 1.934313 5.598741 5.320578 1.694612S x x x x =-+-+
3. 计算2
1,(110,9099)100
k x k k =-+
= 和相应的()
,k k y f x =()n k N x ，()n k S x
20n =时
()110()0.04000.04160.04330.04510.04710.04910.05130.05360.05610.0588f x =
()
9099()0.05880.05610.05360.05130.04910.04710.04510.04330.04160.0400f x =
()
110()58.238150.864428.662310.33450.0471 3.9541 4.0691 2.6744 1.13530.0588N x =----
()
9099()0.0588 1.1353 2.6744 4.0691 3.95410.047110.334528.662650.864458.2381N x =----
()110()0.03900.04030.04230.04660.04710.04940.05160.05390.05630.0588S x =
()9099()0.05880.05640.05420.05200.04970.04710.04410.04130.03910.0380S x =
5n =时
()110()0.04000.04160.04330.04510.04710.04910.05130.05360.05610.0588f x =
()
9099()0.05880.05610.05360.05130.04910.04710.04510.04330.04160.0400f x =
()
110()0.01370.00690.02360.03660.04600.05220.05530.05550.05300.0481N x =---------
()
9099()0.04810.05300.05550.05530.05220.04600.03660.02360.00690.0137N x =---------
()
110()0.03830.03780.03710.03630.03550.03480.03440.03430.03470.0356S x =
()9099()0.03100.02980.02930.02940.03000.03100.03220.03770.03530.0369S x =
4.
计算()()max n n k
k k
E N y N x =-，()()
max n
n k k
k
E S y
S x =-如下：
(1) 20n =时 ()2058.278125E N =，()200.002512E S = (2) 5n =时 ()50.109194E N =，()50.027848E S =
这是由于已知数据的误差对Newton 插值多项式计算的结果带来了很大的误差。具体
说，当函数值123,,,n y y y y 有误差而变为 123,,,n y y y y 且误差限为ε。 由
于
N e 插值与Lagrange 插值本质上是一样的，则Newton 插值与Lagrange 插值的误差
估计是一样的。由Lagrange 插值多项式()n L x 的计算值 ()n L
x 也会有误差，于是可以导出Newton 插值多项式的误差
() ()
()()()()()()
n
n
n n n n
i
i
i
i
i i n
n
i i i i i i N x N x L x L x l x y l x y l x y y l x ε====-=-=-≤-≤∑∑∑∑
当n 很大时0
()1n
i i l x ε
=∑
， ()n N
x 的误差就很大，于是就出现了前面观察到的现象。
实验五 数值积分
一、实验目的
掌握数值积分的基本原理与方法。
二、实验内容
用Romberg 方法计算定积分：1
sin()x I dx x =
?，(0)1f =使误差不超过5
10-。 三、问题分析
在变步长积分的基础上，可以证明变步长积分的结果
()221
()3
n n n I f T T T =+-
与复化Simpson 求积公式的关系为：
()221
()41
n n n n S I f T T T ==+
-- 同理可以证明
()222
1
41n n n n C S S S =+
-- ()2231
41
n n n n R C C C =+--
如此进行下去即为Romberg 积分方法，以上称为Romberg 递推公式。
四、matlab 求解
算法组织如下：
(1) 输入被积函数及积分区间[,]a b ，设置Romberg 积分表格S 的大小N 和计算精度ε (2) 调用求解变步长积分的函数variestepn.m 计算1T ，即计算Romberg 积分表S 的元素
(1,1)S
(3) 调用求解变步长积分的函数variestepn.m 计算12n T -，作为Romberg 积分表S 的元素
(,1)S n ，然后根据Romberg 递推公式计算元素(,)S n k ，其计算公式如下：
1
(,1)(1,1)
(,)(,1),231
k S n k S n k S n k S n k k n -----=-+
=-
(4) 如果(,1)(1,1)S n n S n n ε----<，则将(,)S n n 作为积分结果，停止；否则1n n =+，
跳转到步骤()3
其中函数variestepn.m 的调用格式为
[F,n]=variestepn(f,a,b,N)
参数说明：
f 为被积函数，其定义时要先用命令syms 定义自变量，例如
syms x;
[F,n]=variestepn(sin(x)/x,0,1,4)
a,b 为积分区间的下、上限；N 表示将积分区间划分为12N -等分；输出F 和n 即为(,1)F S n =。
编写Romberg 积分函数romberg.m ，其调用格式为：
[F,n]=romberg(f,a,b,N,epsilon)
参数说明：
f 为被积函数；N 为Romber
g 积分表格的大小；a,b 为积分区间的下、
上限；epsilon 为计算精度；输出F 和n 即为(,)F S n n =。
编写主程序homework5.m 如下：
close all clear all clc syms x;
epsilon=1e-3; %设置积分区间 a=0;b=1; f=sin(x)/x;
% Romberg 积分
[F,n]=romberg(f,a,b,20,epsilon)
五、实验结果
运行主程序homework5.m ，结果如下：
0.946084673818906,7F =n =
《计算方法》课内实验报告
《计算方法》实验报告 姓名： 班级： 学号： 实验日期: 2011年10月26日
一、实验题目： 数值积分 二、实验目的： 1．熟悉matlab 编写及运行数值计算程序的方法。 2．进一步理解数值积分的基础理论。 3．进一步掌握应用不同的数值积分方法求解给定的积分并给出数据结果及误差分析。 三、实验内容： 1．分别用复合梯形求积公式及复合辛普森求积公式计算积分xdx x ln 10 ? , 要求计算精度达到410-,给出计算结果并比较两种方法的计算节点数. 2．用龙贝格求积方法计算积分dx x x ?+3 021,使误差不超过510-. 3．用3=n 的高斯-勒让德公式计算积分?3 1 sin x e x ,给出计算结果. 4．用辛普森公式(取2==M N ) 计算二重积分.5 .00 5 .00 dydx e x y ? ? - 四、实验结果： 1.（1）复合梯形法: 将区间[a,b]划分为n 等份，分点n k n a b h kh a x k ,2,1,0,,=-=+=在每个区间[1,+k k x x ](k=0,1,2,···n-1)上采用梯形公式，则得 )()]()([2)()(1 11 1 f R x f x f h dx x f dx x f I n n k k k b a n k x x k k ++===∑?∑? -=+-=+ 故)]()(2)([21 1 b f x f a f h T n k k n ++=∑-=称为复合梯形公式 计算步长和划分的区间 Eps=1E-4 h1=sqrt(Eps/abs(-(1-0)/12*1/(2+1))) h1 =0.0600 N1=ceil(1/h1) N1 =17 用复合梯形需要计算17个结点。 复合梯形： function T=trap(f,a,b,n) h=(b-a)/n;
统计西安交大期末考试试题(含答案)
西安交大统计学考试试卷 一、单项选择题（每小题2 分，共20 分） 1.在企业统计中，下列统计标志中属于数量标志的是（C） A、文化程度 B、职业 C、月工资 D、行业 2.下列属于相对数的综合指标有（B ） A、国民收入 B、人均国民收入 C、国内生产净值 D、设备台数 3.有三个企业的年利润额分别是5000 万元、8000 万元和3900 万元，则这句话中有（B）个变量？ A、0 个 B、两个 C、1 个 D、3 个 4.下列变量中属于连续型变量的是（A ） A、身高 B、产品件数 C、企业人数 D、产品品种 5.下列各项中，属于时点指标的有（A ） A、库存额 B、总收入 C、平均收入 D、人均收入 6.典型调查是（B ）确定调查单位的 A、随机 B、主观 C、随意 D 盲目 7.总体标准差未知时总体均值的假设检验要用到（A ）： A、Z 统计量 B、t 统计量 C、统计量 D、X 统计量 8.把样本总体中全部单位数的集合称为（A ） A、样本 B、小总体 C、样本容量 D、总体容量 9.概率的取值范围是p（D ） A、大于1 B、大于－1 C、小于1 D、在0 与1 之间 10.算术平均数的离差之和等于（A ） A、零 B、1 C、－1 D、2 二、多项选择题（每小题2 分，共10 分。每题全部答对才给分，否则不计分） 1.数据的计量尺度包括（ABCD ）： A、定类尺度 B、定序尺度 C、定距尺度 D、定比尺度 E、测量尺度 2.下列属于连续型变量的有（BE ）： A、工人人数 B、商品销售额 C、商品库存额 D、商品库存量 E、总产值 3.测量变量离中趋势的指标有（ABE ） A、极差 B、平均差 C、几何平均数 D、众数 E、标准差 4.在工业企业的设备调查中（BDE ） A、工业企业是调查对象 B、工业企业的所有设备是调查对象 C、每台设备是 填报单位D、每台设备是调查单位E、每个工业企业是填报单位 5.下列平均数中，容易受数列中极端值影响的平均数有（ABC ） A、算术平均数 B、调和平均数 C、几何平均数 D、中位数 E、众数 三、判断题（在正确答案后写“对”，在错误答案后写“错”。每小题1 分，共10 分） 1、“性别”是品质标志。（对）
西安交通大学计算方法B上机试题
1.计算以下和式：01421181 84858616n n S n n n n ∞ =?? =--- ?++++??∑ ，要求： (1)若保留11个有效数字，给出计算结果，并评价计算的算法； (2)若要保留30个有效数字，则又将如何进行计算。 （1）题目分析 该题是对无穷级数求和，因此在使用matlab 进行累加时需要一个累加的终止条件。这里令?? ? ??+-+-+-+= 681581482184161n n n n a n n ，则 ()()1.016 1 6855844864816114851384128698161 681581482184161148113811282984161111<< ? ??? ????? ??++++++???? ????? ??++++++=??? ????? ??+-+-+-+??? ????? ??+-+-+-+=+++n n n n n n n n n n n n n n n n a a n n n n n n 故近似取其误差为1+≈k a ε，并且有m -1m -111021 21 ?=?=≈+βεk a ， （2）算法依据 使用matlab 编程时用digits 函数和vpa 函数来控制位数。 （3）Matlab 运行程序 %%保留11位有效数字 k1=11; s1=0;%用于存储这一步计算值 for n=0:50 a=(1/16^n)*(4/(8*n+1)-2/(8*n+4)-1/(8*n+5)-1/(8*n+6)); n1=n-1; if a<=0.5*10^(1-k1) break end end; for i=0:1:n1 t=(1/16^i)*(4/(8*i+1)-2/(8*i+4)-1/(8*i+5)-1/(8*i+6)); s1=s1+t; end s11=vpa(s1,k1); disp('保留11位有效数字的结果为：');disp(s11); disp('此时n 值为：');disp(n1); %%保留30位有效数字 clear all; k2=30;
计算方法上机实验报告
. / 《计算方法》上机实验报告 班级：XXXXXX 小组成员：XXXXXXX XXXXXXX XXXXXXX XXXXXXX 任课教师：XXX 二〇一八年五月二十五日
前言 通过进行多次的上机实验，我们结合课本上的内容以及老师对我们的指导，能够较为熟练地掌握Newton 迭代法、Jacobi 迭代法、Gauss-Seidel 迭代法、Newton 插值法、Lagrange 插值法和Gauss 求积公式等六种算法的原理和使用方法，并参考课本例题进行了MATLAB 程序的编写。 以下为本次上机实验报告，按照实验内容共分为六部分。 实验一： 一、实验名称及题目： Newton 迭代法 例2.7(P38):应用Newton 迭代法求在附近的数 值解，并使其满足. 二、解题思路： 设'x 是0)(=x f 的根,选取0x 作为'x 初始近似值,过点())(,00x f x 做曲线)(x f y =的切线L ,L 的方程为))((')(000x x x f x f y -+=,求出L 与x 轴交
点的横坐标) (') (0001x f x f x x - =,称1x 为'x 的一次近似值,过点))(,(11x f x 做曲线)(x f y =的切线,求该切线与x 轴的横坐标) (') (1112x f x f x x - =称2x 为'x 的二次近似值,重复以上过程,得'x 的近似值序列{}n x ,把) (') (1n n n n x f x f x x - =+称为'x 的1+n 次近似值，这种求解方法就是牛顿迭代法。 三、Matlab 程序代码： function newton_iteration(x0,tol) syms z %定义自变量 format long %定义精度 f=z*z*z-z-1; f1=diff(f);%求导 y=subs(f,z,x0); y1=subs(f1,z,x0);%向函数中代值 x1=x0-y/y1; k=1; while abs(x1-x0)>=tol x0=x1; y=subs(f,z,x0); y1=subs(f1,z,x0); x1=x0-y/y1;k=k+1; end x=double(x1) K 四、运行结果：
太原理工大学数值计算方法实验报告
本科实验报告 课程名称：计算机数值方法 实验项目：方程求根、线性方程组的直接解 法、线性方程组的迭代解法、代数插值和最 小二乘拟合多项式 实验地点：行勉楼 专业班级： ******** 学号： ********* 学生姓名： ******** 指导教师：李誌，崔冬华 2016年 4 月 8 日
y = x*x*x + 4 * x*x - 10; return y; } float Calculate(float a,float b) { c = (a + b) / 2; n++; if (GetY(c) == 0 || ((b - a) / 2) < 0.000005) { cout << c <<"为方程的解"<< endl; return 0; } if (GetY(a)*GetY(c) < 0) { return Calculate(a,c); } if (GetY(c)*GetY(b)< 0) { return Calculate(c,b); } } }; int main() { cout << "方程组为：f(x)=x^3+4x^2-10=0" << endl; float a, b; Text text; text.Getab(); a = text.a; b = text.b; text.Calculate(a, b); return 0; } 2.割线法： // 方程求根（割线法）.cpp : 定义控制台应用程序的入口点。// #include "stdafx.h" #include"iostream"
心得体会 使用不同的方法，可以不同程度的求得方程的解，通过二分法计算的程序实现更加了解二分法的特点，二分法过程简单，程序容易实现，但该方法收敛比较慢一般用于求根的初始近似值，不同的方法速度不同。面对一个复杂的问题，要学会简化处理步骤，分步骤一点一点的循序处理，只有这样，才能高效的解决一个复杂问题。
数值计算实验报告
(此文档为word格式，下载后您可任意编辑修改！) 2012级6班###(学号)计算机数值方法 实验报告成绩册 姓名:宋元台 学号： 成绩：
数值计算方法与算法实验报告 学期： 2014 至 2015 第 1 学期 2014年 12月1日课程名称: 数值计算方法与算法专业:信息与计算科学班级 12级5班 实验编号： 1实验项目Neton插值多项式指导教师：孙峪怀 姓名：宋元台学号：实验成绩： 一、实验目的及要求 实验目的： 掌握Newton插值多项式的算法，理解Newton插值多项式构造过程中基函数的继承特点，掌握差商表的计算特点。 实验要求： 1. 给出Newton插值算法 2. 用C语言实现算法 二、实验内容 三、实验步骤(该部分不够填写.请填写附页)
1.算法分析： 下面用伪码描述Newton插值多项式的算法： Step1 输入插值节点数n，插值点序列{x(i),f(i)},i=1,2,……,n,要计算的插值点x. Step2 形成差商表 for i=0 to n for j=n to i f(j)=((f(j)-f(j-1)(x(j)-x(j-1-i)); Step3 置初始值temp=1,newton=f(0) Step4 for i=1 to n temp=(x-x(i-1))*temp*由temp(k)=(x-x(k-1))*temp(k-1)形成 (x-x(0).....(x-x(i-1)* Newton=newton+temp*f(i); Step5 输出f（x）的近似数值newton(x)=newton. 2.用C语言实现算法的程序代码 #include
c 计算器实验报告
简单计算器 姓名: 周吉祥 实验目的：模仿日常生活中所用的计算器，自行设计一个简单的计算器程序，实现简单的计算功能。 实验内容： (1)体系设计： 程序是一个简单的计算器，能正确输入数据，能实现加、减、乘、除等算术运算，运算结果能正确显示，可以清楚数据等。 (2)设计思路： 1)先在Visual C++ 6.0中建立一个MFC工程文件，名为 calculator. 2)在对话框中添加适当的编辑框、按钮、静态文件、复选框和单 选框 3)设计按钮，并修改其相应的ID与Caption. 4)选择和设置各控件的单击鼠标事件。 5)为编辑框添加double类型的关联变量m_edit1. 6)在calculatorDlg.h中添加math.h头文件，然后添加public成 员。 7)打开calculatorDlg.cpp文件，在构造函数中，进行成员初始 化和完善各控件的响应函数代码。 (3)程序清单：
●添加的public成员： double tempvalue; //存储中间变量 double result; //存储显示结果的值 int sort; //判断后面是何种运算：1.加法2.减法3. 乘法 4.除法 int append; //判断后面是否添加数字 ●成员初始化： CCalculatorDlg::CCalculatorDlg(CWnd* pParent /*=NULL*/) : CDialog(CCalculatorDlg::IDD, pParent) { //{{AFX_DATA_INIT(CCalculatorDlg) m_edit1 = 0.0; //}}AFX_DATA_INIT // Note that LoadIcon does not require a subsequent DestroyIcon in Win32 m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME); tempvalue=0; result=0; sort=0; append=0; }
成本会计-学习指南 西交大考试题库及答案
成本会计-学习指南 一、单项选择题 1．产品成本计算分类法的成本计算对象是（A ） A．产品类别B．产品品种 C．产品规格D．产品加工步骤 2．生产经营费用按费用的（B ）分类形成要素费用。 A．经济内容B．经济性质 C．经济用途D．经济作用 3．对大量大批生产的产品，应当以（A ）作为产品成本计算对象。 A．产品的品种B．产品的批次 C．产品的生产步骤D．产品的类别 4．最基本的产品成本计算方法是（C ） A．分批法B．分步法 C．品种法D．分类法 5．李某本月生产甲零件2000只，其中合格品1950只，工废品30只，料废品20只。本月李某计算计件工资的甲零件数量是（ C ） A．2000 B．1980 C．1970 D．1950 6．成本会计的对象是（D ） A．产品生产成本的形成 B．各项期间费用的支出和归集 C. 生产费用和期间费用 D．各行业企业生产经营业务的成本和有关的期间费用 7．下列制造费用分配方法中，使制造费用账户可能出现余额的是（D ）A．工时比例法B．工资比例法 C．机时比例法D．年度计划分配率法 8．成本会计的最基本职能是（C ） A．成本预测B．成本决策 C．成本核算D．成本分析 9．下列企业中，适合运用品种法计算产品成本的是（A ） A．发电厂B．纺织厂 C．拖拉机厂D．造船厂 10．王某去年8月参加工作（病假扣发比例为40%），月标准工资418元，本月日历天数为31天，出勤19天，双休日8天，病假4天（合双休日1天）。若按月薪制计算，月工作天数为20.9天，则本月应付王某的计时工资是（ B ）A．386元B．394元C．396元D．418元 11．下列报表中不属于产品成本报表的是（D ） A．主要产品单位成本表B．制造费用明细表 C．营业费用明细表D．主营业务收支明细表 12．甲、乙两种产品的重量不同、材料单位消耗量基本相同、企业没有制定材料单位消耗定额、材料领用时未能区分每种材料的消耗量，则对甲、乙产品共同消耗的材料费用，可以用作为分配标准的是（ B ）
计算方法实验报告
计算方法实验报告（四） 方程和方程组的迭代解法 一、实验问题 利用简单迭代法，两种加速技术，牛顿法，改进牛顿法，弦割法求解习题5-1，5-2，5-3中的一题，并尽可能准确。 选取5-3：求在x=1.5附近的根。 二、问题的分析（描述算法的步骤等） （1）简单迭代法算法： 给定初始近似值,求的解。 Step 1 令i=0； Step 2 令(计算)； Step 3 如果,则迭代终止，否则重复Step 2。 （2）Aitken加速法算法 Step 1 令k=0，利用简单迭代算法得到迭代序列； Step 2 令-(计算得到一个新的序列，其中k=0,1,2…)；Step 3 如果,则迭代终止，否则重复Step 2。 (3)插值加速法算法 Step 1 令k=0，利用简单迭代算法得到迭代序列； Step 2 令+(计算得到一个新的序列，其中k=1,2,3…)； Step 3 如果,则迭代终止，否则重复Step 2。 (4)牛顿法算法
Step 1给定初始近似值； Step 2令,其中k计算得到的序列； Step 3如果,则迭代终止，否则重复Step 2。 （5）改进牛顿法的算法 Step 1给定初始近似值； Step 2令,其中k迭代计算得到的序列； Step 3如果,则迭代终止，否则重复Step 2。 （6）弦割法算法（双点弦割法） Step 1给定初始近似值，； Step 2令其中k计算得到的序列； Step 3如果,则迭代终止，否则重复Step 2。 三、程序设计 (1)简单迭代法 利用迭代公式进行迭代运算。 #include
计算方法实验报告格式
计算方法实验报告格式 小组名称: 组长姓名(班号)： 小组成员姓名(班号): 按贡献排序情况: 指导教师评语: 小组所得分数: 一个完整的实验，应包括数据准备、理论基础、实验内容及方法，最终对实验结果进行分析，以达到对理论知识的感性认识，进一步加深对相关算法的理解，数值实验以实验报告形式完成，实验报告格式如下： 一、实验名称 实验者可根据报告形式需要适当写出． 二、实验目的及要求 首先要求做实验者明确，为什么要做某个实验，实验目的是什么，做完该实验应达到什么结果，在实验过程中的注意事项，实验方法对结果的影响也可以以实验目的的形式列出． 三、算法描述（实验原理与基础理论） 数值实验本身就是为了加深对基础理论及方法的理解而设置的，所以要求将实验涉及到的理论基础，算法原理详尽列出． 四、实验内容 实验内容主要包括实验的实施方案、步骤、实验数据准备、实验的算法以及可能用到的仪器设备． 五、程序流程图 画出程序实现过程的流程图，以便更好的对程序执行的过程有清楚的认识，在程序调试过程中更容易发现问题． 六、实验结果 实验结果应包括实验的原始数据、中间结果及实验的最终结果，复杂的结果可以用表格
形式列出，较为简单的结果可以与实验结果分析合并出现． 七、实验结果分析 实验结果分析包括对对算法的理解与分析、改进与建议． 数值实验报告范例 为了更好地做好数值实验并写出规范的数值实验报告，下面给出一简单范例供读者参考． 数值实验报告 小组名称: 小组成员(班号): 按贡献排序情况: 指导教师评语: 小组所得分数: 一、实验名称 误差传播与算法稳定性． 二、实验目的 1．理解数值计算稳定性的概念． 2．了解数值计算方法的必要性． 3．体会数值计算的收敛性与收敛速度． 三、实验内容 计算dx x x I n n ? += 1 10 ，1,2,,10n = ． 四、算法描述 由 dx x x I n n ? += 1 10 ，知 dx x x I n n ?+=--101110，则
西安交通大学计算方法B大作业
计算方法上机报告 姓名： 学号： 班级：
目录 题目一------------------------------------------------------------------------------------------ - 4 - 1.1题目内容 ---------------------------------------------------------------------------- - 4 - 1.2算法思想 ---------------------------------------------------------------------------- - 4 - 1.3Matlab源程序----------------------------------------------------------------------- - 5 - 1.4计算结果及总结 ------------------------------------------------------------------- - 5 - 题目二------------------------------------------------------------------------------------------ - 7 - 2.1题目内容 ---------------------------------------------------------------------------- - 7 - 2.2算法思想 ---------------------------------------------------------------------------- - 7 - 2.3 Matlab源程序---------------------------------------------------------------------- - 8 - 2.4计算结果及总结 ------------------------------------------------------------------- - 9 - 题目三----------------------------------------------------------------------------------------- - 11 - 3.1题目内容 --------------------------------------------------------------------------- - 11 - 3.2算法思想 --------------------------------------------------------------------------- - 11 - 3.3Matlab源程序---------------------------------------------------------------------- - 13 - 3.4计算结果及总结 ------------------------------------------------------------------ - 14 - 题目四----------------------------------------------------------------------------------------- - 15 - 4.1题目内容 --------------------------------------------------------------------------- - 15 - 4.2算法思想 --------------------------------------------------------------------------- - 15 - 4.3Matlab源程序---------------------------------------------------------------------- - 15 - 4.4计算结果及总结 ------------------------------------------------------------------ - 16 - 题目五----------------------------------------------------------------------------------------- - 18 -
数学计算方法实验报告
数学计算方法实验报告 习题二 2.估计用二分法求方程f(x)=x3+4x2-10=0在区间[1,2]内根的近似值，为使方程不超过10时所需的二分次数。f(x k) 程序过程： function two (tolerance) a=1;b=2;counter=0; while (abs(b-a)>tolerance) c=(a+b)/2; fa=a^3+4*a^2-10;
fb=b^3+4*b^2-10; fc=c^3+4*c^2-10; if ((fa==0|fb==0)) disp(counter); elseif (fa*fc<0) b=c;counter=counter+1; elseif (fb*fc<0) a=c;counter=counter+1; elseif (fb==0) disp(counter); end end solution=(a+b)/2; disp(solution); disp(counter); 实验结果： 6.取x0=1.5,用牛顿迭代法求第三中的方程根.f(x)=x3+4x2-10=0的近似值（精确到||x k+1-x k|≦10-5,并将迭代次数与3题比较。 程序过程： function six (g) a=1.5; fa=a^3+4*a^2-10;
ga=3*a^2+8*a; b=a-fa/ga; k=1; while(abs(b-a)>g) a=b; fa=a^3+4*a^2-10; ga=3*a^2+8*a; b=a-fa/ga; k=k+1; end format long; disp(a); disp(k); 实验结果：程序结果计算结果 8.用弦割法求方程f(x)=x3-3x2-x+9=0在区间[-2,-1]内的一个实根近似值x k,|f(x k)|≦10-5. 程序过程： function eight (t) a=-2; b=-1; fa=a^3-3*a^2-a+9; fb=b^3-3*b^2-b+9; c=b-fb*(b-a)/(fb-fa); k=1; while(abs(c-b)>t) a=b; b=c; fa=a^3-3*a^2-a+9; fb=b^3-3*b^2-b+9; c=b-fb*(b-a)/(fb-fa); k=k+1; end
计算方法实验报告 拟合
南京信息工程大学实验（实习）报告 一、实验目的: 用最小二乘法将给定的十个点拟合成三次多项式。 二、实验步骤： 用matlab编制以函数为基的多项式最小二乘拟合程序，并用于对下列数据作三次多项式最小二乘拟合(取权函数wi=1) x -2.5 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 y -2.30 -1 -0.14 -0.25 0.61 1.03 1.75 2.75 4.42 6.94 给定直线方程为：y=1/4*x3+1/2*x2+x+1 三、实验结论： 最小二乘法：通过最小化误差的平方和寻找数据的最佳函数匹配。利用最小二乘法可以简便地求得未知的数据，并使得这些求得的数据与实际数据之间误差的平方和为最小。最小二乘法还可用于曲线拟合。 一般地。当测量数据的散布图无明显的规律时，习惯上取n次代数多项式。 程序运行结果为： a = 0.9731 1.1023 0.4862 0.2238 即拟合的三次方程为：y=0.9731+1.1023x+0.4862*x2+0.2238*x3
-2.5 -2-1.5-1-0.5 00.51 1.52 2.5 -4-20246 81012 x 轴 y 轴 拟合图 离散点 y=a(1)+a(2)*x+a(3)*x.2+a(4)*x.3 结论： 一般情况下，拟合函数使得所有的残差为零是不可能的。由图形可以看出最小二乘解决了残差的正负相互抵消的问题，使得拟合函数更加密合实验数据。 优点：曲线拟合是使拟合函数和一系列的离散点与观测值的偏差平方和达到最小。 缺点：由于计算方法简单，若要保证数据的精确度，需要大量的数据代入计算。
西交大计算方法上机报告
计算方法（B）实验报告 姓名： 学号： 学院： 专业：
实验一 三对角方程组Tx f =的求解 一、 实验目的 掌握三对角方程组Tx f =求解的方法。 二、 实验内容 求三对角方程组Tx f =的解，其中： 4 -1 -1 4 -1 -1 4 1 -1 4T ????????=?? ?? ???? ， 3223f ?? ? ? ?= ? ? ??? 三、 算法组织 设系数矩阵为三对角矩阵 11222333111 b c a b c a b c a b c b n n n n T ---???????? =?????? ?????? 则方程组Tx f =称为三对角方程组。 设矩阵T 非奇异，T 可分解为T=LU ，其中L 为下三角矩阵，U 为单位上三角矩阵，记 1 1 212 313 1 1 1111 ,11n n n n n r l r l r L U l r l μμμμμ---???? ? ? ? ? ? ?== ? ? ? ? ? ? ? ? ? ?? ? ? ? 可先依次求出,L U 中的元素后，令Ux y =，先求解下三角方程组Ly f =得出 y ，再求解上三角方程组Ux y =。 追赶法的算法组织如下： 1.输入三对角矩阵T 和右端向量f ；
2.将Tx f =压缩为四个一维数组{}{}{}{}i i i i a b c d 、、、，{}{}{}i i i a b c 、、是T 的三对角线性方程组的三个对角，{}i d 是右端向量。将分解矩阵压缩为三个一维数组 {}{}{}i i i l r μ、、。 3.对T 做Crout 分解（也可以用Doolittle 分解）导出追赶法的计算步骤如下： 1111,b r c μ== for 2i n = 111, , ,i i i i i i i i i i i i i l a b a r r c y d l y μμ---==-==- end 4.回代求解x /n n n x y μ= for 11i n =- 1()/i i i i i x y c x μ+=- end 5. 停止，输出结果。 四、 MATLAB 程序 MATLAB 程序见附件1. 五、 结果及分析 实验结果为： (1.0000 1.0000 1.0000 1.0000)T x =
计算方法实验报告
实验报告 一、求方程f(x)=x^3-sinx-12x+1的全部根, ε=1e -6 1、 用一般迭代法; 2、 用牛顿迭代法; 并比较两种迭代的收敛速度。 一、首先，由题可求得：12cos 3)(2 ' --=x x x f . 其次，分析得到其根所在的区间。 ① 令()0=x f ，可得到x x x sin 1123 =+-. ② 用一阶导数分析得到1123 +-x x 和x sin 两个函数的增减区间；再用二阶导数分析得到 两个函数的拐点以及凹凸区间. ③ 在直角坐标轴上描摹出01123 =+-x x 和0sin =x 的图，在图上可以看到他们的交点，然后估计交点所在的区间，即是所要求的根的区间。经过估计，得到根所在的区间为 []3,4--，[]1,0和[]4,3. 1、 一般迭代法 （1）算法步骤： 设ε为给定的允许精度，迭代法的计算步骤为： ① 选定初值0x .由()0=x f 确定函数()x g ，得等价形式()x g x =. ② 计算()0x g .由迭代公式得()01x g x =. ③ 如果ε≤-01x x ，则迭代结束，取1x 为解的近似值；否则，用1x 代替0x ，重复步骤②和步骤③. （2）程序代码： ① 在区间[]3,4--内， 代码： clc
x0=-3.5; %初值0x iter_max=100; %迭代的最大次数 ep=1e-6; %允许精度 ε k=0; while k<=iter_max %k 从0开始到iter_max 循环 x1=(sin(x0)+12*x0-1).^(1/3); %代入0x ，算出1x 的值 if abs(x1-x0) 西安电子科技大学出版社《计算方法》任传祥等编著第九章计算方法上机参考答案 实验一，算法一 #include #include 计算方法A 上机大作业 1. 共轭梯度法求解线性方程组 算法原理：由定理3.4.1可知系数矩阵A 是对称正定矩阵的线性方程组Ax=b 的解与求解二次函数1()2 T T f x x Ax b x =-极小点具有等价性，所以可以利用共轭梯度法求解1()2 T T f x x Ax b x = -的极小点来达到求解Ax=b 的目的。 共轭梯度法在形式上具有迭代法的特征，在给定初始值情况下，根据迭代公式： (1)()()k k k k x x d α+=+ 产生的迭代序列(1)(2)(3)x x x ，，，... 在无舍入误差假定下，最多经过n 次迭代，就可求得()f x 的最小值，也就是方程Ax=b 的解。 首先导出最佳步长k α的计算式。 假设迭代点()k x 和搜索方向()k d 已经给定，便可以通过()()()() k k f x d φαα=+的极小化 ()()min ()()k k f x d φαα=+ 来求得，根据多元复合函数的求导法则得: ()()()'()()k k T k f x d d φαα=?+ 令'()0φα=，得到: ()() ()()k T k k k T k r d d Ad α=，其中()()k k r b Ax =- 然后确定搜索方向()k d 。给定初始向量(0)x 后，由于负梯度方向是函数下降最快的方向，故第一次迭代取搜索方向(0) (0)(0)(0)()d r f x b Ax ==-?=-。令 (1)(0)00x x d α=+ 其中(0)(0)0(0)(0) T T r d d Ad α=。第二次迭代时，从(1) x 出发的搜索方向不再取(1)r ，而是选取(1) (1)(0)0d r d β=+，使得(1)d 与(0)d 是关于矩阵A 的共轭向量，由此可 求得参数0β: 《数据分析方法与技术》上机实验——实验1描述性统计方法 学号： 姓名： 日期： 实验项目（一）：描述性统计方法 一、实验内容 1．实验目的 掌握常用的描述性图表展示方法的原理及操作，包括：频数分布表、分组频数表、列联表、茎叶图、箱线图、误差图、散点图等； 掌握常用的描述性统计方法的原理及操作，包括：算术平均值、中位数、众数、四分位数、极差、平均差、方差、标准差、标准分数、离散系数等。 2. 实验内容和要求 实验内容：基于标准数据集，属性描述性图表展示方法（数分布表、分组频数表、列联表、茎叶图、箱线图、误差图、散点图等），对统计指标（算术平均值、中位数、众数、极差、平均差、方差、标准差、标准分数、离散系数、偏态峰态）进行计算。 实验要求：掌握各种描述性统计指标的计算思路及其在SPSS或EXCEL环境下的操作方法，掌握输出结果的解释。 二、实验过程 1、数据集介绍 1.数据库标题:鲍鱼数据 2.该数据库共计4177行数据 3.该数据有八个属性（包含性别共有九项） 4.以下是关于属性的描述，包括属性的名称，数据类型，测量单元和一个简短的描述： Name Data TypeMeas.Description ---- --------- ----- ----------- Sex nominal M, F, and I (infant)鲍鱼宝宝 Length continuousmm Longest shell measurement最长壳 Diameter continuousmm perpendicular to length垂直长度 Height continuousmm with meat in shell有肉的壳高度 Whole weightcontinuousgramswhole abalone整个鲍鱼 Shucked weightcontinuousgramsweight of meat肉的重量 Viscera weightcontinuousgramsgut weight (after bleeding)放血后内脏重 Shell weightcontinuousgramsafter being dried弄干后重量 Rings integer +1.5 gives the age in years +1.5=年龄 5.数据的值域 第一次作业 三、主观题(共14道小题) 21. 管理心理学的研究重点是组织管理中具体的社会、心理现象，以及（）、群体、（）、组织中的具体心理活动的规律性。 答：个体，领导 22. 现场实验又称为（） 答：自然实验法 23. 霍桑实验发现并证实了（）的存在 答：“非正式组织” 24. 超Y理论是由（）和（）提出来的。 答：莫尔斯,洛希 25. 投射是一种通过（）的方法而达到（）的目的。 答案：以己度人,心理防御 26. 挫折是人们在有目的的活动中遇到了无法克服或自以为是无法克服的障碍和干扰，其（）和（）不能满足时所产生的消极的情绪反应。 参考答案：需要,动机 27. 名词解释：观察法--- 答案： 观察法，是在未受控制的日常生活中，了解和分析人的言行、表情等，借此来判断被观察者心理活动的一种研究方法。 28. ：复杂人假设--- 答案：复杂人假设是指人是很复杂的，人们的需要与潜在的欲望多种多样，而且这些需要的模式也是随着年龄与发展阶段的变迁，随着所扮演的角色的变化，随着所处境遇及人际关系的演变而不断变化的。 29. 名词解释：知觉防御--- 答案：知觉防御是指人们对不利于自己的信息会视而不见或加以歪曲，以达到防御的目的。 30. 名词解释：角色知觉--- 答案：角色知觉是指人对于自己所处的特定的社会与组织中的地位的知觉。 31. 名词解释：心理疏导--- 答案：心理疏导是指运用一定的心理诱导的策略和方法使受挫者在别人引导下发挥内在潜力，达到消除心理障碍、明确前进方向、排除不良情绪和行为的目的。 32. 麦格雷戈关于人性假定的论述是什么？ 答案：（1）管理的理论与管理者的观念是第一位的，而管理的政策与具体措施是第二位的，不能本末倒置，也不能简单混同、不加区分。 （2）强调在管理中要着重开发人力资源，发觉人的“潜在力量”。 （3）管理人员采取哪种理论假定要看具体情况，但是所持理论的观点要旗帜鲜明。 33. 一个完整的角色知觉过程应该包括哪些成分？ 答案：一个完整的角色知觉过程应该包括以下四个成分：角色认知、角色行为、角色期望、角色评价。 角色认知是指一个人对自己应该在社会与组织中所处地位的认识。 角色行为是指一个人按照特定的社会与组织所赋予角色的特定的行为模式而进行的行为。角色期望是指他人对一个人所应承担角色的希望与寄托。 角色评价是指他人对一个人的角色扮演的评论与估价。 其中，角色认知与角色行为是角色扮演者主观方面的因素；而角色期望与角色评价是指他人对角色扮演者的反馈信息，属于客观方面的因素。角色知觉作为复杂的社会认知与社会知觉中的一个方面，只有在主客观因素相互作用的条件下，才能最后完整、正确地形成。 34. 生活压力源具体包括哪些方面？ 答案：生活压力源指应激起源于与员工个人生活有关的因素，具体包括四个方面： （1）重要人员的影响。包括员工家庭成员、师长、邻里或亲朋好友的期望与态度。 （2）个人生活事件的影响。包括结婚、离婚，家庭成员的生产、死亡等个人生活经历中的突发事件、重大变化，这些事件足以扰乱人们的生理与心理稳定。 （3）生活方式的变化。主要体现为现代生活的节奏加快，使人们产生不适感，以及消费导向的迷惘感的压力、对生活质量的高期望值与实际生活之间的差异造成的失望感和压力等。（4）经济收入压力。一方面，收入低会产生生活中入不敷出的压力；另一方面，收入高的人则可能有请客、救助，甚至道德等方面的压力。 第二次作业 三、主观题(共14道小题) 21. 人们对不利于自己的信息会视而不见或加以歪曲，以达到防御目的是指（）. 答案：知觉防御 22. 自我认识的内容包括以下三个方面：物质自我、社会自我和（）。 答案：精神自我 23. 形成个性的原因基本上可以归结为两个方面：（）和（）。 答：遗传因素，环境因素 24. 马斯洛的需要层次理论将人的需要分为了五个层次：（）、安全需要、爱的需要、（）和（）。 答案：生理需要，尊重需要，自我实现需要西安电子科技大学出版社计算方法上机答案
西交计算方法A上机大作业
数据分析方法与技术- 实验报告模板
西安交大《管理心理学B》作业高分答案