跳转到内容

STEP2 2004 -- Pure Mathematics

此内容尚不支持你的语言。

STEP2 2004 — Section A (Pure Mathematics)

Section titled “STEP2 2004 — Section A (Pure Mathematics)”

Exam: STEP2  |  Year: 2004  |  Questions: Q1—Q8  |  Total marks per question: 20

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

QTopicDifficultyKey Techniques
1代数/方程求解 Algebra/Equation SolvingRoutine平方去根号,因式分解,验根
2代数/不等式 Algebra/InequalitiesStandard配方,绝对值处理,韦达定理
3微积分/曲线分析 Calculus/Curve AnalysisStandard乘积求导法则,二阶导数判别法,曲线作图
4微积分/最优化 Calculus/OptimisationChallenging三角函数求导,约束优化,消元法
5积分/函数方程 Integration/Functional EquationsStandard分部积分,三角恒等式展开,待定系数法
6向量 VectorsChallenging向量投影,点积运算,正交分解
7数值方法 Numerical MethodsStandard二分法,函数符号判断,三角函数近似计算
8微分方程 Differential EquationsChallenging分离变量法,三角换元,不等式证明,图像分析

Topic: 代数/方程求解 Algebra/Equation Solving  |  Difficulty: Routine  |  Marks: 20

1 Find all real values of xx that satisfy:

(i) 3x2+1+x2x1=0\sqrt{3x^2 + 1} + \sqrt{x} - 2x - 1 = 0;

(ii) 3x2+12x+x1=0\sqrt{3x^2 + 1} - 2\sqrt{x} + x - 1 = 0;

(iii) 3x2+12xx+1=0\sqrt{3x^2 + 1} - 2\sqrt{x} - x + 1 = 0.

Hint

Q1(i) Put the terms with radicals on one side and the terms without on the other and square. Repeat this strategy (S) and the equation x^4 - 6x^3 + 9x^2 - 4x = 0 () will be obtained. The roots of () are x = 0, 1, 4.

Squaring may introduce spurious roots, so these numbers must be checked to see that they are roots of the original equation. In fact, they are.

(ii) Application of S again leads to (*). Checking shows that x = 0, x = 1 are roots of the second equation but that x = 4 is not.

(iii) Again application of S leads to (*). Checking shows that x = 1, x = 4 are roots of the third equation but that x = 0 is not.

Model Solution

Part (i)

Rearrange the equation as

3x2+1+x=2x+1\sqrt{3x^2 + 1} + \sqrt{x} = 2x + 1

The domain requires x0x \geq 0 (for x\sqrt{x} to exist). For x0x \geq 0, both sides are non-negative, so squaring preserves equivalence.

Square both sides:

3x2+1+2x(3x2+1)+x=4x2+4x+13x^2 + 1 + 2\sqrt{x(3x^2 + 1)} + x = 4x^2 + 4x + 1

Rearrange to isolate the remaining radical:

2x(3x2+1)=4x2+4x+13x21x=x2+3x2\sqrt{x(3x^2 + 1)} = 4x^2 + 4x + 1 - 3x^2 - 1 - x = x^2 + 3x

So

2x(3x2+1)=x(x+3)(*)2\sqrt{x(3x^2 + 1)} = x(x + 3) \qquad \text{(*)}

For x0x \geq 0, both sides of (*) are non-negative. Square again:

4x(3x2+1)=x2(x+3)24x(3x^2 + 1) = x^2(x + 3)^2

If x=0x = 0: both sides equal 00, so x=0x = 0 satisfies this equation. Check in the original: 1+001=0\sqrt{1} + 0 - 0 - 1 = 0 ✓.

If x>0x > 0, divide both sides by xx:

4(3x2+1)=x(x+3)24(3x^2 + 1) = x(x + 3)^2

12x2+4=x3+6x2+9x12x^2 + 4 = x^3 + 6x^2 + 9x

x36x2+9x4=0x^3 - 6x^2 + 9x - 4 = 0

Test x=1x = 1: 16+94=01 - 6 + 9 - 4 = 0 ✓. Factor out (x1)(x - 1):

x36x2+9x4=(x1)(x25x+4)=(x1)(x1)(x4)=(x1)2(x4)x^3 - 6x^2 + 9x - 4 = (x - 1)(x^2 - 5x + 4) = (x - 1)(x - 1)(x - 4) = (x - 1)^2(x - 4)

So the candidates are x=1x = 1 and x=4x = 4.

Check in the original equation 3x2+1+x2x1=0\sqrt{3x^2 + 1} + \sqrt{x} - 2x - 1 = 0:

  • x=1x = 1: 4+121=2+13=0\sqrt{4} + 1 - 2 - 1 = 2 + 1 - 3 = 0
  • x=4x = 4: 49+281=7+29=0\sqrt{49} + 2 - 8 - 1 = 7 + 2 - 9 = 0

Solutions: x{0,1,4}x \in \{0,\, 1,\, 4\}.

Part (ii)

Rearrange as

3x2+12x=1x\sqrt{3x^2 + 1} - 2\sqrt{x} = 1 - x

Square both sides:

3x2+14x(3x2+1)+4x=12x+x23x^2 + 1 - 4\sqrt{x(3x^2 + 1)} + 4x = 1 - 2x + x^2

Rearrange:

4x(3x2+1)=12x+x23x214x=2x26x-4\sqrt{x(3x^2 + 1)} = 1 - 2x + x^2 - 3x^2 - 1 - 4x = -2x^2 - 6x

4x(3x2+1)=2x2+6x=2x(x+3)4\sqrt{x(3x^2 + 1)} = 2x^2 + 6x = 2x(x + 3)

2x(3x2+1)=x(x+3)2\sqrt{x(3x^2 + 1)} = x(x + 3)

This is identical to equation (*) from part (i). The same polynomial x36x2+9x4=(x1)2(x4)=0x^3 - 6x^2 + 9x - 4 = (x - 1)^2(x - 4) = 0 is obtained, giving candidates x=0,1,4x = 0, 1, 4.

Check in the original equation 3x2+12x+x1=0\sqrt{3x^2 + 1} - 2\sqrt{x} + x - 1 = 0:

  • x=0x = 0: 10+01=01 - 0 + 0 - 1 = 0
  • x=1x = 1: 22+11=02 - 2 + 1 - 1 = 0
  • x=4x = 4: 74+41=607 - 4 + 4 - 1 = 6 \neq 0

Solutions: x{0,1}x \in \{0,\, 1\}.

Part (iii)

Rearrange as

3x2+12x=x1\sqrt{3x^2 + 1} - 2\sqrt{x} = x - 1

Square both sides:

3x2+14x(3x2+1)+4x=x22x+13x^2 + 1 - 4\sqrt{x(3x^2 + 1)} + 4x = x^2 - 2x + 1

Rearrange:

4x(3x2+1)=x22x+13x214x=2x26x-4\sqrt{x(3x^2 + 1)} = x^2 - 2x + 1 - 3x^2 - 1 - 4x = -2x^2 - 6x

2x(3x2+1)=x(x+3)2\sqrt{x(3x^2 + 1)} = x(x + 3)

Again the same equation (*). Same candidates x=0,1,4x = 0, 1, 4.

Check in the original equation 3x2+12xx+1=0\sqrt{3x^2 + 1} - 2\sqrt{x} - x + 1 = 0:

  • x=0x = 0: 100+1=201 - 0 - 0 + 1 = 2 \neq 0
  • x=1x = 1: 221+1=02 - 2 - 1 + 1 = 0
  • x=4x = 4: 744+1=07 - 4 - 4 + 1 = 0

Solutions: x{1,4}x \in \{1,\, 4\}.


Topic: 代数/不等式 Algebra/Inequalities  |  Difficulty: Standard  |  Marks: 20

2 Prove that, if α<22|\alpha| < 2\sqrt{2}, then there is no value of xx for which

x2αx+2<0.(*)x^2 - \alpha|x| + 2 < 0 . \qquad \text{(*)}

Find the solution set of (*) for α=3\alpha = 3.

For α>22\alpha > 2\sqrt{2}, the sum of the lengths of the intervals in which xx satisfies (*) is denoted by SS. Find SS in terms of α\alpha and deduce that S<2αS < 2\alpha.

Sketch the graph of SS against α\alpha.

Hint

Q2 Write Q = x^2 - α|x| + 2 = [|x| - α/2]^2 + 2 - α^2/4.

Thus α < 2√2 => 2 - α^2/4 > 0 => Q > 0 for all x.

It is therefore unnecessary to consider x > 0 and x < 0 separately and even more unnecessary to use calculus methods.

  • if α = 3 then Q = (|x| - 1)(|x| - 2), in which case the solution set of Q < 0 is {x : -2 < x < -1} ∪ {x : 1 < x < 2}.

  • The solutions in x of the equation Q = 0 are of the form -x2, -x1, x1, x2, where 0 < x1 < x2, so that S = 2(x2 - x1). Use of the identity x2 - x1 = sqrt((x2 + x1)^2 - 4x1x2) will lead immediately to S = 2sqrt(α^2 - 8). Thus S < 2sqrt(α^2) = 2α.

  • The graph of S as a function of α is that part of the hyperbola 4α^2 - S^2 = 32 which is in the first quadrant. A sketch of this graph should, therefore, leave the other quadrants empty. It should also show the curve starting at the point (2√2, 0) and asymptotically approaching the line S = 2α.

Model Solution

Proof that no solution exists when α<22|\alpha| < 2\sqrt{2}:

Let u=x0u = |x| \geq 0. Since x2=x2=u2x^2 = |x|^2 = u^2, the inequality becomes

u2αu+2<0u^2 - \alpha u + 2 < 0

Complete the square:

(uα2)2+2α24<0\left(u - \frac{\alpha}{2}\right)^2 + 2 - \frac{\alpha^2}{4} < 0

This requires

(uα2)2<α242\left(u - \frac{\alpha}{2}\right)^2 < \frac{\alpha^2}{4} - 2

If α<22|\alpha| < 2\sqrt{2}, then α2<8\alpha^2 < 8, so α24<2\frac{\alpha^2}{4} < 2, hence α242<0\frac{\alpha^2}{4} - 2 < 0.

Since (uα2)20\left(u - \frac{\alpha}{2}\right)^2 \geq 0 for all real uu, we would need a non-negative quantity to be strictly less than a negative quantity, which is impossible.

Therefore, when α<22|\alpha| < 2\sqrt{2}, there is no value of xx for which x2αx+2<0x^2 - \alpha|x| + 2 < 0.

Solution set for α=3\alpha = 3:

x23x+2=(x1)(x2)x^2 - 3|x| + 2 = (|x| - 1)(|x| - 2)

The inequality (x1)(x2)<0(|x| - 1)(|x| - 2) < 0 holds when one factor is positive and the other is negative, i.e., 1<x<21 < |x| < 2.

This gives 2<x<1-2 < x < -1 or 1<x<21 < x < 2.

Solution set: {x:2<x<1}{x:1<x<2}\{x : -2 < x < -1\} \cup \{x : 1 < x < 2\}.

Finding SS for α>22\alpha > 2\sqrt{2}:

The equation x2αx+2=0x^2 - \alpha|x| + 2 = 0 is equivalent to x2αx+2=0|x|^2 - \alpha|x| + 2 = 0. By the quadratic formula:

x=α±α282|x| = \frac{\alpha \pm \sqrt{\alpha^2 - 8}}{2}

Since α>22\alpha > 2\sqrt{2}, we have α2>8\alpha^2 > 8, so α28\sqrt{\alpha^2 - 8} is real and positive. Define:

x1=αα282,x2=α+α282x_1 = \frac{\alpha - \sqrt{\alpha^2 - 8}}{2}, \qquad x_2 = \frac{\alpha + \sqrt{\alpha^2 - 8}}{2}

Since α28<α\sqrt{\alpha^2 - 8} < \alpha (as 8>08 > 0), both roots are positive with 0<x1<x20 < x_1 < x_2.

The inequality x2αx+2<0x^2 - \alpha|x| + 2 < 0 holds when x1<x<x2x_1 < |x| < x_2, i.e.,

x(x2,x1)(x1,x2)x \in (-x_2,\, -x_1) \cup (x_1,\, x_2)

The total length of these two intervals is

S=(x2x1)+(x2x1)=2(x2x1)S = (x_2 - x_1) + (x_2 - x_1) = 2(x_2 - x_1)

By Vieta’s formulas applied to x2αx+2=0|x|^2 - \alpha|x| + 2 = 0:

x1+x2=α,x1x2=2x_1 + x_2 = \alpha, \qquad x_1 x_2 = 2

Therefore

(x2x1)2=(x1+x2)24x1x2=α28(x_2 - x_1)^2 = (x_1 + x_2)^2 - 4x_1 x_2 = \alpha^2 - 8

x2x1=α28x_2 - x_1 = \sqrt{\alpha^2 - 8}

S=2α28S = 2\sqrt{\alpha^2 - 8}

Showing S<2αS < 2\alpha:

Since α28<α2\alpha^2 - 8 < \alpha^2 (as 8>08 > 0), taking positive square roots:

α28<α2=α\sqrt{\alpha^2 - 8} < \sqrt{\alpha^2} = \alpha

(since α>22>0\alpha > 2\sqrt{2} > 0). Therefore

S=2α28<2αS = 2\sqrt{\alpha^2 - 8} < 2\alpha

Sketch of SS against α\alpha:

From S=2α28S = 2\sqrt{\alpha^2 - 8}, squaring gives S2=4α232S^2 = 4\alpha^2 - 32, i.e.,

4α2S2=324\alpha^2 - S^2 = 32

This is a hyperbola in the (α,S)(\alpha, S)-plane with asymptotes S=±2αS = \pm 2\alpha. For α>22\alpha > 2\sqrt{2} and S>0S > 0, we take the branch in the first quadrant:

  • The curve starts at the point (22,0)(2\sqrt{2},\, 0).
  • As α\alpha \to \infty, the curve approaches the line S=2αS = 2\alpha from below.
  • The curve lies entirely below the line S=2αS = 2\alpha for all α>22\alpha > 2\sqrt{2}, consistent with the result S<2αS < 2\alpha.

A sketch should show only this first-quadrant branch, with the asymptote S=2αS = 2\alpha drawn as a dashed line.


Topic: 微积分/曲线分析 Calculus/Curve Analysis  |  Difficulty: Standard  |  Marks: 20

3 The curve CC has equation

y=x(x+1)(x2)4.y = x(x + 1)(x - 2)^4 .

Show that the gradient of CC is (x2)3(6x2+x2)(x - 2)^3 (6x^2 + x - 2) and find the coordinates of all the stationary points. Determine the nature of each stationary point and sketch CC.

In separate diagrams draw sketches of the curves whose equations are:

(i) y2=x(x+1)(x2)4y^2 = x(x + 1)(x - 2)^4;

(ii) y=x2(x2+1)(x22)4y = x^2(x^2 + 1)(x^2 - 2)^4.

In each case, you should pay particular attention to the points where the curve meets the xx axis.

Hint

Q3 The obtaining of dy/dx in the form required is a routine exercise in differentiation followed by some algebra.

Setting dy/dx = 0 shows that there are stationary points where x = -2/3, 1/2, 2. Moreover d^2y/dx^2 = (x - 2)^3(12x + 1) + a term which is necessarily zero when x = -2/3, 1/2, 2. Thus d^2y/dx^2 is positive when x = -2/3 and negative when x = 1/2, so that C has a minimum at (-2/3, -8192/729) and a maximum at (1/2, 243/64). (Note that it is unnecessary to determine a simplified version of d^2y/dx^2 before inserting values of x.)

The argument d^2y/dx^2 = 0 at x = 2 => C has a point of inflexion at (2, 0) is false. In fact, in the neighbourhood of this point, y ≈ 6(x - 2)^4, so that it is obvious that C has a minimum there.

The sketch of C must have correct overall shape, location and orientation, and also show correct forms at (0,0), (2,0) and at ∞.

(i) This sketch may be deduced from that of C. It has symmetry about the x-axis and no part of it appears in the region -1 < x < 0.

(ii) This sketch may also be deduced from that of C. It has symmetry about the y-axis and no part of it appears in the region y < 0.

Model Solution

Finding dydx\dfrac{dy}{dx}:

Let y=x(x+1)(x2)4y = x(x+1)(x-2)^4. We use the product rule with u=x(x+1)=x2+xu = x(x+1) = x^2 + x and v=(x2)4v = (x-2)^4:

dudx=2x+1,dvdx=4(x2)3\frac{du}{dx} = 2x + 1, \qquad \frac{dv}{dx} = 4(x-2)^3

dydx=(2x+1)(x2)4+(x2+x)4(x2)3\frac{dy}{dx} = (2x+1)(x-2)^4 + (x^2+x) \cdot 4(x-2)^3

Factor out (x2)3(x-2)^3:

dydx=(x2)3[(2x+1)(x2)+4(x2+x)]\frac{dy}{dx} = (x-2)^3 \left[ (2x+1)(x-2) + 4(x^2+x) \right]

Expand the bracket:

(2x+1)(x2)+4x2+4x=2x24x+x2+4x2+4x=6x2+x2(2x+1)(x-2) + 4x^2 + 4x = 2x^2 - 4x + x - 2 + 4x^2 + 4x = 6x^2 + x - 2

Therefore

dydx=(x2)3(6x2+x2)\frac{dy}{dx} = (x-2)^3(6x^2 + x - 2) \qquad \checkmark

Stationary points:

Set dydx=0\dfrac{dy}{dx} = 0: either (x2)3=0(x-2)^3 = 0 giving x=2x = 2, or 6x2+x2=06x^2 + x - 2 = 0.

For 6x2+x2=06x^2 + x - 2 = 0, factor: 6x2+x2=(2x1)(3x+2)=06x^2 + x - 2 = (2x - 1)(3x + 2) = 0, so x=12x = \tfrac{1}{2} or x=23x = -\tfrac{2}{3}.

Compute yy at each:

  • x=23x = -\tfrac{2}{3}: y=(23)(13)(83)4=(23)(13)(409681)=8192729y = (-\tfrac{2}{3})(\tfrac{1}{3})(-\tfrac{8}{3})^4 = (-\tfrac{2}{3})(\tfrac{1}{3})(\tfrac{4096}{81}) = -\dfrac{8192}{729}

  • x=12x = \tfrac{1}{2}: y=(12)(32)(32)4=(34)(8116)=24364y = (\tfrac{1}{2})(\tfrac{3}{2})(-\tfrac{3}{2})^4 = (\tfrac{3}{4})(\tfrac{81}{16}) = \dfrac{243}{64}

  • x=2x = 2: y=2304=0y = 2 \cdot 3 \cdot 0^4 = 0

Stationary points: (23,8192729)\left(-\tfrac{2}{3},\, -\tfrac{8192}{729}\right), (12,24364)\left(\tfrac{1}{2},\, \tfrac{243}{64}\right), (2,0)(2,\, 0).

Nature of each stationary point:

Using the product rule on dydx=(x2)3(6x2+x2)\dfrac{dy}{dx} = (x-2)^3(6x^2 + x - 2):

d2ydx2=3(x2)2(6x2+x2)+(x2)3(12x+1)\frac{d^2y}{dx^2} = 3(x-2)^2(6x^2+x-2) + (x-2)^3(12x+1)

At x=23x = -\tfrac{2}{3}:

d2ydx2=3 ⁣(83) ⁣2 ⁣0+(83) ⁣3 ⁣(7)=51227×(7)=358427>0\frac{d^2y}{dx^2} = 3\!\left(-\tfrac{8}{3}\right)^{\!2}\!\cdot 0 + \left(-\tfrac{8}{3}\right)^{\!3}\!\left(-7\right) = -\frac{512}{27} \times (-7) = \frac{3584}{27} > 0

So (23,8192729)\left(-\tfrac{2}{3},\, -\tfrac{8192}{729}\right) is a minimum.

At x=12x = \tfrac{1}{2}:

d2ydx2=3 ⁣(32) ⁣2 ⁣0+(32) ⁣3 ⁣(7)=278×7=1898<0\frac{d^2y}{dx^2} = 3\!\left(-\tfrac{3}{2}\right)^{\!2}\!\cdot 0 + \left(-\tfrac{3}{2}\right)^{\!3}\!\left(7\right) = -\frac{27}{8} \times 7 = -\frac{189}{8} < 0

So (12,24364)\left(\tfrac{1}{2},\, \tfrac{243}{64}\right) is a maximum.

At x=2x = 2: both terms vanish, so y=0y'' = 0. We need a different approach. Near x=2x = 2:

y=x(x+1)(x2)423(x2)4=6(x2)4y = x(x+1)(x-2)^4 \approx 2 \cdot 3 \cdot (x-2)^4 = 6(x-2)^4

Since 6(x2)406(x-2)^4 \geq 0 with equality only at x=2x = 2, the point (2,0)(2, 0) is a minimum (a flat minimum, where the curve touches the xx-axis).

Sketch of CC:

Key features: y=0y = 0 at x=0,1,2x = 0, -1, 2. The curve has a minimum at (23,8192729)(0.67,11.24)\left(-\tfrac{2}{3},\, -\tfrac{8192}{729}\right) \approx (-0.67,\, -11.24), a maximum at (12,24364)(0.5,3.80)\left(\tfrac{1}{2},\, \tfrac{243}{64}\right) \approx (0.5,\, 3.80), and a flat minimum at (2,0)(2, 0). As x+x \to +\infty, y+y \to +\infty; as xx \to -\infty, y+y \to +\infty (since x(x+1)+x(x+1) \to +\infty and (x2)4+(x-2)^4 \to +\infty).

Sketch (i): y2=x(x+1)(x2)4y^2 = x(x+1)(x-2)^4

This curve has symmetry about the xx-axis (replacing yy by y-y leaves the equation unchanged).

It exists only where x(x+1)(x2)40x(x+1)(x-2)^4 \geq 0. Since (x2)40(x-2)^4 \geq 0 always, we need x(x+1)0x(x+1) \geq 0, i.e., x1x \leq -1 or x0x \geq 0. The region 1<x<0-1 < x < 0 has no curve.

The curve meets the xx-axis at x=1,0,2x = -1, 0, 2. Away from these points, y=±x(x+1)(x2)4y = \pm\sqrt{x(x+1)(x-2)^4}, so the curve is a symmetric pair of branches above and below the xx-axis, shaped like the original curve CC but reflected.

Sketch (ii): y=x2(x2+1)(x22)4y = x^2(x^2+1)(x^2-2)^4

This function is even: f(x)=f(x)f(-x) = f(x), so the curve is symmetric about the yy-axis.

Since x20x^2 \geq 0, (x2+1)>0(x^2+1) > 0, and (x22)40(x^2-2)^4 \geq 0, we have y0y \geq 0 for all xx. The curve lies entirely in the upper half-plane.

Zeros: x=0x = 0 (from x2x^2, double root — curve touches axis and turns back) and x=±2x = \pm\sqrt{2} (from (x22)4(x^2-2)^4, even multiplicity — flat minimum touching axis).

The shape in the right half-plane (x0x \geq 0) mirrors that of CC with xx replaced by x2x^2: a local maximum between x=0x = 0 and x=2x = \sqrt{2}, then a flat minimum at x=2x = \sqrt{2}, then yy \to \infty as xx \to \infty. By symmetry, the left half-plane is the mirror image.


Topic: 微积分/最优化 Calculus/Optimisation  |  Difficulty: Challenging  |  Marks: 20

4

Figure 1 Figure 2

(i) An attempt is made to move a rod of length LL from a corridor of width aa into a corridor of width bb, where aba \neq b. The corridors meet at right angles, as shown in Figure 1 and the rod remains horizontal. Show that if the attempt is to be successful then

Lacosecα+bsecα,L \le a \operatorname{cosec} \alpha + b \sec \alpha ,

where α\alpha satisfies

tan3α=ab.\tan^3 \alpha = \frac{a}{b} .

(ii) An attempt is made to move a rectangular table-top, of width ww and length ll, from one corridor to the other, as shown in the Figure 2. The table-top remains horizontal. Show that if the attempt is to be successful then

lacosecβ+bsecβ2wcosec2β,l \le a \operatorname{cosec} \beta + b \sec \beta - 2w \operatorname{cosec} 2\beta,

where β\beta satisfies

w=(abtan3β1tan2β)cosβ.w = \left( \frac{a - b \tan^3 \beta}{1 - \tan^2 \beta} \right) \cos \beta .

Hint

Q4 It is important to realise at the outset that α is a constant defined by a and b and that β is a constant defined by a, b and w. Variable angles θ/φ are needed to define the orientation of the rod/table in the general situation.

(i) Clearly, for all θ ∈ (0, π/2), it is necessary that f(θ) ≥ L, where f(θ) = a csc θ + b sec θ. Setting f’(θ) = 0 will then lead to the required result.

(ii) Here, for all φ ∈ (0, π/2), it is necessary that y ≥ l, where y is such that b = (y-x) cos φ + w sin φ and x is such that a = x sin φ + w cos φ. (Other formulations are possible.) Elimination of x leads to y = a csc φ + b sec φ - 2w csc 2φ

Setting y’(φ) = 0 plus some further working will then produce the required result.

Model Solution

Part (i).

Let θ\theta be the variable angle the rod makes with the corridor of width bb. The rod extends a distance asinθ=acosecθ\dfrac{a}{\sin\theta} = a\operatorname{cosec}\theta in the aa-corridor and bcosθ=bsecθ\dfrac{b}{\cos\theta} = b\sec\theta in the bb-corridor.

For the rod to fit through the corner:

Lacosecθ+bsecθfor all θ(0,π/2)(*)L \leq a\operatorname{cosec}\theta + b\sec\theta \quad \text{for all } \theta \in (0,\, \pi/2) \qquad \text{(*)}

Define f(θ)=acosecθ+bsecθf(\theta) = a\operatorname{cosec}\theta + b\sec\theta. The tightest constraint comes at the minimum of ff.

f(θ)=acosecθcotθ+bsecθtanθf'(\theta) = -a\operatorname{cosec}\theta\cot\theta + b\sec\theta\tan\theta

Setting f(θ)=0f'(\theta) = 0:

bsecθtanθ=acosecθcotθb\sec\theta\tan\theta = a\operatorname{cosec}\theta\cot\theta

bsinθcos2θ=acosθsin2θ\frac{b\sin\theta}{\cos^2\theta} = \frac{a\cos\theta}{\sin^2\theta}

Cross-multiply:

bsin3θ=acos3θb\sin^3\theta = a\cos^3\theta

tan3θ=ab\tan^3\theta = \frac{a}{b}

The critical angle α\alpha satisfies tan3α=ab\tan^3\alpha = \dfrac{a}{b}. Since f(θ)f(\theta) \to \infty as θ0+\theta \to 0^+ (because cosecθ\operatorname{cosec}\theta \to \infty) and as θπ2\theta \to \frac{\pi}{2}^- (because secθ\sec\theta \to \infty), this critical point must be a minimum.

From (*), evaluated at θ=α\theta = \alpha:

Lacosecα+bsecαL \leq a\operatorname{cosec}\alpha + b\sec\alpha

Part (ii).

Let ϕ\phi be the variable angle the table-top makes with the bb-corridor. Let xx be the distance from the corner along the table’s inner edge to the point where it meets the aa-corridor wall, and let yxy - x be the distance from the corner to the point where it meets the bb-corridor wall.

The table’s width ww is perpendicular to its length. Projecting onto each corridor direction:

a=xsinϕ+wcosϕ(*)a = x\sin\phi + w\cos\phi \qquad \text{(*)}

b=(yx)cosϕ+wsinϕ(**)b = (y - x)\cos\phi + w\sin\phi \qquad \text{(**)}

From (*): x=awcosϕsinϕx = \dfrac{a - w\cos\phi}{\sin\phi}.

From (**): yx=bwsinϕcosϕy - x = \dfrac{b - w\sin\phi}{\cos\phi}.

Adding:

y=awcosϕsinϕ+bwsinϕcosϕy = \frac{a - w\cos\phi}{\sin\phi} + \frac{b - w\sin\phi}{\cos\phi}

=asinϕwcosϕsinϕ+bcosϕwsinϕcosϕ= \frac{a}{\sin\phi} - \frac{w\cos\phi}{\sin\phi} + \frac{b}{\cos\phi} - \frac{w\sin\phi}{\cos\phi}

=acosecϕ+bsecϕw(cosϕsinϕ+sinϕcosϕ)= a\operatorname{cosec}\phi + b\sec\phi - w\left(\frac{\cos\phi}{\sin\phi} + \frac{\sin\phi}{\cos\phi}\right)

=acosecϕ+bsecϕwcos2ϕ+sin2ϕsinϕcosϕ= a\operatorname{cosec}\phi + b\sec\phi - w \cdot \frac{\cos^2\phi + \sin^2\phi}{\sin\phi\cos\phi}

=acosecϕ+bsecϕwsinϕcosϕ= a\operatorname{cosec}\phi + b\sec\phi - \frac{w}{\sin\phi\cos\phi}

Since sinϕcosϕ=12sin2ϕ\sin\phi\cos\phi = \frac{1}{2}\sin 2\phi, we have 1sinϕcosϕ=2cosec2ϕ\dfrac{1}{\sin\phi\cos\phi} = 2\operatorname{cosec}2\phi, so

y=acosecϕ+bsecϕ2wcosec2ϕy = a\operatorname{cosec}\phi + b\sec\phi - 2w\operatorname{cosec}2\phi

For the table to fit: lyl \leq y for all ϕ(0,π/2)\phi \in (0,\, \pi/2).

The critical angle β\beta minimizes yy. Setting dydϕ=0\dfrac{dy}{d\phi} = 0:

dydϕ=acosecϕcotϕ+bsecϕtanϕ+4wcosec2ϕcot2ϕ=0\frac{dy}{d\phi} = -a\operatorname{cosec}\phi\cot\phi + b\sec\phi\tan\phi + 4w\operatorname{cosec}2\phi\cot 2\phi = 0

Converting to sin\sin and cos\cos and using sin22ϕ=4sin2ϕcos2ϕ\sin^2 2\phi = 4\sin^2\phi\cos^2\phi:

acosϕsin2ϕ+bsinϕcos2ϕ+wcos2ϕsin2ϕcos2ϕ=0-\frac{a\cos\phi}{\sin^2\phi} + \frac{b\sin\phi}{\cos^2\phi} + \frac{w\cos 2\phi}{\sin^2\phi\cos^2\phi} = 0

Multiply through by sin2ϕcos2ϕ\sin^2\phi\cos^2\phi:

acos3ϕ+bsin3ϕ+wcos2ϕ=0-a\cos^3\phi + b\sin^3\phi + w\cos 2\phi = 0

Using cos2ϕ=cos2ϕsin2ϕ\cos 2\phi = \cos^2\phi - \sin^2\phi:

acos3ϕ+bsin3ϕ+w(cos2ϕsin2ϕ)=0-a\cos^3\phi + b\sin^3\phi + w(\cos^2\phi - \sin^2\phi) = 0

Divide by cos3ϕ\cos^3\phi:

a+btan3ϕ+w(1cosϕsin2ϕcos3ϕ)=0-a + b\tan^3\phi + w\left(\frac{1}{\cos\phi} - \frac{\sin^2\phi}{\cos^3\phi}\right) = 0

a+btan3ϕ+w(cos2ϕsin2ϕ)cos3ϕ=0-a + b\tan^3\phi + \frac{w(\cos^2\phi - \sin^2\phi)}{\cos^3\phi} = 0

a+btan3ϕ+w(1tan2ϕ)cosϕ=0-a + b\tan^3\phi + \frac{w(1 - \tan^2\phi)}{\cos\phi} = 0

w(1tan2ϕ)cosϕ=abtan3ϕ\frac{w(1 - \tan^2\phi)}{\cos\phi} = a - b\tan^3\phi

w=(abtan3ϕ)cosϕ1tan2ϕw = \frac{(a - b\tan^3\phi)\cos\phi}{1 - \tan^2\phi}

At ϕ=β\phi = \beta:

w=(abtan3β)cosβ1tan2βw = \frac{(a - b\tan^3\beta)\cos\beta}{1 - \tan^2\beta}


Topic: 积分/函数方程 Integration/Functional Equations  |  Difficulty: Standard  |  Marks: 20

5 Evaluate 0πxsinxdx\int_0^\pi x \sin x \, dx and 0πxcosxdx\int_0^\pi x \cos x \, dx.

The function ff satisfies the equation

f(t)=t+0πf(x)sin(x+t)dx.(*)f(t) = t + \int_0^\pi f(x) \sin(x + t) \, dx . \qquad \text{(*)}

Show that

f(t)=t+Asint+Bcost,(**)f(t) = t + A \sin t + B \cos t , \qquad \text{(**)}

where A=0πf(x)cosxdxA = \int_0^\pi f(x) \cos x \, dx and B=0πf(x)sinxdxB = \int_0^\pi f(x) \sin x \, dx.

Use the expression (**) to find AA and BB by substituting for f(t)f(t) and f(x)f(x) in (*) and equating coefficients of sint\sin t and cost\cos t.

Hint

Q5 Using the integration by parts rule it is easy to establish the results ∫₀^π x sin x dx = π and ∫₀^π x cos x dx = -2.

  • Write sin(x + t) = sin x cos t + sin t cos x and the result f(t) = t + A sin t + B cos t, where A and B are as defined in the question, follows immediately.

  • Hence write t + A sin t + B cos t = t + ∫₀^π (x + A sin x + B cos x) sin(x + t) dx (***) so that as

∫₀^π x sin(x + t) dx = … = π cos t - 2 sin t,

∫₀^π sin x sin(x + t) dx = … = (π/2) cos t,

∫₀^π cos x sin(x + t) dx = … = (π/2) sin t,

then, by considering the coefficients of cos t and sin t on both sides of (***), it follows that

A = -2 + (π/2)B, B = π + (π/2)A => A = -2, B = 0.

Alternatively, equations for A and B can be obtained by putting t = 0 and t = π/2 in (***).

Model Solution

Part 1: Evaluating the integrals

Integral 1: I1=0πxsinxdx\displaystyle I_1 = \int_0^\pi x \sin x \, dx

Integration by parts with u=xu = x, dv=sinxdxdv = \sin x \, dx (so du=dxdu = dx, v=cosxv = -\cos x):

I1=[xcosx]0π+0πcosxdx=πcosπ+0+[sinx]0π=π+0=πI_1 = \left[-x\cos x\right]_0^\pi + \int_0^\pi \cos x \, dx = -\pi\cos\pi + 0 + \left[\sin x\right]_0^\pi = \pi + 0 = \pi

0πxsinxdx=π(I)\int_0^\pi x\sin x \, dx = \pi \qquad \qquad \text{(I)}

Integral 2: I2=0πxcosxdx\displaystyle I_2 = \int_0^\pi x \cos x \, dx

Integration by parts with u=xu = x, dv=cosxdxdv = \cos x \, dx (so du=dxdu = dx, v=sinxv = \sin x):

I2=[xsinx]0π0πsinxdx=0[cosx]0π=cosπcos0=2I_2 = \left[x\sin x\right]_0^\pi - \int_0^\pi \sin x \, dx = 0 - \left[-\cos x\right]_0^\pi = \cos\pi - \cos 0 = -2

0πxcosxdx=2(II)\int_0^\pi x\cos x \, dx = -2 \qquad \qquad \text{(II)}

Part 2: Showing f(t)=t+Asint+Bcostf(t) = t + A\sin t + B\cos t

Expand sin(x+t)=sinxcost+cosxsint\sin(x + t) = \sin x \cos t + \cos x \sin t in equation (* *):

f(t)=t+0πf(x)(sinxcost+cosxsint)dxf(t) = t + \int_0^\pi f(x)(\sin x\cos t + \cos x\sin t) \, dx

=t+cost0πf(x)sinxdx+sint0πf(x)cosxdx= t + \cos t \int_0^\pi f(x)\sin x \, dx + \sin t \int_0^\pi f(x)\cos x \, dx

=t+Bcost+Asint= t + B\cos t + A\sin t \qquad \checkmark

Part 3: Finding AA and BB

Substitute f(x)=x+Asinx+Bcosxf(x) = x + A\sin x + B\cos x into A=0πf(x)cosxdxA = \int_0^\pi f(x)\cos x \, dx:

A=0π(x+Asinx+Bcosx)cosxdxA = \int_0^\pi (x + A\sin x + B\cos x)\cos x \, dx

=0πxcosxdx+A0πsinxcosxdx+B0πcos2xdx= \int_0^\pi x\cos x \, dx + A\int_0^\pi \sin x\cos x \, dx + B\int_0^\pi \cos^2 x \, dx

Evaluate each piece. We already know 0πxcosxdx=2\int_0^\pi x\cos x \, dx = -2 from (II).

For the second integral: 0πsinxcosxdx=120πsin2xdx=12[12cos2x]0π=14(cos2πcos0)=0\int_0^\pi \sin x\cos x \, dx = \tfrac{1}{2}\int_0^\pi \sin 2x \, dx = \tfrac{1}{2}\left[-\tfrac{1}{2}\cos 2x\right]_0^\pi = -\tfrac{1}{4}(\cos 2\pi - \cos 0) = 0.

For the third integral: 0πcos2xdx=120π(1+cos2x)dx=12[x+12sin2x]0π=π2\int_0^\pi \cos^2 x \, dx = \tfrac{1}{2}\int_0^\pi (1 + \cos 2x) \, dx = \tfrac{1}{2}\left[x + \tfrac{1}{2}\sin 2x\right]_0^\pi = \tfrac{\pi}{2}.

Therefore

A=2+0+π2B(III)A = -2 + 0 + \frac{\pi}{2}B \qquad \qquad \text{(III)}

Substitute into B=0πf(x)sinxdxB = \int_0^\pi f(x)\sin x \, dx:

B=0π(x+Asinx+Bcosx)sinxdxB = \int_0^\pi (x + A\sin x + B\cos x)\sin x \, dx

=0πxsinxdx+A0πsin2xdx+B0πcosxsinxdx= \int_0^\pi x\sin x \, dx + A\int_0^\pi \sin^2 x \, dx + B\int_0^\pi \cos x\sin x \, dx

We know 0πxsinxdx=π\int_0^\pi x\sin x \, dx = \pi from (I), 0πsin2xdx=π2\int_0^\pi \sin^2 x \, dx = \tfrac{\pi}{2} (by symmetry with cos2\cos^2), and 0πsinxcosxdx=0\int_0^\pi \sin x\cos x \, dx = 0.

Therefore

B=π+π2A+0(IV)B = \pi + \frac{\pi}{2}A + 0 \qquad \qquad \text{(IV)}

From (III): A=2+π2BA = -2 + \frac{\pi}{2}B.

Substitute into (IV):

B=π+π2(2+π2B)=ππ+π24B=π24BB = \pi + \frac{\pi}{2}\left(-2 + \frac{\pi}{2}B\right) = \pi - \pi + \frac{\pi^2}{4}B = \frac{\pi^2}{4}B

Bπ24B=0B(1π24)=0B - \frac{\pi^2}{4}B = 0 \qquad \Longrightarrow \qquad B\left(1 - \frac{\pi^2}{4}\right) = 0

Since π24\pi^2 \neq 4, we must have B=0B = 0.

Then from (III): A=2+0=2A = -2 + 0 = -2.

A=2,B=0\boxed{A = -2, \quad B = 0}

so f(t)=t2sintf(t) = t - 2\sin t.


Topic: 向量 Vectors  |  Difficulty: Challenging  |  Marks: 20

6 The vectors a and b lie in the plane Π\Pi. Given that a=1|\mathbf{a}| = 1 and ab=3\mathbf{a} \cdot \mathbf{b} = 3, find, in terms of a and b, a vector p parallel to a and a vector q perpendicular to a, both lying in the plane Π\Pi, such that

p+q=a+b.\mathbf{p} + \mathbf{q} = \mathbf{a} + \mathbf{b} .

The vector c is not parallel to the plane Π\Pi and is such that ac=2\mathbf{a} \cdot \mathbf{c} = -2 and bc=2\mathbf{b} \cdot \mathbf{c} = 2. Given that b=5|\mathbf{b}| = 5, find, in terms of a, b and c, vectors P, Q and R such that P and Q are parallel to p and q, respectively, R is perpendicular to the plane Π\Pi and

P+Q+R=a+b+c.\mathbf{P} + \mathbf{Q} + \mathbf{R} = \mathbf{a} + \mathbf{b} + \mathbf{c} .

Hint

From the data it follows that the component of b in the direction of a is 3a3\mathbf{a}.

Hence p=4a\mathbf{p} = 4\mathbf{a} and q=b3a\mathbf{q} = \mathbf{b} - 3\mathbf{a}.

  • Again from the data, it follows that (c.a)a=2a(\mathbf{c.a})\mathbf{a} = -2\mathbf{a} and

q2=b.b6a.b+9a.a=2518+9=16q=4|\mathbf{q}|^2 = \mathbf{b.b} - 6\mathbf{a.b} + 9\mathbf{a.a} = 25 - 18 + 9 = 16 \Rightarrow |\mathbf{q}| = 4, so that

\left[(\mathbf{c.q})/|\mathbf{q}|^2 ight] \mathbf{q} = (1/2)\mathbf{b} - (3/2)\mathbf{a}.

Thus P=2aP = 2a, Q=(9/2)a+(3/2)bQ = -(9/2)a + (3/2)b, R=(7/2)a(1/2)b+cR = (7/2)a - (1/2)b + c.

Model Solution

Part 1: Finding p\mathbf{p} and q\mathbf{q}

We need p\mathbf{p} parallel to a\mathbf{a} and q\mathbf{q} perpendicular to a\mathbf{a}, both in Π\Pi, with p+q=a+b\mathbf{p} + \mathbf{q} = \mathbf{a} + \mathbf{b}.

Since p\mathbf{p} is parallel to a\mathbf{a}, write p=λa\mathbf{p} = \lambda\mathbf{a} for some scalar λ\lambda. Then q=a+bλa\mathbf{q} = \mathbf{a} + \mathbf{b} - \lambda\mathbf{a}.

The condition qa\mathbf{q} \perp \mathbf{a} gives qa=0\mathbf{q} \cdot \mathbf{a} = 0:

(a+bλa)a=0(\mathbf{a} + \mathbf{b} - \lambda\mathbf{a}) \cdot \mathbf{a} = 0

aa+baλaa=0\mathbf{a} \cdot \mathbf{a} + \mathbf{b} \cdot \mathbf{a} - \lambda\,\mathbf{a} \cdot \mathbf{a} = 0

1+3λ1=0λ=41 + 3 - \lambda \cdot 1 = 0 \qquad \Longrightarrow \qquad \lambda = 4

Therefore

p=4a,q=b3a\mathbf{p} = 4\mathbf{a}, \qquad \mathbf{q} = \mathbf{b} - 3\mathbf{a}

Part 2: Finding P\mathbf{P}, Q\mathbf{Q} and R\mathbf{R}

We need P\mathbf{P} parallel to p\mathbf{p} (i.e., parallel to a\mathbf{a}), Q\mathbf{Q} parallel to q\mathbf{q}, R\mathbf{R} perpendicular to Π\Pi, and P+Q+R=a+b+c\mathbf{P} + \mathbf{Q} + \mathbf{R} = \mathbf{a} + \mathbf{b} + \mathbf{c}.

Write P=μa\mathbf{P} = \mu\mathbf{a} and Q=νq=ν(b3a)\mathbf{Q} = \nu\mathbf{q} = \nu(\mathbf{b} - 3\mathbf{a}) for scalars μ,ν\mu, \nu.

Then R=a+b+cμaν(b3a)\mathbf{R} = \mathbf{a} + \mathbf{b} + \mathbf{c} - \mu\mathbf{a} - \nu(\mathbf{b} - 3\mathbf{a}).

Since RΠ\mathbf{R} \perp \Pi and a,bΠ\mathbf{a}, \mathbf{b} \in \Pi, we need Ra=0\mathbf{R} \cdot \mathbf{a} = 0 and Rq=0\mathbf{R} \cdot \mathbf{q} = 0 (where q\mathbf{q} also lies in Π\Pi and is perpendicular to a\mathbf{a}, so {a,q}\{\mathbf{a}, \mathbf{q}\} spans Π\Pi).

Condition Ra=0\mathbf{R} \cdot \mathbf{a} = 0:

Ra=(a+b+cμaνb+3νa)a=0\mathbf{R} \cdot \mathbf{a} = (\mathbf{a} + \mathbf{b} + \mathbf{c} - \mu\mathbf{a} - \nu\mathbf{b} + 3\nu\mathbf{a}) \cdot \mathbf{a} = 0

=aa+ba+caμaaνba+3νaa=0= \mathbf{a}\cdot\mathbf{a} + \mathbf{b}\cdot\mathbf{a} + \mathbf{c}\cdot\mathbf{a} - \mu\,\mathbf{a}\cdot\mathbf{a} - \nu\,\mathbf{b}\cdot\mathbf{a} + 3\nu\,\mathbf{a}\cdot\mathbf{a} = 0

=1+3+(2)μ3ν+3ν=0= 1 + 3 + (-2) - \mu - 3\nu + 3\nu = 0

2μ=0μ=22 - \mu = 0 \qquad \Longrightarrow \qquad \mu = 2

So P=2a\mathbf{P} = 2\mathbf{a}.

Condition Rq=0\mathbf{R} \cdot \mathbf{q} = 0:

First compute qq=q2\mathbf{q} \cdot \mathbf{q} = |\mathbf{q}|^2:

q2=(b3a)(b3a)=b26ab+9a2=2518+9=16|\mathbf{q}|^2 = (\mathbf{b} - 3\mathbf{a}) \cdot (\mathbf{b} - 3\mathbf{a}) = |\mathbf{b}|^2 - 6\,\mathbf{a}\cdot\mathbf{b} + 9|\mathbf{a}|^2 = 25 - 18 + 9 = 16

So q=4|\mathbf{q}| = 4.

Now compute cq=c(b3a)=bc3ac=23(2)=8\mathbf{c} \cdot \mathbf{q} = \mathbf{c} \cdot (\mathbf{b} - 3\mathbf{a}) = \mathbf{b}\cdot\mathbf{c} - 3\,\mathbf{a}\cdot\mathbf{c} = 2 - 3(-2) = 8.

Also aq=a(b3a)=33=0\mathbf{a} \cdot \mathbf{q} = \mathbf{a} \cdot (\mathbf{b} - 3\mathbf{a}) = 3 - 3 = 0 (as expected, since qa\mathbf{q} \perp \mathbf{a}).

Now Rq=0\mathbf{R} \cdot \mathbf{q} = 0:

Rq=(a+b+c2aνb+3νa)q=0\mathbf{R} \cdot \mathbf{q} = (\mathbf{a} + \mathbf{b} + \mathbf{c} - 2\mathbf{a} - \nu\mathbf{b} + 3\nu\mathbf{a}) \cdot \mathbf{q} = 0

=(a+(1ν)b+c+3νa)q=0= (-\mathbf{a} + (1-\nu)\mathbf{b} + \mathbf{c} + 3\nu\mathbf{a}) \cdot \mathbf{q} = 0

=((1+3ν)a+(1ν)b+c)q=0= ((-1+3\nu)\mathbf{a} + (1-\nu)\mathbf{b} + \mathbf{c}) \cdot \mathbf{q} = 0

Since aq=0\mathbf{a} \cdot \mathbf{q} = 0:

(1ν)bq+cq=0(1-\nu)\,\mathbf{b}\cdot\mathbf{q} + \mathbf{c}\cdot\mathbf{q} = 0

We have bq=b(b3a)=259=16\mathbf{b}\cdot\mathbf{q} = \mathbf{b}\cdot(\mathbf{b}-3\mathbf{a}) = 25 - 9 = 16.

So: (1ν)(16)+8=0(1-\nu)(16) + 8 = 0, giving 1ν=121 - \nu = -\tfrac{1}{2}, hence ν=32\nu = \tfrac{3}{2}.

Therefore

Q=32(b3a)=32b92a\mathbf{Q} = \tfrac{3}{2}(\mathbf{b} - 3\mathbf{a}) = \tfrac{3}{2}\mathbf{b} - \tfrac{9}{2}\mathbf{a}

Finding R\mathbf{R}:

R=a+b+cPQ=a+b+c2a32b+92a\mathbf{R} = \mathbf{a} + \mathbf{b} + \mathbf{c} - \mathbf{P} - \mathbf{Q} = \mathbf{a} + \mathbf{b} + \mathbf{c} - 2\mathbf{a} - \tfrac{3}{2}\mathbf{b} + \tfrac{9}{2}\mathbf{a}

=(12+92)a+(132)b+c=72a12b+c= \left(1 - 2 + \tfrac{9}{2}\right)\mathbf{a} + \left(1 - \tfrac{3}{2}\right)\mathbf{b} + \mathbf{c} = \tfrac{7}{2}\mathbf{a} - \tfrac{1}{2}\mathbf{b} + \mathbf{c}

P=2a,Q=92a+32b,R=72a12b+c\boxed{\mathbf{P} = 2\mathbf{a}, \quad \mathbf{Q} = -\tfrac{9}{2}\mathbf{a} + \tfrac{3}{2}\mathbf{b}, \quad \mathbf{R} = \tfrac{7}{2}\mathbf{a} - \tfrac{1}{2}\mathbf{b} + \mathbf{c}}

Verification: P+Q+R=(292+72)a+(0+3212)b+c=a+b+c\mathbf{P} + \mathbf{Q} + \mathbf{R} = (2 - \tfrac{9}{2} + \tfrac{7}{2})\mathbf{a} + (0 + \tfrac{3}{2} - \tfrac{1}{2})\mathbf{b} + \mathbf{c} = \mathbf{a} + \mathbf{b} + \mathbf{c}


Topic: 数值方法 Numerical Methods  |  Difficulty: Standard  |  Marks: 20

7 The function ff is defined by

f(x)=2sinxx.f(x) = 2 \sin x - x .

Show graphically that the equation f(x)=0f(x) = 0 has exactly one root in the interval [12π,π][\frac{1}{2}\pi, \pi]. This interval is denoted I0I_0.

In order to determine the root, a sequence of intervals I1,I2,I_1, I_2, \dots is generated in the following way. If the interval In=[an,bn]I_n = [a_n, b_n], and cn=(an+bn)/2c_n = (a_n + b_n)/2, then

In+1={[an,cn]if f(an)f(cn)<0;[cn,bn]if f(cn)f(bn)<0.I_{n+1} = \begin{cases} [a_n, c_n] & \text{if } f(a_n)f(c_n) < 0 ; \\ [c_n, b_n] & \text{if } f(c_n)f(b_n) < 0 . \end{cases}

By using the approximations 120.7\frac{1}{\sqrt{2}} \approx 0.7 and π10\pi \approx \sqrt{10}, show that I2=[12π,58π]I_2 = [\frac{1}{2}\pi, \frac{5}{8}\pi] and find I3I_3.

Hint

Good sketch graphs of y=xy = x and y=2sinxy = 2 \sin x, in the same diagram and over the interval [0,π][0, \pi], will readily show that the equation f(x)=0f(x) = 0 has exactly one root in the interval [π/2,π][\pi/2, \pi].

  • f(3π/4)=23π/4f(3\pi/4) = \sqrt{2} - 3\pi/4 has the same sign as 2 - 9\pi^2/16 pprox 2 - 45/8 = -29/8 < 0. Hence as f(π/2)=2π/2>0f(\pi/2) = 2 - \pi/2 > 0 and f(π)=π<0f(\pi) = -\pi < 0, then I1=[π/2,3π/4]I_1 = [\pi/2, 3\pi/4].

  • x=sin5π/82x1x2=sin3π/4=1/28x48x2+1=0x = \sin 5\pi/8 \Rightarrow 2x\sqrt{1 - x^2} = \sin 3\pi/4 = 1/\sqrt{2} \Rightarrow 8x^4 - 8x^2 + 1 = 0 (*) \Rightarrow x^2 = 1/2 + 1/(2\sqrt{2}) pprox 0.85. (**). The sign of f(5π/8)f(5\pi/8) is the same as that of 4x^2 - 25\pi^2/64 pprox 17/5 - 125/32 = -81/625 < 0. Hence I2=[π/2,5π/8]I_2 = [\pi/2, 5\pi/8].

  • A good approximation to x=sin9π/16x = \sin 9\pi/16 may also be obtained in a similar way. In fact, it will be found that f(9π/16)>0f(9\pi/16) > 0 so that I3=[9π/16,5π/8]I_3 = [9\pi/16, 5\pi/8].

Model Solution

Graphical argument for exactly one root in [π2,π]\left[\tfrac{\pi}{2},\, \pi\right]:

Consider f(x)=2sinxxf(x) = 2\sin x - x on [0,π][0, \pi]. At the endpoints of the interval:

f ⁣(π2)=2sinπ2π2=2π221.57=0.43>0f\!\left(\tfrac{\pi}{2}\right) = 2\sin\tfrac{\pi}{2} - \tfrac{\pi}{2} = 2 - \tfrac{\pi}{2} \approx 2 - 1.57 = 0.43 > 0

f(π)=2sinππ=π<0f(\pi) = 2\sin\pi - \pi = -\pi < 0

Since ff is continuous and changes sign on [π2,π]\left[\tfrac{\pi}{2},\, \pi\right], there is at least one root. Moreover, f(x)=2cosx1f'(x) = 2\cos x - 1, which is negative on [π2,π]\left[\tfrac{\pi}{2},\, \pi\right] (since cosx0\cos x \leq 0 there), so ff is strictly decreasing. Hence there is exactly one root in this interval.

Bisection method:

We use the approximations 120.7\dfrac{1}{\sqrt{2}} \approx 0.7 and π10\pi \approx \sqrt{10}.

Step 1: I0=[π2,π]I_0 = \left[\tfrac{\pi}{2},\, \pi\right]. Midpoint c0=3π4c_0 = \tfrac{3\pi}{4}.

Evaluate f ⁣(3π4)=2sin3π43π4=2123π4=23π4f\!\left(\tfrac{3\pi}{4}\right) = 2\sin\tfrac{3\pi}{4} - \tfrac{3\pi}{4} = 2 \cdot \tfrac{1}{\sqrt{2}} - \tfrac{3\pi}{4} = \sqrt{2} - \tfrac{3\pi}{4}.

To determine the sign, note (2)2=2\left(\sqrt{2}\right)^2 = 2 and (3π4)2=9π2169016=458=5.625\left(\tfrac{3\pi}{4}\right)^2 = \tfrac{9\pi^2}{16} \approx \tfrac{90}{16} = \tfrac{45}{8} = 5.625. Since 2<5.6252 < 5.625, we have 2<3π4\sqrt{2} < \tfrac{3\pi}{4}, so f ⁣(3π4)<0f\!\left(\tfrac{3\pi}{4}\right) < 0.

Since f ⁣(π2)>0f\!\left(\tfrac{\pi}{2}\right) > 0 and f ⁣(3π4)<0f\!\left(\tfrac{3\pi}{4}\right) < 0, the root lies in I1=[π2,3π4]I_1 = \left[\tfrac{\pi}{2},\, \tfrac{3\pi}{4}\right].

Step 2: Midpoint c1=12 ⁣(π2+3π4)=5π8c_1 = \tfrac{1}{2}\!\left(\tfrac{\pi}{2} + \tfrac{3\pi}{4}\right) = \tfrac{5\pi}{8}.

Evaluate f ⁣(5π8)=2sin5π85π8f\!\left(\tfrac{5\pi}{8}\right) = 2\sin\tfrac{5\pi}{8} - \tfrac{5\pi}{8}.

We need sin5π8\sin\tfrac{5\pi}{8}. Using the double angle formula with 5π8=π2+π8\tfrac{5\pi}{8} = \tfrac{\pi}{2} + \tfrac{\pi}{8}:

sin5π8=cosπ8\sin\tfrac{5\pi}{8} = \cos\tfrac{\pi}{8}

By the half-angle formula: cosπ8=1+cos(π/4)2=1+1/22=1+0.72=0.850.92\cos\tfrac{\pi}{8} = \sqrt{\tfrac{1+\cos(\pi/4)}{2}} = \sqrt{\tfrac{1 + 1/\sqrt{2}}{2}} = \sqrt{\tfrac{1 + 0.7}{2}} = \sqrt{0.85} \approx 0.92.

So f ⁣(5π8)2(0.92)51081.845×3.1681.841.98=0.14<0f\!\left(\tfrac{5\pi}{8}\right) \approx 2(0.92) - \tfrac{5\sqrt{10}}{8} \approx 1.84 - \tfrac{5 \times 3.16}{8} \approx 1.84 - 1.98 = -0.14 < 0.

More precisely: f ⁣(5π8)f\!\left(\tfrac{5\pi}{8}\right) has the same sign as (2sin5π8)2(5π8)2=4sin25π825π264\left(2\sin\tfrac{5\pi}{8}\right)^2 - \left(\tfrac{5\pi}{8}\right)^2 = 4\sin^2\tfrac{5\pi}{8} - \tfrac{25\pi^2}{64}.

From sin5π8=cosπ8\sin\tfrac{5\pi}{8} = \cos\tfrac{\pi}{8}: 4cos2π8=41+1/22=2+22+1.41=3.414\cos^2\tfrac{\pi}{8} = 4 \cdot \tfrac{1+1/\sqrt{2}}{2} = 2 + \sqrt{2} \approx 2 + 1.41 = 3.41.

And 25π26425064=3.91\tfrac{25\pi^2}{64} \approx \tfrac{250}{64} = 3.91. Since 3.41<3.913.41 < 3.91, f ⁣(5π8)<0f\!\left(\tfrac{5\pi}{8}\right) < 0.

Since f ⁣(π2)>0f\!\left(\tfrac{\pi}{2}\right) > 0 and f ⁣(5π8)<0f\!\left(\tfrac{5\pi}{8}\right) < 0:

I2=[π2,5π8]I_2 = \left[\tfrac{\pi}{2},\, \tfrac{5\pi}{8}\right] \qquad \checkmark

Step 3: Midpoint c2=12 ⁣(π2+5π8)=9π16c_2 = \tfrac{1}{2}\!\left(\tfrac{\pi}{2} + \tfrac{5\pi}{8}\right) = \tfrac{9\pi}{16}.

Evaluate f ⁣(9π16)=2sin9π169π16f\!\left(\tfrac{9\pi}{16}\right) = 2\sin\tfrac{9\pi}{16} - \tfrac{9\pi}{16}.

We need sin9π16\sin\tfrac{9\pi}{16}. Since 9π16=π2+π16\tfrac{9\pi}{16} = \tfrac{\pi}{2} + \tfrac{\pi}{16}:

sin9π16=cosπ16\sin\tfrac{9\pi}{16} = \cos\tfrac{\pi}{16}

By the half-angle formula: cosπ16=1+cos(π/8)2\cos\tfrac{\pi}{16} = \sqrt{\tfrac{1+\cos(\pi/8)}{2}}.

We already found cosπ80.92\cos\tfrac{\pi}{8} \approx 0.92, so cosπ161.922=0.960.98\cos\tfrac{\pi}{16} \approx \sqrt{\tfrac{1.92}{2}} = \sqrt{0.96} \approx 0.98.

Then f ⁣(9π16)2(0.98)910161.9628.46161.961.78=0.18>0f\!\left(\tfrac{9\pi}{16}\right) \approx 2(0.98) - \tfrac{9\sqrt{10}}{16} \approx 1.96 - \tfrac{28.46}{16} \approx 1.96 - 1.78 = 0.18 > 0.

Check precisely: 4cos2π16=41+cos(π/8)2=2+2cosπ8=2+21+1/222+2(0.92)=3.844\cos^2\tfrac{\pi}{16} = 4 \cdot \tfrac{1 + \cos(\pi/8)}{2} = 2 + 2\cos\tfrac{\pi}{8} = 2 + 2\sqrt{\tfrac{1+1/\sqrt{2}}{2}} \approx 2 + 2(0.92) = 3.84.

And (9π16)2=81π2256810256=3.16\left(\tfrac{9\pi}{16}\right)^2 = \tfrac{81\pi^2}{256} \approx \tfrac{810}{256} = 3.16. Since 3.84>3.163.84 > 3.16, f ⁣(9π16)>0f\!\left(\tfrac{9\pi}{16}\right) > 0.

Since f ⁣(9π16)>0f\!\left(\tfrac{9\pi}{16}\right) > 0 and f ⁣(5π8)<0f\!\left(\tfrac{5\pi}{8}\right) < 0:

I3=[9π16,5π8]I_3 = \left[\tfrac{9\pi}{16},\, \tfrac{5\pi}{8}\right]


Topic: 微分方程 Differential Equations  |  Difficulty: Challenging  |  Marks: 20

8 Let xx satisfy the differential equation

dxdt=(1xn)1/n\frac{dx}{dt} = (1 - x^n)^{1/n}

and the condition x=0x = 0 when t=0t = 0.

(i) Solve the equation in the case n=1n = 1 and sketch the graph of the solution for t>0t > 0.

(ii) Prove that 1x<(1x2)1/21 - x < (1 - x^2)^{1/2} for 0<x<10 < x < 1.

Use this result to sketch the graph of the solution in the case n=2n = 2 for 0<t<12π0 < t < \frac{1}{2}\pi, using the same axes as your previous sketch.

By setting x=sinyx = \sin y, solve the equation in this case.

(iii) Use the result (which you need not prove)

(1x2)1/2<(1x3)1/3for0<x<1,(1 - x^2)^{1/2} < (1 - x^3)^{1/3} \quad \text{for} \quad 0 < x < 1 ,

to sketch, without solving the equation, the graph of the solution of the equation in the case n=3n = 3 using the same axes as your previous sketches. Use your sketch to show that x=1x = 1 at a value of tt less than 12π\frac{1}{2}\pi.

Hint

Q8(i) Integration leads to the general solution t=Aln(1x)t = A - \ln(1 - x) and x(0)=0A=0x(0) = 0 \Rightarrow A = 0. Thus x=1etx = 1 - e^{-t}.

(ii) Obviously, (1x)1/2<(1+x)1/2(1 - x)^{1/2} < (1 + x)^{1/2} for all x(0,1]x \in (0, 1]. Hence multiplying this inequality through by (1x)1/2(1 - x)^{1/2} leads immediately to the required result.

Arguments which go in the wrong direction, e.g., 1x<(1x2)1/2xx2>01 - x < (1 - x^2)^{1/2} \Rightarrow \dots \Rightarrow x - x^2 > 0, etc., are invalid. It may be possible to salvage them by replacing ’\Rightarrow' by ’\Leftarrow'.

In the case n=2n = 2, the substitution x=sinyx = \sin y will lead to t=y+Bt = y + B and hence to t=sin1(x)+Bt = \sin^{-1}(x) + B as the general solution. In particular, x(0)=0B=0x=sintx(0) = 0 \Rightarrow B = 0 \Rightarrow x = \sin t.

Note that the question does not allow the use of the standard form (1x2)1/2dx=sin1(x)+\int (1 - x^2)^{-1/2} dx = \sin^{-1}(x) + an arbitrary constant, without proof.

(iii) If GnG_n is the graph of xx for 0x10 \le x \le 1, then the given inequality shows that the gradient of G3G_3 is greater than the gradient of G2G_2 for each xx in this interval. (The inequality of (ii) shows that the same is true of G2G_2 in relation to G1G_1.) These considerations will help to clarify ideas when drawing sketches of GnG_n for n=1,2,3n = 1, 2, 3 in the same diagram. In particular, the sketch of G3G_3 should make it clear that once xx reaches the value 1 it remains there.

Model Solution

Part (i): n=1n = 1

The equation becomes dxdt=1x\dfrac{dx}{dt} = 1 - x, with x(0)=0x(0) = 0.

Separate variables:

dx1x=dt(x1)\frac{dx}{1 - x} = dt \qquad (x \neq 1)

ln1x=t+C-\ln|1 - x| = t + C

At t=0t = 0, x=0x = 0: ln1=C-\ln 1 = C, so C=0C = 0.

ln(1x)=t1x=etx=1et-\ln(1 - x) = t \qquad \Longrightarrow \qquad 1 - x = e^{-t} \qquad \Longrightarrow \qquad x = 1 - e^{-t}

(valid for x<1x < 1, which holds for all finite t>0t > 0).

Sketch for n=1n = 1: The curve starts at the origin, rises steeply at first (gradient 1 at t=0t = 0), and approaches x=1x = 1 asymptotically as tt \to \infty. It never actually reaches x=1x = 1.

Part (ii): Proving 1x<(1x2)1/21 - x < (1 - x^2)^{1/2} for 0<x<10 < x < 1

Both sides are positive for 0<x<10 < x < 1, so squaring preserves the inequality. We need to show:

(1x)2<1x2(1 - x)^2 < 1 - x^2

Expand the left side: 12x+x21 - 2x + x^2.

So we need: 12x+x2<1x21 - 2x + x^2 < 1 - x^2, i.e., 2x22x<02x^2 - 2x < 0, i.e., 2x(x1)<02x(x - 1) < 0.

Since 0<x<10 < x < 1, we have x>0x > 0 and x1<0x - 1 < 0, so 2x(x1)<02x(x-1) < 0 ✓.

This is a valid proof: we started from 2x(x1)<02x(x-1) < 0 (true for 0<x<10 < x < 1) and derived (1x)2<1x2(1-x)^2 < 1 - x^2 by reversible algebraic steps, from which 1x<(1x2)1/21 - x < (1-x^2)^{1/2} follows since both sides are positive.

Solving for n=2n = 2:

The equation is dxdt=(1x2)1/2\dfrac{dx}{dt} = (1 - x^2)^{1/2} with x(0)=0x(0) = 0.

Set x=sinyx = \sin y, so dxdt=cosydydt\dfrac{dx}{dt} = \cos y \cdot \dfrac{dy}{dt}. Then:

cosydydt=(1sin2y)1/2=cosy\cos y \cdot \frac{dy}{dt} = (1 - \sin^2 y)^{1/2} = \cos y

For 0<y<π20 < y < \tfrac{\pi}{2}, cosy>0\cos y > 0, so:

dydt=1y=t+C\frac{dy}{dt} = 1 \qquad \Longrightarrow \qquad y = t + C

At t=0t = 0, x=0x = 0, so siny=0\sin y = 0, giving y=0y = 0 and C=0C = 0.

Therefore y=ty = t, and x=sintx = \sin t.

Sketch for n=2n = 2: x=sintx = \sin t starts at the origin, rises more steeply than x=1etx = 1 - e^{-t} (since (1x2)1/2>1x(1-x^2)^{1/2} > 1 - x for 0<x<10 < x < 1), and reaches x=1x = 1 at t=π2t = \tfrac{\pi}{2}. For t(0,π2)t \in (0, \tfrac{\pi}{2}), the curve x=sintx = \sin t lies above x=1etx = 1 - e^{-t}.

Part (iii): Sketch for n=3n = 3

We are given that (1x2)1/2<(1x3)1/3(1 - x^2)^{1/2} < (1 - x^3)^{1/3} for 0<x<10 < x < 1.

This means dxdtn=3>dxdtn=2\dfrac{dx}{dt}\Big|_{n=3} > \dfrac{dx}{dt}\Big|_{n=2} for each x(0,1)x \in (0, 1). So at any given xx, the solution for n=3n = 3 has a larger gradient than for n=2n = 2.

Combined with part (ii), for 0<x<10 < x < 1:

1xgradient for n=1<(1x2)1/2gradient for n=2<(1x3)1/3gradient for n=3\underbrace{1 - x}_{\text{gradient for } n=1} < \underbrace{(1-x^2)^{1/2}}_{\text{gradient for } n=2} < \underbrace{(1-x^3)^{1/3}}_{\text{gradient for } n=3}

Since the gradient is larger at every xx, the solution x(t)x(t) for n=3n = 3 rises more steeply than for n=2n = 2, which rises more steeply than for n=1n = 1. Starting from the origin, the three curves are ordered: G1G_1 (lowest), G2G_2 (middle), G3G_3 (highest) for t>0t > 0.

Since G2G_2 reaches x=1x = 1 at t=π2t = \tfrac{\pi}{2}, and G3G_3 rises more steeply, G3G_3 must reach x=1x = 1 at some value t<π2t^* < \tfrac{\pi}{2}.

Once x=1x = 1, we have dxdt=(11n)1/n=0\dfrac{dx}{dt} = (1 - 1^n)^{1/n} = 0, so xx stays at 11 for all ttt \geq t^*. The graph of G3G_3 is therefore a curve from the origin, rising steeply, reaching x=1x = 1 at t=t<π2t = t^* < \tfrac{\pi}{2}, and then remaining at x=1x = 1 (a horizontal line) for ttt \geq t^*.

Showing x=1x = 1 at a value less than π2\tfrac{\pi}{2} more formally:

The time to reach x=1x = 1 is t=01dx(1x3)1/3t^* = \int_0^1 \dfrac{dx}{(1-x^3)^{1/3}}. We need to show this is less than π2=01dx(1x2)1/2\tfrac{\pi}{2} = \int_0^1 \dfrac{dx}{(1-x^2)^{1/2}}.

Since (1x3)1/3>(1x2)1/2(1-x^3)^{1/3} > (1-x^2)^{1/2} for 0<x<10 < x < 1, we have 1(1x3)1/3<1(1x2)1/2\dfrac{1}{(1-x^3)^{1/3}} < \dfrac{1}{(1-x^2)^{1/2}} on this interval, so

t=01dx(1x3)1/3<01dx(1x2)1/2=π2t^* = \int_0^1 \frac{dx}{(1-x^3)^{1/3}} < \int_0^1 \frac{dx}{(1-x^2)^{1/2}} = \frac{\pi}{2}

as required.