西安交通大学计算方法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;
digits(k2+2);
s2=vpa(0);%用于存储这一步计算值
for n=0:50
a=vpa((1/16^n)*(4/(8*n+1)-2/(8*n+4)-1/(8*n+5)-1/(8*n+6)));
n2=n-1;
if a<=0.5*10^(1-k2)
break
end
end;
for i=0:1:n2
t=vpa((1/16^i)*(4/(8*i+1)-2/(8*i+4)-1/(8*i+5)-1/(8*i+6)));
s2=vpa(s2+t);
end
s30=vpa(s2,k2);
disp('保留30位有效数字的结果为：');disp(s30);
disp('此时n值为：');disp(n2);
2.某通信公司在一次施工中，需要在水面宽度为20米的河沟底部沿直线走向铺设一条沟底光缆。在铺设光缆之前需要对沟底的地形进行初步探测，从而估计所需光缆的长度，为工程预算提供依据。已探测到一组等分点位置的深度数据(单位:米)如下表所示：
分点0 1 2 3 4 5 6
深度9.01 8.96 7.96 7.97 8.02 9.05 10.13
分点7 8 9 10 11 12 13
深度11.18 12.26 13.28 13.32 12.61 11.29 10.22
分点14 15 16 17 18 19 20
深度9.15 7.90 7.95 8.86 9.81 10.80 10.93
(1)请用合适的曲线拟合所测数据点；
(2)估算所需光缆长度的近似值，并作出铺设河底光缆的曲线图；（2）算法选择
由于多项式插值只适用于数据点较少的插值案例，此时可以避免产生较大的误差；当数据点过多时，多项式插值中的龙格现象会使误差变大。因此当所用数据点过多时，我们采用分段插值方式，即用分段多项式代替单个多项式作插值。分段多项式是由一些在相互连接的区间上的不同多项式连接而成的一条连续曲线，其中三次样条插值方法是一种具有较好“光滑性”的分段插值方法。
（3）matlab算法
clear all;
x=0:1:20; %产生水面宽度取值点的数组
y=[9.01,8.96,7.96,7.97,8.02,9.05,10.13,11.18,12.26,13.28,13.32,12.61, 11.29,10.22,9.15,7.90,7.95,8.86,9.81,10.80,10.93];%对应水面宽度处的水深
n=length(x);%取样点数目
M=y;
for k=2:3;%计算二阶差商
for i=n:-1:k;
M(i)=(M(i)-M(i-1))/(x(i)-x(i-k+1));
end
end
h(1)=x(2)-x(1);
for i=2:n-1;
h(i)=x(i+1)-x(i);
c(i)=h(i)/(h(i)+h(i-1));
a(i)=1-c(i);
b(i)=2;
d(i)=6*M(i+1);
end
M(1)=0; %选择自然边界条件
M(n)=0;
b(1)=2;
b(n)=2;
c(1)=0;
a(n)=0;
d(1)=0;
d(n)=0;
u(1)=b(1); %对三对角矩阵进行LU分解
y1(1)=d(1);
for k=2:n;
l(k)=a(k)/u(k-1);
u(k)=b(k)-l(k)*c(k-1);
y1(k)=d(k)-l(k)*y1(k-1);
end
M(n)=y1(n)/u(n); %追赶法求解样条参数M(i)
for k=n-1:-1:1;
M(k)=(y1(k)-c(k)*M(k+1))/u(k);
end
x1=1:0.1:20; %确定取值点
n1=length(x1); %确定取值点个数
s=zeros(1,n1);
for m=1:n1; k=1;
for i=2:n-1
if x1(m)<=x(i); k=i-1;break; else k=i; end end
h1=x(k+1)-x(k);%在各区间用三次样条插值函数计算X 点处的值 x11=x(k+1)-x1(m);x12=x1(m)-x(k);
sum(m)=(M(k)*(x11^3)/6+M(k+1)*(x12^3)/6+(y(k)-(M(k)*(h1^2)/6))*x11+(y (k+1)-(M(k+1)*(h1^2)/6))*x12)/h1; end
l=0; %计算所需光缆长度 for i=2:n1
l=l+sqrt((x1(i)-x1(i-1))^2+(sum(i)-sum(i-1))^2); end
disp('所需光缆长度为'); disp(l); figure(1)
plot(x,y,'*',x1,sum,'-') xlabel('位置'); ylabel('深度/m');
title('铺设河底光缆的曲线图'); grid on;
（4）结果分析
铺设海底光缆的曲线图如下图：
2468
101214161820
78
9
10
11
12
13
14
位置
深度/m
铺设河底光缆的曲线图
使用matlab 仿真计算出所需光缆的长度为 25.4525m 。 3.假定某天的气温变化记录如下表所示，试用数据拟合的方法找出这一天的气温变化的规律;试计算这一天的平均气温,并试估计误差。
时刻
1
2
3
4
5
6
7
8
9 10 11 12
平均气温 15 14 14 14 14 15 16 18 20 20 23 25 28 时刻
13 14 15 16 17 18 19 20 21 22 23 24
平均气温 31 34 31 29 27 25 24 22 20 18 17 16
（1） 选择算法
在本题中，采用“最小二乘法”来近似这一天的气温变化规律，分别使用二次函数、三次函数、五次函数以及指数型函数来对数据点进行拟合，并计算相应的误差。 （2）算法框架
a ）使用多项式进行最小二乘近似时，首先生成矩阵G ，然后利用算法LSS 求解多项式的系数和误差22E 。
b ）使用指数函数进行最小二乘近似时，设其指数形式为
2
321d x d x d e
y +=，再对函数形式进行变换
332211a ,a ,ln a ,ln d d d y z ====，即
2321a a a x x z ++=。然后生成矩阵G ，再利用算法LSS 求解多
项式的系数和误差22E 。 （3）matlab 程序
%%多项式的最小二乘近似 clear all; x=0:24;
y=[15,14,14,14,14,15,16,18,20,20,23,25,28,31,34,31,29,27,25,24,22,20,18,17,16];
m=length(x); %将数据点存入m
s=0; %求一天的平均气温
for i=1:1:24
s=s+y(i);
end
T=s/m;
disp('一天的平均气温为');
disp(T);
n=h+1; %h次函数的最小二乘近似，h分别取2、3、5
for j=1:n; %生成矩阵G
for i=1:25
G(i,j)=x(i).^(j-1);
end
end
for i=1:m; %将y作为G的最后一列存放
G(i,j+1)=y(i);
end
for k=1:n
if G(k,k)>0 %形成矩阵Q（k)
sgn=1;
elseif G(k,k)==0;
sgn=0;
else sgn=-1;
end
a1=-sgn*sqrt(sum(G(k:m,k).^2));
w=zeros(1,n);
w(k)=G(k,k)-a1;
for j=k+1:m
w(j)=G(j,k);
end
b1=a1*w(k);
G(k,k)=a1; %变换Gk-1到Gk
for j=k+1:n+1
t=sum(w(k:m)*G(k:m,j))/b1;
for i=k:m
G(i,j)=G(i,j)+t*w(i);
end
end
end
a(n)=G(n,n+1)/G(n,n); %求解方程
for i=n-1:-1:1
a(i)=(G(i,n+1)-sum(G(i,i+1:n).*a(i+1:n)))/G(i,i); end
E=sum(G(n+1:m,n+1).^2); %计算误差
t=0:0.01:24;
a=fliplr(a); %将系数数组左右翻转
y1=poly2sym(a); %将系数数组转化为多项式
subs(y1,'x',t);
y1=double(ans);
figure(1)
plot(x,y,'*',t,y1,'r-');
xlabel('时刻');
plot(x,y,'*',t,y1,'r-');
xlabel('时刻');
ylabel('温度');
title('5次函数的最小二乘曲线');
grid on;
disp('误差为');disp(E);
%%指数形式的最小二乘近似
clear all;
x=0:24;
y=[15,14,14,14,14,15,16,18,20,20,23,25,28,31,34,31,29,27,25,24,22,20, 18,17,16];
m=length(x);
n=3;
z=zeros(1,m);
for i=1:m
z(i)=log(y(i));
end
for j=1:n; %生成矩阵G
for i=1:m;
G(i,j)=x(i).^(j-1);
end
end
for i=1:m; %将y作为G的最后一列存放
G(i,j+1)=z(i);
end
for k=1:n
if G(k,k)>0; %形成矩阵Q(k)
sgn=1;
elseif G(k,k)==0;
sgn=0;
else sgn=-1;
end
a1=-sgn*sqrt(sum(G(k:m,k).^2));
w=zeros(1,n);
w(k)=G(k,k)-a1;
for j=k+1:m
w(j)=G(j,k);
end
b1=a1*w(k);
G(k,k)=a1; %变换Gk-1到Gk
for j=k+1:n+1
t=sum(w(k:m)*G(k:m,j))/b1;
for i=k:m;
G(i,j)=G(i,j)+t*w(i);
end
end
end
a(n)=G(n,n+1)/G(n,n); %解方程组
for i=n-1:-1:1
a(i)=(G(i,n+1)-sum(G(i,i+1:n).*a(i+1:n)))/G(i,i);
end
d3=a(3);
d2=a(2);
d1=exp(a(1));
E=0;
for i=1:25
y1(i)=d1*exp(d2*x(i)+d3*x(i)*x(i));
E=(y1(i)-y(i)).^2+E; %计算误差
end
t=0:0.01:24;
y1=d1.*exp(d2.*t+d3.*t.*t);
fprintf('指数函数拟合的系数是：d1=%f,d2=%f,d3=%f',d1,d2,d3); figure(1)
plot(x,y,'*',t,y1,'r-');
xlabel('时刻');
ylabel('温度');
title('指数函数的最小二乘曲线');
grid on;
fprintf('\n指数函数拟合的系数是：%f',E);
（4）结果分析
a）二次多项式的最小二乘近似
一天的平均气温为：20.5600
误差为：280.3395
拟合曲线如下图所示：
0510
152025
5
10
15
20
25
30
35
时刻
温度
2次函数的最小二乘曲线
b ） 三次多项式的最小二乘近似 误差为：131.0618 近似函数曲线如下图所示：
510
152025
10152025
30
35
时刻
温度
3次函数的最小二乘曲线
c ）五次多项式的最小二乘近似 误差为：33.1446
近似函数曲线如下图所示：
0510
152025
10
15
20
2530
35
时刻
温度
5次函数的最小二乘曲线
d ）指数形式的最小二乘近似
指数函数拟合的系数是：d1=10.842462,d2=0.125093,d3=-0.004442 指数函数拟合的误差E 2
是：215.700187 近似函数曲线如下图所示：
510
152025
10152025
30
35
时刻
温度
指数函数的最小二乘曲线
从上述几种拟合可以得出，
多项式的次数越高，拟合的效果越好，误差越小；指数函数的二次多项式拟合相对于同阶的多项式最小二乘近似其误差相对较小，但不如高阶多项式最小二乘近似精确。
4.设计算法，求出非线性方程52645200x x -+=的所有实根，并使误差不超过410-。 （1） 问题分析
首先对函数进行分析，并确定函数值域的范围。令
()()x x x f x x x f 9030 ,20456425-='+-=则。由此可以确定函数()x f 出现
极值的x 的范围在()21-，
，故先选取()1010-，∈x 对函数作图如下：
-10
-8-6-4-2
0246810
-8-6-4-20
2
4
6
x 10
5
x
y
再进一步缩小取值区间()22-，∈x ，如下图：
-2
-1.5-1-0.5
00.51 1.52
-400
-350-300-250-200-150
-100-50050x
y
上图已经粗略的标识出了解的取值范围，我们可以参考上图确定合适的取值范围来对问题进行求解。 （2） 选择算法
这里采用二分法的基本原理进行求解，找到满足精度要求的解。二分法属于区间法的一种，区间法有一个重要的优点：这种方法总是收敛的。
二分法首先产生一串区间，并缩短旧区间的长度来得到新的区间，下一步区间的长度是上一步区间长度的一半，且有一个端点是上一步的端点。其取值方法是：
()
()()[
]122
1+++=k k k x x x
若()()()()02<+k k x f x f ，则取()()()[]11,++=k k k x x I ，否则取()()()[]121,+++=k k k x x I ，并且重复上述过程。显然，每次迭代使区间的长度减小一半，因此二分法总是收敛的。 （4） matlab 程序
clear all; x=-10:0.5:10;
y=6*(x.^5)-45*(x.^2)+20; g=[];
for i=-10:0.5:10 %确定解的区间 k=i+0.5;
if (y(x==i)*y(x==k) syms x; f=6*x^5-45*x^2+20; n=length(g); %确定根的个数 for j=1:n x0=g(j);%求根区间左端点 x1=g(j)+1;%求根区间右端点 while (x1-x0)>=10^(-4) if subs(f,x,x0)*subs(f,x,(x0+x1)/2)>eps x0=(x0+x1)/2; else x1=(x0+x1)/2; end end root=x0 %输出方程的根 end （5）结果分析 使用matlab运行程序得到该非线性方程组的三个实根，分别为-0.6546，0.6812，1.8708。 西安交大统计学考试试卷 一、单项选择题（每小题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、“性别”是品质标志。（对） 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; 数值分析上机题 第一章 17.（上机题）舍入误差与有效数 设∑=-= N j N j S 2 2 11 ，其精确值为）111-23（21+-N N 。 （1）编制按从大到小的顺序1 -1 ···1-311-21222N S N +++=，计算N S 的通用 程序； （2）编制按从小到大的顺序1 21 ···1)1(111 222-++--+ -=N N S N ，计算N S 的通用程序； （3）按两种顺序分别计算210S ,410S ,610S ,并指出有效位数（编制程序时用单精度）； （4）通过本上机题，你明白了什么？ 解： 程序： （1）从大到小的顺序计算1 -1 ···1-311-21222N S N +++= ： function sn1=fromlarge(n) %从大到小计算sn1 format long ; sn1=single(0); for m=2:1:n sn1=sn1+1/(m^2-1); end end （2）从小到大计算1 21 ···1)1(111 2 22 -++--+-= N N S N function sn2=fromsmall(n) %从小到大计算sn2 format long ; sn2=single(0); for m=n:-1:2 sn2=sn2+1/(m^2-1); end end （3） 总的编程程序为： function p203() clear all format long; n=input('please enter a number as the n:') sn=1/2*(3/2-1/n-1/(n+1));%精确值为sn fprintf('精确值为%f\n',sn); sn1=fromlarge(n); fprintf('从大到小计算的值为%f\n',sn1); sn2=fromsmall(n); fprintf('从小到大计算的值为%f\n',sn2); function sn1=fromlarge(n) %从大到小计算sn1 format long; sn1=single(0); for m=2:1:n sn1=sn1+1/(m^2-1); end end function sn2=fromsmall(n) %从小到大计算sn2 format long; sn2=single(0); for m=n:-1:2 sn2=sn2+1/(m^2-1); end end end 运行结果： 数值分析上机题目 1、 分别用不动点迭代与Newton 法求解方程250x x e -+=的正根与负根。 2、 Use each of the following methods to find a solution in [0.1,1] accurate to within 10^-4 for 4326005502002010x x x x -+--= a. Bisection method b. Newton’s method c. Secant method d. Method of False Position e. Muller’s method 3、 应用Newton 法求f （x ）的零点，e=10^-6,这里f （x ）=x-sin （x ）。 再用求重根的两种方法求f （x ）的零点。 4、 应用Newton 法求f （x ）的零点，e=10^-6,f(x)=x-sin(x) 再用Steffensen’s method 加速其收敛。 5、 用Neville’s 迭代差值算法，对于函数2 1 (),11125f x x x = -≤≤+进行lagrange 插值。取不同的等分数n=5,10，将区间[-1,1]n 等分，取等距节点。把f(x)和插值多项式的曲线画在同一张图上进行比较。 6、 画狗的轮廓图 7、 Use Romberg integration to compute the following approximations to ? a 、 Determine R1,1,R2,1,R3,1,R4,1and R5,1,and use these approximations to predict the value of the integral. b 、 Determine R2,2 ,R3,3 ,R4,4 ,and R5,5,and modify your prediction. c 、 Determine R6,1 ,R6,2 ,R6,3 ,R6,4 ,R6,5 and R6,6,and modify your prediction. 西安交大成本会计在线 作业答案 文件编码（008-TTIG-UTITD-GKBTT-PUUTI-WYTUI-8256） 《成本会计》3（2017） 试卷总分：100 测试时间：-- 一、单选题（共25道试题，共50分。） 1.如果同一时期内，在几张定单中规定有相同的产品，则计算成本时 可以(D )。 A. 按定单分批组织生产 B. 按品种分批组织生产 C. 按产品的组成部分分批组织生产 D. 将相同产品合为一批组织生产 满分：2分 2.不在“财务费用”账户核算的项目是(A )。 A. 业务招待费 B. 利息费用 C. 汇兑损失 D. 金融机构结算手续费 满分：2分 3.“基本生产成本”月末借方余额表示(B )。 A. 本期发生的生产费用 B. 完工产品成本 C. 月末在产品成本 D. 累计发生的生产费用 满分：2分 4.下列不属于成本计算基本方法的是(C )。 A. 品种法 B. 分批法 C. 分类法 满分：2分 5.成本还原的对象是（D ）。 A. 产成品成本 B. 各步骤半成品成本 C. 最后步骤产成品成本 D. 产成品成本中所耗上步骤半成品成本费用 满分：2分 6.采用计划成本分配法分配辅助生产费用，辅助生产的实际成本是 (B )。 A. 按计划成本分配前的实际费用 B. 按计划成本分配前的实际费用加上按计划成本分配转入的费用 C. 按计划成本分配前的实际费用减去按计划成本分配转出的费用 D. 按计划成本分配前实际费用加上按计划成本分配转入的费用， 减去按计划成本分配转出的费用 满分：2分 7.成本会计最基本的任务和中心环节是（ C）。 A. 进行成本预测，编制成本计划 B. 审核和控制各项费用的支出 C. 进行成本核算，提供实际成本的核算资料 D. 参与企业的生产经营决策 满分：2分 8.下列各项属于产品成本项目的有（C ）。 A. 财务费用 B. 管理费用 成本会计-学习指南 一、单项选择题 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 ） 第一章 一、题目 设∑ =-= N N j S 2 j 2 1 1，其精确值为)11 123(21+--N N 。 1) 编制按从大到小的顺序1 1 13112122 2-+??+-+-=N S N ，计算S N 的通用程序。 2) 编制按从小到大的顺序1 21 1)1(111222-+ ??+--+-= N N S N ，计算S N 的通用程序。 3) 按两种顺序分别计算64210,10,10S S S ,并指出有效位数。（编制程序时用单精度） 4) 通过本次上机题，你明白了什么？ 二、通用程序 N=input('Please Input an N (N>1):'); AccurateValue=single((0-1/(N+1)-1/N+3/2)/2); Sn1=single(0); for a=2:N; Sn1=Sn1+1/(a^2-1); end Sn2=single(0); for a=2:N; Sn2=Sn2+1/((N-a+2)^2-1); end fprintf('The value of Sn (N=%d)\n',N); fprintf('Accurate Calculation %f\n',AccurateValue); fprintf('Caculate from large to small %f\n',Sn1); fprintf('Caculate from small to large %f\n',Sn2); disp('____________________________________________________') 三、结果 从结果可以看出有效位数是6位。 感想：可以得出，算法对误差的传播有一定的影响，在计算时选一种好的算法可以使结果更为精确。从以上的结果可以看到从大到小的顺序导致大数吃小数的现象，容易产生较大的误差，求和运算从小数到大数所得到的结果才比较准确。 数值分析上机题1 设2 21 1N N j S j ==-∑ ，其精确值为1311221N N ??-- ?+?? 。 （1）编制按从大到小的顺序222 111 21311 N S N = +++---，计算N S 的通用程序。 （2）编制按从小到大的顺序22 21111(1)121 N S N N =+++----，计算N S 的通用程序。 （3）按两种顺序分别计算210S ，410S ，610S ，并指出有效位数。（编制程序时用单精度） （4）通过本上机题，你明白了什么？ 程序代码（matlab 编程）： clc clear a=single(1./([2:10^7].^2-1)); S1(1)=single(0); S1(2)=1/(2^2-1); for N=3:10^2 S1(N)=a(1); for i=2:N-1 S1(N)=S1(N)+a(i); end end S2(1)=single(0); S2(2)=1/(2^2-1); for N=3:10^2 S2(N)=a(N-1); for i=linspace(N-2,1,N-2) S2(N)=S2(N)+a(i); end end S1表示按从大到小的顺序的S N S2表示按从小到大的顺序的S N 计算结果 通过本上机题，看出按两种不同的顺序计算的结果是不相同的，按从大到小的顺序计算的值与精确值有较大的误差，而按从小到大的顺序计算的值与精确值吻合。从大到小的顺序计算得到的结果的有效位数少。计算机在进行数值计算时会出现“大数吃小数”的现象，导致计算结果的精度有所降低，我们在计算机中进行同号数的加法时，采用绝对值较小者先加的算法，其结果的相对误差较小。 计算方法上机报告 姓名： 学号： 班级： 目录 题目一------------------------------------------------------------------------------------------ - 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 - 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; 计算方法（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 = 如有帮助欢迎下载支持 数值分析上机题 姓名：陈作添 学号： 040816 习题 1 20．（上机题）舍入误差与有效数 N 1 1 3 1 1 设 S N ，其精确值为 。 2 2 2 N N 1 j 2 j 1 （1）编制按从大到小的顺序 1 1 1 ，计算 S 的通用程序。 S N 1 32 1 N 2 1 N 2 2 （2）编制按从小到大的顺序 1 1 1 ，计算 S 的通用程序。 S N 1 (N 1)2 1 22 1 N N 2 （3）按两种顺序分别计算 S 102 ， S 104 ， S 106 ，并指出有效位数。 （编制程序时用单精度） （4）通过本上机题，你明白了什么？ 按从大到小的顺序计算 S N 的通用程序为： 按从小到大的顺序计算 S N 的通用程序为： #include 西安电子科技大学出版社《计算方法》任传祥等编著第九章计算方法上机参考答案 实验一，算法一 #include #include 2015.1.9 上机作业题报告 JONMMX 2000 1.Chapter 1 1.1题目 设S N =∑1j 2?1 N j=2 ，其精确值为 )1 1 123(21+--N N 。 （1）编制按从大到小的顺序1 1 131121222-+ ??+-+-=N S N ，计算S N 的通用程序。 （2）编制按从小到大的顺序1 21 1)1(111222-+ ??+--+-= N N S N ，计算S N 的通用程序。 （3）按两种顺序分别计算64210,10,10S S S ,并指出有效位数。（编制程序时用单精度） （4）通过本次上机题，你明白了什么？ 1.2程序 1.3运行结果 1.4结果分析 按从大到小的顺序，有效位数分别为：6，4，3。 按从小到大的顺序，有效位数分别为：5，6，6。 可以看出，不同的算法造成的误差限是不同的，好的算法可以让结果更加精确。当采用从大到小的顺序累加的算法时，误差限随着N 的增大而增大，可见在累加的过程中，误差在放大，造成结果的误差较大。因此，采取从小到大的顺序累加得到的结果更加精确。 2.Chapter 2 2.1题目 （1）给定初值0x 及容许误差ε，编制牛顿法解方程f(x)=0的通用程序。 （2）给定方程03 )(3 =-=x x x f ,易知其有三个根3,0,3321= *=*-=*x x x ○1由牛顿方法的局部收敛性可知存在,0>δ当),(0δδ+-∈x 时，Newton 迭代序列收敛于根x2*。试确定尽可能大的δ。 ○2试取若干初始值，观察当),1(),1,(),,(),,1(),1,(0+∞+-----∞∈δδδδx 时Newton 序列的收敛性以及收敛于哪一个根。 （3）通过本上机题，你明白了什么？ 2.2程序 计算方法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β: 数值分析上机题(matlab版)(东南大学) 数值分析上机报告 第一章 一、题目 精确值为)1 1 123(21+--N N 。 1) 编制按从大到小的顺序 1 1 131121222-+??+-+-= N S N ，计算S N 的通用程序。 2) 编制按从小到大的顺序 1 21 1)1(111222-+??+--+-= N N S N ，计算S N 的通用程序。 3) 按两种顺序分别计算6 42 10,10, 10S S S ,并指出有效位 数。（编制程序时用单精度） 4) 通过本次上机题，你明白了什么？ 二、通用程序 clear N=input('Please Input an N (N>1):'); AccurateValue=single((0-1/(N+1)-1/N+3/2)/2); Sn1=single(0); for a=2:N; Sn1=Sn1+1/(a^2-1); end Sn2=single(0); for a=2:N; Sn2=Sn2+1/((N-a+2)^2-1); end fprintf('The value of Sn using different algorithms (N=%d)\n',N); disp('____________________________________________________') fprintf('Accurate Calculation %f\n',AccurateValue); fprintf('Caculate from large to small %f\n',Sn1); fprintf('Caculate from small to large %f\n',Sn2); 第一次作业 三、主观题(共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. 马斯洛的需要层次理论将人的需要分为了五个层次：（）、安全需要、爱的需要、（）和（）。 答案：生理需要，尊重需要，自我实现需要 如有帮助欢迎下载支持 数值分析上机题 姓名：陈作添 学号：040816 习题1 20．（上机题）舍入误差与有效数 设2 21 1N N j S j ==-∑ ，其精确值为1311221N N ??-- ?+?? 。 （1）编制按从大到小的顺序2 22 111 21311 N S N = +++---，计算N S 的通用程序。 （2）编制按从小到大的顺序2221111(1)121 N S N N =+++----，计算N S 的通用程序。 （3）按两种顺序分别计算210S ，410S ，610S ，并指出有效位数。（编制程序时用单精度） （4）通过本上机题，你明白了什么？ 按从大到小的顺序计算N S 的通用程序为： #include 实验一 矩阵的分解 一、实验目的 掌握矩阵的分解原理和一般方法，学会利用矩阵分解直接求解线性方程组。 二、实验内容 求矩阵() 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统计西安交大期末考试试题(含答案)
西安交通大学计算方法B上机试题
东南大学数值分析上机题答案
数值分析上机题目
西安交大成本会计在线作业答案精编版
成本会计-学习指南 西交大考试题库及答案
数值分析上机题目详解
东南大学《数值分析》-上机题
西安交通大学计算方法B大作业
(完整版)数值计算方法上机实习题答案
西交大计算方法上机报告
数值分析上机题参考答案.docx
西安电子科技大学出版社计算方法上机答案
东南大学-数值分析上机题作业-MATLAB版
西交计算方法A上机大作业
数值分析上机题(matlab版)(东南大学)
西安交大《管理心理学B》作业高分答案
数值分析上机题参考答案
西安交通大学计算方法A实验报告
- 西安交通大学有限元分析word版第一章
- 西安交通大学《计算方法》课件-第一章
- 西安交通大学计算方法A2012考试题
- 西安交通大学 计算方法上机报告
- 西安交通大学计算方法B上机试题
- 2010西安交通大学计算方法考题B(附答案)
- 西安交通大学计算方法B大作业资料
- 西安交通大学计算方法B上机报告
- 西安交通大学计算方法上机作业
- 西安交通大学计算方法A上机作业
- 西安交通大学计算方法10年考试题
- 西安交通大学2004年计算方法A试题
- 西安交大并行计算作业
- 西安交通大学公用房分配定额面积与核算办法
- 统计西安交大期末考试试题(含答案)
- 西交计算方法A上机大作业
- 西安交大计算方法B2017大作业
- 2010西安交通大学计算方法考题B(附答案)
- 西安交通大学计算方法C讲义
- 西交计算方法总结