《高等工程数学》习题卷(一)
第1大题(共8小题,每小题5分,共40分)
1.1 给定线性空间 V = s p a n { ( − 1 , 1 , 0 ) , ( 0 , 1 , 1 ) , ( 1 , 1 , 1 ) } V=\mathrm{span}\{(-1,1,0),(0,1,1),(1,1,1)\} V=span{(−1,1,0),(0,1,1),(1,1,1)},向量 v = ( 1 , 2 , 3 ) \mathbf{v}=(1,2,3) v=(1,2,3)在该基下的坐标为多少?
解答过程:
设
v
=
a
(
−
1
,
1
,
0
)
+
b
(
0
,
1
,
1
)
+
c
(
1
,
1
,
1
)
\mathbf{v}=a(-1,1,0)+b(0,1,1)+c(1,1,1)
v=a(−1,1,0)+b(0,1,1)+c(1,1,1),则有
\begin{align*}
a-b+c&=1\
b+c&=2\
c&=3
\end{align*}
解得
a
=
−
2
,
b
=
−
1
,
c
=
3
a=-2,b=-1,c=3
a=−2,b=−1,c=3,因此
v
\mathbf{v}
v在该基下的坐标为
(
−
2
,
−
1
,
3
)
(-2,-1,3)
(−2,−1,3)。
答案: ( − 2 , − 1 , 3 ) (-2,-1,3) (−2,−1,3)。
1.2 给定欧氏空间 R 3 \mathbb{R}^3 R3的标准正交基 { e 1 , e 2 , e 3 } \{\mathbf{e}_1,\mathbf{e}_2,\mathbf{e}_3\} {e1,e2,e3},向量 v = ( 1 , 2 , 2 ) \mathbf{v}=(1,2,2) v=(1,2,2)的长度为多少?
解答过程:
v \mathbf{v} v的长度为 ∥ v ∥ = v ⋅ v = 1 2 + 2 2 + 2 2 = 9 = 3 \|\mathbf{v}\|=\sqrt{\mathbf{v}\cdot\mathbf{v}}=\sqrt{1^2+2^2+2^2}=\sqrt{9}=3 ∥v∥=v⋅v=12+22+22=9=3。
答案: 3 3 3。
1.3 给定矩阵 A = ( 1 2 3 4 5 6 7 8 9 ) A=\begin{pmatrix}1&2&3\\4&5&6\\7&8&9\end{pmatrix} A=⎝⎛147258369⎠⎞,求其2-范数、无穷范数、1-范数。
解答过程:
A A A的2-范数为 ∥ A ∥ 2 = λ max ( A T A ) = λ max ( 66 78 90 78 93 108 90 108 126 ) = 294 \|A\|_2=\sqrt{\lambda_{\max}(A^TA)}=\sqrt{\lambda_{\max}\begin{pmatrix}66&78&90\\78&93&108\\90&108&126\end{pmatrix}}=\sqrt{294} ∥A∥2=λmax(ATA)=λmax⎝⎛667890789310890108126⎠⎞=294。
A A A的无穷范数为 ∥ A ∥ ∞ = max 1 ≤ i ≤ 3 ∑ j = 1 3 ∣ a i j ∣ = max { 6 , 15 , 24 } = 24 \|A\|_{\infty}=\max_{1\leq i\leq 3}\sum_{j=1}^3|a_{ij}|=\max\{6,15,24\}=24 ∥A∥∞=max1≤i≤3∑j=13∣aij∣=max{6,15,24}=24。
A A A的1-范数为 ∥ A ∥ 1 = max 1 ≤ j ≤ 3 ∑ i = 1 3 ∣ a i j ∣ = max { 12 , 15 , 18 } = 18 \|A\|_1=\max_{1\leq j\leq 3}\sum_{i=1}^3|a_{ij}|=\max\{12,15,18\}=18 ∥A∥1=max1≤j≤3∑i=13∣aij∣=max{12,15,18}=18。
答案: 294 , 24 , 18 \sqrt{294},24,18 294,24,18。
1.4 确定方阵 A A A幂级数 ∑ k = 0 ∞ A k \sum_{k=0}^{\infty}A^k ∑k=0∞Ak收敛的条件。
解答过程:
当 ρ ( A ) < 1 \rho(A)<1 ρ(A)<1时,幂级数 ∑ k = 0 ∞ A k \sum_{k=0}^{\infty}A^k ∑k=0∞Ak收敛,其中 ρ ( A ) \rho(A) ρ(A)为 A A A的谱半径。
答案: ρ ( A ) < 1 \rho(A)<1 ρ(A)<1。
1.5 给定样本 { x 1 , x 2 , … , x n } \{x_1,x_2,\ldots,x_n\} {x1,x2,…,xn},计算其平均数、众数、中位数。
解答过程:
样本的平均数为 x ˉ = 1 n ∑ i = 1 n x i \bar{x}=\frac{1}{n}\sum_{i=1}^nx_i xˉ=n1∑i=1nxi。
样本的众数为出现次数最多的数。
样本的中位数为将样本从小到大排序后,位于中间位置的数,若样本个数为偶数,则中位数为中间两个数的平均数。
答案:平均数为 x ˉ \bar{x} xˉ,众数为出现次数最多的数,中位数为样本的中位数。
1.6 给定来自正态总体 N ( μ , σ 2 ) N(\mu,\sigma^2) N(μ,σ2)的样本 { x 1 , x 2 , … , x n } \{x_1,x_2,\ldots,x_n\} {x1,x2,…,xn},确定统计量 x ˉ − μ s / n \frac{\bar{x}-\mu}{s/\sqrt{n}} s/nxˉ−μ的分布以及常数。
解答过程:
由中心极限定理可知,当样本容量 n n n足够大时, x ˉ − μ s / n \frac{\bar{x}-\mu}{s/\sqrt{n}} s/nxˉ−μ的分布近似于标准正态分布 N ( 0 , 1 ) N(0,1) N(0,1)。
其中 s s s为样本标准差,常数为 n \sqrt{n} n。
答案: N ( 0 , 1 ) N(0,1) N(0,1),常数为 n \sqrt{n} n。
1.7 给定样本 { x 1 , x 2 , … , x n } \{x_1,x_2,\ldots,x_n\} {x1,x2,…,xn},其概率密度函数为 f ( x ; θ ) = 1 θ e − x θ f(x;\theta)=\frac{1}{\theta}\mathrm{e}^{-\frac{x}{\theta}} f(x;θ)=θ1e−θx,其中 θ > 0 \theta>0 θ>0为未知参数。设 θ ^ \hat{\theta} θ^为 θ \theta θ的矩估计,计算 θ ^ \hat{\theta} θ^的值。
解答过程:
θ ^ \hat{\theta} θ^为 θ \theta θ的一阶矩,即 θ ^ = 1 n ∑ i = 1 n x i \hat{\theta}=\frac{1}{n}\sum_{i=1}^nx_i θ^=n1∑i=1nxi。
答案: θ ^ = 1 n ∑ i = 1 n x i \hat{\theta}=\frac{1}{n}\sum_{i=1}^nx_i θ^=n1∑i=1nxi。
1.8 给定正态总体 N ( μ , σ 2 ) N(\mu,\sigma^2) N(μ,σ2)的样本 { x 1 , x 2 , … , x n } \{x_1,x_2,\ldots,x_n\} {x1,x2,…,xn},求 μ \mu μ的置信区间,置信水平为 1 − α 1-\alpha 1−α。
解答过程:
由于
x
ˉ
−
μ
s
/
n
\frac{\bar{x}-\mu}{s/\sqrt{n}}
s/nxˉ−μ的分布近似于标准正态分布
N
(
0
,
1
)
N(0,1)
N(0,1),因此有
P
(
−
z
α
/
2
<
x
ˉ
−
μ
s
/
n
<
z
α
/
2
)
=
1
−
α
P\left(-z_{\alpha/2}<\frac{\bar{x}-\mu}{s/\sqrt{n}}<z_{\alpha/2}\right)=1-\alpha
P(−zα/2<s/nxˉ−μ<zα/2)=1−α
其中
z
α
/
2
z_{\alpha/2}
zα/2为标准正态分布的上
α
/
2
\alpha/2
α/2分位数。
移项得
P
(
x
ˉ
−
z
α
/
2
s
n
<
μ
<
x
ˉ
+
z
α
/
2
s
n
)
=
1
−
α
P\left(\bar{x}-\frac{z_{\alpha/2}s}{\sqrt{n}}<\mu<\bar{x}+\frac{z_{\alpha/2}s}{\sqrt{n}}\right)=1-\alpha
P(xˉ−nzα/2s<μ<xˉ+nzα/2s)=1−α
因此
μ
\mu
μ的置信区间为
(
x
ˉ
−
z
α
/
2
s
n
,
x
ˉ
+
z
α
/
2
s
n
)
\left(\bar{x}-\frac{z_{\alpha/2}s}{\sqrt{n}},\bar{x}+\frac{z_{\alpha/2}s}{\sqrt{n}}\right)
(xˉ−nzα/2s,xˉ+nzα/2s)。
答案: ( x ˉ − z α / 2 s n , x ˉ + z α / 2 s n ) \left(\bar{x}-\frac{z_{\alpha/2}s}{\sqrt{n}},\bar{x}+\frac{z_{\alpha/2}s}{\sqrt{n}}\right) (xˉ−nzα/2s,xˉ+nzα/2s)。
第2大题(共2小题,每小题5分,共10分)
2.1 给定线性变换
T
:
R
3
→
R
3
T:\mathbb{R}^3\rightarrow\mathbb{R}^3
T:R3→R3,使得
T
(
1
0
0
)
=
(
1
1
0
)
T\begin{pmatrix}1\\0\\0\end{pmatrix}=\begin{pmatrix}1\\1\\0\end{pmatrix}
T⎝⎛100⎠⎞=⎝⎛110⎠⎞,
T
(
0
1
0
)
=
(
0
1
1
)
T\begin{pmatrix}0\\1\\0\end{pmatrix}=\begin{pmatrix}0\\1\\1\end{pmatrix}
T⎝⎛010⎠⎞=⎝⎛011⎠⎞,
T
(
0
0
1
)
=
(
1
1
1
)
T\begin{pmatrix}0\\0\\1\end{pmatrix}=\begin{pmatrix}1\\1\\1\end{pmatrix}
T⎝⎛001⎠⎞=⎝⎛111⎠⎞。求
(
2
−
1
3
)
\begin{pmatrix}2\\-1\\3\end{pmatrix}
⎝⎛2−13⎠⎞在
T
T
T下的像。
解:设 ( 2 − 1 3 ) \begin{pmatrix}2\\-1\\3\end{pmatrix} ⎝⎛2−13⎠⎞在基 { ( 1 0 0 ) , ( 0 1 0 ) , ( 0 0 1 ) } \{\begin{pmatrix}1\\0\\0\end{pmatrix},\begin{pmatrix}0\\1\\0\end{pmatrix},\begin{pmatrix}0\\0\\1\end{pmatrix}\} {⎝⎛100⎠⎞,⎝⎛010⎠⎞,⎝⎛001⎠⎞}下的坐标为 ( a b c ) \begin{pmatrix}a\\b\\c\end{pmatrix} ⎝⎛abc⎠⎞,则有
( 2 − 1 3 ) = a ( 1 0 0 ) + b ( 0 1 0 ) + c ( 0 0 1 ) \begin{pmatrix}2\\-1\\3\end{pmatrix}=a\begin{pmatrix}1\\0\\0\end{pmatrix}+b\begin{pmatrix}0\\1\\0\end{pmatrix}+c\begin{pmatrix}0\\0\\1\end{pmatrix} ⎝⎛2−13⎠⎞=a⎝⎛100⎠⎞+b⎝⎛010⎠⎞+c⎝⎛001⎠⎞
解得 a = 2 , b = − 1 , c = 3 a=2,b=-1,c=3 a=2,b=−1,c=3,因此 ( 2 − 1 3 ) \begin{pmatrix}2\\-1\\3\end{pmatrix} ⎝⎛2−13⎠⎞在基 { ( 1 0 0 ) , ( 0 1 0 ) , ( 0 0 1 ) } \{\begin{pmatrix}1\\0\\0\end{pmatrix},\begin{pmatrix}0\\1\\0\end{pmatrix},\begin{pmatrix}0\\0\\1\end{pmatrix}\} {⎝⎛100⎠⎞,⎝⎛010⎠⎞,⎝⎛001⎠⎞}下的坐标为 ( 2 − 1 3 ) \begin{pmatrix}2\\-1\\3\end{pmatrix} ⎝⎛2−13⎠⎞。
由于 ( 1 0 0 ) \begin{pmatrix}1\\0\\0\end{pmatrix} ⎝⎛100⎠⎞在 T T T下的像为 ( 1 1 0 ) \begin{pmatrix}1\\1\\0\end{pmatrix} ⎝⎛110⎠⎞, ( 0 1 0 ) \begin{pmatrix}0\\1\\0\end{pmatrix} ⎝⎛010⎠⎞在 T T T下的像为 ( 0 1 1 ) \begin{pmatrix}0\\1\\1\end{pmatrix} ⎝⎛011⎠⎞, ( 0 0 1 ) \begin{pmatrix}0\\0\\1\end{pmatrix} ⎝⎛001⎠⎞在 T T T下的像为 ( 1 1 1 ) \begin{pmatrix}1\\1\\1\end{pmatrix} ⎝⎛111⎠⎞,因此 T T T在基 { ( 1 0 0 ) , ( 0 1 0 ) , ( 0 0 1 ) } \{\begin{pmatrix}1\\0\\0\end{pmatrix},\begin{pmatrix}0\\1\\0\end{pmatrix},\begin{pmatrix}0\\0\\1\end{pmatrix}\} {⎝⎛100⎠⎞,⎝⎛010⎠⎞,⎝⎛001⎠⎞}下的矩阵为
[ T ] = ( 1 0 1 1 1 1 0 1 1 ) [T]=\begin{pmatrix}1&0&1\\1&1&1\\0&1&1\end{pmatrix} [T]=⎝⎛110011111⎠⎞
因此, ( 2 − 1 3 ) \begin{pmatrix}2\\-1\\3\end{pmatrix} ⎝⎛2−13⎠⎞在 T T T下的像为
( 1 0 1 1 1 1 0 1 1 ) ( 2 − 1 3 ) = ( 5 4 2 ) \begin{pmatrix}1&0&1\\1&1&1\\0&1&1\end{pmatrix}\begin{pmatrix}2\\-1\\3\end{pmatrix}=\begin{pmatrix}5\\4\\2\end{pmatrix} ⎝⎛110011111⎠⎞⎝⎛2−13⎠⎞=⎝⎛542⎠⎞
答案: ( 5 4 2 ) \begin{pmatrix}5\\4\\2\end{pmatrix} ⎝⎛542⎠⎞。
2.2 给定线性变换 T : R 3 → R 3 T:\mathbb{R}^3\rightarrow\mathbb{R}^3 T:R3→R3,使得 T ( 1 0 0 ) = ( 1 1 0 ) T\begin{pmatrix}1\\0\\0\end{pmatrix}=\begin{pmatrix}1\\1\\0\end{pmatrix} T⎝⎛100⎠⎞=⎝⎛110⎠⎞, T ( 0 1 0 ) = ( 0 1 1 ) T\begin{pmatrix}0\\1\\0\end{pmatrix}=\begin{pmatrix}0\\1\\1\end{pmatrix} T⎝⎛010⎠⎞=⎝⎛011⎠⎞, T ( 0 0 1 ) = ( 1 1 1 ) T\begin{pmatrix}0\\0\\1\end{pmatrix}=\begin{pmatrix}1\\1\\1\end{pmatrix} T⎝⎛001⎠⎞=⎝⎛111⎠⎞。求 T T T在基 { ( 1 0 0 ) , ( 0 1 0 ) , ( 0 0 1 ) } \{\begin{pmatrix}1\\0\\0\end{pmatrix},\begin{pmatrix}0\\1\\0\end{pmatrix},\begin{pmatrix}0\\0\\1\end{pmatrix}\} {⎝⎛100⎠⎞,⎝⎛010⎠⎞,⎝⎛001⎠⎞}下的矩阵。
解:由于 ( 1 0 0 ) \begin{pmatrix}1\\0\\0\end{pmatrix} ⎝⎛100⎠⎞在 T T T下的像为 ( 1 1 0 ) \begin{pmatrix}1\\1\\0\end{pmatrix} ⎝⎛110⎠⎞, ( 0 1 0 ) \begin{pmatrix}0\\1\\0\end{pmatrix} ⎝⎛010⎠⎞在 T T T下的像为 ( 0 1 1 ) \begin{pmatrix}0\\1\\1\end{pmatrix} ⎝⎛011⎠⎞, ( 0 0 1 ) \begin{pmatrix}0\\0\\1\end{pmatrix} ⎝⎛001⎠⎞在 T T T下的像为 ( 1 1 1 ) \begin{pmatrix}1\\1\\1\end{pmatrix} ⎝⎛111⎠⎞,因此 T T T在基 { ( 1 0 0 ) , ( 0 1 0 ) , ( 0 0 1 ) } \{\begin{pmatrix}1\\0\\0\end{pmatrix},\begin{pmatrix}0\\1\\0\end{pmatrix},\begin{pmatrix}0\\0\\1\end{pmatrix}\} {⎝⎛100⎠⎞,⎝⎛010⎠⎞,⎝⎛001⎠⎞}下的矩阵为
[ T ] = ( 1 0 1 1 1 1 0 1 1 ) [T]=\begin{pmatrix}1&0&1\\1&1&1\\0&1&1\end{pmatrix} [T]=⎝⎛110011111⎠⎞
答案: ( 1 0 1 1 1 1 0 1 1 ) \begin{pmatrix}1&0&1\\1&1&1\\0&1&1\end{pmatrix} ⎝⎛110011111⎠⎞。
试卷
第3大题(共10分)
已知线性变换 T : R 3 → R 3 T:\mathbb{R}^3\rightarrow\mathbb{R}^3 T:R3→R3的矩阵为 A = ( 1 2 1 0 1 1 0 0 2 ) A=\begin{pmatrix}1&2&1\\0&1&1\\0&0&2\end{pmatrix} A=⎝⎛100210112⎠⎞,求:
(1)线性变换 T T T的零空间和值空间;
(2)判断线性变换 T T T是否可以相似对角化,并说明理由。
解:
(1)线性变换 T T T的零空间和值空间分别为:
ker ( T ) = Nul ( A ) = { ( − 2 t − s − t − s t ) ∣ t , s ∈ R } , im ( T ) = Col ( A ) = Span { ( 1 0 0 ) , ( 2 1 0 ) , ( 1 1 2 ) } . \begin{aligned} \text{ker}(T)&=\text{Nul}(A)=\left\{\begin{pmatrix}-2t-s\\-t-s\\t\end{pmatrix}\Bigg|t,s\in\mathbb{R}\right\},\\ \text{im}(T)&=\text{Col}(A)=\text{Span}\left\{\begin{pmatrix}1\\0\\0\end{pmatrix},\begin{pmatrix}2\\1\\0\end{pmatrix},\begin{pmatrix}1\\1\\2\end{pmatrix}\right\}. \end{aligned} ker(T)im(T)=Nul(A)=⎩⎨⎧⎝⎛−2t−s−t−st⎠⎞∣∣∣∣∣t,s∈R⎭⎬⎫,=Col(A)=Span⎩⎨⎧⎝⎛100⎠⎞,⎝⎛210⎠⎞,⎝⎛112⎠⎞⎭⎬⎫.
(2)线性变换 T T T的特征多项式为:
det ( λ I − A ) = ∣ λ − 1 − 2 − 1 0 λ − 1 − 1 0 0 λ − 2 ∣ = ( λ − 1 ) 2 ( λ − 2 ) . \begin{aligned} \text{det}(\lambda I-A)&=\begin{vmatrix}\lambda-1&-2&-1\\0&\lambda-1&-1\\0&0&\lambda-2\end{vmatrix}\\ &= (\lambda-1)^2(\lambda-2). \end{aligned} det(λI−A)=∣∣∣∣∣∣λ−100−2λ−10−1−1λ−2∣∣∣∣∣∣=(λ−1)2(λ−2).
由于特征多项式的根不重复,所以线性变换 T T T可以相似对角化。
设 P P P为 A A A的对角化矩阵,即 A = P D P − 1 A=PDP^{-1} A=PDP−1,其中 D D D为 A A A的对角化矩阵,即
D = ( 1 0 0 0 1 0 0 0 2 ) . D=\begin{pmatrix}1&0&0\\0&1&0\\0&0&2\end{pmatrix}. D=⎝⎛100010002⎠⎞.
因此, T T T可以相似对角化,且对角化矩阵为 D D D。
答案:(1) ker ( T ) = { ( − 2 t − s − t − s t ) ∣ t , s ∈ R } \text{ker}(T)=\left\{\begin{pmatrix}-2t-s\\-t-s\\t\end{pmatrix}\Bigg|t,s\in\mathbb{R}\right\} ker(T)=⎩⎨⎧⎝⎛−2t−s−t−st⎠⎞∣∣∣∣∣t,s∈R⎭⎬⎫, im ( T ) = Span { ( 1 0 0 ) , ( 2 1 0 ) , ( 1 1 2 ) } \text{im}(T)=\text{Span}\left\{\begin{pmatrix}1\\0\\0\end{pmatrix},\begin{pmatrix}2\\1\\0\end{pmatrix},\begin{pmatrix}1\\1\\2\end{pmatrix}\right\} im(T)=Span⎩⎨⎧⎝⎛100⎠⎞,⎝⎛210⎠⎞,⎝⎛112⎠⎞⎭⎬⎫;(2)可以相似对角化。
第4大题(共10分)
已知矩阵 A = ( 1 1 0 0 1 1 0 0 1 ) A=\begin{pmatrix}1&1&0\\0&1&1\\0&0&1\end{pmatrix} A=⎝⎛100110011⎠⎞,求:
(1)矩阵 A A A的最小多项式和Jordan标准形;
(2)设 f ( x ) = ( x 1 0 0 x 1 0 0 x ) f(x)=\begin{pmatrix}x&1&0\\0&x&1\\0&0&x\end{pmatrix} f(x)=⎝⎛x001x001x⎠⎞,求 f ( A ) f(A) f(A)的Jordan标准形。
解:
(1)矩阵 A A A的特征多项式为:
det ( λ I − A ) = ∣ λ − 1 − 1 0 0 λ − 1 − 1 0 0 λ − 1 ∣ = ( λ − 1 ) 3 . \begin{aligned} \text{det}(\lambda I-A)&=\begin{vmatrix}\lambda-1&-1&0\\0&\lambda-1&-1\\0&0&\lambda-1\end{vmatrix}\\ &= (\lambda-1)^3. \end{aligned} det(λI−A)=∣∣∣∣∣∣λ−100−1λ−100−1λ−1∣∣∣∣∣∣=(λ−1)3.
因此,矩阵 A A A的特征值为 λ = 1 \lambda=1 λ=1,且代数重数为 3 3 3。
将 λ = 1 \lambda=1 λ=1代入 ( λ I − A ) (\lambda I-A) (λI−A),得到
( λ I − A ) = ( 0 − 1 0 0 0 − 1 0 0 0 ) . (\lambda I-A)=\begin{pmatrix}0&-1&0\\0&0&-1\\0&0&0\end{pmatrix}. (λI−A)=⎝⎛000−1000−10⎠⎞.
因此, ( λ I − A ) 3 = 0 (\lambda I-A)^3=0 (λI−A)3=0,即 ( A − I ) 3 = 0 (A-I)^3=0 (A−I)3=0。
所以,矩阵 A A A的最小多项式为 m ( x ) = ( x − 1 ) 2 m(x)=(x-1)^2 m(x)=(x−1)2,Jordan标准形为
J = ( 1 1 0 0 1 1 0 0 1 ) . J=\begin{pmatrix}1&1&0\\0&1&1\\0&0&1\end{pmatrix}. J=⎝⎛100110011⎠⎞.
(2)矩阵 f ( A ) f(A) f(A)的特征多项式为:
det ( λ I − f ( A ) ) = ∣ λ − A 11 − A 12 − A 13 − A 21 λ − A 22 − A 23 − A 31 − A 32 λ − A 33 ∣ = ∣ λ − 1 − 1 0 0 λ − 1 − 1 0 0 λ ∣ = λ ( λ − 1 ) 2 . \begin{aligned} \text{det}(\lambda I-f(A))&=\begin{vmatrix}\lambda-A_{11}&-A_{12}&-A_{13}\\-A_{21}&\lambda-A_{22}&-A_{23}\\-A_{31}&-A_{32}&\lambda-A_{33}\end{vmatrix}\\ &=\begin{vmatrix}\lambda-1&-1&0\\0&\lambda-1&-1\\0&0&\lambda\end{vmatrix}\\ &=\lambda(\lambda-1)^2. \end{aligned} det(λI−f(A))=∣∣∣∣∣∣λ−A11−A21−A31−A12λ−A22−A32−A13−A23λ−A33∣∣∣∣∣∣=∣∣∣∣∣∣λ−100−1λ−100−1λ∣∣∣∣∣∣=λ(λ−1)2.
因此,矩阵 f ( A ) f(A) f(A)的特征值为 λ = 0 , 1 \lambda=0,1 λ=0,1,且代数重数分别为 1 , 2 1,2 1,2。
将 λ = 0 \lambda=0 λ=0代入 ( λ I − f ( A ) ) (\lambda I-f(A)) (λI−f(A)),得到
( λ I − f ( A ) ) = ( − 1 − 1 0 0 − 1 − 1 0 0 0 ) . (\lambda I-f(A))=\begin{pmatrix}-1&-1&0\\0&-1&-1\\0&0&0\end{pmatrix}. (λI−f(A))=⎝⎛−100−1−100−10⎠⎞.
因此, ( λ I − f ( A ) ) 2 = ( 1 2 1 0 1 2 0 0 0 ) (\lambda I-f(A))^2=\begin{pmatrix}1&2&1\\0&1&2\\0&0&0\end{pmatrix} (λI−f(A))2=⎝⎛100210120⎠⎞。
将 λ = 1 \lambda=1 λ=1代入 ( λ I − f ( A ) ) (\lambda I-f(A)) (λI−f(A)),得到
( λ I − f ( A ) ) = ( 0 − 1 0 0 0 − 1 0 0 − 1 ) . (\lambda I-f(A))=\begin{pmatrix}0&-1&0\\0&0&-1\\0&0&-1\end{pmatrix}. (λI−f(A))=⎝⎛000−1000−1−1⎠⎞.
因此, ( λ I − f ( A ) ) 2 = 0 (\lambda I-f(A))^2=0 (λI−f(A))2=0。
所以,矩阵 f ( A ) f(A) f(A)的Jordan标准形为
J = ( 1 1 0 0 1 1 0 0 0 ) . J=\begin{pmatrix}1&1&0\\0&1&1\\0&0&0\end{pmatrix}. J=⎝⎛100110010⎠⎞.
答案:(1) m ( x ) = ( x − 1 ) 2 m(x)=(x-1)^2 m(x)=(x−1)2, J = ( 1 1 0 0 1 1 0 0 1 ) J=\begin{pmatrix}1&1&0\\0&1&1\\0&0&1\end{pmatrix} J=⎝⎛100110011⎠⎞;(2) J = ( 1 1 0 0 1 1 0 0 0 ) J=\begin{pmatrix}1&1&0\\0&1&1\\0&0&0\end{pmatrix} J=⎝⎛100110010⎠⎞。
第5大题(共10分)
假设总体 X X X 服从参数为 θ \theta θ 的指数分布,即 X ∼ E x p ( θ ) X\sim Exp(\theta) X∼Exp(θ),样本容量为 n n n,样本值为 x 1 , x 2 , ⋯ , x n x_1,x_2,\cdots,x_n x1,x2,⋯,xn。
(1) 求 θ \theta θ 的极大似然估计量。
解:总体 X X X 的概率密度函数为 f ( x ; θ ) = 1 θ e − x θ , x > 0 f(x;\theta)=\frac{1}{\theta}e^{-\frac{x}{\theta}},x>0 f(x;θ)=θ1e−θx,x>0。
样本的似然函数为 L ( θ ; x 1 , x 2 , ⋯ , x n ) = ∏ i = 1 n f ( x i ; θ ) = ∏ i = 1 n 1 θ e − x i θ = 1 θ n e − ∑ i = 1 n x i θ L(\theta;x_1,x_2,\cdots,x_n)=\prod_{i=1}^n f(x_i;\theta)=\prod_{i=1}^n\frac{1}{\theta}e^{-\frac{x_i}{\theta}}=\frac{1}{\theta^n}e^{-\frac{\sum_{i=1}^n x_i}{\theta}} L(θ;x1,x2,⋯,xn)=∏i=1nf(xi;θ)=∏i=1nθ1e−θxi=θn1e−θ∑i=1nxi。
取对数得到 ln L ( θ ; x 1 , x 2 , ⋯ , x n ) = − n ln θ − ∑ i = 1 n x i θ \ln L(\theta;x_1,x_2,\cdots,x_n)=-n\ln\theta-\frac{\sum_{i=1}^n x_i}{\theta} lnL(θ;x1,x2,⋯,xn)=−nlnθ−θ∑i=1nxi。
对 θ \theta θ 求导得到 d ln L ( θ ; x 1 , x 2 , ⋯ , x n ) d θ = − n θ + ∑ i = 1 n x i θ 2 \frac{\mathrm{d}\ln L(\theta;x_1,x_2,\cdots,x_n)}{\mathrm{d}\theta}=-\frac{n}{\theta}+\frac{\sum_{i=1}^n x_i}{\theta^2} dθdlnL(θ;x1,x2,⋯,xn)=−θn+θ2∑i=1nxi。
令导数等于零,解得 θ ^ = ∑ i = 1 n x i n \hat{\theta}=\frac{\sum_{i=1}^n x_i}{n} θ^=n∑i=1nxi。
(2) 证明 θ ^ \hat{\theta} θ^ 是 θ \theta θ 的最小方差无偏估计。
解:由于 θ ^ = ∑ i = 1 n x i n \hat{\theta}=\frac{\sum_{i=1}^n x_i}{n} θ^=n∑i=1nxi,因此 θ ^ \hat{\theta} θ^ 是 θ \theta θ 的无偏估计。
又因为 V a r ( θ ^ ) = V a r ( ∑ i = 1 n x i n ) = 1 n 2 ∑ i = 1 n V a r ( x i ) = θ 2 n Var(\hat{\theta})=Var(\frac{\sum_{i=1}^n x_i}{n})=\frac{1}{n^2}\sum_{i=1}^n Var(x_i)=\frac{\theta^2}{n} Var(θ^)=Var(n∑i=1nxi)=n21∑i=1nVar(xi)=nθ2。
因此, θ ^ \hat{\theta} θ^ 的方差为 θ 2 n \frac{\theta^2}{n} nθ2,是 θ \theta θ 的最小方差无偏估计。
第6大题(共10分)
假设总体 X X X 服从正态分布 N ( μ , σ 2 ) N(\mu,\sigma^2) N(μ,σ2),样本容量为 n n n,样本均值为 x ‾ \overline{x} x,样本标准差为 s s s。
(1) 假设检验的原假设和备择假设分别为 H 0 : μ = μ 0 H_0:\mu=\mu_0 H0:μ=μ0 和 H 1 : μ ≠ μ 0 H_1:\mu\neq\mu_0 H1:μ=μ0,其中 μ 0 \mu_0 μ0 是已知的常数。给出检验统计量的表达式。
解:检验统计量为 t = x ‾ − μ 0 s / n t=\frac{\overline{x}-\mu_0}{s/\sqrt{n}} t=s/nx−μ0,其中 x ‾ \overline{x} x 是样本均值, s s s 是样本标准差, n n n 是样本容量。
(2) 假设检验的拒绝域为 ∣ t ∣ > t α / 2 ( n − 1 ) |t|>t_{\alpha/2}(n-1) ∣t∣>tα/2(n−1),其中 α \alpha α 是显著性水平, t α / 2 ( n − 1 ) t_{\alpha/2}(n-1) tα/2(n−1) 是自由度为 n − 1 n-1 n−1 的 t t t 分布上侧 α / 2 \alpha/2 α/2 分位数。给出样本容量的确定方法,使得第Ⅱ类错误概率不超过 β \beta β。
解:第Ⅱ类错误概率为 β \beta β,即 P ( ∣ t ∣ ≤ t α / 2 ( n − 1 ) ∣ μ = μ 1 ) = 1 − β P(|t|\leq t_{\alpha/2}(n-1)|\mu=\mu_1)=1-\beta P(∣t∣≤tα/2(n−1)∣μ=μ1)=1−β,其中 μ 1 \mu_1 μ1 是真实的总体均值。
由于 t t t 统计量服从自由度为 n − 1 n-1 n−1 的 t t t 分布,因此 t α / 2 ( n − 1 ) t_{\alpha/2}(n-1) tα/2(n−1) 是常数。为了使第Ⅱ类错误概率不超过 β \beta β,需要满足 P ( ∣ t ∣ ≤ t α / 2 ( n − 1 ) ∣ μ = μ 1 ) = 1 − β P(|t|\leq t_{\alpha/2}(n-1)|\mu=\mu_1)=1-\beta P(∣t∣≤tα/2(n−1)∣μ=μ1)=1−β,即 P ( − t α / 2 ( n − 1 ) ≤ x ‾ − μ 1 s / n ≤ t α / 2 ( n − 1 ) ) = 1 − β P(-t_{\alpha/2}(n-1)\leq\frac{\overline{x}-\mu_1}{s/\sqrt{n}}\leq t_{\alpha/2}(n-1))=1-\beta P(−tα/2(n−1)≤s/nx−μ1≤tα/2(n−1))=1−β。
根据标准正态分布的性质,可得 P ( − z α / 2 ≤ x ‾ − μ 1 s / n ≤ z α / 2 ) = 1 − β P(-z_{\alpha/2}\leq\frac{\overline{x}-\mu_1}{s/\sqrt{n}}\leq z_{\alpha/2})=1-\beta P(−zα/2≤s/nx−μ1≤zα/2)=1−β,其中 z α / 2 z_{\alpha/2} zα/2 是标准正态分布上侧 α / 2 \alpha/2 α/2 分位数。
因此,需要满足 z α / 2 s n = ∣ μ 1 − μ 0 ∣ n ≤ t α / 2 ( n − 1 ) \frac{z_{\alpha/2}s}{\sqrt{n}}=\frac{|\mu_1-\mu_0|}{\sqrt{n}}\leq t_{\alpha/2}(n-1) nzα/2s=n∣μ1−μ0∣≤tα/2(n−1),即 n ≥ ( z α / 2 s ∣ μ 1 − μ 0 ∣ ) 2 n\geq(\frac{z_{\alpha/2}s}{|\mu_1-\mu_0|})^2 n≥(∣μ1−μ0∣zα/2s)2。
第7大题(共10分)
已知线性正态回归模型
y
i
=
β
1
x
i
1
+
β
2
x
i
2
+
⋯
+
β
p
x
i
p
+
ϵ
i
y_i = \beta_1 x_{i1} + \beta_2 x_{i2} + \cdots + \beta_p x_{ip} + \epsilon_i
yi=β1xi1+β2xi2+⋯+βpxip+ϵi,其中
i
=
1
,
2
,
⋯
,
n
i=1,2,\cdots,n
i=1,2,⋯,n,
ϵ
i
∼
N
(
0
,
σ
2
)
\epsilon_i\sim N(0,\sigma^2)
ϵi∼N(0,σ2),
x
i
1
,
x
i
2
,
⋯
,
x
i
p
x_{i1},x_{i2},\cdots,x_{ip}
xi1,xi2,⋯,xip 是已知的常数,
y
1
,
y
2
,
⋯
,
y
n
y_1,y_2,\cdots,y_n
y1,y2,⋯,yn 是观测值。现在要求该模型参数的最小二乘估计,并判断这些估计量是否独立。
解答过程:
定义观测数据的矩阵表示 Y = ( y 1 , y 2 , ⋯ , y n ) ′ Y=(y_1,y_2,\cdots,y_n)' Y=(y1,y2,⋯,yn)′,协变量的矩阵表示 X = ( x i 1 , x i 2 , ⋯ , x i p ) X=(x_{i1},x_{i2},\cdots,x_{ip}) X=(xi1,xi2,⋯,xip),将 β \beta β 的最小二乘估计量表示为 β ^ = ( β 1 ^ , β 2 ^ , ⋯ , β p ^ ) ′ \hat{\beta}=(\hat{\beta_1},\hat{\beta_2},\cdots,\hat{\beta_p})' β^=(β1^,β2^,⋯,βp^)′。
首先,计算 β ^ \hat{\beta} β^:
β ^ = ( X ′ X ) − 1 X ′ Y \hat{\beta}=(X'X)^{-1}X'Y β^=(X′X)−1X′Y
其中, ( X ′ X ) − 1 (X'X)^{-1} (X′X)−1 表示矩阵 ( X ′ X ) (X'X) (X′X) 的逆矩阵。这个公式就是最小二乘估计量的标准表达式。下面需要证明这个公式的正确性。
构造最小二乘估计量的误差函数:
S ( β ) = ∑ i = 1 n ( y i − β 1 x i 1 − β 2 x i 2 − ⋯ − β p x i p ) 2 S(\beta)=\sum_{i=1}^n (y_i-\beta_1 x_{i1}-\beta_2 x_{i2}-\cdots-\beta_p x_{ip})^2 S(β)=i=1∑n(yi−β1xi1−β2xi2−⋯−βpxip)2
对 β \beta β 求导,令导数为零:
∂ S ∂ β = 2 ∑ i = 1 n x i ( y i − β 1 x i 1 − β 2 x i 2 − ⋯ − β p x i p ) = 0 \frac{\partial S}{\partial \beta}=2\sum_{i=1}^n x_i(y_i-\beta_1 x_{i1}-\beta_2 x_{i2}-\cdots-\beta_p x_{ip})=0 ∂β∂S=2i=1∑nxi(yi−β1xi1−β2xi2−⋯−βpxip)=0
其中, x i x_i xi 表示第 i i i 个观测对应的协变量向量。
化简得:
X ′ ( Y − X β ) = 0 X'(Y-X\beta)=0 X′(Y−Xβ)=0
解得:
β ^ = ( X ′ X ) − 1 X ′ Y \hat{\beta}=(X'X)^{-1}X'Y β^=(X′X)−1X′Y
因此,最小二乘估计量的标准公式得到证明。
接下来,需要判断这些估计量是否独立。
由于 ϵ i \epsilon_i ϵi 是正态分布的,因此
β ^ ∼ N p ( β , ( X ′ X ) − 1 σ 2 ) \hat{\beta} \sim N_p(\beta,(X'X)^{-1}\sigma^2) β^∼Np(β,(X′X)−1σ2)
其中, N p N_p Np 表示 p p p 维正态分布, ( X ′ X ) − 1 σ 2 (X'X)^{-1}\sigma^2 (X′X)−1σ2 是 β ^ \hat{\beta} β^ 的协方差矩阵。
可以发现, β ^ \hat{\beta} β^ 的分布不依赖于 σ 2 \sigma^2 σ2,因此, β ^ \hat{\beta} β^ 与 σ 2 \sigma^2 σ2 独立。
而这意味着, β ^ \hat{\beta} β^ 的每个分量都与 σ 2 \sigma^2 σ2 独立,即最小二乘估计量的参数是独立的。
答案:
线性正态回归模型参数的最小二乘估计是 β ^ = ( X ′ X ) − 1 X ′ Y \hat{\beta}=(X'X)^{-1}X'Y β^=(X′X)−1X′Y,这些估计量独立。