2018年秋北航工科数分期末考题(含参考解答)

北京航空航天大学

2018－2019 学年第一学期期末考试

《工科数学分析（Ⅰ）》

（A卷）

班号学号姓名

主讲教师考场成绩

2019年01月 11日

一、 选择题（每题4分，满20分）

1. 曲线

sin cos ,0cos sin 2x t t t t y t t t

π

=-?≤≤?

=+?的弧长为（ D ）.

A. 2

;π B. 2;2π C. 2;4π D. 2

.8

π 2. 设()f x 满足等式1201()()d 1x

f x e f x x x

-=++?，则()f x =（ D ）. A. 121;12x e x -++ B. 12

1;14

x e x -++ C.

121;12x e x π-++ D. 12

1.14

x

e x π-++ 3. 设函数()

f x 可导，则0d ()d x e xf t dt x ??= ???

?（ C ）.

A.

();x xf e B. ();x x xe f e

C. 0

()();x

e x

x

xe f e f t dt +

?

D. 0

()().x

e x

xf e f t dt +?.

4. 下列广义积分中，收敛的是（ B ）.

A.

1

dx +∞?

； B.

3+arctan d 2x x

x x ∞+?；

C.

2

1

ln dx x x

+∞?

； D. +1

1

sin d .x x

∞

?

5. 设函数

()f x 在区间[,]a b 上可积，则下列4个结论中正确的是 （ A ）.

(1) 若()()F x f x '=

，则()()dF x F x C =+?；

(2)若[,][,]a b c d ?，则()()b

d

a

c

f x dx f x dx ≥?

?；

(3)若

()f x 为奇函数, 则其原函数()F x 必为偶函数； (4) 若

()f x 为周期函数, 则其原函数()F x 必为周期函数；

A. (1)(3)；

B. (2)(4)；

C. (1)(2)(3)；

D. (1)(2)(3)(4).

二、 计算题（每题6分，满分30分）

1.

2

2+3

d 310

x x x x +-?

解：

2

2+32+3310(5)(2)52

x x A B

x x x x x x ==++-+-+-, 计算可得1,1A B ==

22+3ln |5|ln |2|31052

x A B

dx dx dx x x C x x x x =+=++-++-+-?

?? 2.

2

arcsin d x x ?

解：22

arcsin arcsin xdx x x =-?

222arcsin 2arcsin (arcsin 2arcsin 2x x xd x x x dx x x x x C

=-=-=-+??

3. 2201921

21

sin cos d 1x x x x x

--+? 解： 由对称性，2201922

1

12211sin cos d d 11x x x x x x x

x ---=++?? 21

11

221

111d (1)d (arctan )2112x x x x x x x π---=-=-=-++?? 4.

220

2

d lim

cos(2)sin x t x e t

x x

→?

解：

22220

2

2

d d lim

lim

cos(2)sin x x t

t x x e t

e t x x

x

→→=?

?2

202lim

12x x xe

x

→==

5. 计算瑕积分

20

d .x ?

解： 2

是瑕点，220

00d lim arcsin

22

x

x ε

επ

+-→==

?

三、 （本题10分，每题5分） (1)利用定积分定义，求极限22

23

1lim (123)n n n →∞++++.

解

22

222223

11123

lim

(123)lim [()()()()]n n n

n n n n n n n

→∞→∞++++=+++

+1

2

01

=3x dx =? (2)求极限1

0lim d 1n

n x x x

→∞+?. 解法1：1100

100111

n n

n n

x x x dx x dx x x n ≤≤?≤≤=+++??，夹逼定理得到极限为0 解法2：

1

1

00

111

lim lim lim

01111n n n n n x dx x dx x n

ξ

ξ→∞→∞→∞==?=++++??

解法3：111001,1111n n

n

x x x dx dx dx x x x

εεε--=++++???为任意小于的正数 10

0(1)(1)11n n n x x ε

ξεεξ

-≤=-≤-++?

，夹逼定理得到此部分极限为0 1

1110,1n

x dx dx x εεεε--≤≤≤+??由的任意性可得此极限为0。故10lim d =01n

n x x x

→∞+?

四、 (本题8分)

(1)求二阶线性非齐次常微分方程332x

y y y e '''-+=的通解； (2)求上述方程满足(0)0,(0)1y y '==的特解 .

解： （1）对应齐次方程的特征方程为2

320λ

λ-+=, 特征根为1,2

故对应齐次方程的通解为212x

x y c e c e =+

设非齐次方程的特解为*

3x y

Ae =，代入原方程得1

2

A =

, 故此非齐次方程通解为231212

x

x

x

y c e c e e =++

（2）将初始条件代入方程通解得

121213

0,2122

c c c c ++

=++=，解得121,02c c =-=，故特解为31122

x x

y e e =-+

五、 (本题10分)

过点(1,1)作抛物线22y

x =-的切线，求由此切线，抛物线以及x 轴所围成的公共区域的面

积以及此区域绕x 轴旋转一周所得旋转体的体积.

解：(1)1

122x x y x =='=-=-，切线方程为23y x =-+，切线与x 轴交于点3,02

（），抛物线与x

轴交于点（）

±

所围成的曲边三角形面积为

(

)3

322

21

1

(23))3232123123x x dx x dx x x x S ?=-+--=-+- -=??

(2) 132221342

211(23))(23321)4)6V x dx x dx x x x dx ππ-??

????

=-+--=--++-???????????

5314=4653143+6159130x x x πππ??-+-? ????=???=-??

六、 (本题10分)

讨论无穷广义积分1

)

(1arctan )d 1

＋x x x x ∞

++?

的敛散性，若收敛，说明是绝对还是条

件收敛 .

解：积分

1

)

d 1

＋x x x ∞

+?

中，1()sin(2)d [1,)在上有界A F A x x =+∞?

0,1

1

当x

x

'

??

≤≥

?

+

??

，故

+

[1,)lim=0

11

在上单调且

x

x x

→∞

+∞

++

，由Dirichlet

判别法可知

1

)

d

1

＋x x

x

∞

+

?收敛

又由于1arctan在[1,+）上单调有界

x

+∞，故由Abel判别法可知

1

)

(1arctan)d

1

＋x x x

x

∞

+

+

?收敛

2 111

sin(2(2)

arctan)d d

1

＋＋＋

x

x

x x x x x x ∞∞∞

+≥≥

+

???

11

)

=d d

2(1)2(1)

＋＋

x

x x

x x

∞∞

-

++

??

2(1)

时，

x

x

→+∞

+

，比较判别法得出积分

1

d

2(1)

＋

x

x

∞

+

?发散，

类似于

1

)

d

1

＋x x

x

∞

+

?的收敛判别法得到

1

)

d

2(1)

＋x x

x

∞

+

?收敛

故

1

arctan)d

＋

x x

+

?发散，原积分条件收敛

七、(本题12分,每题6分)

(1)设()

f x在[0,1]上有一阶导数，且满足(0)0,0()1

'

=≤≤

f f x.

证明：

2

11

3

00

()d()d

??

≤??

??

??

f x x f x x.

(2) 设()

f x在[,]

a b上连续，在(,)

a b内可导，且满足()=

f a a，

22

1

()d()

2

=-

?b a f x x b a. 证明：至少存在一点(,)a b

ξ∈，使得()()1

f f

ξξξ

'=-+.

证明：（1）设

2

3

00

()()()

??

=-??

??

??

t t

F t f x dx f x dx，则

32

00

()()2()()()(()2())

'=-=-

??

t t

F t f t f t f x dx f t f t f x dx

设2

()()2()

=-?t

G t f t f x dx，则

()2()()2()2()(()1)0'''=-=-≤G t f t f t f t f t f t ，即()G t 单调减， ()(0)0≤=G t G

()()()0'∴=≤F t f t G t ，()()(0)0∴≤=单调减，F t F t F ，原不等式成立

（2）设()

(())-=-x F x e f x x ，则()0=F a

22

1()d ()d (())d 02

=-=?-=?

??b

b b a

a a f x x

b a x x f x x x

由微分中值定理，存在一点[,],η∈a b 使得()0ηη-=f .

若上式只对η

=a

成立，则

[,]()-≥≤在上恒(或）0,且不恒为0，

a b f x x 故(())d 00->

a

f x x x ，这与(())d =0-?b

a

f x x x 矛盾，故存在η≠a ，使

得

(),()0ηηη==即f F . 在[,]ηa 上对()F x 应用Rolle 定理得到结论

2018年北航数学专业考研(0701)考试科目、招生人数、参考书目、复习指导---新祥旭考研

2018年北航数学专业考研（0701）考试科目、招生人数、参考书目、 复习指导 一、招生信息 所属学院：数学与系统科学学院 招生人数：53 所属门类代码、名称：理学[07] 所属一级学科代码、名称：数学[0701] 二、研究方向： 01 代数学及其应用 02 复分析及其应用 03 泛函分析及其应用 04 偏微分方程及其应用 05 微分方程与动力系统 06 信息数学与科学计算 07 概率与数理统计 08 运筹学与控制论 三、初试考试科目： ①101思想政治理论 ②201英语一 ③609数学专业基础 ④891数学专业综合 四、参考书目 《高等代数》第三版高等教育出版社北京大学数学系编 《数学分析》(上册、下册) 高等教育出版社陈纪修等

五、复习指导 1、参考书的阅读方法 （1）目录法：先通读各本参考书的目录，对于知识体系有着初步了解，了解书的内在逻辑结构，然后再去深入研读书的内容。 （2）体系法：为自己所学的知识建立起框架，否则知识内容浩繁，容易遗忘，最好能够闭上眼睛的时候，眼前出现完整的知识体系。 （3）问题法：将自己所学的知识总结成问题写出来，每章的主标题和副标题都是很好的出题素材。尽可能把所有的知识要点都能够整理成问题。 2、学习笔记的整理方法 （1）第一遍学习教材的时候，做笔记主要是归纳主要内容，最好可以整理出知识框架记到笔记本上，同时记下重要知识点，如假设条件，公式，结论，缺陷等。记笔记的过程可以强迫自己对所学内容进行整理，并用自己的语言表达出来，有效地加深印象。第一遍学习记笔记的工作量较大可能影响复习进度，但是切记第一遍学习要夯实基础，不能一味地追求速度。第一遍要以稳、细为主，而记笔记能够帮助考生有效地达到以上两个要求。并且在后期逐步脱离教材以后，笔记是一个很方便携带的知识宝典，可以方便随时查阅相关的知识点。 （2）第一遍的学习笔记和书本知识比较相近，且以基本知识点为主。第二遍学习的时候可以结合第一遍的笔记查漏补缺，记下自己生疏的或者是任何觉得重要的知识点。再到后期做题的时候注意记下典型题目和错题。 （3）做笔记要注意分类和编排，便于查询。可以在不同的阶段使用大小合适的不同的笔记本。也可以使用统一的笔记本但是要注意各项内容不要混杂在以前，不利于以后的查阅。同时注意编好页码等序号。另外注意每隔一定时间对于在此期间自己所做的笔记进行相应的复印备份，以防原件丢失。统一的参考书书店可以买到，但是笔记是独一无二的，笔记是整个复习过程的心血所得，一定要好好保管。

北航2010-2011年研究生数值分析期末模拟试卷1-3

数值分析模拟试卷1 一、填空（共30分，每空3分） 1 设??? ? ??-=1511A ，则A 的谱半径=)(a ρ______，A 的条件数)(1A cond =________. 2 设 ,2,1,0,,53)(2==+=k kh x x x f k ，则],,[21++n n n x x x f ＝________, ],,[321+++n n n n x x x x f ，＝________. 3 设?????≤≤-++≤≤+=2 1,121 0,)(2 323x cx bx x x x x x S ，是以0，1，2为节点的三次样条函数，则b=________,c=________. 4 设∞=0)]([k k x q 是区间［0，1］上权函数为x x =)(ρ的最高项系数为1的正交多项式族，其中1)(0=x q ，则 ?=1 )(dx x xq k ________，=)(2 x q ________. 5 设???? ??????=11001a a a a A ，当∈a ________时，必有分解式，其中L 为下三角阵，当 其对角线元素)3,2,1(=i L ii 满足条件________时，这种分解是唯一的. 二、（14分）设4 9,1,41,)(2102 3 === =x x x x x f , （1）试求)(x f 在]4 9,41[上的三次Hermite 插值多项式)(x H 使满足 2,1,0),()(==i x f x H i i ，)()(11x f x H '='. （2）写出余项)()()(x H x f x R -=的表达式. 三、（14分）设有解方程0cos 2312=+-x x 的迭代公式为n n x x cos 3 2 41+ =+， （1） 证明R x ∈?0均有? ∞ →=x x n x lim （? x 为方程的根）； （2） 取40=x ，用此迭代法求方程根的近似值，误差不超过，列出各次迭代值； （3）此迭代的收敛阶是多少？证明你的结论. 四、(16分) 试确定常数A ，B ，C 和，使得数值积分公式 有尽可能高的代数精度. 试问所得的数值积分公式代数精度是多少？它是否为Gauss 型的？

北航数值分析大作业一

《数值分析B》大作业一 SY1103120 朱舜杰 一．算法设计方案： 1.矩阵A的存储与检索 将带状线性矩阵A[501][501]转存为一个矩阵MatrixC[5][501] . 由于C语言中数组角标都是从0开始的，所以在数组MatrixC[5][501]中检索A的带内元素a ij的方法是： A的带内元素a ij=C中的元素c i-j+2,j 2.求解λ1，λ501，λs ①首先分别使用幂法和反幂法迭代求出矩阵按摸最大和最小的特征值λmax和λmin。λmin即为λs； 如果λmax>0,则λ501=λmax；如果λmax<0,则λ1=λmax。 ②使用带原点平移的幂法（mifa（）函数），令平移量p=λmax，求 出对应的按摸最大的特征值λ，max， 如果λmax>0,则λ1=λ，max+p；如果λmax<0,则λ501=λ，max+p。 3.求解A的与数μk=λ1+k（λ501-λ1）/40的最接近的特征值λik （k=1,2，…，39）。 使用带原点平移的反幂法，令平移量p=μk，即可求出与μk最接近的特征值λik。 4.求解A的（谱范数）条件数cond（A）2和行列式d etA。 ①cond（A）2=|λ1/λn|，其中λ1和λn分别是矩阵A的模最大和 最小特征值。

②矩阵A的行列式可先对矩阵A进行LU分解后，detA等于U所有对角线上元素的乘积。 二．源程序 #include #include #include #include #include #include #include #define E 1.0e-12 /*定义全局变量相对误差限*/ int max2(int a,int b) /*求两个整型数最大值的子程序*/ { if(a>b) return a; else return b; } int min2(int a,int b) /*求两个整型数最小值的子程序*/ { if(a>b) return b; else return a; } int max3(int a,int b,int c) /*求三整型数最大值的子程序*/ { int t; if(a>b) t=a; else t=b; if(t

理工科大学物理知识点总结及典型例题解析

理工科大学物理知识点总结及典型例题解析

第一章 质点运动学 本章提要 1、 参照系：描述物体运动时作参考的其他物体。 2、 运动函数：表示质点位置随时间变化的函数。 位置矢量： k t z j t y i t x t r r )()()()( 位置矢量：)()(t r t t r r 一般情况下：r r 3、速度和加速度： dt r d v ； 2 2dt r d dt v d a 4、匀加速运动： a 常矢量 ； t a v v 0 2 2 10 t a t v r 5、一维匀加速运动： at v v 0 ； 2 210at t v x ax v v 22 02 6、抛体运动： x a ； g a y cos 0v v x ； gt v v y sin 0 t v x cos 0 ； 2 210sin gt t v y 7、圆周运动：t n a a a 法向加速度： 2 2 R R v a n 切向加速度：dt dv a t 8、伽利略速度变换式：u v v 【典型例题分析与解 答】

m j t i t j t i t r r ]2)310[(2322220 (2)由以上可得质点的运动方程的分量式x=x(t) y=y(t) 即 x=10+3t 2 y=2t 2 消去参数t, 3y=2x-20 这是一个直线方程.由 m i r 100 知 x 0=10m,y 0=0.而直线斜率 3 2 tga dy/dx k , 则1433 a 轨迹方程如图所示 3. 质点的运动方程为2 3010t t -x 和2 2015t t-y ,(SI)试求:(1) 初速度的大小和方向;(2)加速度的大小和方向. 解.(1)速度的分量式为 t -dx/dt v x 6010 t -dy/dt v y 4015 当t=0时,v 0x =-10m/s,v 0y =15m/s,则初速度的大小为0182 02 00 .v v v y x m/s 而v 0与x 轴夹角为 1412300 x y v v arctg a (2)加速度的分量式为 260-x x ms dt dv a 2 40-y y ms dt dv a 则其加速度的大小为 17222 . a a a y x ms -2 a 与x 轴的夹角为 1433 -a a arctg x y (或91326 ) X 10

北航数学分析期末考试卷

A 一、填空题（每题5分，共30分） 1． 设向量场),,(222xyz z xy yz x A =，求=divA =rotA 2．求=+?→x x dx ααcos 12100lim 3．设),(y x f 在原点领域连续， 求极限=??≤+→dxdy y x f y x ),(12222 0lim ρρπρ 4．设为自然数，n z y x z y x D },10,10,10|),,{(≤≤≤≤≤≤= 求=+++???dxdydz z y x y x n n n n n D 5．设,)(2)1(cos sin dt e x f t x x +?= 求=)('x f 6．)为右半单位圆 设L (,sin cos :???==θ θy x L 求=?ds y L || 二、（本题满分１0分） 设Ω为椭球体,1222222≤++c z b y a x 计算dxdydz xy z I )2(2+=???Ω

三（本题满分１0分） 计算曲面积分,)(dS z y x ++??∑ 其中∑是平面5=+z y 被柱面2522=+y x 所 截得的部分。 四（本题满分3０分，每题10分） 1． 计算曲线积分 ?-+-+-=L dz y x dy x z dx z y I ,)()()(02222=++=++z y x a z y x L 与平面是球面其中取逆时针方向。轴正向看去的交线，从L z

2．计算曲面积分.zdxdy ydzdx xdydz ++??∑ 其中)0(:22h y z x y ≤≤+=∑，方 向取左侧。 3．计算,4)4()(.22y x dy y x dx y x L +++-?其中L 为单位圆周,.122=+y x 方向为逆时针方向。

北航数理统计回归分析大作业

应用数理统计第一次大作业 学号： 姓名： 班级： 2013年12月

国家财政收入的多元线性回归模型 摘 要 本文以多元线性回归为出发点，选取我国自1990至2008年连续19年的财政收入为因变量，初步选取了7个影响因素,并利用统计软件PASW Statistics 17.0对各影响因素进行了筛选，最终确定了能反映财政收入与各因素之间关系的“最优”回归方程： 46?578.4790.1990.733y x x =++ 从而得出了结论，最后我们用2009年的数据进行了验证，得出的结果在误差范围内，表明这个模型可以正确反映影响财政收入的各因素的情况。 关键词：多元线性回归，逐步回归法，财政收入，SPSS 0符号说明 变 量 符号 财政收入 Y 工 业 X 1 农 业 X 2 受灾面积 X 3 建 筑 业 X 4 人 口 X 5 商品销售额 X 6

进出口总额X7

1 引言 中国作为世界第一大发展中国家，要实现中华民族的伟大复兴，必须把发展放在第一位。近年来，随着国家经济水平的飞速进步，人民生活水平日益提高，综合国力日渐强大。经济上的飞速发展并带动了国家财政收入的飞速增加，国家财政的状况对整个社会的发展影响巨大。政府有了强有力的财政保证才能够对全局进行把握和调控，对于整个国家和社会的健康快速发展有着重要的意义。所以对国家财政的收入状况进行研究是十分必要的。 国家财政收入的增长，宏观上必然与整个国家的经济有着必然的关系，但是具体到各个方面的影响因素又有着十分复杂的相关原因。为了研究影响国家财政收入的因素，我们就很有必要对其财政收入和影响财政收入的因素作必要的认识，如果能对他们之间的关系作一下回归，并利用我们所知道的数据建立起回归模型这对我们很有作用。而影响财政收入的因素有很多，如人口状况、引进的外资总额，第一产业的发展情况，第二产业的发展情况，第三产业的发展情况等等。本文从国家统计信息网上选取了1990-2009年这20年间的年度财政收入及主要影响因素的数据，包括工业，农业，建筑业，批发和零售贸易餐饮业，人口总数等。文中主要应用逐步回归的统计方法，对数据进行分析处理，最终得出能够反映各个因素对财政收入影响的最“优”模型。 2解决问题的方法和计算结果 2.1 样本数据的选取与整理 本文在进行统计时，查阅《中国统计年鉴2010》中收录的1990年至2009年连续20年的全国财政收入为因变量，考虑一些与能源消耗关系密切并且直观上

北航大一上工科数分期中考试试卷

北京航空航天大学2011-2012学年第一学期期中考试 工科数学分析试卷（2011.12.25） 一、计算(5’*8=40’) 1) 用Stolz 定理计算极限41233122123lim n n n n n +→∞++++L . 2) 设32()(1)x f x x x x =++,求()f x '. 3) 求极限1 0(1)e lim x x x x →+-. 4) 求函数2()(4)f x x x = -的拐点。 5) 设(cos sin )()=(sin cos )x a t t f x y a t t t =+??=-?,求d d y x . 6) 求函数()ln f x x x =在(0,)+∞上的最值. 7) 判断函数21 1()=e x n f x x -?间断点的类型. 8) 求函数2()=ln(1)f x x x ++在0x =处直到四阶的Taylor 展开（Peano 余项形式）. 二、证明(15’) 1) 3 sin (0)6 x x x x >-> 2) 设函数1()=ln ()n f x x x n -+∈￠,证明()(1)!n n y x -=. 三、(10’) 设1110,0,(2),1,2,n n n A x x x Ax n A +><<=-=L ,证明不等式11n n x x A +<<对任意

n +∈￠成立，并求出极限lim n n x →∞ . 四、(10’) 用Cauchy 收敛原理证明数列2sin (sin )n n k kx x k k kx == +∑收敛. 五、(15’) 设()f x 在0x 处二次可导，且()0f x ''≠，由Lagrange 中值定理知存在0()1h θ<<,使得式子000(+)()(())f x h f x f x h h h θ'=++成立，计算或者证明下列结论： 1) 写出()f x 和()f x '在0x x =处的Taylor 公式； 2) 证明01lim ()2 h h θ→=. 六、(10’) 设()f x '在(0,]a 连续，且极限lim ()x x →'存在，证明()f x 在(0,]a 上一致连续. [附加题] 七、(10’) 以下题目任选其一： 1) 设()[01]f x ∈￡，，且()0f x >,令0()max (),[0,1]t x M x f t x ≤≤=∈, 证明：函数()()lim ()n n f x Q x M x →∞??=???? 连续的充要条件是()f x 单调递增. 2) 证明开区间套定理 1. 设开区间序列(,),n n n I a b n +=∈￥ 满足12121n n n a a a b b b b -<<<<<<<

北航数理统计第二次大作业-数据分析模板

数理统计第二次大作业材料行业股票的聚类分析与判别分析 2015年12月26日

材料行业股票的聚类分析与判别分析摘要

1 引言 2 数据采集及标准化处理 2.1 数据采集 本文选取的数据来自大智慧软件的股票基本资料分析数据，从材料行业的股票中选取了30支股票2015年1月至9月的7项财务指标作为分类的自变量，分别是每股收益（单位：元）、净资产收益率（单位：%）、每股经营现金流（单位：元）、主营业务收入同比增长率（单位：%）、净利润同比增长率（单位：%）、流通股本（单位：万股）、每股净资产（单位：元）。各变量的符号说明见表2.1，整理后的数据如表2.2。 表2.1 各变量的符号说明 自变量符号 每股收益（单位：元）X1 净资产收益率（单位：%）X2 每股经营现金流（单位：元）X3 主营业务收入同比增长率（单位：%）X4 净利润同比增长率（单位：%）X5 流通股本（单位：万股）X6 每股净资产（单位：元）X7 表2.2 30支股票的财务指标 股票代码X1 X2 X3 X4 X5 X6 X7 武钢股份600005-0.0990-2.81-0.0237-35.21-200.231009377.98 3.4444宝钢股份6000190.1400 1.980.9351-14.90-55.011642427.88 6.9197山东钢铁600022-0.11650.060.0938-20.5421.76643629.58 1.8734北方稀土6001110.0830 3.640.652218.33-24.02221920.48 2.2856

杭钢股份600126-0.4900-13.190.4184-36.59-8191.0283893.88 3.4497抚顺特钢6003990.219310.080.1703-14.26714.18112962.28 1.4667盛和资源6003920.0247 1.84-0.2141-5.96-19.3739150.00 1.2796宁夏建材6004490.04000.510.3795-22.15-92.3447818.108.7321宝钛股份600456-0.2090-2.53-0.3313-14.81-6070.2043026.578.1497山东药玻6005290.4404 5.26 1.2013 6.5016.7825738.018.5230国睿科技6005620.410011.53-0.2949 3.3018.9416817.86 3.6765海螺水泥600585 1.15169.05 1.1960-13.06-25.33399970.2612.9100华建集团6006290.224012.75-0.57877.90-6.4034799.98 1.8421福耀玻璃6006600.790014.250.9015 3.6017.27200298.63 6.2419宁波富邦600768-0.2200-35.02-0.5129 3.1217.8813374.720.5188马钢股份600808-0.3344-11.710.3939-21.85-689.22596775.12 2.6854亚泰集团6008810.02000.600.1400-23.63-68.16189473.21 4.5127博闻科技6008830.503516.71-0.1010-10.992612.8023608.80 3.0126新疆众和6008880.0523 1.04-0.910662.64162.0464122.59 5.0385西部黄金6010690.0969 3.940.115115.5125.5712600.00 2.4965中国铝业601600-0.0700-2.920.2066-9.0882.79958052.19 2.3811明泰铝业6016770.2688 4.66-1.09040.8227.8640770.247.4850金隅股份6019920.1989 3.390.3310-10.05-39.01311140.26 6.7772松发股份6032680.35007.00-0.3195-4.43-9.622200.00 6.0244方大集团0000550.0950 5.66-0.480939.2920.6742017.94 1.6961铜陵有色0006300.0200 1.220.6132 3.23-30.74956045.21 1.5443鞍钢股份000898-0.1230-1.870.7067-27.32-196.21614893.17 6.4932中钢国际0009280.572714.45-0.4048-14.33410.2441286.57 4.2449中材科技0020800.684610.27 1.219547.69282.1740000.00 6.8936中南重工0024450.1100 4.300.340518.8445.0950155.00 2.7030 2.2 数据的标准化处理 由于不同的变量之间存在着较大的数量级的差别，因此要对数据变量进行标准化处理。本文采用Z得分值法标准化的方法进行标准化，用x的值减去x的均值再除以样本的方差。也就是把个案转换为样本均值为0、标准差为1的样本。如果不同变量的变量值数值相差太大，会导致计算个案间距离时，由于绝对值较小的数值权数较小，个案距离的大小几乎由大数值决定，标准化过程可以解决此类问题，使不同变量的数值具有同等的重要性。经Z标准化输出结果见表 2.2。 表2.2 经Z标准化后的数据 ZX1ZX2ZX3ZX4ZX5ZX6ZX7

北航数值分析报告第三次大作业

数值分析第三次大作业 一、算法的设计方案： （一）、总体方案设计： x y当作已知量代入题目给定的非线性方程组，求（1）解非线性方程组。将给定的(,) i i

得与(,)i i x y 相对应的数组t[i][j],u[i][j]。 （2）分片二次代数插值。通过分片二次代数插值运算，得到与数组t[11][21],u[11][21]]对应的数组z[11][21]，得到二元函数z=(,)i i f x y 。 （3）曲面拟合。利用x[i],y[j],z[11][21]建立二维函数表，再根据精度的要求选择适当k 值，并得到曲面拟合的系数矩阵C[r][s]。 （4）观察和(,)i i p x y 的逼近效果。观察逼近效果只需要重复上面（1）和（2）的过程，得到与新的插值节点(,)i i x y 对应的(,)i i f x y ，再与对应的(,)i i p x y 比较即可，这里求解 (,)i i p x y 可以直接使用（3）中的C[r][s]和k 。 （二）具体算法设计： （1）解非线性方程组 牛顿法解方程组()0F x =的解* x ，可采用如下算法： 1）在* x 附近选取(0) x D ∈，给定精度水平0ε>和最大迭代次数M 。 2）对于0,1, k M =执行 ① 计算() ()k F x 和()()k F x '。 ② 求解关于() k x ?的线性方程组 () ()()()()k k k F x x F x '?=- ③ 若() () k k x x ε∞∞ ?≤，则取*()k x x ≈，并停止计算；否则转④。 ④ 计算(1) ()()k k k x x x +=+?。 ⑤ 若k M <，则继续，否则，输出M 次迭代不成功的信息，并停止计算。 （2）分片双二次插值 给定已知数表以及需要插值的节点，进行分片二次插值的算法： 设已知数表中的点为： 00(0,1,,) (0,1,,)i j x x ih i n y y j j m τ=+=???=+=?? ，需要插值的节点为(,)x y 。 1) 根据(,)x y 选择插值节点(,)i j x y ： 若12h x x ≤+ 或12 n h x x ->-，插值节点对应取1i =或1i n =-，

北航2015-2016年工科数分(1)期末_A卷_答案

北京航空航天大学 2015－2016 学年第一学期期末考试 《工科数学分析（Ⅰ）》 （A卷） 班号学号姓名 主讲教师考场成绩 2016年01月20日

1. 下列命题中错误的是 （ D ） A. 若()f x 在区间(,)a b 内的原函数是常数，则()f x 在(,)a b 内恒为0; B. 若],[)(b a x f 在上可积, 则],[)(b a x f 在上必有界 ； C. 若],[)(b a x f 在上可积, 则()f x 在区间[,]a b 上也可积 ； D. 若],[)(b a x f 在上不连续，则],[)(b a x f 在上必不可积 . 2. 设 ()f x 满足等式1 2 0()2()d f x x f x x =-?，则1 ()d f x x ?=（ B ） A. 1; B. 1;9 C. 1;- D. 1 .3 - 3. 设函数()f x 可导，则（ C ） A. ()d ();f x x f x =? B. ()d ();f x x f x '=? C. () d ()d ();d f x x f x x =? D. () d ()d ().d f x x f x C x =+? 4. 下列广义积分中，发散的是（ C ） A. 1 dx +∞ ? ； B. 21 1 dx x +∞? ； C. 1 1sin d x x x +∞ +? ； D. 1 sin d .x e x x +∞ -? 5. 瑕积分 3 1 ln dx x x =? （ C ） A. l n l n 3; B. 0; C. ;+∞ D. 1.

北航数值分析大作业第二题

数值分析第二次大作业 史立峰 SY1505327

一、 方案 （1）利用循环结构将sin(0.50.2)() 1.5cos( 1.2)() {i j i j ij i j i j a +≠+==(i,j=1,2,……,10)进行赋值，得到需要变换的 矩阵A ； （2）然后，对矩阵A 利用Householder 矩阵进行相似变换，把A 化为上三角矩阵A (n-1)。 对A 拟上三角化，得到拟上三角矩阵A (n-1)，具体算法如下： 记A(1)=A ，并记A(r)的第r 列至第n 列的元素为()n r r j n i a r ij ,,1,;,,2,1) ( +==。 对于2,,2,1-=n r 执行 1. 若 ()n r r i a r ir ,,3,2) ( ++=全为零，则令A(r+1) =A(r),转5；否则转2。 2. 计算 () ∑+== n r i r ir r a d 1 2 )( ()( )r r r r r r r r r r d c a d a c ==-=++则取,0sgn ) (,1)(,1若 )(,12r r r r r r a c c h +-= 3. 令 () n T r nr r r r r r r r r R a a c a u ∈-=++) ()(,2)(,1,,,,0,,0 。 4. 计算 r r T r r h u A p /)(= r r r r h u A q /)(= r r T r r h u p t /= r r r r u t q -=ω T r r T r r r r p u u A A --=+ω)()1( 5. 继续。 （3）使用带双步位移的QR 方法计算矩阵A (n-1)的全部特征值，也是A 的全部特征值，具体算法如下： 1. 给定精度水平0>ε和迭代最大次数L 。 2. 记n n ij n a A A ?-==][) 1()1()1(，令n m k ==,1。

北航数值分析报告大作业第八题

北京航空航天大学 数值分析大作业八 学院名称自动化 专业方向控制工程 学号 学生姓名许阳 教师孙玉泉 日期2014 年11月26 日

一．题目 关于x , y , t , u , v , w 的方程组(A.3) ???? ?? ?=-+++=-+++=-+++=-+++79 .0sin 5.074.3cos 5.007.1cos sin 5.067.2cos 5.0y w v u t x w v u t y w v u t x w v u t (A.3) 以及关于z , t , u 的二维数表（见表A-1）确定了一个二元函数z =f (x , y )。 表A-1 二维数表 t z u 0 0.4 0.8 1.2 1.6 2 0 -0.5 -0.34 0.14 0.94 2.06 3.5 0.2 -0.42 -0.5 -0.26 0.3 1.18 2.38 0.4 -0.18 -0.5 -0.5 -0.18 0.46 1.42 0.6 0.22 -0.34 -0.58 -0.5 -0.1 0.62 0.8 0.78 -0.02 -0.5 -0.66 -0.5 -0.02 1.0 1.5 0.46 -0.26 -0.66 -0.74 -0.5 1. 试用数值方法求出f (x , y ) 在区域}5.15.0,8.00|), {≤≤≤≤=y x y x D （上的近似表达式 ∑∑===k i k j s r rs y x c y x p 00 ),( 要求p (x , y )以最小的k 值达到以下的精度 ∑∑==-≤-=10020 7210)],(),([i j i i i i y x p y x f σ 其中j y i x i i 05.05.0,08.0+==。 2. 计算),(),,(* ***j i j i y x p y x f (i =1,2,…,8 ; j =1,2,…,5) 的值，以观察p (x , y ) 逼 近f (x , y )的效果，其中j y i x j i 2.05.0,1.0**+==。

北航数理统计回归分析大作业

数理统计(课程大作业1) 逐步回归分析 学院：机械工程学院 专业：材料加工工程 日期：2014年12月7日

摘要：本文介绍多元线性回归分析方法以及逐步回归法，然后结合实际，以我国1995-2012年的财政收入为因变量，选取了8个可能的影响因素，选用逐步回归法对各影响因素进行了筛选分析，最终确定了其“最优”回归方程。 关键字：多元线性回归 逐步回归法 财政收入 SPSS 1 引言 自然界中任何事物都是普遍联系的，客观事物之间往往都存在着某种程度的关联关系。为了研究变量之间的相关关系，人们常用回归分析的方法，而回归分析是数理统计中一种常用方法。数理统计作为一种实用有效的工具，广泛应用于国民经济的各个方面，在解决实际问题中发挥了巨大的作用，是一种理论联系实践、指导实践的科学方法。 财政收入，是指政府为履行其职能、实施公共政策和提供公共物品与服务需要而筹集的一切资金的总和。财政收入表现为政府部门在一定时期内（一般为一个财政年度）所取得的货币收入。财政收入是衡量一国政府财力的重要指标，政府在社会经济活动中提供公共物品和服务的范围和数量，在很大程度上决定于财政收入的充裕状况。 本文将以回归分析为方法，运用数理统计工具探求财政收入与各种统计指标之间的关系，总结主要影响因素，并对其作用、前景进行分析和展望。 2 多元线性回归 2.1 多元线性回归简介 在实际问题中，某一因素的变化往往受到许多因素的影响，多元回归分析的任务就是要找出这些因素之间的某种联系。由于许多非线性的情形都可以通过变换转化为线性回归来处理，因此，一般的实际问题都是基于多元线性回归问题进行处理的。对多元线性回归模型简要介绍如下： 如果随机变量y 与m )2(≥m 个普通变量m x x x 21,有关，且满足关系式： εββββ++++=m m x x x y 22110 2,0σεε==D E （2.1） 其中，2210,,,σββββm 是与m x x x 21,无关的未知参数，ε是不可观测的随机变量，),0(~2N I N σε。

数值分析

习 题 1. 指出有效数49×102，0.0490，490.00的绝对误差限、相对误差限和有效数字位数． 2. 将 3.142作为π的近似值，它有几位有效数字，相对误差限和绝对误差限各为多少？ 3. 要使101的近似值x * 的相对误差限不超过4102 1?×，问查开方表时x * 需要保留几位有效数字？ 4. 已知近似数x * 有两位有效数字，试估计其相对误差限． 5. 设x * 为x 的近似数, 证明n x * 的相对误差大约为x * 相对误差的n 1倍． 6. 某矩形的长和宽大约为100cm 和50cm, 应该选用最小刻度为多少cm 的测量工具, 才能保证计算出的面积误差（绝对值）不超过0.15cm 2． 7. 已知三角形面积c ab S sin 2 1=，测量a , b , c 时产生的相对误差为)(*a e r ，)(*b e r ，)(*c e r ，其中2 ,0*π<>2时的情形．用所设计的算法以及二次方程求根公式计算05.240=p ，00.1=q 时方程根的近似值(计算过程保留2位小数)，并给出它们的相对误差限（根的准确值为L 0916683.4801?=x , L 002082935.02?=x ）．

北航数值分析第二次大作业--QR分解

《数值分析A》

一、算法设计方案 整个程序主要分为四个函数，主函数，拟上三角化函数，QR分解函数以及使用双步位移求解矩阵特征值、特征向量的函数。因为在最后一个函数中也存在QR分解，所以我没有采用参考书上把矩阵M进行的QR分解与矩阵Ak的迭代合并的方法，而是在该函数中调用了QR分解函数，这样增强了代码的复用性，减少了程序长度；但由于时间关系，对阵中方法的运算速度没有进行深入研究。 1.为了减少QR分解法应用时的迭代次数，首先对给定矩阵进行拟上三角化处理。 2.对经过拟上三角化处理的矩阵进行QR分解。 3.注意到计算特征值与特征向量的过程首先要应用前面两个函数，于是在拟上三角化矩阵的基础上对QR分解函数进行了调用。计算过程中，没有采用goto语句，而是根据流程图采用其他循环方式完成了设计，通过对迭代过程的合并，简化了程序的循环次数，最后在计算特征向量的时候采用了列主元高斯消去法。

二、源程序代码 #include #include #include int i,j,k,l,m; //定义外部变量double d,h,b,c,t,s; double A[10][10],AA[10][10],R[10][10],Q[10][10],RQ[10][10]; double X[10][10],Y[10][10],Qt[10][10],M[10][10]; double U[10],P[10],T[10],W[10],Re[10]={0},Im[10]={0}; double epsilon=1e-12; void main() { void Quasiuppertriangular(double A[][10]); void QRdecomposition(double A[][10]); void DoublestepsQR(double A[][10]); int i,j; for(i=0;i<10;i++) { for(j=0;j<10;j++) { A[i][j]=sin(0.5*(i+1)+0.2*(j+1)); Q[i][j]=0; AA[i][j]=A[i][j]; } A[i][i]=1.5*cos(2.2*(i+1)); AA[i][i]=A[i][i];

北航数值分析大作业第二次

《数值分析》计算实习作业 （第二题）

算法设计方案： 1、对矩阵A 赋值，取计算精度ε=1×10-12； 2、对矩阵A 进行拟上三角化，得到A (n-1)，并输出A (n-1)； 对矩阵A 的拟上三角化，通过直接调用子函数inftrianglize(A)来实现；拟上三角化得到的矩阵A (n-1)输出至文件solution.txt 中。 3、对A (n-1)进行QR 分解并输出Q 、R 及RQ 矩阵； QR 分解通过直接调用子函数QRdescom(A,Q,R, n)实现。 4、运用QR 方法求所有的特征值，并输出； (1)初始时令m=n ，在m>2的条件下执行； (2)判断如果|A mm-1|<ε，则得到一个特征值，m=m-1，转(4)；否则转(3)； (3)判断如果|A m-1m-2|<ε，则得到两个特征值，m=m-2，转(4)； (4)判断如果m ≤2，转(6)；否则转(5)； (5)执行相似迭代，转(2)； k k T k k k k k k k k k k Q A Q A R Q M I D A D tr A M ==+-=+1)2)det(( (6)求出最后的一个或两个特征值； (7)输出全部的特征值至文件solution.txt 中。 5、输出QR 分解法迭代结束之后的A (n-1)至文件solution.txt 中； 6、通过反幂法求出所有实特征值的特征向量并输出。 首先令B=(A-λi I)，其中λi 是实特征值；反幂法通过调用子函数Bpowmethod(B,x1)实现，最终λi 对应的特征向量就是x1；最后将所有的实特征值的特征向量输出。

北航数值分析课程第一次大作业讲解

《数值分析A》计算实习题目第一题 一．算法设计方案： 1.矩阵A的存储与检索 将带状线性矩阵A[501][501]转存为一个矩阵MatrixC[5][501] . 由于C语言中数组角标都是从0开始的，所以在数组MatrixC[5][501]中检索A的带内元素a ij的方法是： A的带内元素a ij=C中的元素c i-j+2,j 2.求解λ1，λ501，λs ①首先分别使用幂法和反幂法迭代求出矩阵按摸最大和最小的特征值λmax和λmin。λmin即为λs； 如果λmax>0,则λ501=λmax；如果λmax<0,则λ1=λmax。 ②使用带原点平移的幂法（mifa（）函数），令平移量p=λmax，求出对应的按摸最大的特征值λ，max， 如果λmax>0,则λ1=λ，max+p；如果λmax<0,则λ501=λ，max+p。 3.求解A的与数μk=λ1+k（λ501-λ1）/40的最接近的特征值λik （k=1,2，…，39）。 使用带原点平移的反幂法，令平移量p=μk，即可求出与μk最接近的特征值λik。 4.求解A的（谱范数）条件数cond（A）2和行列式d etA。 ①cond（A）2=|λ1/λn|，其中λ1和λn分别是矩阵A的模最大和最小特征值。 ②矩阵A的行列式可先对矩阵A进行LU分解后，detA等于U所有

对角线上元素的乘积。 二．源程序(VS2010环境下，C++语言) #include #include #include #include #include #include #include #define E 1.0e-12 /*定义全局变量相对误差限*/ int max2(int a,int b) /*求两个整型数最大值的子程序*/ { if(a>b) return a; else return b; } int min2(int a,int b) /*求两个整型数最小值的子程序*/ { if(a>b) return b; else return a; } int max3(int a,int b,int c) /*求三整型数最大值的子程序*/ { int t; if(a>b) t=a; else t=b; if(t

北航-数理统计大作业

对中国各地财政收入情况的聚类分析和判 别分析 应用数理统计第二次大作业 学院名称 学号 学生姓名 摘要 我国幅员辽阔，由于人才、地理位置、自然资源等条件的不同，各地区的财政收入类型各自呈现出不一样的发展趋势，通过准确定位中国各地区财政收入情况对于正确认识我国财政收入具有重要的意义。本文以中国各地财政收入情况为研究对象，从《中国统计年鉴》中选取2011年期间中国各地财政收入情况为因

变量，选取国内增值税、营业税、企业所得税、个人所得税、城市维护建设税、土地增值税、契税、专项收入、行政事业性收费收入、国有资本经营收入和国有资源（资产）有偿使用收入11个可能影响中国各地财政收入的因素为自变量，利用统计软件SPSS，对27个地区的财政收入进行了聚类分析，并对另外4个地区的财政收入进行了判别分析，并最终确定了中国各地区根据财政收入类型的分类情况。 关键词：聚类分析，判别分析，SPSS，中国各地财政收入类型 1、引言 财政收入，是指政府为履行其职能、实施公共政策和提供公共物品与服务需要而筹集的一切资金的总和。财政收入表现为政府部门在一定时期内（一般为一个财政年度）所取得的货币收入。财政收入是衡量一国政府财力的重要指标，政府在社会经济活动中提供公共物品和服务的范围和数量，在很大程度上决定于财政收入的充裕状况。通过准确定位中国各地区财政收入情况对于正确认识我国财政收入具有重要的意义。 本文利用统计软件SPSS，根据各地区的财政收入情况，对北京、天津、河北等27个地区进行聚类分析，并对青海、重庆、四川、贵州4个省市进行判别分析，判断属于聚类分析结果中的哪种财政收入类型。 1.1 聚类分析 聚类分析是根据研究对象的特征对研究对象进行分类的多元统计分析技术的总称，它直接比较各事物之间的性质，将性质相近的归为一类，将性质差别较大的归入不同的类。本文采用的是系统聚类分析，它又称集群分析，是聚类分析中应用最广的一种方法，其基本思想是：首先将每个聚类对象看作一类，然后根据对象间的相似程度，将相似程度最高的两类进行合并，并计算合并后的类与其他类之间的距离，再选择相近者进行合并，每合并一次减少一类，直至所有的对象都并为一类为止。 系统聚类分为Q型聚类和R型聚类两种：Q型聚类是对样本进行聚类，它使具有相似特征的样本聚集在一起，使差异性大的样本分离开来；R型聚类是对变量进行聚类，它使差异性大的变量分离开来，相似的变量聚集在一起，这样就