Skip to content

STEP3 1998 -- Pure Mathematics

STEP3 1998 — Section A (Pure Mathematics)

Section titled “STEP3 1998 — Section A (Pure Mathematics)”

Exam: STEP3  |  Year: 1998  |  Questions: Q1—Q8  |  Total marks per question: 20

All questions are pure mathematics. Solutions and examiner commentary are included below.

QTopicDifficultyKey Techniques
1三角函数与函数变换 Trigonometric Functions and TransformationsChallenging三角恒等变换,求导求驻点,复合函数求导,分式函数分析
2积分与特殊函数 Integration and Special FunctionsChallenging分部积分,递推关系,对称性论证,Beta函数
3递推关系与差分方程 Recurrence RelationsStandard一阶线性递推,特解与齐次解,极限过渡
4极坐标与面积计算 Polar Coordinates and AreaChallenging极坐标转直角坐标,极坐标面积公式,三角函数周期性,对称性
5矩阵指数与几何变换 Matrix Exponential and Geometric TransformationsChallenging矩阵幂级数,泰勒展开,M^2=-I的利用,旋转与剪切变换
6几何 GeometryChallenging坐标几何,体积公式,对称性分析,向量运算
7微积分 CalculusHard函数图像分析,求导找极值,换元积分,量纲分析,近似估算
8向量与几何 Vectors and GeometryChallenging向量方程,点积运算,球面方程,正交条件,切面分析

Topic: 三角函数与函数变换 Trigonometric Functions and Transformations  |  Difficulty: Challenging  |  Marks: 20

1 Let f(x)=sin2x+2cosx+1f(x) = \sin^2 x + 2 \cos x + 1 for 0x2π0 \leqslant x \leqslant 2\pi. Sketch the curve y=f(x)y = f(x), giving the coordinates of the stationary points. Now let g(x)=af(x)+bcf(x)+dadbc, d3c, dc.g(x) = \frac{af(x) + b}{cf(x) + d} \qquad ad \neq bc \, , \ d \neq -3c \, , \ d \neq c \, . Show that the stationary points of y=g(x)y = g(x) occur at the same values of xx as those of y=f(x)y = f(x), and find the corresponding values of g(x)g(x).

Explain why, if d/c<3d/c < -3 or d/c>1d/c > 1, g(x)|g(x)| cannot be arbitrarily large.

Model Solution

Stationary points of f(x)f(x)

Using sin2x=1cos2x\sin^2 x = 1 - \cos^2 x, we rewrite:

f(x)=1cos2x+2cosx+1=2+2cosxcos2x=3(cosx1)2.f(x) = 1 - \cos^2 x + 2\cos x + 1 = 2 + 2\cos x - \cos^2 x = 3 - (\cos x - 1)^2.

Differentiating:

f(x)=2sinxcosx2sinx=2sinx(cosx1).f'(x) = 2\sin x \cos x - 2\sin x = 2\sin x(\cos x - 1).

Setting f(x)=0f'(x) = 0: either sinx=0\sin x = 0 or cosx=1\cos x = 1.

In [0,2π][0, 2\pi]: sinx=0\sin x = 0 at x=0,π,2πx = 0, \pi, 2\pi, and cosx=1\cos x = 1 at x=0,2πx = 0, 2\pi. So the stationary points are x=0,π,2πx = 0, \pi, 2\pi.

Classifying the stationary points:

f(x)=2cosx(cosx1)+2sinx(sinx)=2cos2x2cosx2sin2x.f''(x) = 2\cos x(\cos x - 1) + 2\sin x(-\sin x) = 2\cos^2 x - 2\cos x - 2\sin^2 x.

Using sin2x=1cos2x\sin^2 x = 1 - \cos^2 x:

f(x)=2cos2x2cosx2(1cos2x)=4cos2x2cosx2.f''(x) = 2\cos^2 x - 2\cos x - 2(1 - \cos^2 x) = 4\cos^2 x - 2\cos x - 2.

At x=πx = \pi: f(π)=4(1)2(1)2=4>0f''(\pi) = 4(1) - 2(-1) - 2 = 4 > 0, so x=πx = \pi is a minimum with f(π)=3(11)2=34=1f(\pi) = 3 - (-1-1)^2 = 3 - 4 = -1.

At x=0x = 0: f(0)=4(1)2(1)2=0f''(0) = 4(1) - 2(1) - 2 = 0, so the second derivative test is inconclusive. But since f(x)=3(cosx1)2f(x) = 3 - (\cos x - 1)^2 and (cosx1)20(\cos x - 1)^2 \geqslant 0 with equality at x=0,2πx = 0, 2\pi, these are maxima with f=3f = 3.

Sketch: ff starts at f(0)=3f(0) = 3, decreases monotonically to f(π)=1f(\pi) = -1, then increases monotonically back to f(2π)=3f(2\pi) = 3. The curve is symmetric about x=πx = \pi.


Stationary points of g(x)g(x)

Let u=f(x)u = f(x). Then g(x)=au+bcu+dg(x) = \dfrac{au + b}{cu + d}, and by the chain rule:

g(x)=ddu ⁣(au+bcu+d)f(x)=(adbc)(cu+d)2f(x).g'(x) = \frac{d}{du}\!\left(\frac{au + b}{cu + d}\right) \cdot f'(x) = \frac{(ad - bc)}{(cu + d)^2} \cdot f'(x).

Since adbcad \neq bc, the factor adbc(cu+d)20\dfrac{ad - bc}{(cu + d)^2} \neq 0 wherever it is defined. Thus g(x)=0g'(x) = 0 if and only if f(x)=0f'(x) = 0, provided the denominator cf(x)+d0cf(x) + d \neq 0.

The stationary values of ff are f=3f = 3 and f=1f = -1. The constraints d3cd \neq -3c (i.e. 3c+d03c + d \neq 0) and dcd \neq c (i.e. c+d0-c + d \neq 0) ensure cf(x)+d0cf(x) + d \neq 0 at both stationary values. Therefore g(x)=0g'(x) = 0 at exactly the same xx-values as f(x)=0f'(x) = 0. \blacksquare

Corresponding values of gg:

At x=0x = 0 and x=2πx = 2\pi (where f=3f = 3):

g=3a+b3c+d.g = \frac{3a + b}{3c + d}.

At x=πx = \pi (where f=1f = -1):

g=a+bc+d=badc.g = \frac{-a + b}{-c + d} = \frac{b - a}{d - c}.


Why g(x)|g(x)| cannot be arbitrarily large when d/c<3d/c < -3 or d/c>1d/c > 1

We write cf(x)+d=c(f(x)+d/c)cf(x) + d = c(f(x) + d/c). For g(x)|g(x)| to be arbitrarily large, we need cf(x)+dcf(x) + d to approach zero, which requires f(x)=d/cf(x) = -d/c for some x[0,2π]x \in [0, 2\pi].

The range of ff on [0,2π][0, 2\pi] is [1,3][-1, 3].

  • If d/c>1d/c > 1: then d/c<1-d/c < -1, which is below the range of ff. So f(x)+d/c>0f(x) + d/c > 0 for all xx, and in fact f(x)+d/c1+d/c>0f(x) + d/c \geqslant -1 + d/c > 0.
  • If d/c<3d/c < -3: then d/c>3-d/c > 3, which is above the range of ff. So f(x)+d/c<0f(x) + d/c < 0 for all xx, and in fact f(x)+d/c3+d/c<0f(x) + d/c \leqslant 3 + d/c < 0.

In both cases, cf(x)+d=cf(x)+d/c|cf(x) + d| = |c| \cdot |f(x) + d/c| is bounded below by a positive constant. Since af(x)+b|af(x) + b| is bounded (as ff is continuous on the closed interval [0,2π][0, 2\pi]), it follows that g(x)|g(x)| is bounded and cannot be arbitrarily large. \blacksquare


Topic: 积分与特殊函数 Integration and Special Functions  |  Difficulty: Challenging  |  Marks: 20

2 Let I(a,b)=01ta(1t)bdt(a0, b0).I(a, b) = \int_0^1 t^a (1 - t)^b \, dt \qquad (a \geqslant 0, \ b \geqslant 0). (i) Show that I(a,b)=I(b,a)I(a, b) = I(b, a),

(ii) Show that I(a,b)=I(a+1,b)+I(a,b+1)I(a, b) = I(a + 1, b) + I(a, b + 1).

(iii) Show that (a+1)I(a,b)=bI(a+1,b1)(a+1)I(a, b) = bI(a+1, b-1) when aa and bb are positive and hence calculate I(a,b)I(a, b) when aa and bb are positive integers.

Model Solution

Part (i): Show that I(a,b)=I(b,a)I(a, b) = I(b, a)

Substitute u=1tu = 1 - t, so du=dtdu = -dt. When t=0t = 0, u=1u = 1; when t=1t = 1, u=0u = 0.

I(a,b)=01ta(1t)bdt=10(1u)aub(du)=01ub(1u)adu=I(b,a).I(a, b) = \int_0^1 t^a (1-t)^b \, dt = \int_1^0 (1-u)^a u^b (-du) = \int_0^1 u^b (1-u)^a \, du = I(b, a). \qquad \blacksquare


Part (ii): Show that I(a,b)=I(a+1,b)+I(a,b+1)I(a, b) = I(a+1, b) + I(a, b+1)

Since 1=t+(1t)1 = t + (1 - t), we multiply the integrand by 11:

I(a,b)=01ta(1t)b1dt=01ta(1t)b[t+(1t)]dtI(a, b) = \int_0^1 t^a(1-t)^b \cdot 1 \, dt = \int_0^1 t^a(1-t)^b \bigl[t + (1-t)\bigr] \, dt

=01ta+1(1t)bdt+01ta(1t)b+1dt=I(a+1,b)+I(a,b+1).= \int_0^1 t^{a+1}(1-t)^b \, dt + \int_0^1 t^a(1-t)^{b+1} \, dt = I(a+1, b) + I(a, b+1). \qquad \blacksquare


Part (iii): Show (a+1)I(a,b)=bI(a+1,b1)(a+1)I(a, b) = bI(a+1, b-1) and calculate I(a,b)I(a,b) for positive integers

Deriving the recurrence:

Use integration by parts on I(a,b)=01ta(1t)bdtI(a, b) = \int_0^1 t^a(1-t)^b \, dt with:

u=ta,dv=(1t)bdtdu=ata1dt,v=(1t)b+1b+1.u = t^a, \quad dv = (1-t)^b \, dt \qquad \Longrightarrow \qquad du = at^{a-1} \, dt, \quad v = -\frac{(1-t)^{b+1}}{b+1}.

I(a,b)=[ta(1t)b+1b+1]01+ab+101ta1(1t)b+1dt.I(a,b) = \left[-\frac{t^a(1-t)^{b+1}}{b+1}\right]_0^1 + \frac{a}{b+1}\int_0^1 t^{a-1}(1-t)^{b+1} \, dt.

The boundary term vanishes (at t=0t=0 the factor ta=0t^a = 0 for a>0a > 0; at t=1t=1 the factor (1t)b+1=0(1-t)^{b+1} = 0). So:

I(a,b)=ab+1I(a1,b+1).()I(a,b) = \frac{a}{b+1}\,I(a-1, b+1). \qquad (\star)

Now apply the symmetry from Part (i) to ()(\star) with (a,b)(a,b) replaced by (b,a)(b,a):

I(b,a)=ba+1I(b1,a+1).I(b,a) = \frac{b}{a+1}\,I(b-1, a+1).

Using I(b,a)=I(a,b)I(b,a) = I(a,b) and I(b1,a+1)=I(a+1,b1)I(b-1, a+1) = I(a+1, b-1) (by symmetry again):

I(a,b)=ba+1I(a+1,b1),I(a,b) = \frac{b}{a+1}\,I(a+1, b-1),

which gives (a+1)I(a,b)=bI(a+1,b1)(a+1)I(a,b) = bI(a+1, b-1). \qquad \blacksquare


Calculating I(a,b)I(a, b) for positive integers:

The recurrence (a+1)I(a,b)=bI(a+1,b1)(a+1)I(a,b) = bI(a+1, b-1) can be rearranged as:

I(a+1,b1)=a+1bI(a,b).I(a+1, b-1) = \frac{a+1}{b}\,I(a,b).

Applying this repeatedly to reduce the second argument:

I(a,b)=a+1bI(a+1,b1)=(a+1)(a+2)b(b1)I(a+2,b2)=I(a,b) = \frac{a+1}{b}\,I(a+1, b-1) = \frac{(a+1)(a+2)}{b(b-1)}\,I(a+2, b-2) = \cdots

After bb applications:

I(a,b)=(a+1)(a+2)(a+b)b!I(a+b,0).I(a,b) = \frac{(a+1)(a+2)\cdots(a+b)}{b!}\,I(a+b, 0).

We compute the base case directly:

I(a+b,0)=01ta+bdt=1a+b+1.I(a+b, 0) = \int_0^1 t^{a+b} \, dt = \frac{1}{a+b+1}.

Therefore:

I(a,b)=(a+1)(a+2)(a+b)b!1a+b+1=(a+b)!a!b!1a+b+1=a!b!(a+b+1)!.I(a,b) = \frac{(a+1)(a+2)\cdots(a+b)}{b!} \cdot \frac{1}{a+b+1} = \frac{(a+b)!}{a!\,b!} \cdot \frac{1}{a+b+1} = \frac{a!\,b!}{(a+b+1)!}.

Hence for positive integers aa and bb:

I(a,b)=a!b!(a+b+1)!.I(a,b) = \frac{a!\,b!}{(a+b+1)!}. \qquad \blacksquare

Verification: I(1,2)=01t(1t)2dt=01(t2t2+t3)dt=1223+14=112I(1,2) = \int_0^1 t(1-t)^2\,dt = \int_0^1 (t - 2t^2 + t^3)\,dt = \frac{1}{2} - \frac{2}{3} + \frac{1}{4} = \frac{1}{12}, and 1!2!4!=224=112\frac{1!\,2!}{4!} = \frac{2}{24} = \frac{1}{12}. \checkmark


Topic: 递推关系与差分方程 Recurrence Relations  |  Difficulty: Standard  |  Marks: 20

3 The value VNV_N of a bond after NN days is determined by the equation VN+1=(1+c)VNd(c>0, d>0),V_{N+1} = (1 + c)V_N - d \qquad (c > 0, \ d > 0), where cc and dd are given constants. By looking for solutions of the form VT=AkT+BV_T = Ak^T + B for some constants A,BA, B and kk, or otherwise, find VNV_N in terms of V0V_0.

What is the solution for c=0c = 0? Show that this is the limit (for fixed NN) as c0c \to 0 of your solution for c>0c > 0.

Model Solution

We seek a solution of the form VT=AkT+BV_T = Ak^T + B for constants AA, BB, and kk.

Part (i): Finding VNV_N for c>0c > 0

Substituting VT=AkT+BV_T = Ak^T + B into the recurrence VT+1=(1+c)VTdV_{T+1} = (1+c)V_T - d:

AkT+1+B=(1+c)(AkT+B)d=(1+c)AkT+(1+c)Bd.Ak^{T+1} + B = (1+c)(Ak^T + B) - d = (1+c)Ak^T + (1+c)B - d.

Comparing coefficients of kTk^T:

Ak=(1+c)A    k=1+c(A0).Ak = (1+c)A \implies k = 1 + c \qquad (A \neq 0).

Comparing constant terms:

B=(1+c)Bd    cB=d    B=dc.B = (1+c)B - d \implies cB = d \implies B = \frac{d}{c}.

So the general solution is VN=A(1+c)N+dcV_N = A(1+c)^N + \dfrac{d}{c} for arbitrary constant AA.

Setting N=0N = 0: V0=A+dcV_0 = A + \dfrac{d}{c}, so A=V0dcA = V_0 - \dfrac{d}{c}.

Therefore:

VN=(V0dc)(1+c)N+dc(c>0).V_N = \left(V_0 - \frac{d}{c}\right)(1+c)^N + \frac{d}{c} \qquad (c > 0).

Part (ii): Solution for c=0c = 0

When c=0c = 0, the recurrence becomes VN+1=VNdV_{N+1} = V_N - d, which is a simple arithmetic progression:

VN=V0Nd.V_N = V_0 - Nd.

Part (iii): Showing the limit as c0c \to 0

We need to show that for fixed NN:

limc0[(V0dc)(1+c)N+dc]=V0Nd.\lim_{c \to 0} \left[\left(V_0 - \frac{d}{c}\right)(1+c)^N + \frac{d}{c}\right] = V_0 - Nd.

Using the binomial expansion (valid for fixed NN and small cc):

(1+c)N=1+Nc+(N2)c2+(1+c)^N = 1 + Nc + \binom{N}{2}c^2 + \cdots

Substituting:

(V0dc) ⁣(1+Nc+O(c2))+dc\left(V_0 - \frac{d}{c}\right)\!\left(1 + Nc + O(c^2)\right) + \frac{d}{c}

=V0+V0Nc+V0O(c2)dcdNdO(c)+dc= V_0 + V_0 Nc + V_0 \cdot O(c^2) - \frac{d}{c} - dN - d \cdot O(c) + \frac{d}{c}

The dc\dfrac{d}{c} terms cancel:

=V0+V0Nc+O(c2)dNdO(c)= V_0 + V_0 Nc + O(c^2) - dN - d \cdot O(c)

=V0dN+c(V0NdO(1))+O(c2).= V_0 - dN + c\bigl(V_0 N - d \cdot O(1)\bigr) + O(c^2).

As c0c \to 0:

limc0VN=V0Nd,\lim_{c \to 0} V_N = V_0 - Nd,

which is exactly the solution for c=0c = 0. \blacksquare


Topic: 极坐标与面积计算 Polar Coordinates and Area  |  Difficulty: Challenging  |  Marks: 20

4 Show that the equation (in plane polar coordinates) r=cosθr = \cos \theta, for π2θπ2-\frac{\pi}{2} \leqslant \theta \leqslant \frac{\pi}{2}, represents a circle.

Sketch the curve r=cos2θr = \cos 2\theta for 0θ2π0 \leqslant \theta \leqslant 2\pi, and describe the curves r=cos2nθr = \cos 2n\theta, where nn is an integer. Show that the area enclosed by such a curve is independent of nn.

Sketch also the curve r=cos3θr = \cos 3\theta for 0θ2π0 \leqslant \theta \leqslant 2\pi.

Model Solution

Part (i): r=cosθr = \cos\theta is a circle

Multiply both sides by rr:

r2=rcosθ.r^2 = r\cos\theta.

Using r2=x2+y2r^2 = x^2 + y^2 and rcosθ=xr\cos\theta = x:

x2+y2=x.x^2 + y^2 = x.

Completing the square:

(x12)2+y2=14.\left(x - \frac{1}{2}\right)^2 + y^2 = \frac{1}{4}.

This is a circle with centre (12,0)\left(\dfrac{1}{2}, 0\right) and radius 12\dfrac{1}{2}.

For π2θπ2-\dfrac{\pi}{2} \leqslant \theta \leqslant \dfrac{\pi}{2}, we have r=cosθ0r = \cos\theta \geqslant 0, so no extra points arise from negative rr. The curve traces out the full circle exactly once. \blacksquare


Part (ii): r=cos2θr = \cos 2\theta and the general case r=cos2nθr = \cos 2n\theta

Sketch of r=cos2θr = \cos 2\theta:

The function cos2θ\cos 2\theta has period π\pi. Over [0,2π][0, 2\pi] it completes two full periods. The curve has r>0r > 0 in the intervals [0,π4]\left[0, \frac{\pi}{4}\right], [3π4,5π4]\left[\frac{3\pi}{4}, \frac{5\pi}{4}\right], [7π4,2π]\left[\frac{7\pi}{4}, 2\pi\right] and r<0r < 0 in [π4,3π4]\left[\frac{\pi}{4}, \frac{3\pi}{4}\right], [5π4,7π4]\left[\frac{5\pi}{4}, \frac{7\pi}{4}\right]. When r<0r < 0, the point is plotted in the opposite direction, which traces the same petals again. The result is a four-petalled rose with petals along the xx- and yy-axes.

General curves r=cos2nθr = \cos 2n\theta:

For integer n1n \geqslant 1, cos2nθ\cos 2n\theta has period πn\dfrac{\pi}{n}. Over [0,2π][0, 2\pi] the cosine completes 2n2n full periods, producing 4n4n petals (each period gives two intervals where r>0r > 0 and two where r<0r < 0, but the negative-rr intervals trace the same petals). The petals are equally spaced, with angular spacing π2n\dfrac{\pi}{2n} between adjacent petals.

Area enclosed by r=cos2nθr = \cos 2n\theta:

Using the polar area formula and the 4n4n-fold symmetry (each petal is identical):

A=4n120π/(4n)r2dθ=2n0π/(4n)cos22nθdθ.A = 4n \cdot \frac{1}{2}\int_0^{\pi/(4n)} r^2 \, d\theta = 2n\int_0^{\pi/(4n)} \cos^2 2n\theta \, d\theta.

Using cos22nθ=1+cos4nθ2\cos^2 2n\theta = \dfrac{1 + \cos 4n\theta}{2}:

A=2n0π/(4n)1+cos4nθ2dθ=n[θ+sin4nθ4n]0π/(4n)A = 2n \int_0^{\pi/(4n)} \frac{1 + \cos 4n\theta}{2} \, d\theta = n\left[\theta + \frac{\sin 4n\theta}{4n}\right]_0^{\pi/(4n)}

=n[π4n+sinπ4n0]=nπ4n=π4.= n\left[\frac{\pi}{4n} + \frac{\sin\pi}{4n} - 0\right] = n \cdot \frac{\pi}{4n} = \frac{\pi}{4}.

Since the nn cancels, the area π4\dfrac{\pi}{4} is independent of nn. \blacksquare


Part (iii): Sketch of r=cos3θr = \cos 3\theta

The function cos3θ\cos 3\theta has period 2π3\dfrac{2\pi}{3}. Over [0,2π][0, 2\pi] it completes three full periods. In each period, cos3θ>0\cos 3\theta > 0 on an interval of length π3\dfrac{\pi}{3} and cos3θ<0\cos 3\theta < 0 on another interval of length π3\dfrac{\pi}{3}. The positive-rr intervals produce petals in one set of directions, and the negative-rr intervals produce petals in between. The result is a three-petalled rose (since 33 is odd, r=cos3θr = \cos 3\theta produces exactly 33 petals, not 66).

The petals are centred at θ=0\theta = 0, 2π3\dfrac{2\pi}{3}, and 4π3\dfrac{4\pi}{3}, each with maximum radius r=1r = 1. The petals are equally spaced at 120°120° intervals.


Topic: 矩阵指数与几何变换 Matrix Exponential and Geometric Transformations  |  Difficulty: Challenging  |  Marks: 20

5 The exponential of a square matrix A\mathbf{A} is defined to be

exp(A)=r=01r!Ar,\exp(\mathbf{A}) = \sum_{r=0}^{\infty} \frac{1}{r!} \mathbf{A}^r ,

where A0=I\mathbf{A}^0 = \mathbf{I} and I\mathbf{I} is the identity matrix.

Let

M=(0110).\mathbf{M} = \begin{pmatrix} 0 & -1 \\ 1 & 0 \end{pmatrix} .

Show that M2=I\mathbf{M}^2 = -\mathbf{I} and hence express exp(θM)\exp(\theta\mathbf{M}) as a single 2×22 \times 2 matrix, where θ\theta is a real number. Explain the geometrical significance of exp(θM)\exp(\theta\mathbf{M}).

Let

N=(0100).\mathbf{N} = \begin{pmatrix} 0 & 1 \\ 0 & 0 \end{pmatrix} .

Express similarly exp(sN)\exp(s\mathbf{N}), where ss is a real number, and explain the geometrical significance of exp(sN)\exp(s\mathbf{N}).

For which values of θ\theta does

exp(sN)exp(θM)=exp(θM)exp(sN)\exp(s\mathbf{N}) \exp(\theta\mathbf{M}) = \exp(\theta\mathbf{M}) \exp(s\mathbf{N})

for all ss? Interpret this fact geometrically.

Model Solution

Showing M2=I\mathbf{M}^2 = -\mathbf{I}:

M2=(0110)(0110)=((0)(0)+(1)(1)(0)(1)+(1)(0)(1)(0)+(0)(1)(1)(1)+(0)(0))=(1001)=I\mathbf{M}^2 = \begin{pmatrix} 0 & -1 \\ 1 & 0 \end{pmatrix}\begin{pmatrix} 0 & -1 \\ 1 & 0 \end{pmatrix} = \begin{pmatrix} (0)(0)+(-1)(1) & (0)(-1)+(-1)(0) \\ (1)(0)+(0)(1) & (1)(-1)+(0)(0) \end{pmatrix} = \begin{pmatrix} -1 & 0 \\ 0 & -1 \end{pmatrix} = -\mathbf{I}

Computing exp(θM)\exp(\theta\mathbf{M}):

From M2=I\mathbf{M}^2 = -\mathbf{I}, we can compute all powers of M\mathbf{M}:

M0=I,M1=M,M2=I,M3=M,M4=I,\mathbf{M}^0 = \mathbf{I}, \quad \mathbf{M}^1 = \mathbf{M}, \quad \mathbf{M}^2 = -\mathbf{I}, \quad \mathbf{M}^3 = -\mathbf{M}, \quad \mathbf{M}^4 = \mathbf{I}, \quad \ldots

In general, for even powers M2k=(1)kI\mathbf{M}^{2k} = (-1)^k \mathbf{I} and for odd powers M2k+1=(1)kM\mathbf{M}^{2k+1} = (-1)^k \mathbf{M}.

Substituting into the definition:

exp(θM)=r=0θrr!Mr=k=0θ2k(2k)!M2k+k=0θ2k+1(2k+1)!M2k+1\exp(\theta\mathbf{M}) = \sum_{r=0}^{\infty} \frac{\theta^r}{r!} \mathbf{M}^r = \sum_{k=0}^{\infty} \frac{\theta^{2k}}{(2k)!} \mathbf{M}^{2k} + \sum_{k=0}^{\infty} \frac{\theta^{2k+1}}{(2k+1)!} \mathbf{M}^{2k+1}

=k=0θ2k(2k)!(1)kI+k=0θ2k+1(2k+1)!(1)kM= \sum_{k=0}^{\infty} \frac{\theta^{2k}}{(2k)!} (-1)^k \mathbf{I} + \sum_{k=0}^{\infty} \frac{\theta^{2k+1}}{(2k+1)!} (-1)^k \mathbf{M}

=(k=0(1)kθ2k(2k)!)I+(k=0(1)kθ2k+1(2k+1)!)M= \left(\sum_{k=0}^{\infty} \frac{(-1)^k \theta^{2k}}{(2k)!}\right) \mathbf{I} + \left(\sum_{k=0}^{\infty} \frac{(-1)^k \theta^{2k+1}}{(2k+1)!}\right) \mathbf{M}

Recognising the Taylor series for cosine and sine:

exp(θM)=cosθI+sinθM=cosθ(1001)+sinθ(0110)\exp(\theta\mathbf{M}) = \cos\theta \cdot \mathbf{I} + \sin\theta \cdot \mathbf{M} = \cos\theta \begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix} + \sin\theta \begin{pmatrix} 0 & -1 \\ 1 & 0 \end{pmatrix}

exp(θM)=(cosθsinθsinθcosθ)\exp(\theta\mathbf{M}) = \begin{pmatrix} \cos\theta & -\sin\theta \\ \sin\theta & \cos\theta \end{pmatrix}

Geometrical significance: This is the matrix for anticlockwise rotation by angle θ\theta about the origin. The matrix M\mathbf{M} is the generator of rotations in two dimensions.


Computing exp(sN)\exp(s\mathbf{N}):

First we compute the powers of N\mathbf{N}:

N2=(0100)(0100)=(0000)=0\mathbf{N}^2 = \begin{pmatrix} 0 & 1 \\ 0 & 0 \end{pmatrix}\begin{pmatrix} 0 & 1 \\ 0 & 0 \end{pmatrix} = \begin{pmatrix} 0 & 0 \\ 0 & 0 \end{pmatrix} = \mathbf{0}

Since N2=0\mathbf{N}^2 = \mathbf{0}, all higher powers are also zero: Nr=0\mathbf{N}^r = \mathbf{0} for r2r \geqslant 2. The exponential series therefore terminates after two terms:

exp(sN)=I+sN+s22!N2+=I+sN\exp(s\mathbf{N}) = \mathbf{I} + s\mathbf{N} + \frac{s^2}{2!}\mathbf{N}^2 + \cdots = \mathbf{I} + s\mathbf{N}

exp(sN)=(1001)+(0s00)=(1s01)\exp(s\mathbf{N}) = \begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix} + \begin{pmatrix} 0 & s \\ 0 & 0 \end{pmatrix} = \begin{pmatrix} 1 & s \\ 0 & 1 \end{pmatrix}

Geometrical significance: The matrix (1s01)\begin{pmatrix} 1 & s \\ 0 & 1 \end{pmatrix} represents a horizontal shear. It maps (x,y)(x+sy,y)(x, y) \mapsto (x + sy, y): each point is shifted horizontally by ss times its yy-coordinate, while the yy-coordinate is unchanged. Horizontal lines remain fixed, and vertical lines are tilted.


Finding θ\theta for commutativity:

We compute both products.

exp(sN)exp(θM)=(1s01)(cosθsinθsinθcosθ)=(cosθ+ssinθsinθ+scosθsinθcosθ)\exp(s\mathbf{N})\exp(\theta\mathbf{M}) = \begin{pmatrix} 1 & s \\ 0 & 1 \end{pmatrix}\begin{pmatrix} \cos\theta & -\sin\theta \\ \sin\theta & \cos\theta \end{pmatrix} = \begin{pmatrix} \cos\theta + s\sin\theta & -\sin\theta + s\cos\theta \\ \sin\theta & \cos\theta \end{pmatrix}

exp(θM)exp(sN)=(cosθsinθsinθcosθ)(1s01)=(cosθscosθsinθsinθssinθ+cosθ)\exp(\theta\mathbf{M})\exp(s\mathbf{N}) = \begin{pmatrix} \cos\theta & -\sin\theta \\ \sin\theta & \cos\theta \end{pmatrix}\begin{pmatrix} 1 & s \\ 0 & 1 \end{pmatrix} = \begin{pmatrix} \cos\theta & s\cos\theta - \sin\theta \\ \sin\theta & s\sin\theta + \cos\theta \end{pmatrix}

Comparing entries:

  • Entry (1,1)(1,1): cosθ+ssinθ=cosθ    ssinθ=0\cos\theta + s\sin\theta = \cos\theta \implies s\sin\theta = 0
  • Entry (2,2)(2,2): cosθ=ssinθ+cosθ    ssinθ=0\cos\theta = s\sin\theta + \cos\theta \implies s\sin\theta = 0

For the equation to hold for all ss, we need sinθ=0\sin\theta = 0, i.e. θ=nπ\theta = n\pi for integer nn.

(We can verify the remaining entries are automatically satisfied: when sinθ=0\sin\theta = 0 and cosθ=(1)n\cos\theta = (-1)^n, entry (1,2)(1,2) gives s(1)n=s(1)ns(-1)^n = s(-1)^n on both sides, and entry (2,1)(2,1) gives 0=00 = 0.)

Geometrical interpretation: θ=nπ\theta = n\pi corresponds to rotations by multiples of 180°180° (including the identity at θ=0\theta = 0). A horizontal shear followed by a 180°180° rotation is the same as the rotation followed by the shear: the rotation flips both the shear direction and the yy-axis, and these two effects cancel. For any other rotation angle, the shear direction changes non-trivially, so the operations do not commute.


Topic: 几何 Geometry  |  Difficulty: Challenging  |  Marks: 20

6 (i) Show that four vertices of a cube, no two of which are adjacent, form the vertices of a regular tetrahedron. Hence, or otherwise, find the volume of a regular tetrahedron whose edges are of unit length.

(ii) Find the volume of a regular octahedron whose edges are of unit length.

(iii) Show that the centres of the faces of a cube form the vertices of a regular octahedron. Show that its volume is half that of the tetrahedron whose vertices are the vertices of the cube.

[A regular tetrahedron (octahedron) has four (eight) faces, all equilateral triangles.]

Model Solution

Part (i): Four non-adjacent vertices of a cube form a regular tetrahedron

Take a unit cube with vertices at all points (x,y,z)(x,y,z) where each coordinate is 00 or 11. A vertex is adjacent to another if they differ in exactly one coordinate. Four non-adjacent vertices (differing in all three coordinates pairwise) are:

A=(0,0,0),B=(1,1,0),C=(1,0,1),D=(0,1,1).A = (0,0,0), \quad B = (1,1,0), \quad C = (1,0,1), \quad D = (0,1,1).

We compute all six pairwise distances:

AB=1+1+0=2,AC=1+0+1=2,AD=0+1+1=2,|AB| = \sqrt{1+1+0} = \sqrt{2}, \quad |AC| = \sqrt{1+0+1} = \sqrt{2}, \quad |AD| = \sqrt{0+1+1} = \sqrt{2},

BC=0+1+1=2,BD=1+0+1=2,CD=1+1+0=2.|BC| = \sqrt{0+1+1} = \sqrt{2}, \quad |BD| = \sqrt{1+0+1} = \sqrt{2}, \quad |CD| = \sqrt{1+1+0} = \sqrt{2}.

All six edges have equal length 2\sqrt{2}, so these four vertices form a regular tetrahedron. \blacksquare

Volume of a unit regular tetrahedron:

Scale the tetrahedron so that each edge has length 11. The scaling factor from edge-length 2\sqrt{2} to edge-length 11 is 12\frac{1}{\sqrt{2}}, so we multiply all coordinates by 12\frac{1}{\sqrt{2}}:

A=(0,0,0),B=(12,12,0),C=(12,0,12),D=(0,12,12).A' = (0,0,0), \quad B' = \left(\frac{1}{\sqrt{2}}, \frac{1}{\sqrt{2}}, 0\right), \quad C' = \left(\frac{1}{\sqrt{2}}, 0, \frac{1}{\sqrt{2}}\right), \quad D' = \left(0, \frac{1}{\sqrt{2}}, \frac{1}{\sqrt{2}}\right).

The volume is:

V=16AB(AC×AD).V = \frac{1}{6}\left|\overrightarrow{A'B'} \cdot (\overrightarrow{A'C'} \times \overrightarrow{A'D'})\right|.

AB=(12,12,0),AC=(12,0,12),AD=(0,12,12).\overrightarrow{A'B'} = \left(\frac{1}{\sqrt{2}}, \frac{1}{\sqrt{2}}, 0\right), \quad \overrightarrow{A'C'} = \left(\frac{1}{\sqrt{2}}, 0, \frac{1}{\sqrt{2}}\right), \quad \overrightarrow{A'D'} = \left(0, \frac{1}{\sqrt{2}}, \frac{1}{\sqrt{2}}\right).

Computing the cross product:

AC×AD=ijk1201201212=(0121212)i(1212120)j+(121200)k\overrightarrow{A'C'} \times \overrightarrow{A'D'} = \begin{vmatrix} \mathbf{i} & \mathbf{j} & \mathbf{k} \\ \frac{1}{\sqrt{2}} & 0 & \frac{1}{\sqrt{2}} \\ 0 & \frac{1}{\sqrt{2}} & \frac{1}{\sqrt{2}} \end{vmatrix} = \left(0 \cdot \frac{1}{\sqrt{2}} - \frac{1}{\sqrt{2}} \cdot \frac{1}{\sqrt{2}}\right)\mathbf{i} - \left(\frac{1}{\sqrt{2}} \cdot \frac{1}{\sqrt{2}} - \frac{1}{\sqrt{2}} \cdot 0\right)\mathbf{j} + \left(\frac{1}{\sqrt{2}} \cdot \frac{1}{\sqrt{2}} - 0 \cdot 0\right)\mathbf{k}

=(12,12,12).= \left(-\frac{1}{2}, -\frac{1}{2}, \frac{1}{2}\right).

Dot product with AB\overrightarrow{A'B'}:

(12,12,0)(12,12,12)=122122+0=12.\left(\frac{1}{\sqrt{2}}, \frac{1}{\sqrt{2}}, 0\right) \cdot \left(-\frac{1}{2}, -\frac{1}{2}, \frac{1}{2}\right) = -\frac{1}{2\sqrt{2}} - \frac{1}{2\sqrt{2}} + 0 = -\frac{1}{\sqrt{2}}.

V=1612=162=212.V = \frac{1}{6} \cdot \frac{1}{\sqrt{2}} = \frac{1}{6\sqrt{2}} = \frac{\sqrt{2}}{12}.

So the volume of a regular tetrahedron with unit edges is 212\dfrac{\sqrt{2}}{12}. \blacksquare


Part (ii): Volume of a regular octahedron with unit edges

Place the octahedron with vertices at (±a,0,0)(\pm a, 0, 0), (0,±a,0)(0, \pm a, 0), (0,0,±a)(0, 0, \pm a). Adjacent vertices (e.g. (a,0,0)(a,0,0) and (0,a,0)(0,a,0)) have distance a2a\sqrt{2}. Setting a2=1a\sqrt{2} = 1 gives a=12a = \dfrac{1}{\sqrt{2}}.

The octahedron consists of 88 congruent tetrahedra, each with one vertex at the origin and three vertices on the positive or negative axes. Consider the tetrahedron with vertices O=(0,0,0)O = (0,0,0), P1=(a,0,0)P_1 = (a,0,0), P2=(0,a,0)P_2 = (0,a,0), P3=(0,0,a)P_3 = (0,0,a):

Vtet=16det(a000a000a)=a36.V_{\text{tet}} = \frac{1}{6} \left|\det\begin{pmatrix} a & 0 & 0 \\ 0 & a & 0 \\ 0 & 0 & a \end{pmatrix}\right| = \frac{a^3}{6}.

Total volume:

Voct=8a36=4a33=43122=462=232=23.V_{\text{oct}} = 8 \cdot \frac{a^3}{6} = \frac{4a^3}{3} = \frac{4}{3} \cdot \frac{1}{2\sqrt{2}} = \frac{4}{6\sqrt{2}} = \frac{2}{3\sqrt{2}} = \frac{\sqrt{2}}{3}.

So the volume of a regular octahedron with unit edges is 23\dfrac{\sqrt{2}}{3}. \blacksquare


Part (iii): Face centres of a cube form a regular octahedron

Take a cube of side aa with vertices at all (x,y,z)(x,y,z) where each coordinate is 00 or aa. The six face centres are:

(a2,a2,0),(a2,a2,a),(a2,0,a2),(a2,a,a2),(0,a2,a2),(a,a2,a2).\left(\frac{a}{2}, \frac{a}{2}, 0\right), \quad \left(\frac{a}{2}, \frac{a}{2}, a\right), \quad \left(\frac{a}{2}, 0, \frac{a}{2}\right), \quad \left(\frac{a}{2}, a, \frac{a}{2}\right), \quad \left(0, \frac{a}{2}, \frac{a}{2}\right), \quad \left(a, \frac{a}{2}, \frac{a}{2}\right).

Shifting the origin to the centre of the cube (a2,a2,a2)\left(\frac{a}{2}, \frac{a}{2}, \frac{a}{2}\right), these become:

(0,0,a2),(0,0,a2),(0,a2,0),(0,a2,0),(a2,0,0),(a2,0,0).\left(0, 0, -\frac{a}{2}\right), \quad \left(0, 0, \frac{a}{2}\right), \quad \left(0, -\frac{a}{2}, 0\right), \quad \left(0, \frac{a}{2}, 0\right), \quad \left(-\frac{a}{2}, 0, 0\right), \quad \left(\frac{a}{2}, 0, 0\right).

These are the six vertices (±a2,0,0)(\pm \frac{a}{2}, 0, 0), (0,±a2,0)(0, \pm \frac{a}{2}, 0), (0,0,±a2)(0, 0, \pm \frac{a}{2}), which is a regular octahedron (as shown in Part (ii) with edge length a2\frac{a}{\sqrt{2}}). \blacksquare

Comparing volumes:

The tetrahedron from Part (i) inscribed in this cube (side aa) has edge length a2a\sqrt{2} and volume:

Vtet=(a2)3621212=a321222=a33.V_{\text{tet}} = \left(\frac{a}{\sqrt{2}}\right)^3 \cdot \frac{6\sqrt{2}}{1} \cdot \frac{\sqrt{2}}{12} = a^3 \cdot \frac{\sqrt{2}}{12} \cdot 2\sqrt{2} = \frac{a^3}{3}.

(Using the unit-edge formula scaled by (a2)3(a\sqrt{2})^3: Vtet=212(a2)3=21222a3=a33V_{\text{tet}} = \frac{\sqrt{2}}{12}(a\sqrt{2})^3 = \frac{\sqrt{2}}{12} \cdot 2\sqrt{2}\,a^3 = \frac{a^3}{3}.)

The octahedron formed by face centres has edge length a2\frac{a}{\sqrt{2}} and volume:

Voct=23(a2)3=23a322=a36.V_{\text{oct}} = \frac{\sqrt{2}}{3} \cdot \left(\frac{a}{\sqrt{2}}\right)^3 = \frac{\sqrt{2}}{3} \cdot \frac{a^3}{2\sqrt{2}} = \frac{a^3}{6}.

Therefore Voct=12VtetV_{\text{oct}} = \dfrac{1}{2}\,V_{\text{tet}}. The volume of the octahedron is half that of the tetrahedron. \blacksquare


Topic: 微积分 Calculus  |  Difficulty: Hard  |  Marks: 20

7 Sketch the graph of f(s)=es(s3)+3f(s) = e^s(s - 3) + 3 for 0s<0 \leqslant s < \infty. Taking e2.7e \approx 2.7, find the smallest positive integer, mm, such that f(m)>0f(m) > 0.

Now let

b(x)=x3ex/T1b(x) = \frac{x^3}{e^{x/T} - 1}

where TT is a positive constant. Show that b(x)b(x) has a single turning point in 0<x<0 < x < \infty. By considering the behaviour for small xx and for large xx, sketch b(x)b(x) for 0x<0 \leqslant x < \infty.

Let

0b(x)dx=B,\int_0^\infty b(x) \, dx = B,

which may be assumed to be finite. Show that B=KTnB = KT^n where KK is a constant, and nn is an integer which you should determine.

Given that B20Tmb(x)dxB \approx 2 \int\limits_0^{Tm} b(x) dx, use your graph of b(x)b(x) to find a rough estimate for KK.

Model Solution

Part (i): Sketch of f(s)=es(s3)+3f(s) = e^s(s-3) + 3 and finding mm

Derivative:

f(s)=es(s3)+es=es(s2).f'(s) = e^s(s-3) + e^s = e^s(s-2).

Since es>0e^s > 0 always: f(s)<0f'(s) < 0 for s<2s < 2 and f(s)>0f'(s) > 0 for s>2s > 2. So ff has a minimum at s=2s = 2.

Key values:

f(0)=e0(03)+3=3+3=0.f(0) = e^0(0-3) + 3 = -3 + 3 = 0.

f(2)=e2(23)+3=3e237.29=4.29.f(2) = e^2(2-3) + 3 = 3 - e^2 \approx 3 - 7.29 = -4.29.

f(3)=e3(33)+3=3.f(3) = e^3(3-3) + 3 = 3.

Behaviour: ff starts at 00, decreases to a minimum of 3e24.293 - e^2 \approx -4.29 at s=2s = 2, then increases. At s=3s = 3 we have f(3)=3>0f(3) = 3 > 0. For large ss, f(s)sesf(s) \approx se^s \to \infty.

Sketch: The curve starts at the origin, dips down to a minimum near s=2s = 2, then rises steeply through f(3)=3f(3) = 3 and continues to infinity.

Finding mm: We need the smallest positive integer with f(m)>0f(m) > 0:

  • f(1)=e1(13)+3=2e+35.4+3=2.4<0f(1) = e^1(1-3) + 3 = -2e + 3 \approx -5.4 + 3 = -2.4 < 0
  • f(2)=3e24.3<0f(2) = 3 - e^2 \approx -4.3 < 0
  • f(3)=3>0f(3) = 3 > 0

So m=3m = 3.


Part (ii): b(x)b(x) has a single turning point

b(x)=x3ex/T1.b(x) = \frac{x^3}{e^{x/T} - 1}.

Using the quotient rule:

b(x)=3x2(ex/T1)x31Tex/T(ex/T1)2=x2[3(ex/T1)xTex/T](ex/T1)2.b'(x) = \frac{3x^2(e^{x/T} - 1) - x^3 \cdot \frac{1}{T}e^{x/T}}{(e^{x/T} - 1)^2} = \frac{x^2\left[3(e^{x/T} - 1) - \frac{x}{T}e^{x/T}\right]}{(e^{x/T} - 1)^2}.

Setting the numerator to zero (for x>0x > 0):

3(ex/T1)xTex/T=0ex/T(3xT)=3.()3(e^{x/T} - 1) - \frac{x}{T}e^{x/T} = 0 \qquad \Longleftrightarrow \qquad e^{x/T}\left(3 - \frac{x}{T}\right) = 3. \qquad (\star)

Let u=x/T>0u = x/T > 0 and define h(u)=eu(3u)3h(u) = e^u(3 - u) - 3. Then:

h(u)=eu(3u)eu=eu(2u).h'(u) = e^u(3-u) - e^u = e^u(2 - u).

So h(u)>0h'(u) > 0 for u<2u < 2 and h(u)<0h'(u) < 0 for u>2u > 2: hh increases on (0,2)(0, 2) and decreases on (2,)(2, \infty).

At the boundaries: h(0)=33=0h(0) = 3 - 3 = 0 and h(3)=e303=3<0h(3) = e^3 \cdot 0 - 3 = -3 < 0. Also h(2)=e213=e23>0h(2) = e^2 \cdot 1 - 3 = e^2 - 3 > 0.

Since h(0)=0h(0) = 0, hh increases from 00 to a positive maximum at u=2u = 2, then decreases to negative values. It crosses zero exactly once for u>0u > 0 (at some u0(2,3)u_0 \in (2, 3)). This gives a unique turning point of b(x)b(x) at x0=Tu0x_0 = Tu_0. \blacksquare

Behaviour of b(x)b(x):

  • Near x=0x = 0: Using ex/T1xTe^{x/T} - 1 \approx \frac{x}{T}: b(x)x3x/T=Tx20\quad b(x) \approx \frac{x^3}{x/T} = Tx^2 \to 0.
  • At x=0x = 0: By L’Hopital (or the above), b(0)=0b(0) = 0.
  • For large xx: ex/Te^{x/T} grows exponentially, so b(x)x3ex/T0b(x) \approx x^3 e^{-x/T} \to 0.
  • At the turning point x0x_0: bb attains its maximum.

The curve starts at the origin, rises to a single maximum, then decays back toward zero.


Part (iii): Show B=KT4B = KT^4

Substituting u=x/Tu = x/T (so x=Tux = Tu, dx=Tdudx = T\,du):

B=0x3ex/T1dx=0T3u3eu1Tdu=T40u3eu1du.B = \int_0^\infty \frac{x^3}{e^{x/T} - 1}\,dx = \int_0^\infty \frac{T^3 u^3}{e^u - 1}\,T\,du = T^4 \int_0^\infty \frac{u^3}{e^u - 1}\,du.

The integral K=0u3eu1du\displaystyle K = \int_0^\infty \frac{u^3}{e^u - 1}\,du is a dimensionless constant (it equals π415\dfrac{\pi^4}{15}, though this is not required here). Therefore B=KT4B = KT^4, so n=4n = 4. \blacksquare


Part (iv): Rough estimate for KK

Using m=3m = 3 and B20Tmb(x)dx=203Tb(x)dxB \approx 2\int_0^{Tm} b(x)\,dx = 2\int_0^{3T} b(x)\,dx, substitute u=x/Tu = x/T:

KT42T403u3eu1duK203u3eu1du.KT^4 \approx 2T^4 \int_0^3 \frac{u^3}{e^u - 1}\,du \qquad \Longrightarrow \qquad K \approx 2\int_0^3 \frac{u^3}{e^u - 1}\,du.

We estimate the integral using Simpson’s rule with h=0.5h = 0.5 and values at u=0,0.5,1,1.5,2,2.5,3u = 0, 0.5, 1, 1.5, 2, 2.5, 3:

uu000.50.5111.51.5222.52.533
f(u)=u3eu1f(u) = \frac{u^3}{e^u - 1}110.2380.2380.5820.5820.7530.7530.7670.7670.6920.6920.5880.588

(Near u=0u = 0: u3eu1u2\frac{u^3}{e^u-1} \approx u^2, so f(0)=0f(0) = 0; at u=0.5u = 0.5: 0.125e0.51=0.1250.6487=0.193\frac{0.125}{e^{0.5}-1} = \frac{0.125}{0.6487} = 0.193.)

Let me recompute more carefully. At u=0u = 0: f(0)=0f(0) = 0 (by L’Hopital, limu0u3eu1=0\lim_{u\to 0} \frac{u^3}{e^u-1} = 0).

uu000.50.5111.51.5222.52.533
f(u)f(u)000.1930.1930.5820.5820.7530.7530.7670.7670.6920.6920.5880.588

Simpson’s rule with h=0.5h = 0.5:

03f(u)du0.53[f(0)+4f(0.5)+2f(1)+4f(1.5)+2f(2)+4f(2.5)+f(3)]\int_0^3 f(u)\,du \approx \frac{0.5}{3}\left[f(0) + 4f(0.5) + 2f(1) + 4f(1.5) + 2f(2) + 4f(2.5) + f(3)\right]

=16[0+4(0.193)+2(0.582)+4(0.753)+2(0.767)+4(0.692)+0.588]= \frac{1}{6}\left[0 + 4(0.193) + 2(0.582) + 4(0.753) + 2(0.767) + 4(0.692) + 0.588\right]

=16[0+0.772+1.164+3.012+1.534+2.768+0.588]=9.83861.64.= \frac{1}{6}\left[0 + 0.772 + 1.164 + 3.012 + 1.534 + 2.768 + 0.588\right] = \frac{9.838}{6} \approx 1.64.

Therefore K2×1.64=3.3K \approx 2 \times 1.64 = 3.3.

(For comparison, the exact value is K=π4156.49K = \frac{\pi^4}{15} \approx 6.49; our approximation underestimates because we are missing the tail integral from 33 to \infty, which contributes about half the total.)


Topic: 向量与几何 Vectors and Geometry  |  Difficulty: Challenging  |  Marks: 20

8 (i) Show that the line r=b+λm\mathbf{r} = \mathbf{b} + \lambda\mathbf{m}, where m\mathbf{m} is a unit vector, intersects the sphere r.r=a2\mathbf{r}.\mathbf{r} = a^2 at two points if

a2>b.b(b.m)2.a^2 > \mathbf{b}.\mathbf{b} - (\mathbf{b}.\mathbf{m})^2.

Write down the corresponding condition for there to be precisely one point of intersection. If this point has position vector p\mathbf{p}, show that m.p=0\mathbf{m}.\mathbf{p} = 0.

(ii) Now consider a second sphere of radius aa and a plane perpendicular to a unit vector n\mathbf{n}. The centre of the sphere has position vector d\mathbf{d} and the minimum distance from the origin to the plane is ll. What is the condition for the plane to be tangential to this second sphere?

(iii) Show that the first and second spheres intersect at right angles (i.e. the two radii to each point of intersection are perpendicular) if

d.d=2a2.\mathbf{d}.\mathbf{d} = 2a^2.

Model Solution

Part (i): Line-sphere intersection

The line is r=b+λm\mathbf{r} = \mathbf{b} + \lambda\mathbf{m} where m=1|\mathbf{m}| = 1. The sphere is rr=a2\mathbf{r} \cdot \mathbf{r} = a^2. Substituting:

(b+λm)(b+λm)=a2(\mathbf{b} + \lambda\mathbf{m}) \cdot (\mathbf{b} + \lambda\mathbf{m}) = a^2

bb+2λ(bm)+λ2(mm)=a2.\mathbf{b} \cdot \mathbf{b} + 2\lambda(\mathbf{b} \cdot \mathbf{m}) + \lambda^2(\mathbf{m} \cdot \mathbf{m}) = a^2.

Since m=1|\mathbf{m}| = 1:

λ2+2(bm)λ+(bba2)=0.()\lambda^2 + 2(\mathbf{b} \cdot \mathbf{m})\lambda + (\mathbf{b} \cdot \mathbf{b} - a^2) = 0. \qquad (\star)

This is a quadratic in λ\lambda with discriminant:

Δ=4(bm)24(bba2)=4[a2bb+(bm)2].\Delta = 4(\mathbf{b} \cdot \mathbf{m})^2 - 4(\mathbf{b} \cdot \mathbf{b} - a^2) = 4\left[a^2 - \mathbf{b} \cdot \mathbf{b} + (\mathbf{b} \cdot \mathbf{m})^2\right].

Two distinct real roots (two intersection points) exist if and only if Δ>0\Delta > 0:

a2>bb(bm)2.a^2 > \mathbf{b} \cdot \mathbf{b} - (\mathbf{b} \cdot \mathbf{m})^2. \qquad \blacksquare

Condition for exactly one intersection (tangency):

a2=bb(bm)2.a^2 = \mathbf{b} \cdot \mathbf{b} - (\mathbf{b} \cdot \mathbf{m})^2.

When Δ=0\Delta = 0, the unique solution of ()(\star) is λ=bm\lambda = -\mathbf{b} \cdot \mathbf{m}, giving:

p=b(bm)m.\mathbf{p} = \mathbf{b} - (\mathbf{b} \cdot \mathbf{m})\mathbf{m}.

We verify mp=0\mathbf{m} \cdot \mathbf{p} = 0:

mp=mb(bm)(mm)=bmbm=0.\mathbf{m} \cdot \mathbf{p} = \mathbf{m} \cdot \mathbf{b} - (\mathbf{b} \cdot \mathbf{m})(\mathbf{m} \cdot \mathbf{m}) = \mathbf{b} \cdot \mathbf{m} - \mathbf{b} \cdot \mathbf{m} = 0. \qquad \blacksquare

(Geometrically, p\mathbf{p} is the foot of the perpendicular from the centre of the sphere to the line, and tangency occurs when this perpendicular has length aa.)


Part (ii): Tangency condition for a plane and a sphere

The plane has equation rn=l\mathbf{r} \cdot \mathbf{n} = l (where n=1|\mathbf{n}| = 1 and l>0l > 0 is the minimum distance from the origin to the plane). The second sphere has centre d\mathbf{d} and radius aa, so its equation is:

(rd)(rd)=a2.(\mathbf{r} - \mathbf{d}) \cdot (\mathbf{r} - \mathbf{d}) = a^2.

The distance from the centre d\mathbf{d} to the plane rn=l\mathbf{r} \cdot \mathbf{n} = l is dnl|\mathbf{d} \cdot \mathbf{n} - l|. The plane is tangential to the sphere if and only if this distance equals the radius:

(dnl)2=a2dnl=a.(\mathbf{d} \cdot \mathbf{n} - l)^2 = a^2 \qquad \Longleftrightarrow \qquad |\mathbf{d} \cdot \mathbf{n} - l| = a. \qquad \blacksquare


Part (iii): Two spheres intersect at right angles

The first sphere (centre OO, radius aa): rr=a2\mathbf{r} \cdot \mathbf{r} = a^2.

The second sphere (centre d\mathbf{d}, radius aa): (rd)(rd)=a2(\mathbf{r} - \mathbf{d}) \cdot (\mathbf{r} - \mathbf{d}) = a^2.

Finding the intersection: Expanding the second equation:

rr2rd+dd=a2.\mathbf{r} \cdot \mathbf{r} - 2\mathbf{r} \cdot \mathbf{d} + \mathbf{d} \cdot \mathbf{d} = a^2.

Substituting rr=a2\mathbf{r} \cdot \mathbf{r} = a^2 from the first sphere:

a22rd+dd=a2rd=dd2.()a^2 - 2\mathbf{r} \cdot \mathbf{d} + \mathbf{d} \cdot \mathbf{d} = a^2 \qquad \Longrightarrow \qquad \mathbf{r} \cdot \mathbf{d} = \frac{\mathbf{d} \cdot \mathbf{d}}{2}. \qquad (\star)

So the intersection lies in the plane rd=d22\mathbf{r} \cdot \mathbf{d} = \frac{|\mathbf{d}|^2}{2}.

Orthogonality condition: Two spheres intersect at right angles if, at every point of intersection r\mathbf{r}, the radius vectors from the two centres are perpendicular. The radius vectors are r\mathbf{r} (from OO) and rd\mathbf{r} - \mathbf{d} (from d\mathbf{d}). We need:

r(rd)=0rrrd=0.\mathbf{r} \cdot (\mathbf{r} - \mathbf{d}) = 0 \qquad \Longleftrightarrow \qquad \mathbf{r} \cdot \mathbf{r} - \mathbf{r} \cdot \mathbf{d} = 0.

Using rr=a2\mathbf{r} \cdot \mathbf{r} = a^2 and ()(\star):

a2dd2=0dd=2a2.a^2 - \frac{\mathbf{d} \cdot \mathbf{d}}{2} = 0 \qquad \Longleftrightarrow \qquad \mathbf{d} \cdot \mathbf{d} = 2a^2. \qquad \blacksquare