Loading [MathJax]/jax/output/HTML-CSS/jax.js

2020年8月31日 星期一

國立臺灣大學一百學年度研究所碩士班入學考試試題:微積分乙(不含線性代數)

  1. (30%) Suppose a random value y follows the g distribution, i.e., y can be expressed as a function of a normally distributed variable as follows.

    y=δ+λexp(gz)1g,

    where δ and g are real numbers, λ is a positive real number, and z follows the normal distribution with mean to be μ and standard deviation to be σ. Given the constraint [g(yδ)/λ+1]>0, solve the following problems.
    1. (4%) Express z in terms of y, i.e., find the function h(y)=z.
    2. (4%) With the function h(y), derive the density function of y, f(y), according to the definition f(y)=1σ2π|h(y)|e12σ2(h(y)μ)2, where μ and σ are the mean and the standard deviation of z.
    3. (8%) Derive the expectation of y, E(y).
    4. (8%) Derive the variance of y, var(y).
    5. (6%) Show that when g is zero, y follows a normal distribution. [Hint: Consider the result of limg0exp(gz)1g.]
  2. 訣竅前兩小題根據題目給定的資訊整理即可,而後三小題則分別按照期望值、變異數與常態分配的定義處理。
    解法
    1. 移項便有

      g(yδ)λ+1=exp(gz)

      由於 [g(yδ)λ+1]>0,故可同取自然對數便有 gz=ln(1+g(yδ)λ),故得

      z=1gln(1+g(yδ)λ):=h(y)

    2. 先求 h 的導函數有

      h(y)=1ggλ+g(yδ)=1λ+g(yδ)

      由於 g(yδ)/λ+1>0,故 λ+g(yδ)>0,故 |h(y)|=h(y)。從而 y 的機率密度函數為

      f(y)=1σ2π1λ+g(yδ)exp[12σ2(ln(λ+g(yδ))lnλgμ)2]

      其中 y(y,)y=λ/g+δ
    3. 那麼由期望值的定義可知

      E(y)=yyf(y)dy

      t=h(y)μσ2,那
      • yy 時有 t
      • y 時有 t
      • 求導有 dt=h(y)σ2dy。再者整理有 y=h1(μ+tσ2)
      如此所求為

      E(y)=1πh1(μ+tσ2)et2dt=1π[δλg+λgexp(gμ+tgσ2)]et2dt=1π(δλg)et2dt+λegμgπexp(t2+tgσ2)dt=δλg+λegμ+g2σ2/2g

    4. 按照變異數的定義有 var(y)=E(y2)(E(y))2,故我們僅需計算 E(y2),即計算

      E(y2)=yy2f(y)dy

      沿用前一小題的代換可知

      E(y2)=1π[δλg+λgexp(gμ+tgσ2)]2et2dt=1π(δλg)2et2dt+2λegμgπ(δλg)exp(t2+tgσ2)dt+λ2e2gμg2πexp(t2+2tgσ2)dt=(δλg)2+2λegμ+g2σ2/2g(δλg)+λ2e2gμ+2g2σ2g2

      因此所求為

      var(y)=[(δλg)2+2λegμ+g2σ2/2g(δλg)+λ2e2gμ+2g2σ2g2](δλg+λegμ+g2σ2/2g)2=λ2g2e2gμ+g2σ2(eg2σ21)

    5. 由於當 g 趨於零時可以發現 limg0exp(gz)1g=limg0zexp(gz)1=z。故當 g=0 時有 y=δ+λz,此表明經由線性轉換能使 y 成為平均為 δ+λμ,而標準差則為 λσ。由 (c) 與 (d) 能知 limg0E(y)=δ+λμlimg0var(y)=λ2σ2

  3. (8%) Find the maximum and minimum values of f(x,y,z)=x+3yz subject to 4x2+2y2+z2=4.
  4. 訣竅運用初等不等式求極值;亦可使用拉格朗日乘子法求解。
    解法一由柯西不等式可知

    23=[(2x)2+(2y)2+z2][(12)2+(32)2+(1)2](x+3yz)2

    因此函數 f 在給定的限制條件下的最大最小值分別為 2323,此時的等號成立條件為 4x=2y3=z,即 (123,623,423)(123,623,423)
    解法二設定拉格朗日乘子函數如下

    F(x,y,z,λ)=x+3yz+λ(4x2+2y2+z24)

    據此解下列的聯立方程組

    {Fx(x,y,z,λ)=1+8λx=0Fy(x,y,z,λ)=3+4λy=0Fz(x,y,z,λ)=1+2λz=0Fλ(x,y,z,λ)=4x2+2y2+z24=0

    由於明顯 λ0,故 x=18λy=34λz=12λ,由此代入第四式為

    464λ2+2916λ2+14λ2=4

    即得 λ2=2364,故有 λ=±238。如此得兩座標

    (x,y,z)=(1,6,4)23, (x,y,z)=(1,6,4)23

    而最大最小值分別為 2323

  5. (12%) A 10,000-cubic-foot-room has 500 smoke particles per cubic foot. A ventilation system is turned on that each minute brings in 500 cubic feet of smoke-free air, while an equal volume of air leaves the room. Also, during each minute, smoker in the room add a total of 10,000 particles of smoke to the room. Assume that the air in the room mixes thoroughly.
    1. (4%) Find a differential equation and initial condition that govern the total number y(t) of smoke particles in the room after t minutes.
    2. (4%) Solve this differential equation.
    3. (4%) Find how soon the smoke level will fall to 100 smoke particles per cubic foot. (If the answers are with decimal numbers, please round to the nearest hundredth.)
  6. 訣竅按照濃度的變化率寫出出微分方程並求解。
    解法
    1. 由於 y(t) 表在時刻 t 分鐘後的煙霧粒子的個數,那麼 y(t) 便表示在 t 分鐘後增加煙霧粒子的瞬間增加量,按照題意可知

      y(t)=10000y(t)10000500

      再者一開始房間內有 500104 個煙霧粒子,因此 y(0)=5106
    2. 經由移項可知

      dy200000y=dt20

      同取積分有

      ln|200000y|=t20+C

      透過初始條件可知 C=ln4800000,從而所求為

      y(t)=200000+4800000et/20

    3. 要達到濃度為每立方英尺有 100 顆煙霧粒子,則該房間要有總數 106 顆煙霧粒子。因此解方程式 y(t)=106,即 2+48et/20=10,故有 t=20ln6

  7. Evaluate:
    1. (5%) ϕ(α)=0αeαxdx for α>0.
    2. (10%) Evaluate 10ex2dx, approximately by using Taylor's theorem of the mean and estimate the maximum error.
  8. 訣竅直接演算瑕積分並透過帶餘項的泰勒定理估算瑕積分並給出誤差估計。
    解法
    1. 直接計算瑕積分可知

      ϕ(α)=eαx|0=1

    2. 回憶自然指數的泰勒展開式為

      ex=nk=0xkk!+eξ(n+1)!xn+1

      此處 ξ 介於 0x 之間。現用 x2 取代 x

      ex2=nk=0x2kk!+eη2(n+1)!x2n+2

      而此時的 η 介於 0x 之間。由此計算定積分可知其近似值為

      10ex2dx10nk=0x2kk!dx=nk=01k!(2k+1)

      而其最大的誤差估計為

      |10ex2dxnk=01k!(2k+1)|10eη2x2n+2(n+1)!dx=eη2(n+1)!(2n+3)e(n+1)!(2n+3)

      這是因為 η[0,x][0,1],因此 eη2e

  9. Let φ be a continuous real function on [0,α]. Assume that φ(x)>0 if 0<xα and that φ(x)Axr as x0 (A>0, r0). For all t>0 set

    F(t)=α0dxt+φ(x).

    1. (10%) Show that if r>1, then as t0

      F(t)πrA1/rsin(π/r)1t11/r.

    2. (10%) Show that if r=1, then as t0

      F(t)1Alog1t.

  10. 訣竅由於題目所提供的資訊為在原點附近的漸近行為,故適當的將給定的積分範圍分割後以便計算。
    解法
    1. 給定正數 ε<1,那麼存在 δ>0使得 0<x<δ 蘊含 Axr(1ε)<φ(x)<Axr(1+ε)。據此我們可寫 F 如下

      F(t)=δ0dxt+φ(x)+αδdxt+φ(x)

      對於第二項可得如下的估算

      0<αδdxt+φ(x)αδdxmin[δ,α]φ=αδmin[δ,α]φ<

      故有 t1 使得 0<t<t1 蘊含

      αδ1t+φ(x)<ε2rA1/r(1+ε)1/rt11/r

      至於第一項則容易知道

      1tδ0dx1+Axr(1+ε)/t=δ0dxt+Axr(1+ε)<δ0dxt+φ(x)<δ0dxt+Axr(1ε)=1tδ0dx1+Axr(1ε)/t

      首先我們處理積分 δ0dx1+Axr(1+ε)/t,令 y=1+Axr(1+ε)/t,那麼有
      • x=0 時有 y=1
      • x=δ 時有 y=1+Aδr(1+ε)/t
      • 整理有 x=(t(y1)A(1+ε))1/r,求導有 dx=1r(tA(1+ε))1/r(y1)(1r)/rdy
      如此可知

      δ0dx1+Axr(1+ε)/t=t1/rrA1/r(1+ε)1/r1+Aδr(1+ε)/t11y(y1)11/rdy=t1/rrA1/r(1+ε)1/r[11y(y1)11/rdy1+Aδr(1+ε)/t1y(y1)11/rdy]

      由於 limy1y(y1)11/r1y=0,故瑕積分 1+Aδr(1+ε)/t1y(y1)11/rdy 有限,從而 limt0+1+Aδr(1+ε)/t1y(y1)11/rdy=0,即當 t<t2 時有

      1+Aδr(1+ε)/t1y(y1)11/rdy<ε2

      另一方面,透過變換 u=y1 可知

      11y(y1)11/rdy=0u1/r11+udu

      運用複變函數論中的留數定理與路徑積分(可參見「複變數-原理及題解」,中央圖書出版社,吳新生編著,第 403 頁)可知

      0u1/r11+udu=πsin(π/r)

      綜合上述的推導可得

      1rA1/r(1+ε)1/r(πsin(π/r)ε)1t11/r<F(t)

      類似地,我們也有

      F(t)<1rA1/r(1ε)1/r(πsin(π/r)+ε)1t11/r

      將兩者結合可知

      1rA1/r(1+ε)1/r(πsin(π/r)ε)<t11/rF(t)<1rA1/r(1ε)1/r(πsin(π/r)+ε)

      故取由 t0+ε 的任意性可知

      limt0+t11/rF(t)=πrA1/rsin(π/r)

      這就給出了

      F(t)πrA1/rsin(π/r)1t11/r

    2. 承前述步驟可知

      δ0dxt+Ax(1+ε)<δ0dxt+φ(x)<δ0dxt+Ax(1ε)

      直接計算可給出

      1A(1+ε)log(1+Aδ(1+ε)t)<δ0dxt+φ(x)<1A(1ε)log(1+Aδ(1ε)t)

      同乘以 Alog(1/t)

      11+εlog(1+Aδ(1+ε)/t)log(1/t)<Alog(1/t)δ0dxt+φ(x)<11εlog(1+Aδ(1ε)/t)log(1/t)

      t 趨於零且由 ε 的任意性可知 limt0+Alog(1/t)δ0dxt+φ(x)=1,從而有

      limt0+Alog(1/t)F(t)=1

      這就給出了

      F(t)1Alog1t


  11. Suppose that anc as n and that {αi}i=1 is a sequence of positive terms for which ni=1ai as n.
    1. (10%) Show that

      ni=1αiaini=1αic as n.

      In particular, if αi=1 for all i, then

      1nni=1aic as n.

    2. (5%) Show that the converse of the special case in (a) does not always hold by giving a counterexample of a sequence {an}n=1 that does not converges, yet 1nni=1ai converges as n.
  12. 訣竅運用Stolz–Cesàro定理(一種數列版本的羅必達法則);對於反例可考慮經典的發散數列。
    解法
    1. bn=ni=1αicn=ni=1αiai,按題設可知 limnbn=。再者,有

      limncn+1cnbn+1bn=limnαn+1an+1αn+1=limnan+1=c

      那麼根據Stolz-Cesàro定理可知 limncnbn=c,此即所欲證。而當 αi=1,那麼容易知道 bn=ncn=ni=1ai
    2. 考慮 an=(1)n,容易直接計算

      ni=1ai={1,if n is odd,0,if n is even.

      由此可知 limn1nni=1ai=0,但數列 {ai}i=1 發散。

沒有留言:

張貼留言