2021 National Graduate Entrance Examination Mathematics I Past Paper

I. Multiple Choice

I. Multiple Choice: Questions 1-10, 5 points each, 50 points in total. For each of the following questions, only one of the four options given is correct; write the letter in front of the chosen option on the designated position of the answer sheet.

1. The function

f(x)={ex−1x,x≠0,1,x=0f(x)= \begin{cases} \dfrac{e^x-1}{x}, & x\neq 0,\\ 1, & x=0 \end{cases}

is ( ) at x=0x=0

(A) continuous and attains a maximum
(B) continuous and attains a minimum
(C) differentiable with derivative equal to 00
(D) differentiable with derivative not equal to 00

limitdifferentiability

Approach:

This question tests the limit and the differentiability of a piecewise function at its breakpoint. First check continuity by comparing the limit with the value at the point, and then compute the derivative at the point from the difference quotient; the two results together determine the correct option.

Steps:

  1. Continuity at x=0x=0: lim⁡x→0ex−1x=1=f(0)\lim\limits_{x\to0} \frac{e^x-1}{x}=1=f(0), so ff is continuous at 00; this excludes the options claiming a maximum or a minimum, since the function is given by one and the same formula on both sides of 00 and therefore has no extremum there.
  2. Differentiability at x=0x=0: lim⁡x→0f(x)−f(0)x−0=lim⁡x→0ex−1−xx2=12\lim\limits_{x\to0}\dfrac{f(x)-f(0)}{x-0}=\lim\limits_{x\to0}\dfrac{e^x-1 - x}{x^2}=\tfrac{1}{2}, so the derivative exists and is not 00, which rules out (C) and leaves (D).

Answer:

D.

2. Let the function f(x,y)f(x,y) be differentiable and satisfy

fx(x+1,ey)+eyfy(x+1,ey)=(x+1)2+2x(x+1)(1)f_x(x+1,e^y)+e^y f_y(x+1,e^y)=(x+1)^2+2x(x+1) \quad (1) fx(x,x2)+2xfy(x,x2)=4xln⁡x+2x(2)f_x(x,x^2)+2x f_y(x,x^2)=4x\ln x+2x \quad (2)

then df(1,1)=df(1,1)= ( )

(A) dx+dydx+dy
(B) dx−dydx-dy
(C) dydy
(D) −dy-dy

partial derivativetotal differential

Approach:

This question tests the use of the chain rule along two different curves in order to determine the partial derivatives of a function of two variables, and then the construction of the total differential from them. The two identities hold on two curves; evaluating them at suitable points and solving the resulting linear system gives fx(1,1)f_x(1,1) and fy(1,1)f_y(1,1).

Steps:

  1. In (1) take x=0,y=0x=0,y=0: the arguments become x+1=1x+1=1 and ey=1e^y=1, and the right-hand side is (0+1)2+2⋅0⋅(0+1)=1(0+1)^2+2\cdot0\cdot(0+1)=1, so fx(1,1)+fy(1,1)=1f_x(1,1)+f_y(1,1)=1.
  2. In (2) take x=1x=1: the arguments become x=1x=1 and x2=1x^2=1, and the right-hand side is 4⋅1ln⁡1+2⋅1=24\cdot1\ln1+2\cdot1=2, so fx(1,1)+2fy(1,1)=2f_x(1,1)+2f_y(1,1)=2.
  3. Solving these two linear equations gives fx(1,1)=0, fy(1,1)=1f_x(1,1)=0,\ f_y(1,1)=1.
  4. Hence df(1,1)=fx(1,1)dx+fy(1,1)dy=dydf(1,1)=f_x(1,1)dx+f_y(1,1)dy=dy, which is option (C).

Answer:

C.

3. Let the cubic Taylor polynomial of the function f(x)=sin⁡x1+x2f(x)=\dfrac{\sin x}{1+x^2} at x=0x=0 be ax+bx2+cx3ax+bx^2+cx^3, then ( )

(A) a=1, b=0, c=−76a=1,\ b=0,\ c=-\dfrac{7}{6}
(B) a=1, b=0, c=76a=1,\ b=0,\ c=\dfrac{7}{6}
(C) a=1, b=−1, c=−76a=1,\ b=-1,\ c=-\dfrac{7}{6}
(D) a=−1, b=−1, c=−76a=-1,\ b=-1,\ c=-\dfrac{7}{6}

Taylor expansioninfinitesimal

Approach:

This question tests the composition of Taylor expansions. The function is a product of two factors, so both factors are expanded to the required order and the product is then formed, keeping the terms up to x3x^3.

Steps:

  1. Expand the first factor: sin⁡x=x−x36+o(x3)\sin x = x-\dfrac{x^3}{6}+o(x^3).
  2. Expand the second factor: 11+x2=1−x2+o(x2)\dfrac{1}{1+x^2}=1 - x^2 + o(x^2).
  3. Multiply and keep the terms up to x3x^3: f(x)=(x−x36+o(x3))(1−x2+o(x2))=x−x36−x3+o(x3)=x−76x3+o(x3)f(x)=\left(x-\dfrac{x^3}{6}+o(x^3)\right)\left(1-x^2+o(x^2)\right)=x-\dfrac{x^3}{6}-x^3+o(x^3)=x-\dfrac{7}{6}x^3+o(x^3).
  4. Comparing this with ax+bx2+cx3ax+bx^2+cx^3 gives a=1, b=0, c=−76a=1,\ b=0,\ c=-\dfrac{7}{6}, so (A) is the correct choice.

Answer:

A.

4. Let the function f(x)f(x) be continuous on the interval [0,1][0,1], then ∫01f(x) dx=\int_0^1 f(x)\,dx= ( )

(A) lim⁡n→∞∑k=1nf ⁣(2k−12n)12n\displaystyle \lim_{n\to\infty}\sum_{k=1}^n f\!\left(\frac{2k-1}{2n}\right)\frac{1}{2n}
(B) lim⁡n→∞∑k=1nf ⁣(2k−12n)1n\displaystyle \lim_{n\to\infty}\sum_{k=1}^n f\!\left(\frac{2k-1}{2n}\right)\frac{1}{n}
(C) lim⁡n→∞∑k=12nf ⁣(k−12n)1n\displaystyle \lim_{n\to\infty}\sum_{k=1}^{2n} f\!\left(\frac{k-1}{2n}\right)\frac{1}{n}
(D) lim⁡n→∞∑k=12nf ⁣(k2n)2n\displaystyle \lim_{n\to\infty}\sum_{k=1}^{2n} f\!\left(\frac{k}{2n}\right)\frac{2}{n}

definite integralRiemann sum

Approach:

This question tests the definition of the definite integral as a limit of Riemann sums. Divide [0,1][0,1] into nn equal subintervals and check, option by option, whether the number of sample points, the length of each subinterval and the step size are consistent with each other.

Steps:

  1. The definition of the definite integral as a Riemann sum: divide [0,1][0,1] into nn equal subintervals, each of length 1n\dfrac1n, and take the midpoint ξk=2k−12n\xi_k=\dfrac{2k-1}{2n} of each subinterval; then ∫01f(x)dx=lim⁡n→∞∑k=1nf ⁣(2k−12n)1n\displaystyle\int_0^1f(x)dx=\lim_{n\to\infty}\sum_{k=1}^nf\!\left(\dfrac{2k-1}{2n}\right)\dfrac1n (this holds when ff is continuous on [0,1][0,1]).
  2. Check the options one by one: in (A) there are nn sample points but the step size is 12n\dfrac{1}{2n}, so the total length is 12\dfrac12 and it is not ∫01f(x)dx\int_0^1f(x)dx; in (C) there are 2n2n sample points but the step size is 1n\dfrac1n, so the total length is 22, which is likewise wrong.
  3. In option (D) the sample points k2n\dfrac{k}{2n} are the right endpoints rather than the midpoints (the first one is 00, an endpoint), and the step size is 2n\dfrac2n, so the total length is 22; this is not the required integral either.
  4. Only in option (B) are the sample points the midpoints, the number of them is nn and the step size is 1n\dfrac1n, so the three match the Riemann sum of ∫01f(x)dx\int_0^1f(x)dx exactly, and (B) is the correct choice.

Answer:

B.

5. The positive index and the negative index of inertia of the quadratic form f(x1,x2,x3)=(x1+x2)2+(x2+x3)2−(x3−x1)2f(x_1,x_2,x_3)=(x_1+x_2)^2+(x_2+x_3)^2-(x_3-x_1)^2 are respectively ( )

(A) 2, 0
(B) 1, 1
(C) 2, 1
(D) 1, 2

quadratic formindex of inertia

Approach:

This question tests the index of inertia of a quadratic form. First expand the quadratic form and write it in matrix form, then compute the eigenvalues; the signs of the eigenvalues give the positive and the negative index of inertia.

Steps:

  1. Expand the quadratic form: f(x1,x2,x3)=(x1+x2)2+(x2+x3)2−(x3−x1)2=2x22+2x1x2+2x2x3+2x1x3f(x_1,x_2,x_3)=(x_1+x_2)^2+(x_2+x_3)^2-(x_3-x_1)^2=2x_2^2+2x_1x_2+2x_2x_3+2x_1x_3, so the matrix of the quadratic form is A=(011121110)A=\begin{pmatrix}0&1&1\\1&2&1\\1&1&0\end{pmatrix}.
  2. Compute the characteristic polynomial: ∣λE−A∣=λ(λ−3)(λ+1)|\lambda E-A|=\lambda(\lambda-3)(\lambda+1), so the eigenvalues are λ1=3, λ2=0, λ3=−1\lambda_1=3,\ \lambda_2=0,\ \lambda_3=-1.
  3. The values above are confirmed by tr⁡A=2\operatorname{tr}A=2, ∣A∣=0|A|=0 and tr⁡A2=10\operatorname{tr}A^2=10.
  4. There is only one positive eigenvalue, namely 33, and only one negative eigenvalue, namely −1-1, so the positive index of inertia is 1 and the negative index of inertia is 1, and (B) is the correct choice.

Answer:

B.

6. Let α1=(100)\alpha_1=\begin{pmatrix}1\\0\\0\end{pmatrix}, α2=(121)\alpha_2=\begin{pmatrix}1\\2\\1\end{pmatrix}, α3=(321)\alpha_3=\begin{pmatrix}3\\2\\1\end{pmatrix}, and put β1=α1, β2=α2−kβ1, β3=α3−l1β1−l2β2\beta_1=\alpha_1,\ \beta_2=\alpha_2-k\beta_1,\ \beta_3=\alpha_3-l_1\beta_1-l_2\beta_2. If β1,β2,β3\beta_1,\beta_2,\beta_3 are pairwise orthogonal, then l1,l2l_1,l_2 are respectively ( )

(A) 52, 12\dfrac{5}{2},\ \dfrac{1}{2}
(B) 52, 12⋅2\dfrac{5}{2},\ \dfrac{1}{2\cdot2}
(C) 52, −12\dfrac{5}{2},\ -\dfrac{1}{2}
(D) −52, 12-\dfrac{5}{2},\ \dfrac{1}{2}

orthogonalizationGram-Schmidt

Approach:

This question tests the Gram-Schmidt orthogonalisation of a system of vectors. The orthogonalisation formula is used to remove, in turn, the components of each vector along the directions of the vectors already obtained.

Steps:

  1. First β1=α1=(1,0,0)T\beta_1=\alpha_1=(1,0,0)^T is already normalised, and the coefficient that makes β2\beta_2 orthogonal to it is k=[α2,β1][β1,β1]=2k=\dfrac{[\alpha_2,\beta_1]}{[\beta_1,\beta_1]}=2, so β2=α2−2α1=(0,2,0)T\beta_2=\alpha_2-2\alpha_1=(0,2,0)^T.
  2. The coefficient that makes β3\beta_3 orthogonal to β1\beta_1 is l1=[α3,β1][β1,β1]=52l_1=\dfrac{[\alpha_3,\beta_1]}{[\beta_1,\beta_1]}=\dfrac{5}{2}.
  3. The coefficient that makes β3\beta_3 orthogonal to β2\beta_2 is l2=[α3,β2][β2,β2]=12l_2=\dfrac{[\alpha_3,\beta_2]}{[\beta_2,\beta_2]}=\dfrac{1}{2}.
  4. With these values the three vectors are pairwise orthogonal, so l1=52l_1=\dfrac{5}{2} and l2=12l_2=\dfrac{1}{2}, and (A) is the correct choice.

Answer:

A.

7. Let A,BA,B be real matrices of order nn, then the statement that does not hold is ( )

(A) r(A00ATA)=2r(A)r\begin{pmatrix}A&0\\0&A^T A\end{pmatrix}=2r(A)
(B) r(A0ABAT)=2r(A)r\begin{pmatrix}A&0\\AB&A^T\end{pmatrix}=2r(A)
(C) r(A0BAAT)=2r(A)r\begin{pmatrix}A&0\\BA&A^T\end{pmatrix}=2r(A)
(D) r(AB0AT)=2r(A)r\begin{pmatrix}A&B\\0&A^T\end{pmatrix}=2r(A)

matrix rankblock matrix

Approach:

This question tests the rank of block matrices. The rank of a block matrix (A0CAT)\begin{pmatrix}A&0\\ C&A^T\end{pmatrix} is the dimension obtained by superimposing the row (column) spaces of AA and of CC and then superimposing the row space of ATA^T, so each option has to be examined separately.

Steps:

  1. For option (A): the block matrix is (A00ATA)\begin{pmatrix}A&0\\0&A^TA\end{pmatrix}, the two blocks are independent of each other, and the rank is r(A)+r(ATA)=2r(A)r(A)+r(A^TA)=2r(A), so it holds.
  2. For option (B): (A0ABAT)\begin{pmatrix}A&0\\AB&A^T\end{pmatrix}; the row vectors of ABAB are linear combinations of the row vectors of AA and add no new dimension to the row space, so the rank is r(A)+r(AT)=2r(A)r(A)+r(A^T)=2r(A), and it holds.
  3. For option (C): (A0BAAT)\begin{pmatrix}A&0\\BA&A^T\end{pmatrix}; the rows of the lower block BABA are in general not linear combinations of the rows of AA (they combine the rows of BB with the columns of AA) and may bring an extra dimension to the row space, making the rank greater than 2r(A)2r(A), so the equality does not always hold.
  4. For option (D): (AB0AT)\begin{pmatrix}A&B\\0&A^T\end{pmatrix}; in the same way the rows of BB may bring an extra dimension to the row space, so this equality does not always hold either. The statement that does not hold is option (C).

Answer:

C.

8. Let A,BA,B be random events with 0<P(B)<10<P(B)<1, then the statement that does not hold is ( )

(A) If P(A∣B)=P(A)P(A|B)=P(A), then P(A∣Bˉ)=P(A)P(A|\bar B)=P(A)
(B) If P(A∣B)>P(A)P(A|B)>P(A), then P(A∣Bˉ)>P(A)P(A|\bar B)>P(A)
(C) If P(A∣B)>P(A∣Bˉ)P(A|B)>P(A|\bar B), then P(A∣B)>P(A)P(A|B)>P(A)
(D) If P(A∣A∪B)≥P(Aˉ∣A∪B)P(A|A\cup B)\geq P(\bar A|A\cup B), then P(A)>P(B)P(A)>P(B)

conditional probabilitylaw of total probability

Approach:

This question tests the monotonic relations between conditional probabilities. Each statement is checked by expanding the conditional probabilities with the law of total probability, and a counterexample is looked for where the relation fails.

Steps:

  1. By the law of total probability P(A∣A∪B)=P(A)P(A)+P(B)−P(AB)P(A|A\cup B)=\dfrac{P(A)}{P(A)+P(B)-P(AB)} and P(Aˉ∣A∪B)=P(A∪B)−P(A)P(A)+P(B)−P(AB)P(\bar A|A\cup B)=\dfrac{P(A\cup B)-P(A)}{P(A)+P(B)-P(AB)}, so the condition P(A∣A∪B)≥P(Aˉ∣A∪B)P(A|A\cup B)\geq P(\bar A|A\cup B) only says P(A)≥P(A∪B)−P(A)P(A)\geq P(A\cup B)-P(A), that is 2P(A)≥P(A∪B)−P(A)2P(A)\geq P(A\cup B)-P(A), and it does not force P(A)>P(B)P(A)>P(B).
  2. For instance, taking P(A)=P(B)=12P(A)=P(B)=\tfrac12 and P(AB)=0P(AB)=0 gives P(A∣A∪B)=P(Aˉ∣A∪B)=12P(A|A\cup B)=P(\bar A|A\cup B)=\tfrac12, so the hypothesis of (D) holds while P(A)=P(B)P(A)=P(B), and (D) does not hold.
  3. For (A): P(A∣B)=P(A)P(A|B)=P(A) means that AA and BB are independent, and then P(A∣Bˉ)=P(A)P(A|\bar B)=P(A) as well, so (A) holds.
  4. For (B) and (C): by the law of total probability P(A)=P(A∣B)P(B)+P(A∣Bˉ)(1−P(B))P(A)=P(A|B)P(B)+P(A|\bar B)(1-P(B)), so P(A∣B)>P(A)P(A|B)>P(A) forces P(A∣Bˉ)>P(A)P(A|\bar B)>P(A), and P(A∣B)>P(A∣Bˉ)P(A|B)>P(A|\bar B) forces P(A∣B)>P(A)P(A|B)>P(A); hence (B) and (C) hold, and the statement that does not hold is (D).

Answer:

D.

9. Let (X1,Y1),…,(Xn,Yn)(X_1,Y_1),\ldots,(X_n,Y_n) be a simple random sample from the population N(μ1,μ2;σ12,σ22;ρ)N(\mu_1,\mu_2;\sigma_1^2,\sigma_2^2;\rho), and put θ=μ1−μ2, Xˉ=1n∑i=1nXi, Yˉ=1n∑i=1nYi, θ^=Xˉ−Yˉ\theta=\mu_1-\mu_2,\ \bar X=\frac{1}{n}\sum_{i=1}^n X_i,\ \bar Y=\frac{1}{n}\sum_{i=1}^n Y_i,\ \hat\theta=\bar X-\bar Y, then ( )

(A) θ^\hat\theta is an unbiased estimator of θ\theta and D(θ^)=σ12+σ22nD(\hat\theta)=\dfrac{\sigma_1^2+\sigma_2^2}{n}
(B) θ^\hat\theta is not an unbiased estimator of θ\theta and D(θ^)=σ12+σ22D(\hat\theta)=\sigma_1^2+\sigma_2^2
(C) θ^\hat\theta is an unbiased estimator of θ\theta and D(θ^)=σ12+σ22−2ρσ1σ2nD(\hat\theta)=\dfrac{\sigma_1^2+\sigma_2^2-2\rho\sigma_1\sigma_2}{n}
(D) θ^\hat\theta is an unbiased estimator of θ\theta and D(θ^)=σ12+σ22−2ρσ1σ2nD(\hat\theta)=\dfrac{\sigma_1^2+\sigma_2^2-2\rho\sigma_1\sigma_2}{n}

unbiased estimatorvariancecorrelation coefficient

Approach:

This question tests the unbiasedness of a linear combination of sample means and the computation of its variance. The expectation is obtained from the linearity of the expectation, and the variance is obtained from the covariance between the two sample means.

Steps:

  1. By the linearity of the expectation, E(θ^)=E(Xˉ)−E(Yˉ)=μ1−μ2=θE(\hat\theta)=E(\bar X)-E(\bar Y)=\mu_1-\mu_2=\theta, so θ^\hat\theta is an unbiased estimator of θ\theta.
  2. The variance of a difference of two correlated sample means is D(Xˉ−Yˉ)=D(Xˉ)+D(Yˉ)−2Cov⁡(Xˉ,Yˉ)D(\bar X-\bar Y)=D(\bar X)+D(\bar Y)-2\operatorname{Cov}(\bar X,\bar Y).
  3. Here D(Xˉ)=σ12nD(\bar X)=\dfrac{\sigma_1^2}{n}, D(Yˉ)=σ22nD(\bar Y)=\dfrac{\sigma_2^2}{n}, and Cov⁡(Xˉ,Yˉ)=Cov⁡(X1,Y1)n=ρσ1σ2n\operatorname{Cov}(\bar X,\bar Y)=\dfrac{\operatorname{Cov}(X_1,Y_1)}{n}=\dfrac{\rho\sigma_1\sigma_2}{n}.
  4. Hence D(θ^)=σ12+σ22−2ρσ1σ2nD(\hat\theta)=\dfrac{\sigma_1^2+\sigma_2^2-2\rho\sigma_1\sigma_2}{n}, so (C) is the correct choice.

Answer:

C.

10. Let X1,X2,…,X16X_1,X_2,\ldots,X_{16} be a simple random sample from the population N(μ,4)N(\mu,4), and consider the hypothesis testing problem H0:μ≤10, H1:μ>10H_0:\mu\leq10,\ H_1:\mu>10 with rejection region W={Xˉ≥11}W=\{\bar X\geq11\}, where Xˉ=116∑i=116Xi\bar X=\frac{1}{16}\sum_{i=1}^{16}X_i. Then when μ=11.5\mu=11.5 the probability of a type II error is ( )

(A) 1−Φ(0.5)1-\Phi(0.5)
(B) 1−Φ(1)1-\Phi(1)
(C) 1−Φ(1.5)1-\Phi(1.5)
(D) 1−Φ(2)1-\Phi(2)

hypothesis testingnormal distributiontype II error

Approach:

This question tests the probability of a type II error in a hypothesis test. The sample mean is standardised with the help of the normal distribution, and the probability of not rejecting the null hypothesis when it is false is then read off from the standard normal distribution function.

Steps:

  1. Since X1,…,X16X_1,\dots,X_{16} are independent and follow N(μ,4)N(\mu,4), the sample mean satisfies Xˉ∼N ⁣(μ,14)\bar X\sim N\!\left(\mu,\tfrac{1}{4}\right), so its standard deviation is 12\tfrac12.
  2. When μ=11.5\mu=11.5, the rejection region is {Xˉ≥11}\{\bar X\geq11\}, so a type II error occurs exactly when Xˉ<11\bar X<11.
  3. Therefore P(Xˉ<11)=P ⁣(Z<11−11.51/2)=P(Z<−1)=1−Φ(1)P(\bar X<11)=P\!\left(Z<\dfrac{11-11.5}{1/2}\right)=P(Z<-1)=1-\Phi(1), so (B) is the correct choice.

Answer:

B.

II. Fill in the Blanks

II. Fill in the Blanks: Questions 11-16, 5 points each, 30 points in total. Write the answers on the designated positions of the answer sheet.

11. ∫0+∞dxx2+2x+2=\displaystyle \int_0^{+\infty}\frac{dx}{x^2+2x+2}= ____.

definite integralinverse trigonometric function

Approach:

This question tests the evaluation of an improper integral by completing the square. The quadratic x2+2x+2x^2+2x+2 is rewritten as (x+1)2+1(x+1)^2+1, after which the integral becomes the antiderivative of 11+(x+1)2\dfrac{1}{1+(x+1)^2}, that is, an arctangent, evaluated between the two limits.

Steps:

  1. Complete the square: x2+2x+2=(x+1)2+1x^2+2x+2=(x+1)^2+1.
  2. With the substitution u=x+1u=x+1 the lower limit x=0x=0 becomes u=1u=1, while the improper upper limit stays infinite.
  3. The antiderivative is arctan⁡u\arctan u, so ∫0+∞dx(x+1)2+1=arctan⁡(x+1)∣0+∞=π2−π4=π4\int_0^{+\infty}\dfrac{dx}{(x+1)^2+1}=\left.\arctan(x+1)\right|_0^{+\infty}=\dfrac{\pi}{2}-\dfrac{\pi}{4}=\dfrac{\pi}{4}.

Answer:

π4\dfrac{\pi}{4}.

12. Let the function y=y(x)y=y(x) be determined by the parametric equations

{x=2et+t+1,y=4(t−1)et+t2\begin{cases} x=2e^t+t+1,\\ y=4(t-1)e^t+t^2 \end{cases}

Find y′′∣x=0=y''\big|_{x=0}= ____.

parametric equationsecond derivative

Approach:

This question tests the second derivative of a function given by parametric equations. First dydx\dfrac{dy}{dx} is computed from the ratios of the derivatives with respect to tt, and then the result is differentiated once more with respect to tt and divided by dxdt\dfrac{dx}{dt}.

Steps:

  1. From x=2et+t+1x=2e^t+t+1 and y=4(t−1)et+t2y=4(t-1)e^t+t^2 we get dxdt=2et+1\dfrac{dx}{dt}=2e^t+1 and dydt=4et+2t\dfrac{dy}{dt}=4e^t+2t.
  2. Hence dydx=4et+2t2et+1\dfrac{dy}{dx}=\dfrac{4e^t+2t}{2e^t+1}.
  3. Differentiating once more with respect to tt and dividing by dxdt\dfrac{dx}{dt} gives d2ydx2=(4et+4et+2)(2et+1)−(4et+2t)2et(2et+1)3\dfrac{d^2y}{dx^2}=\dfrac{(4e^t+4e^t+2)(2e^t+1)-(4e^t+2t)2e^t}{(2e^t+1)^3}.
  4. The condition x=0x=0 corresponds to t=0t=0, and substituting t=0t=0 gives 23\dfrac{2}{3}.

Answer:

23\dfrac{2}{3}.

13. The solution of the Euler equation x2y′′+y′−4y=0x^2 y''+y'-4y=0 satisfying y(1)=1, y′(1)=2y(1)=1,\ y'(1)=2 is ____.

Euler equationordinary differential equation

Approach:

This question tests the solution of a Cauchy-Euler equation. Substituting a trial solution y=xλy=x^\lambda gives the characteristic equation, whose roots determine the general solution; the initial conditions then fix the two constants.

Steps:

  1. Put y=xλy=x^\lambda; then y′=λxλ−1y'=\lambda x^{\lambda-1} and y′′=λ(λ−1)xλ−2y''=\lambda(\lambda-1)x^{\lambda-2}, and the equation x2y′′+y′−4y=0x^2y''+y'-4y=0 becomes λ(λ−1)xλ+λxλ−1−4xλ=0\lambda(\lambda-1)x^\lambda+\lambda x^{\lambda-1}-4x^\lambda=0, that is λ2−4=0\lambda^2-4=0.
  2. The roots are λ=2\lambda=2 and λ=−2\lambda=-2, so the general solution is y=C1x2+C2x−2y=C_1 x^2+C_2 x^{-2}.
  3. The initial values give y(1)=C1+C2=1, y′(x)=2C1x−2C2x−3y(1)=C_1+C_2=1,\ y'(x)=2C_1 x-2C_2 x^{-3}, and y′(1)=2C1−2C2=2y'(1)=2C_1-2C_2=2.
  4. Solving the two equations gives C1=1,C2=0C_1=1,C_2=0, so y=x2y=x^2.

Answer:

y=x2y=x^2.

14. Let Σ\Sigma be the outer side of the surface of the region {(x,y,z)∣x2+4y2≤4, 0≤z≤2}\{(x,y,z)\mid x^2+4y^2\leq4,\,0\leq z\leq2\}. Find the surface integral ∬Σx2 dy dz+y2 dz dx+z dx dy=\displaystyle \iint_\Sigma x^2\,dy\,dz + y^2\,dz\,dx + z\,dx\,dy= ____.

surface integraldivergence theorem

Approach:

This question tests the application of the divergence theorem to a surface integral. The original integral is a second-kind surface integral over the whole outer surface of the closed region {(x,y,z)∣x2+4y2≤4, 0≤z≤2}\{(x,y,z)\mid x^2+4y^2\le4,\ 0\le z\le2\}, so it can be turned directly into a triple integral.

Steps:

  1. Write the integral as the flux of a vector field: F=(x2, y2, z)F=(x^2,\ y^2,\ z), and Σ\Sigma is the whole outer surface of the given region (the lateral side together with the upper and the lower cap).
  2. By the divergence theorem, ∬Σx2 dy dz+y2 dz dx+z dx dy=∭Ωdiv⁡F dV\displaystyle\iint_{\Sigma}x^2\,dy\,dz+y^2\,dz\,dx+z\,dx\,dy=\iiint_{\Omega}\operatorname{div}F\,dV, where div⁡F=2x+2y+1\operatorname{div}F=2x+2y+1.
  3. The region Ω={(x,y,z)∣x2+4y2≤4, 0≤z≤2}\Omega=\{(x,y,z)\mid x^2+4y^2\le4,\ 0\le z\le2\} is symmetric about both the xx-axis and the yy-axis, so ∭Ω2x dV=0\iiint_\Omega2x\,dV=0 and ∭Ω2y dV=0\iiint_\Omega2y\,dV=0.
  4. The integral is therefore ∭Ω1 dV\iiint_\Omega1\,dV, that is, the volume. The cross-section is an ellipse with semiaxes 22 and 11, of area π⋅2⋅1=2π\pi\cdot2\cdot1=2\pi, and the height is 22, so the volume is 2π⋅2=4π2\pi\cdot2=4\pi.

Answer:

4π4\pi.

15. Let AA be a 3×33\times3 matrix, let AjA_j be the cofactor of the entry in row jj deleted, and suppose that the entries of every row of AA sum to 22 and that ∣A∣=3|A|=3. Find A11+A21+A31=A_{11}+A_{21}+A_{31}= ____.

determinanteigenvalue

Approach:

This question tests the relation between a determinant and the cofactors of a matrix. Since the entries of every row of AA sum to 22, the vector (1,1,1)T(1,1,1)^T is an eigenvector of AA with eigenvalue 22, and the adjugate can then be written in terms of that eigenvector.

Steps:

  1. The condition on the row sums gives Aα=λαA\alpha=\lambda\alpha with α=(1,1,1)T\alpha=(1,1,1)^T and λ=2\lambda=2.
  2. The adjugate satisfies A∗=∣A∣λααTA^*=\dfrac{|A|}{\lambda}\alpha\alpha^T, and the sum of the entries of the first column of A∗A^* is exactly A11+A21+A31A_{11}+A_{21}+A_{31}.
  3. That sum is ∣A∣λ⋅(1+1+1)⋅1=∣A∣λ\dfrac{|A|}{\lambda}\cdot(1+1+1)\cdot1=\dfrac{|A|}{\lambda}, and with ∣A∣=3|A|=3 and λ=2\lambda=2 this gives A11+A21+A31=∣A∣λ=32A_{11}+A_{21}+A_{31}=\dfrac{|A|}{\lambda}= \dfrac{3}{2}.

Answer:

32\dfrac{3}{2}.

16. Box A and box B each contain 2 red balls and 2 white balls. A ball is first drawn at random from box A, its colour is observed and the ball is then put into box B, after which a ball is drawn at random from box B. Let XX and YY denote the numbers of red balls drawn from box A and from box B respectively. Find the correlation coefficient of XX and YY ____.

discrete distributioncorrelation coefficient

Approach:

This question tests the computation of the correlation coefficient of two discrete random variables. The joint distribution of (X,Y)(X,Y) is listed first, and the covariance and the two variances are then computed from it.

Steps:

  1. The joint distribution is (X,Y)∼{(0,0),(0,1),(1,0),(1,1)}(X,Y)\sim \{(0,0),(0,1),(1,0),(1,1)\} with probabilities 3/10,1/5,1/5,3/103/10,1/5,1/5,3/10 respectively.
  2. Hence EX=EY=12EX=EY=\dfrac12, so DX=DY=14DX=DY=\dfrac14.
  3. Moreover E(XY)=310E(XY)=\dfrac{3}{10} and EX⋅EY=14EX\cdot EY=\dfrac14, so cov(X,Y)=E(XY)−EX⋅EY=310−14=120\mathrm{cov}(X,Y)=E(XY)-EX\cdot EY=\dfrac{3}{10}-\dfrac14=\dfrac1{20}.
  4. Therefore the correlation coefficient is ρXY=1/201/4⋅1/4=15\rho_{XY}=\dfrac{1/20}{\sqrt{1/4\cdot1/4}}=\dfrac{1}{5}.

Answer:

15\dfrac{1}{5}.

III. Problems

III. Problems: Questions 17-22, 70 points in total. Write the solutions on the designated positions of the answer sheet. The solutions should include verbal explanations, proof procedures or computational steps.

17. (10 points) Find the limit

lim⁡x→0[1+∫0xt2etex−1 dt−1sin⁡x]x2.\lim_{x\to0}\frac{\left[1+\int_0^x \frac{t^2 e^t}{e^x-1}\,dt-\dfrac{1}{\sin x}\right]}{x^2}.
L'Hopital's ruleTaylor expansion

Approach:

This question tests the use of a small-order expansion to evaluate a limit whose numerator contains an integral with a variable upper limit. The numerator is expanded up to the order x2x^2, and its two contributions are then combined.

Steps:

  1. For the integral, since ex−1∼xe^x-1\sim x and et∼1e^t\sim1 on [0,x][0,x], we have ∫0xt2etex−1dt∼∫0xt2dt=x33+o(x3)\int_0^x \frac{t^2 e^t}{e^x-1}dt\sim\int_0^x t^2dt=\frac{x^3}{3}+o(x^3).
  2. Also sin⁡x=x−x36+o(x3)\sin x = x-\frac{x^3}{6}+o(x^3), so 1sin⁡x=1x(1+x26+o(x2))\dfrac{1}{\sin x}=\dfrac{1}{x}\left(1+\frac{x^2}{6}+o(x^2)\right).
  3. The numerator is therefore 1+x33+o(x3)−1sin⁡x=12x2+o(x2)1+\frac{x^3}{3}+o(x^3)-\dfrac{1}{\sin x}=\frac{1}{2}x^2+o(x^2).
  4. Dividing by x2x^2 and taking the limit gives 12\dfrac{1}{2}.

Answer:

12\dfrac{1}{2}.

18. (12 points) Let un(x)=e−x 1n(n+1) xn+1 (n=1,2,…)u_n(x)=e^{-x}\,\dfrac{1}{n(n+1)}\,x^{n+1}\ (n=1,2,\ldots). Find the domain of convergence and the sum function of the series ∑n=1∞un(x)\sum_{n=1}^\infty u_n(x).

power seriesuniform convergence

Approach:

The factor e−xe^{-x} does not affect the convergence, so one first studies the power series ∑n=1∞xn+1n(n+1)\sum_{n=1}^{\infty}\dfrac{x^{n+1}}{n(n+1)} and then sums it with the help of 1n(n+1)=1n−1n+1\dfrac{1}{n(n+1)}=\dfrac{1}{n}-\dfrac{1}{n+1}.

Steps:

  1. The ratio test gives the radius R=1R=1.
  2. At x=1x=1 the series ∑1n(n+1)\sum\dfrac{1}{n(n+1)} converges, and at x=−1x=-1 it converges absolutely, so the domain of convergence is [−1,1][-1,1].
  3. For ∣x∣<1|x|<1 we have
S(x)=∑n=1∞e−xxn+1n(n+1)=e−x[(1−x)ln⁡(1−x)+x].S(x)= \sum_{n=1}^{\infty}e^{-x}\dfrac{x^{n+1}}{n(n+1)} =e^{-x}\left[(1-x)\ln(1-x)+x\right].
  1. Substituting the endpoints into the original series gives
S(1)=e−1,S(−1)=e(2ln⁡2−1).S(1)=e^{-1},\qquad S(-1)=e(2\ln2-1).

Answer:

S(x)={e(2ln⁡2−1),x=−1,e−x[(1−x)ln⁡(1−x)+x],−1<x<1,e−1,x=1.S(x)= \begin{cases} e(2\ln2-1),&x=-1,\\[4pt] e^{-x}\left[(1-x)\ln(1-x)+x\right],&-1<x<1,\\[4pt] e^{-1},&x=1. \end{cases}

19. (12 points) Given the curve

{x2+2y2−z=6,4x+2y+z=30\begin{cases} x^2+2y^2 - z = 6,\\ 4x+2y+z = 30 \end{cases}

find the maximum distance from the points of CC to the xOyxOy coordinate plane.

Lagrange multiplier methodspace curve

Approach:

This question tests the maximum distance from the points of a space curve to a coordinate plane. First zz is eliminated by means of the equation of the plane, so that the curve is projected onto the xOyxOy plane as an ellipse, and then the extrema of ∣z∣|z| are found with the Lagrange multiplier method.

Steps:

  1. From the second equation z=30−4x−2yz=30-4x-2y, so the distance to the xOyxOy plane is ∣z∣=∣30−4x−2y∣|z|=|30-4x-2y|.
  2. Substituting z=30−4x−2yz=30-4x-2y into the first equation x2+2y2−z=6x^2+2y^2-z=6 gives x2+4x+2y2+2y−36=0x^2+4x+2y^2+2y-36=0, that is (x+2)2+2(y+12)2=812(x+2)^2+2\left(y+\tfrac12\right)^2=\tfrac{81}{2}, which shows that the projection of the curve onto the xOyxOy plane is an ellipse.
  3. Applying the Lagrange multiplier method to ∣z∣|z| gives the extremal points (4,1,12)(4,1,12) and (−8,−2,66)(-8,-2,66), and substituting back into z=30−4x−2yz=30-4x-2y gives z=12z=12 and z=66z=66 respectively.
  4. Comparing the values of ∣z∣|z| at the two ends, the maximum distance from the points of the curve to the xOyxOy plane is 6666.

Answer:

6666.

20. (12 points) Let D⊂R2D\subset\mathbb{R}^2 be a bounded simply connected closed region, and let I(D)=∬D(4−x2−y2) dx dyI(D)=\displaystyle\iint_D(4-x^2-y^2)\,dx\,dy; denote by D1D_1 the region of integration at which I(D)I(D) attains its maximum.
(1) Find the value of I(D1)I(D_1).
(2) Compute

∮∂D1(xex2+4y2+y) dx+(4yex2+4y2−x) dyx2+4y2,\oint_{\partial D_1}\frac{(xe^{x^2+4y^2}+y)\,dx+(4y e^{x^2+4y^2}-x)\,dy}{x^2+4y^2},

where ∂D1\partial D_1 is the positively oriented boundary of D1D_1.

area integralGreen's theorem

Approach:

The integrand of the area integral is symmetric about the origin, so the maximum corresponds to the disc x2+y2≤4x^2+y^2\le 4; the curve integral is then evaluated with Green’s theorem.

Steps:

  1. The maximum of I(D)I(D) is attained when x2+y2≤4x^2+y^2\le4, and then I(D1)=∫02π∫02(4−r2)r dr dθ=8πI(D_1)=\int_0^{2\pi}\int_0^{2}(4-r^2)r\,dr\,d\theta=8\pi.
  2. Along ∂D1\partial D_1, Green’s theorem turns the curve integral into an area integral, whose value is −π-\pi.

Answer:

I(D1)=8πI(D_1)=8\pi, and the curve integral is −π-\pi.

21. (12 points) Given

A=(a1−11a−1−11a).A=\begin{pmatrix} a&1&-1\\ 1&a&-1\\ -1&1&a \end{pmatrix}.

(1) Find the orthogonal matrix PP such that PTAPP^T A P is diagonal;
(2) Find the positive definite matrix CC such that C2=(a+3)E−AC^2=(a+3)E-A.

eigenvalue decompositionorthogonal similarity

Approach:

The eigenvalues of AA are λ1=a+2, λ2=λ3=a−1\lambda_1=a+2,\ \lambda_2=\lambda_3=a-1, and the corresponding orthogonal eigenvectors are used to build PP; one then sets C=PΛ1/2PTC=P\Lambda^{1/2}P^T.

Steps:

  1. The characteristic polynomial of AA factors as (λ−a−2)(λ−a+1)2(\lambda-a-2)(\lambda-a+1)^2, so the eigenvalues are a+2a+2 and a−1a-1 (double), and the orthogonal eigenvectors give the matrix PP below.
  2. One may take
P=(131216131216−13026),P=\begin{pmatrix} \tfrac{1}{\sqrt3}&\tfrac{1}{\sqrt2}&\tfrac{1}{\sqrt6}\\ \tfrac{1}{\sqrt3}&\tfrac{1}{\sqrt2}&\tfrac{1}{\sqrt6}\\ -\tfrac{1}{\sqrt3}&0&\tfrac{2}{\sqrt6} \end{pmatrix},

and then PTAP=diag⁡(a+2,a−1,a−1)P^T A P=\operatorname{diag}(a+2,a-1,a-1). 3. Using C=PΛ1/2PTC=P\Lambda^{1/2}P^T one obtains

C=(53−1−1−15313−11353).C=\begin{pmatrix} \tfrac{5}{3}&-1&-1\\ -1&\tfrac{5}{3}&\tfrac{1}{3}\\ -1&\tfrac{1}{3}&\tfrac{5}{3} \end{pmatrix}.

Answer:

PP and CC are as above.

22. (12 points) A point is chosen at random on the interval (0,2)(0,2), which is then divided into two parts; let XX be the length of the shorter part and YY the length of the longer part, and put Z=YXZ=\dfrac{Y}{X}.
(1) Find the density of XX;
(2) Find the density of ZZ;
(3) Find E ⁣(XY)E\!\left(\dfrac{X}{Y}\right).

geometric probabilitytransformation of random variables

Approach:

The cutting point is uniformly distributed, so first X∼U(0,1)X\sim U(0,1); then the distribution function of Z=2−XXZ=\dfrac{2-X}{X} is computed and differentiated, and finally the expectation is evaluated from its definition.

Steps:

  1. The cutting point is uniform on (0,2)(0,2), so the shorter part satisfies fX(x)=1, 0<x<1f_X(x)=1,\ 0<x<1, and 00 otherwise.
  2. FZ(z)=P ⁣(2−XX≤z)=P ⁣(X≥2z+1)F_Z(z)=P\!\left(\dfrac{2-X}{X}\le z\right)=P\!\left(X\ge\dfrac{2}{z+1}\right), hence fZ(z)=2(z+1)2, z≥1f_Z(z)=\dfrac{2}{(z+1)^2},\ z\ge1, and 00 otherwise.
  3. E ⁣(XY)=∫01x2−x dx=−1+2ln⁡2E\!\left(\dfrac{X}{Y}\right)=\int_0^1 \dfrac{x}{2-x}\,dx=-1+2\ln2.

Answer:

fX(x)=1 (0<x<1)f_X(x)=1\ (0<x<1); fZ(z)=2(z+1)2 (z≥1)f_Z(z)=\dfrac{2}{(z+1)^2}\ (z\ge1); E ⁣(XY)=−1+2ln⁡2E\!\left(\dfrac{X}{Y}\right)=-1+2\ln2.

Found an issue with this content?

Tell us what could be improved.