矩阵 计算习题及答案
1、选择题
1）下列变量中 A 是合法的。
A. Char_1,i,j
B.x*y,a.1
C. X\y, a1234
D. end, 1bcd
2）下列 C 是合法的常量。
A. 3e10
B. 1e500
C. -1.85e-56
D. 10-2
3）x=uint8(1.2e10)，则x所占的字节是 D 个。
A. 1
B. 2
C. 4
D. 8
4）已知x=0：10，则x有 B 个元素。
A. 9
B. 10
C. 11
D. 12
5）产生对角线元素全为1其余为0的2×3矩阵的命令是 C 。
A. Ones(2,3)
B. Ones(3,2)
C. Eye(2,3)
D. Eye(3,2)
6）a=
123
456
789
??
?
?
?
??
，则a(:,end)是指 C 。
A.所有元素
B. 第一行元素
C. 第三列元素
D. 第三行元素
7）a=
123
456
789
??
?
?
?
??
，则运行a(:,1)=[] 命令后 C 。
A.a变成行向量
B. a数组成2行2列
C. a数组成3行2列
D. a数组没有元素
8）a=
123
456
789
??
?
?
?
??
，则运行命令mean(a)是 B 。
A. 计算a的平均值
B. 计算a每列的平均值
C. 计算a每行的平均值
D.a数组增加一列平均值
9）已知x是一个向量，计算ln(x)的命令是 B 。
A. ln(x)
B. log(x)
C. Ln(x)
D. lg10(x)
10）当a=2.4时，使用取整函数得到3，则该函数名是 C 。
A.fix
B. round
C. ceil
D. floor
11）已知a=0：4，b=1：5，下面的运算表达式出错的是 D 。
A. a+b
B. a./b
C. a'*b
D. a*b
12）已知a=4，b=‘4’，下面说法错误的是 C 。
A. 变量a比变量b占用的空间大
B. 变量a、b可以进行加减乘除运算
C. 变量a、b数据类型相同
D. 变量b可以用eval计算
13）已知s=‘显示“hello”’，则s 元素的个数是 A 。
A. 12
B. 9
C. 7
D. 18
14）运行字符串函数strncmp('s1','s2',2)，则结果为 B 。
A. 1
B. 0
C. true
D. fales
15）命令day（now）是指 C 。
A. 按日期字符串格式提取当前时间
B. 提取当前时间
C. 提取当前时间的日期
D. 按日期字符串格式提取当前日期
16）有一个2行2列的元胞数组c ，则c（2）是指 D 。
A. 第1行第2列元素内容
B. 第2行第1列元素内容
C. 第1行第2列元素 D .第2行第1列元素
17）以下运算中哪个运算级别最高 B 。
A. *
B. ^
C. ~=
D. /
18）运行命令bitand（20，15）的结果是 C 。
A. 15
B. 20
C. 4
D. 5
19）使用检测函数isinteger（15）的结果是 B 。
A. 1
B. 0
C. true
D. fales
20）计算三个多项式s1、s2和s3的乘积，则算式为 C 。
A. conv(s1,s2,s3)
B. s1*s2*s3
C. conv(conv(s1,s2),s3)
D. conv(s1*s2*s3)
以下写出MATLAB命令序列，并给出结果
2.复数向量a=2+3i，b=3-4i，计算a+b,a-b,c=a*b,d=a/b，并计算变量c的实部、虚部、模和相角。
3．用from:step:to的方式和linspace函数分别得到0~4π步长为0.4π的变量x1，0~4π分成10个点的变量x2。
4．输入矩阵a=
123
456
789
??
?
?
?
??
，使用全下标方式提取元素3，使用单下标方式提取元素8，取
出后两行子矩阵块，使用逻辑矩阵提取
13
79
?? ???
。
5．输入a 为3×3的魔方阵，b 为3×3的单位阵，将他们生成3×6的大矩阵c 、6×3的大矩阵d ，将d 的最后一行提取生成小矩阵e 。
6．矩阵a=
123
456
789
??
?
?
?
??
用flipud、fliplr、rot90、diag、triu和tril进行操作。并求其转置、秩、
逆矩阵、矩阵的行列式值及三次幂。
8．解线性方程组
1234
124
1234
1234 2328
36
87 73225
x x x x
x x x
x x x x
x x x x
-++=
?
?++=
?
?
-++=
?
?+-+=
?
。
9．输入字符串变量a 为‘hello ’，将其每个字符后移4个，如‘h ’变为‘l ’，然后再逆序存入变量b 。
10 计算函数2()10sin(4)t f t e t =-，其中t 范围为0到20，步长为0.2，g （t ）为f （t ）大于0的部分，计算g （t ）的值。
11．矩阵a=
123
456
789
??
?
?
?
??
，使用数组信息获取函数求其行列数、元素个数，是否为稀疏矩阵、
是否为字符型。
《数值计算方法》试题集及答案
《数值计算方法》复习试题 一、填空题： 1、????? ?????----=410141014A ，则A 的LU 分解为 A ??? ?????????=? ?????????? ?。 答案： ?? ????????--??????????--=1556141501 4115401411A 2、已知3.1)3(,2.1)2(,0.1)1(===f f f ，则用辛普生（辛卜生）公式计算求得 ?≈3 1 _________ )(dx x f ,用三点式求得≈')1(f 。 答案：， 3、1)3(,2)2(,1)1(==-=f f f ，则过这三点的二次插值多项式中2 x 的系数为 ， 拉格朗日插值多项式为 。 答案：-1， )2)(1(21 )3)(1(2)3)(2(21)(2--------= x x x x x x x L 4、近似值*0.231x =关于真值229.0=x 有( 2 )位有效数字； 5、设)(x f 可微,求方程)(x f x =的牛顿迭代格式是( )； （ 答案 )(1)(1n n n n n x f x f x x x '--- =+ 6、对1)(3 ++=x x x f ,差商=]3,2,1,0[f ( 1 ),=]4,3,2,1,0[f ( 0 )； 7、计算方法主要研究( 截断 )误差和( 舍入 )误差； 8、用二分法求非线性方程 f (x )=0在区间(a ,b )内的根时，二分n 次后的误差限为 ( 1 2+-n a b )； 9、求解一阶常微分方程初值问题y '= f (x ,y )，y (x 0)=y 0的改进的欧拉公式为
( )] ,(),([2111+++++=n n n n n n y x f y x f h y y )； 10、已知f (1)＝2，f (2)＝3，f (4)＝，则二次Newton 插值多项式中x 2系数为( )； 11、 两点式高斯型求积公式?1 d )(x x f ≈( ?++-≈1 )] 321 3()3213([21d )(f f x x f )，代数精 度为( 5 )； 12、 解线性方程组A x =b 的高斯顺序消元法满足的充要条件为(A 的各阶顺序主子式均 不为零)。 13、 为了使计算 32)1(6 )1(41310-- -+-+ =x x x y 的乘除法次数尽量地少，应将该表 达式改写为 11 ,))64(3(10-= -++=x t t t t y ，为了减少舍入误差，应将表达式 19992001-改写为 199920012 + 。 14、 用二分法求方程01)(3 =-+=x x x f 在区间[0,1]内的根,进行一步后根的所在区间 为 ，1 ,进行两步后根的所在区间为 ， 。 15、 、 16、 计算积分?1 5 .0d x x ,取4位有效数字。用梯形公式计算求得的近似值为 ，用辛卜 生公式计算求得的近似值为 ，梯形公式的代数精度为 1 ，辛卜生公式的代数精度为 3 。 17、 求解方程组?? ?=+=+042.01532121x x x x 的高斯—塞德尔迭代格式为 ?????-=-=+++20/3/)51()1(1)1(2)(2)1(1 k k k k x x x x ，该迭 代格式的迭代矩阵的谱半径)(M ρ= 121 。 18、 设46)2(,16)1(,0)0(===f f f ,则=)(1x l )2()(1--=x x x l ，)(x f 的二次牛顿 插值多项式为 )1(716)(2-+=x x x x N 。 19、 求积公式 ?∑=≈b a k n k k x f A x x f )(d )(0 的代数精度以( 高斯型 )求积公式为最高，具 有( 12+n )次代数精度。
《结构力学习题集》(下)-矩阵位移法习题及答案 (2)
第七章 矩阵位移法 一、是非题 1、单元刚度矩阵反映了该单元杆端位移与杆端力之间的关系。 2、单元刚度矩阵均具有对称性和奇异性。 3、局部坐标系与整体坐标系之间的坐标变换矩阵T 是正交矩阵。 4、结构刚度矩阵反映了结构结点位移与荷载之间的关系。 5、用 矩 阵 位 移 法 计 算 连 续 梁 时 无 需 对 单 元 刚 度 矩 阵 作 坐 标 变 换。 6、结 构 刚 度 矩 阵 是 对 称 矩 阵 ，即 有K i j = K j i ，这 可 由 位 移 互 等 定 理 得 到 证 明 。 7、结构刚度方程矩阵形式为：[]{}{}K P ?=，它是整个结构所应满足的变形条件。 ? 8、在直接刚度法的先处理法中，定位向量的物理意义是变形连续条件和位移边界条件。 9、等效结点荷载数值等于汇交于该结点所有固端力的代数和。 10、矩阵位移法中，等效结点荷载的“等效原则”是指与非结点荷载的结点位移相等。 11、矩阵位移法既能计算超静定结构，也能计算静定结构。 二、选择题 1、已知图示刚架各杆EI = 常数，当只考虑弯曲变形，且各杆单元类型相同时，采用先处理法进行结点位移编号，其正确编号是： (0,1,2) (0,0,0) (0,0,0) (0,1,3) (0,0,0)(1,2,0) (0,0,0)(0,0,3) (1,0,2) (0,0,0) (0,0,0)(1,0,3) (0,0,0) (0,1,2) (0,0,0)(0,3,4) A. B. C. D. 2134123412341234 2、平面杆件结构一般情况下的单元刚度矩阵[]k 66?，就其性质而言，是： A ．非对称、奇异矩阵； B ．对称、奇异矩阵； C ．对称、非奇异矩阵； D ．非对称、非奇异矩阵。 — 3、单元i j 在图示两种坐标系中的刚度矩阵相比：
计算方法的课后答案
《计算方法》习题答案 第一章 数值计算中的误差 1．什么是计算方法？（狭义解释） 答：计算方法就是将所求的的数学问题简化为一系列的算术运算和逻辑运算，以便在计算机上编程上机，求出问题的数值解，并对算法的收敛性、稳定性和误差进行分析、计算。 2．一个实际问题利用计算机解决所采取的五个步骤是什么？ 答：一个实际问题当利用计算机来解决时，应采取以下五个步骤： 实际问题→建立数学模型→构造数值算法→编程上机→获得近似结果 4．利用秦九韶算法计算多项式4)(5 3 -+-=x x x x P 在3-=x 处的值，并编程获得解。 解：400)(2 3 4 5 -+?+-?+=x x x x x x P ，从而 所以，多项式4)(5 3 -+-=x x x x P 在3-=x 处的值223)3(-=-P 。 5．叙述误差的种类及来源。 答：误差的种类及来源有如下四个方面： （1）模型误差：数学模型是对实际问题进行抽象，忽略一些次要因素简化得到的，它是原始问题的近似，即使数学模型能求出准确解，也与实际问题的真解不同，我们把数学模型与实际问题之间存在的误差称为模型误差。 （2）观测误差：在建模和具体运算过程中所用的一些原始数据往往都是通过观测、实验得来的，由于仪器的精密性，实验手段的局限性，周围环境的变化以及人们的工作态度和能力等因素，而使数据必然带有误差，这种误差称为观测误差。 （3）截断误差：理论上的精确值往往要求用无限次的运算才能得到，而实际运算时只能用有限次运算的结果来近似，这样引起的误差称为截断误差（或方法误差）。 （4）舍入误差：在数值计算过程中还会用到一些无穷小数，而计算机受机器字长的限制，它所能表示的数据只能是一定的有限数位，需要把数据按四舍五入成一定位数的近似的有理数来代替。这样引起的误差称为舍入误差。 6．掌握绝对误差（限）和相对误差（限）的定义公式。 答：设* x 是某个量的精确值，x 是其近似值，则称差x x e -=* 为近似值x 的绝对误差（简称误差）。若存在一个正数ε使ε≤-=x x e * ，称这个数ε为近似值x 的绝对误差限（简称误差限或精度）。 把绝对误差e 与精确值* x 之比* **x x x x e e r -==称为近似值x 的相对误差，称
《结构力学习题集》-矩阵位移法习题及答案
第八章 矩阵位移法 – 老八校 一、判断题： 1、单元刚度矩阵反映了该单元杆端位移与杆端力之间的关系。 2、单元刚度矩阵均具有对称性和奇异性。 3、局部坐标系与整体坐标系之间的坐标变换矩阵T 是正交矩阵。 4、结构刚度矩阵反映了结构结点位移与荷载之间的关系。 5、结构刚度方程矩阵形式为：[]{}{}K P ?=，它是整个结构所应满足的变形条件。 6、图示结构用矩阵位移法计算时（计轴向变形）未知量数目为8个。 7、在直接刚度法的先处理法中，定位向量的物理意义是变形连续条件和位移边界条件。 8、等效结点荷载数值等于汇交于该结点所有固端力的代数和。 9、矩阵位移法中，等效结点荷载的“等效原则”是指与非结点荷载的结点位移相等。 10、矩阵位移法既能计算超静定结构，也能计算静定结构。 11、已知图示刚架各杆EI = 常数，当只考虑弯曲变形，且各杆单元类型相同时，采用先处理法进行结点位移编号，其正确编号是： (0,1,2) (0,0,0) (0,0,0) (0,1,3) (0,0,0)(1,2,0) (0,0,0)(0,0,3) (1,0,2) (0,0,0) (0,0,0)(1,0,3) (0,0,0) (0,1,2) (0,0,0)(0,3,4) A. B. C. D. 2134123412341234 ( )
二、计算题： 12、用先处理法计算图示结构刚度矩阵的元素133322,,K K K 。 12 3l l 4 l 5EI 2EI EA (0,0,0) (0,0,1) (0,2,3) (0,0,0) (0,2,4)(0,0,0) EI 13、用先处理法计算图示刚架结构刚度矩阵的元素153422,,K K K 。EI ，EA 均为常数。 l 14、计算图示结构整体刚度矩阵的元素665544,,K K K 。E 为常数。 l l 1 3 4 2 A , I A A /222A I , 2A 15、写出图示结构以子矩阵形式表达的结构原始刚度矩阵的子矩阵 [][]K K 22 24 ,。 [][]k k 1112 [][] k k 2122 [] k = i i i i i 单刚分块形式为 ：
2016矩阵论复习题
矩阵论复习题 1. 设+=R V 是正实数集,对于任意的V y x ∈,,定义x 与y 的和为 y x y x ?=⊕ 对于任意的数R k ∈,定义k 与x 的数乘为 k x x k =? 问:对于上述定义加法和数乘运算的集合V ,是否构成线性空间,并说明理由. 2.对任意的2,R y x ∈,),(21x x x =,),(21y y y =定义x 与y 的和为 ),(112211y x y x y x y x +++=⊕ 对于任意的数R k ∈,定义k 与x 的数乘为 )2 )1(,(2121x k k kx kx x k -+=? 问:对于上述定义加法和数乘运算的集合2R ,是否构成线性空间,并说明理由. 3．设},022|),,{(321321R x x x x x x x S i ∈=++=，试证明S 是3R 的子空间，并求S 的一组基和S dim . 4.设)(R P n 表示次数不超过n 的全体多项式构成的线性空间, )}()(,0)0(|)({R P x f f x f S n ∈='= 证明S 是)(R P n 的子空间,并写出S 的一组基和计算S dim . 5. 设33:R R T →是线性变换， ()()321323213212,,2,,x x x x x x x x x x x T -++-+= 求T 的零空间)(T N 和像空间)(T R 的基和维数. 6. 设T 是3R 上的线性变换,对于基},,{k j i 有 k j k j i T -=++)( i k j T =+)( k j i k T 532)(++= 1)确定T 在基},,{k j i 下的矩阵; 2)求T 的像空间的基与维数.
计算方法练习题与答案
练习题与答案 练习题一 练习题二 练习题三 练习题四 练习题五 练习题六 练习题七 练习题八 练习题答案 练习题一 一、是非题 1.–作为x的近似值一定具有6位有效数字，且其误差限。() 2.对两个不同数的近似数，误差越小，有效数位越多。() 3.一个近似数的有效数位愈多，其相对误差限愈小。()
4.用近似表示cos x产生舍入误差。 ( ) 5.和作为的近似值有效数字位数相同。 ( ) 二、填空题 1.为了使计算的乘除法次数尽量少，应将该表达式改写 为； 2.–是x舍入得到的近似值，它有位有效数字，误差限 为，相对误差限为； 3.误差的来源是； 4.截断误差 为； 5.设计算法应遵循的原则 是。 三、选择题 1．–作为x的近似值，它的有效数字位数为( ) 。 (A) 7； (B) 3； (C) 不能确定 (D) 5. 2．舍入误差是( )产生的误差。 (A) 只取有限位数 (B) 模型准确值与用数值方法求得的准确值 (C) 观察与测量 (D) 数学模型准确值与实际值 3．用 1+x近似表示e x所产生的误差是( )误差。 (A). 模型 (B). 观测 (C). 截断 (D). 舍入 4．用s*=g t2表示自由落体运动距离与时间的关系式 (g为重力加速度)，s t是在时间t内的实际距离，则s t s*是（）误差。 (A). 舍入 (B). 观测 (C). 模型 (D). 截断 5．作为的近似值，有( )位有效数字。 (A) 3； (B) 4； (C) 5； (D) 6。
四、计算题 1．，，分别作为的近似值，各有几位有效数字？ 2．设计算球体积允许的相对误差限为1%，问测量球直径的相对误差限最大为多少？ 3．利用等价变换使下列表达式的计算结果比较精确： (1)， (2) (3) ， (4) 4．真空中自由落体运动距离s与时间t的关系式是s=g t2，g为重力加速度。现设g是精确的，而对t有秒的测量误差，证明：当t增加时，距离的绝对误差增加，而相对误差却减少。 5*. 采用迭代法计算，取 k=0,1,…, 若是的具有n位有效数字的近似值，求证是的具有2n位有效数字的近似值。 练习题二 一、是非题 1.单点割线法的收敛阶比双点割线法低。 ( ) 2.牛顿法是二阶收敛的。 ( ) 3.求方程在区间[1, 2]内根的迭代法总是收敛的。( ) 4.迭代法的敛散性与迭代初值的选取无关。 ( ) 5.求非线性方程f (x)=0根的方法均是单步法。 ( ) 二、填空题
《结构力学习题集》下矩阵位移法习题及答案 2
第七章 矩阵位移法 一、就是非题 1、单元刚度矩阵反映了该单元杆端位移与杆端力之间的关系。 2、单元刚度矩阵均具有对称性与奇异性。 3、局部坐标系与整体坐标系之间的坐标变换矩阵T 就是正交矩阵。 4、结构刚度矩阵反映了结构结点位移与荷载之间的关系。 5、用 矩 阵 位 移 法 计 算 连 续 梁 时 无 需 对 单 元 刚 度 矩 阵 作 坐 标 变 换。 6、结 构 刚 度 矩 阵 就是 对 称 矩 阵 ,即 有K i j = K j i ,这 可 由 位 移 互 等 定 理 得 到 证 明 。 7、结构刚度方程矩阵形式为:[]{}{}K P ?=,它就是整个结构所应满足的变形条件。 8、在直接刚度法的先处理法中,定位向量的物理意义就是变形连续条件与位移边界条件。 9、等效结点荷载数值等于汇交于该结点所有固端力的代数与。 10、矩阵位移法中,等效结点荷载的“等效原则”就是指与非结点荷载的结点位移相等。 11、矩阵位移法既能计算超静定结构,也能计算静定结构。 二、选择题 1、已知图示刚架各杆EI = 常数,当只考虑弯曲变形,且各杆单元类型相同时,采用先处理法进行结点位移编号,其正确编号就是: (0,1,2) (0,0,0) (0,0,0) (0,1,3) (0,0,0)(1,2,0) (0,0,0)(0,0,3) (1,0,2) (0,0,0) (0,0,0)(1,0,3) (0,0,0) (0,1,2) (0,0,0)(0,3,4) A. B. C. D. 2134123412341234 2、平面杆件结构一般情况下的单元刚度矩阵[]k 66?,就其性质而言,就是: A.非对称、奇异矩阵; B.对称、奇异矩阵; C.对称、非奇异矩阵; D.非对称、非奇异矩阵。 3、单元i j 在图示两种坐标系中的刚度矩阵相比: A.完全相同; B.第2、3、5、6行(列)等值异号;
上海交大研究生矩阵理论答案
习题 一 1．（1）因 cos sin sin cos nx nx nx nx ?? ? ? -?? cos sin sin cos x x x x ????-??= cos(1) sin(1)sin(1) cos(1)n x n x n x n x ++?? ??-++?? ，故由归纳法知 cos sin sin cos n nx nx A nx nx ?? =??-?? 。 （2）直接计算得4A E =-，故设4(0,1,2,3)n k r r =+=，则4(1)n k r k r A A A A ==-，即只需算出2 3 ,A A 即可。 （3）记J=0 1 0 1 1 0 ?????? ?????????? ，则 ， 1122111 11 () n n n n n n n n n n n n n n i i n i n n i n n n a C a C a C a C a C a A aE J C a J a C a a -----=-?????? ??=+==?? ???????? n ∑。 2．设112 2 (1,0),0 a A P P a A E λλ-??===???? 则由得 2 1112111 1 1 210 0 0 a λλλλλλλ?? ????==?????????????? 1时,不可能。 而由2 112 222 0 0 000 0 0 a λλλλλλ??????==?????????????? 1时,知1i λ=±所以所求矩阵为1 i PB P -， 其中P 为任意满秩矩阵，而 1231 0 1 0 1 0,,0 10 1 0 1B B B -?????? ===?????? --?????? 。 注：2A E =-无实解，n A E =的讨论雷同。 3．设A 为已给矩阵，由条件对任意n 阶方阵X 有AX=XA ，即把X 看作2 n 个未知数时线 性方程AX -XA=0有2 n 个线性无关的解，由线性方程组的理论知其系数矩阵为零矩
矩阵位移法练习题
结构力学自测题（第八单元） 矩阵位移法 姓名 学号 一、是 非 题（将 判 断 结 果 填 入 括 弧 ：以 O 表 示 正 确 ，以 X 表 示 错 误 ） 1、用 矩 阵 位 移 法 计 算 连 续 梁 时 无 需 对 单 元 刚 度 矩 阵 作 坐 标 变 换。 ( ) 2、结 构 刚 度 矩 阵 是 对 称 矩 阵 ，即 有 K ij = K ji ，这 可 由 位 移 互 等 定 理 得 到 证 明 。 () 3、图 示 梁 结 构 刚 度 矩 阵 的 元 素 K EI l 113 24=/ 。 ( ) EI l l EI 212 x y M , θ 附： ????? ?????????? ?????????? ???? ?--- -----l EI l EI l EI l EI l EI l EI l EI l EI l EA l EA l EI l EI l EI l EI l EI l EI l EI l EI l EA l EA 460260612061200000260460 6120612000002 22323222323 4、在 任 意 荷 载 作 用 下 ，刚 架 中 任 一 单 元 由 于 杆 端 位 移 所 引 起 的 杆 端 力 计 算 公 式 为 ：{} [][]{}F T K e e e =δ 。 ( ) 二、选 择 题 （ 将 选 中 答 案 的 字 母 填 入 括 弧 内 ） 1、已 知 图 示 刚 架 各杆 EI = 常 数，当 只 考 虑 弯 曲 变 形 ，且 各 杆 单 元 类 型 相 同 时 ，采 用 先 处 理 法 进 行 结 点 位 移 编 号 ，其 正 确 编 号 是 ： (0,1,2) (0,0,0) (0,0,0) (0,1,3) (0,0,0) (1,2,0) (0,0,0) (0,0,3) (1,0,2) (0,0,0) (0,0,0) (1,0,3) (0,0,0) (0,1,2) (0,0,0) (0,3,4) A. B. C. D. 2 1 3 4 1 2 3 4 1 2 3 4 1 2 3 4 x y M , θ ( ) 2、平 面 杆 件 结 构 一 般 情 况 下 的 单 元 刚 度 矩 阵 []k 66?， 就 其 性 质 而 言 ，是 ： ( ) A ．非 对 称 、奇 异 矩 阵 ； B ．对 称 、奇 异 矩 阵 ； C ．对 称 、非 奇 异 矩 阵 ； D ．非 对 称 、非 奇 异 矩 阵 。 3、单 元 i j 在 图 示 两 种 坐 标 系 中 的 刚 度 矩 阵 相 比 ： A ． 完 全 相 同 ； B ． 第 2、3、5、6 行 (列 ) 等 值 异 号 ； C ． 第 2、5 行 (列 )等 值 异 号 ； D ． 第 3、6 行 (列 ) 等 值 异 号 。 ( ) i j y x i j y x M , θ M , θ 4、矩 阵 位 移 法 中 ，结 构 的 原 始 刚 度 方 程 是 表 示 下 列 两 组 量 值 之 间 的 相 互 关 系 ： ( ) A ．杆 端 力 与 结 点 位 移 ； B ．杆 端 力 与 结 点 力 ； C ．结 点 力 与 结 点 位 移 ； D ．结 点 位 移 与 杆 端 力 。 5、单 元 刚 度 矩 阵 中 元 素 k ij 的 物 理 意 义 是 ： A ．当 且 仅 当 δi =1 时 引 起 的 与 δj 相 应 的 杆 端 力 ； B ．当 且 仅 当 δj =1时 引 起 的 与 δi 相 应 的 杆 端 力 ； C ．当 δj =1时 引 起 的 δi 相 应 的 杆 端 力 ； D ．当 δi =1时 引 起 的 与 δj 相 应 的 杆 端 力。 () 6、用 矩 阵 位 移 法 解 图 示 连 续 梁 时 ，结 点 3 的 综 合 结 点 荷 载 是 ： A ．[]-ql ql 2 12 T 132 ； B ．[]ql ql 2132 12T -； C ．[]--ql ql 2112 12T ； D ．[]ql ql 2112 12T 。 ( ) 123 l /2 l l ql 2 q 4 ql l /2 x y M , θ 7、用 矩 阵 位 移 法 解 图 示 结 构 时 ，已 求 得 1 端 由 杆 端 位 移 引 起 的 杆 端 力 为 {}[] T F 461--=，则 结 点 1 处 的 竖 向 反 力 Y 1 等 于 ： A ．6-； B ．-10； C ．10 ； D ．14 。 ( ) 2m 4m 12 3 M 1 Y 20kN/m 1 x y M , θ 三、填 充 题 （ 将 答 案 写 在 空 格 内） 1、图 示 桁 架 结 构 刚 度 矩 阵 有 个 元 素 ，其 数 值 等 于 。 2m 3m 3m A B C D EA EA EA x y M , θ 2、图 示 刚 架 用 两 种 方 式 进 行 结 点 编 号 ，结 构 刚 度 矩 阵 最 大 带 宽 较 小 的 是 图 。 3 5 641 2 7 1 2345 6 7 (a) (b) 3、图 示 梁 结 构 刚 度 矩 阵 的 主 元 素 K K 1122== , 。 l l 2EI EI 1 2 x y M , θ 四、图 a 、b 所 示 两 结 构 ，各 杆 EI 、l 相 同 ，不 计 轴 向 变 形 ， 已 求 得 图 b 所 示 结 构 的 结 点 位 移 列 阵 为 {}?=-???? ? ?ql EI ql REI ql EI 34396192192 T 。试 求 图 a 所 示 结 构 中 单 元 ① 的 杆 端 力 列 阵。 q 1 2 3 4(a) ql 2 ② ③ ① 1 2 34 (b) ② ③ ① x y M , θ 五、图 a 所 示 结 构 (整 体 坐 标 见 图 b )，图 中 圆 括 号 内 数 码 为 结 点 定 位 向 量 (力 和 位 移 均 按 水 平 、竖 直 、转 动
研究生矩阵论课后习题答案(全)习题二
习题二 1．化下列矩阵为Smith 标准型： （1）222211λλλλ λλλλλ?? -?? -????+-?? ; （2）2222 00 000 00(1)00000λλλλλλ ?? ?? -? ? ??-?? -?? ; （3）2222 232321234353234421λλλλλλλλλλλλλλ?? +--+-??+--+-????+---?? ; (4)23014360220620101003312200λλλλλλλλλλλλλλ????++??????--????---?? . 解：（1）对矩阵作初等变换 23221311(1)100 10 000000(1)00(1)c c c c c c r λλλλλλλλλ+--?-???????????→-???→? ??? ????-++???? ， 则该矩阵为Smith 标准型为 ???? ? ?????+)1(1λλλ； （2）矩阵的各阶行列式因子为 44224321()(1),()(1),()(1),()1D D D D λλλλλλλλλλ=-=-=-=, 从而不变因子为 22 2341234123()()() ()1,()(1),()(1),()(1)()()() D D D d d d d D D D λλλλλλλλλλλλλλλλ== =-==-==-故该矩阵的Smith 标准型为
2210000(1)0000(1)00 00(1)λλλλλλ?? ??-????-?? -??； （3）对矩阵作初等变换 故该矩阵的Smith 标准型为 ?? ?? ??????+--)1()1(112 λλλ; (4)对矩阵作初等变换 在最后的形式中，可求得行列式因子 3254321()(1),()(1),()()()1D D D D D λλλλλλλλλ=-=-===, 于是不变因子为 2541234534()() ()()()1,()(1),()(1)()() D D d d d d d D D λλλλλλλλλλλλλ==== =-==-故该矩阵的Smith 标准形为 2 1 0000 010 0000100000(1)00 00 0(1)λλλλ?????????? -?? ??-?? . 2.求下列λ-矩阵的不变因子： （1） 21 0021002λλλ--????--????-??； （2）100 1000 λαββλα λαββ λα+????-+? ???+??-+?? ；
《数值计算方法》试题集及答案
《数值计算方法》复习试题 一、填空题: 1、????? ?????----=410141014A ，则A 的ＬU 分解为 A ??? ?????????=? ?????????? ?。 答案： ?? ????????--??????????--=1556141501 4115401411A 3、1)3(,2)2(,1)1(==-=f f f ,则过这三点的二次插值多项式中2 x 的系数 为 ，拉格朗日插值多项式为 。 答案:-1， )2)(1(21 )3)(1(2)3)(2(21)(2--------= x x x x x x x L 4、近似值*0.231x =关于真值229.0=x 有( 2 )位有效数字； ５、设)(x f 可微,求方程)(x f x =的牛顿迭代格式是( ）； 答案 )(1)(1n n n n n x f x f x x x '--- =+ ６、对1)(3 ++=x x x f ,差商=]3,2,1,0[f ( 1 )，=]4,3,2,1,0[f ( ０ )； 7、计算方法主要研究（ 截断 ）误差和( 舍入 )误差; 8、用二分法求非线性方程 f (x )=0在区间(a ,b )内的根时,二分n 次后的误差限为 （ 1 2+-n a b ); １０、已知f （１）＝2,f (2)=３,ｆ(４)＝5.９，则二次Ne ｗton 插值多项式中x 2系数为 ( ０.15 )； 11、 解线性方程组A x =b 的高斯顺序消元法满足的充要条件为（A 的各阶顺序主子式均 不为零)。 12、 为了使计算 32)1(6 )1(41310-- -+-+ =x x x y 的乘除法次数尽量地少，应将该
《结构力学习题集》(下)-矩阵位移法习题及答案
第八章 矩阵位移法 1、(O) 2、(X) 3、(O) 4、(X) 5、(X) 6、(O) 7、(O) 8、(X) 9、(O) 10、(O) 11、(A) 一、判断题： 1、单元刚度矩阵反映了该单元杆端位移与杆端力之间的关系。 2、单元刚度矩阵均具有对称性和奇异性。 3、局部坐标系与整体坐标系之间的坐标变换矩阵T 是正交矩阵。 4、结构刚度矩阵反映了结构结点位移与荷载之间的关系。 5、结构刚度方程矩阵形式为：[]{}{}K P ?=，它是整个结构所应满足的变形条件。 6、图示结构用矩阵位移法计算时（计轴向变形）未知量数目为8个。 7、在直接刚度法的先处理法中，定位向量的物理意义是变形连续条件和位移边界条件。 8、等效结点荷载数值等于汇交于该结点所有固端力的代数和。 9、矩阵位移法中，等效结点荷载的“等效原则”是指与非结点荷载的结点位移相等。 10、矩阵位移法既能计算超静定结构，也能计算静定结构。 11、已知图示刚架各杆EI = 常数，当只考虑弯曲变形，且各杆单元类型相同时，采用先处理法进行结点位移编号，其正确编号是： (0,1,2) (0,0,0) (0,0,0) (0,1,3) (0,0,0)(1,2,0) (0,0,0)(0,0,3) (1,0,2) (0,0,0) (0,0,0)(1,0,3) (0,0,0) (0,1,2) (0,0,0)(0,3,4) A. B. C. D. 2134123412341234x y M , θ( )
二、计算题： 12、用先处理法计算图示结构刚度矩阵的元素133322,,K K K 。 123l l 4l l 5EI 2EI EA (0,0,0) (0,0,1) (0,2,3) (0,0,0) (0,2,4)(0,0,0) x y M , θ EI 13、用先处理法计算图示刚架结构刚度矩阵的元素153422,,K K K 。EI ，EA 均为常数。 l (0,0,1) (0,5,0) (2,3,4) l ① ② 123x y M , θ 14、计算图示结构整体刚度矩阵的元素665544,,K K K 。E 为常数。 l l l 1 3 4 2A , I A A /222A I , 2A x y M , θ 15、写出图示结构以子矩阵形式表达的结构原始刚度矩阵的子矩阵 [][]K K 22 24 ,。 3 12① ② ③ [][]k k 1112 [][] k k 2122 [] k = i i i i i 单刚分块形式为 ： 4x y M , θ
研究生矩阵论课后习题答案(全)习题三
习题三 1．证明下列问题： （1）若矩阵序列{}m A 收敛于A ，则{}T m A 收敛于T A ，{} m A 收敛于A ； （2）若方阵级数∑∞ =0m m m A c 收敛，则∑∑∞ =∞==?? ? ??00)(m m T m T m m m A c A c . 证明：（1）设矩阵 ,,2,1,)() (Λ==?m a A n n m ij m 则 ,)()(n n m ji T m a A ?=,)()(n n m ij m a A ?=,,2,1Λ=m 设 ,)(n n ij a A ?= 则 n n ji T a A ?=)(，,)(n n ij a A ?= 若矩阵序列{}m A 收敛于A ，即对任意的n j i ,,2,1,Λ=，有 ij m ij m a a =∞ →) (lim ， 则 ji m ji m a a =∞ →)(lim ，ij m ij m a a =∞ →)(lim ，n j i ,,2,1,Λ=， 故{} T m A 收敛于T A ，{} m A 收敛于A . (2)设方阵级数 ∑∞ =0 m m m A c 的部分和序列为 ΛΛ,,,,21m S S S ， 其中m m m A c A c c S +++=Λ10.
若 ∑∞ =0 m m m A c 收敛，设其和为S ，即 S A c m m m =∑∞ =0 ，或S S m m =∞ →lim ， 则 T T m m S S =∞ →lim . 而级数∑∞ =0 )(m m T m A c 的部分和即为T m S ，故级数∑∞ =0 )(m m T m A c 收敛，且其和为T S ， 即 ∑∑∞ =∞==?? ? ??00)(m m T m T m m m A c A c . 2.已知方阵序列{}m A 收敛于A ，且{} 1-m A ，1 -A 都存在，证明： （1）A A m m =∞ →lim ；（2）{}1 1 lim --∞ →=A A m m . 证明：设矩阵 ,,2,1,)() (Λ==?m a A n n m ij m ,)(n n ij a A ?= 若矩阵序列{}m A 收敛于A ，即对任意的n j i ,,2,1,Λ=，有 ij m ij m a a =∞ →) (lim . (1) 由于对任意的n j j j ,,,21Λ，有 ,lim ) (k k kj m kj m a a =∞ → n k ,,2,1Λ=， 故 ∑-∞ →n n n j j j m nj m j m j j j j m a a a ΛΛΛ2121)()(2)(1) ()1(lim τ ＝ ∑-n n n j j j nj j j j j j a a a ΛΛΛ21212121) ()1(τ ， 而 ∑-= n n n j j j m nj m j m j j j j m a a a A ΛΛΛ2121) ()(2)(1)()1(τ，
第9章 矩阵位移法 例题
第9章 矩阵位移法 习 题 9-1：请给图示结构编号（同时用先处理法和后处理法）及建立坐标。 题9-1图 9-2：求图示连续梁的整体刚度矩阵。 题9-2图 9-3：求图示刚架的整体刚度矩阵。 （c ） （e ）
题9-3图 9-4：求图示组合结构的整体刚度矩阵。 题9-4图 9-5：求图示桁架结构的整体刚度矩阵，所有杆件的EA 均相同。 题9-5图 9-6：求图示排架结构的整体刚度矩阵。 题9-6图 9-7：求图示结构的等效结点荷载，请利用结构的对称性。 1kN/m
题9-7图 9-8：求图示结构的等效结点荷载，请利用结构的对称性。 题9-8图 9-9：求图示结构的等效结点荷载。 题9-9图 9-10：求出图示结构的荷载列阵。 题9-10图 9-11：求出图示结构的荷载列阵，请分别用先处理法和后处理法进行编号。 q q
题9-11图 9-12：求图示结构的荷载列阵，考虑轴向变形。 题9-12图 9-13：求图示结构的荷载列阵。 题9-13图 9-14：图示连续梁中间支座发生了下向的移动a ，请求出其整体刚度方程。 题9-14图 10kN/m q
9-15：请求出图示连续梁的整体刚度方程。 题9-15图 9-16：求图示连续梁的整体刚度矩阵。 题9-16图 9-17：图示结构温度发生了变化，请求出整体刚度方程。杆件的EI 、EA 相同。 题9-17图 9-18：图示结构温度发生了变化，请求出整体刚度方程。 题9-18图 9-19：图示结构发生了支座移动，请画出结构的内力图。 00
数值计算方法习题答案(第二版)(绪论)
数值分析 （p11页） 4 试证：对任给初值x 0, 0)a >的牛顿迭代公式 112(),0,1 ,2,......k a k k x x x k +=+= 恒成立下列关系式： 2112(1)(,0,1,2,.... (2)1,2,...... k k k x k x x k x k +-=≥= 证明： （1 ）(2 1122k k k k k k x a x x x x +-??=+= =? ?? （2） 取初值00>x ，显然有0>k x ，对任意0≥k ， a a x a x x a x x k k k k k ≥+??? ? ??-=???? ??+=+2 12121 6 证明： 若k x 有n 位有效数字，则n k x -?≤ -1102 1 8， 而() k k k k k x x x x x 28882182 1-=-???? ? ?+=-+ n n k k x x 21221102 1 5.22104185 .28--+?=??<-∴>≥ 1k x +∴必有2n 位有效数字。 8 解： 此题的相对误差限通常有两种解法. ①根据本章中所给出的定理： （设x 的近似数* x 可表示为m n a a a x 10......021*?±=，如果* x 具有l 位有效数字，则其相对误差限为 ()11 * *1021 --?≤ -l a x x x ，其中1a 为*x 中第一个非零数）
则7.21=x ，有两位有效数字，相对误差限为 025.0102 21 111=??≤--x x e 71.22=x ，有两位有效数字，相对误差限为 025.0102 21 122=??≤--x x e 3 2.718x =，有两位有效数字，其相对误差限为： 00025.0102 21 333=??≤--x e x ②第二种方法直接根据相对误差限的定义式求解 对于7.21=x ，0183.01<-e x ∴其相对误差限为00678.07 .20183.011≈<-x e x 同理对于71.22=x ，有 003063 .071 .20083 .022≈<-x e x 对于718.23=x ，有 00012.0718 .20003 .033≈<-x e x 备注：（1）两种方法均可得出相对误差限，但第一种是对于所有具有n 位有效数字的近似数都成立的正确结论，故他对误差限的估计偏大，但计算略简单些；而第二种方法给出较好的误差限估计，但计算稍复杂。 （2）采用第二种方法时，分子为绝对误差限，不是单纯的对真实值与近似值差值的四舍五入，绝对误差限大于或等于真实值与近似值的差。 11. 解: ......142857.3722≈，.......1415929.3113 255≈ 2102 1 722-?≤-∴ π，具有3位有效数字
【免费下载】控制中的矩阵理论习题
练习一： 1．设A 、是Hermite 矩阵，证明：AB 是Hermite 矩阵的充分必要条件是n n C B ?∈AB=BA 。2．设，若，则A 为反Hermite 矩阵。试证明：任意一个都n n C A ?∈A A H -=n n C B ?∈可以唯一地表示为一个Hermitet 矩阵与一个反Hermite 矩阵的和。3．证明反Hermite 矩阵的主对角线上的元素或为零，或为纯虚数。4．设是Hermite 矩阵，rank(A)=1，证明：矩阵A 的主对角线上凡不是零的元素n n C A ?∈都是具有同符号的实数；又设是反Hermite 矩阵，rank(B)=1，证明：矩阵B n n C B ?∈的主对角线上凡不是零的元素都是具有同符号的虚部之纯虚数。5．试求一酉矩阵P ，使为对角矩阵，这里AP P AP P H =-1（1）A=； （2）A=。??????????----10001i i i i ??????????-0010010i i 6. 设是Hermite 矩阵。证明A 是Hernite 正定矩阵的充分必要条件是，存在n n C A ?∈Hermite 正定矩阵B ，使得。2 B A =7．设是Hermite 矩阵,则下列条件等价：n n C A ?∈ （1）A 是Hernite 半正定矩阵； （2）A 的特征值全为非负实数； （3）存在矩阵，使得。n n C P ?∈P P A H =练习二：1．用初等变换化下列多项式矩阵为Smith 标准形：（1） ； （2）；()???? ??+-=λλλλλλλ352223A ()??????????-+--=222211λλλλλλλλλλB （3） ；（4）()()220000 001C λλλλλ??+??=????+????。()()??????????????---=00000100000002222λλλλλλλD 2．求下多项式矩阵的不变因子：
计算方法模拟试题及答案
计算方法模拟试题 一、 单项选择题（每小题3分，共15分） 1．近似值210450.0?的误差限为（ ）。 A ． 0.5 B. 0.05 C ． 0.005 D. 0.0005. 2. 求积公式)2(3 1 )1(34)0(31)(2 0f f f dx x f ++≈ ?的代数精确度为（ ）。 A. 1 B. 2 C. 3 D. 4 3. 若实方阵A 满足（ ）时，则存在唯一单位下三角阵L 和上三角阵R ，使LR A =。 A. 0det ≠A B. 某个0 det ≠k A C. )1,1(0det -=≠n k A k D. ),,1(0det n k A k =≠ 4.已知?? ?? ? ?????=531221112A ，则=∞A （ ）。 A. 4 B. 5 C. 6 D 9 5.当实方阵A 满足)2(,221>>-=i i λλλλ，则乘幂法计算公式1e =（ ）。 A. 1+k x B. k k x x 11λ++ C. k x D. k k x x 11λ-+ 二、填空题（每小题3分，共15分） 1. 14159.3=π，具有4位有效数字的近似值为 。 2. 已知近似值21,x x ，则=-?)(21x x 。 3．已知1)(2-=x x f ，则差商=]3,2,1[f 。 4．雅可比法是求实对称阵 的一种变换方法。
5．改进欧拉法的公式为 。 三、计算题（每小题12分 ，共60分） 1. 求矛盾方程组； ??? ??=-=+=+2 42321 2121x x x x x x 的最小二乘解。 2．用列主元法解方程组 ??? ??=++=++=++4 26453426352321 321321x x x x x x x x x 3．已知方程组 ???? ? ?????=????????????????????----131********x x x a a a a （1） 写出雅可比法迭代公式； （2） 证明2