User Tools

Site Tools


math121b:final-sol

Final Solution

\gdef\E{\mathbb E} Due Date : May 10th (Sunday) 11:59PM. Submit online to gradescope.

Policy : You can use textbook and your notes. There should be no discussion or collaborations, since this is suppose to be a exam on your own understanding. If you found some question that is unclear, please let me know via email.


1. Vector spaces and Curvilinear Coordinates (30 pts)

All vectors spaces are finite dimensional over R\R.

1. True or False (10 pts)

  1. (F) Any vector space has a unique basis.
  2. (F) Any vector space has a unique inner product.
  3. (T) Given a basis e1,,ene_1, \cdots, e_n of VV, there exists a basis E1,,EnE^1, \cdots, E^n on VV^*, such that Ei(ej)=δijE^i (e_j) = \delta_{ij}.
  4. (F) If we change e1e_1 in the basis e1,,ene_1, \cdots, e_n, in the dual basis only E1E^1 will change.
    • The correct statement would be, “If we change e1e_1 in the basis e1,,ene_1, \cdots, e_n, in the dual basis, not only E1E_1 will change, all other EiE_i might also change.”
  5. (F) Given a vector space with inner product, there exists a unique orthogonal basis.
  6. (T) Let VV and WW be two vector spaces with inner products and of the same dimension.Then there exists a linear map f:VWf: V \to W, such that for any v1,v2Vv_1, v_2 \in V, $$\la v_1, v_2 \ra = \la f(v_1), f(v_2) \ra.$$
  7. (F) If VV and WW are vector spaces of dimension 33 and 55, then the tensor product VWV \otimes W have dimension 88. (should be 3×5=153 \times 5 = 15)
  8. (T) If VV has dimension 55, then the exterior power 3V\wedge^3 V is a vector space with dimension 1010.
  9. (T) The solution space of equation y(x)+x2y(x)=0y'(x) + x^2 y(x) = 0 forms a vector space.
  10. (F) The solution space of equation y(x)+xy2(x)=0y'(x) + x y^2(x) = 0 forms a vector space.

2. (10 pt) Let VnV_n be the vector space of polynomials whose degree is at most nn. Let f(x)f(x) be any smooth function on [1,1][-1,1]. We fix f(x)f(x) once and for all. Show that there is a unique element fnVnf_n \in V_n (depending on our choice of ff), such that for any gVng \in V_n, we have 11fn(x)g(x)dx=11f(x)g(x)dx. \int_{-1}^1 f_n(x) g(x) dx = \int_{-1}^1 f(x) g(x) dx.

Hint:

  • (1) Equip VnV_n with an inner product $\la g_1, g_2 \ra = \int_{-1}^1 g_1(x) g_2(x) dx$.
  • (2) Show that f(x)f(x) induces an element in VnV_n^*: g11f(x)g(x)dxg \mapsto \int_{-1}^1 f(x) g(x) dx.
  • (3) use inner product to identify VnV_n and VnV_n^*.

A remark: if f(x)f(x) were a polynomial of degree less than nn, then you could just take fn(x)=f(x)f_n(x) = f(x). But, we have limited our choices of fnf_n to be just degree n \leq n polynomial, so we are looking for a 'best approximation' of f(x)f(x) in VnV_n in a sense. Try solve the example case of n=1n=1, f(x)=sin(x)f(x) = \sin(x) if you need some intuition.

Solution: First some remark about what we are trying to prove: say you have an friend called AA who want to challenge you by playing a game:

  • AA provide nn and f(x)f(x) to you,
  • then you need to provide fn(x)f_n(x) to AA,
  • then AA will examine if your fnf_n pass the quality-check, i.e., AA take an arbitrary g(x)Vng(x) \in V_n, and test if 11g(x)f(x)dx=11g(x)fn(x)dx\int_{-1}^1 g(x) f(x) dx = \int_{-1}^1 g(x) f_n(x) dx.

Note that, when you produce fn(x)f_n(x), you have no knowledge of what g(x)g(x) would be.

Here is a solution, that is of a concrete flavor, not quite following the hint. First, we define an inner product $\la g_1, g_2 \ra = \int_{-1}^1 g_1(x) g_2(x) dx$, whenever the integral make sense. Then, we may take an orthonormal basis e0,,ene_0, \cdots, e_n of VnV_n (note dimVn=n+1\dim V_n= n+1), and then $$ f_n(x) = \sum_{i=0}^n \la f, e_i \ra e_i. $$ Then, we have for any gVng \in V_n, $$ \la f_n, g \ra = \sum_{i=0}^n \la f, e_i \ra \la e_i, g \ra = \la f, \sum_{i=0}^n \la e_i, g \ra e_i \ra = \la f, g \ra. $$

Conceptually, if VV denote the space of smooth functions on [0,1][0,1] with inner product, then VVnV \supset V_n, then there is an orthogonal projection Πn:VVn.\Pi_n: V \to V_n. You have seen this orthogonal projection in different guises, for example the least square regression, the truncation of Fourier series expansion of some function, … Here fn=Πn(f)f_n = \Pi_n(f).

3. (10 pt) Let R3\R^3 be equipped with curvilinear coordinate (u,v,w)(u,v,w) where u=x,v=y,w=zx2+y2. u = x, v = y, w = z - x^2 + y^2.

  1. (3pt) Write the vector fields u,v,w\d_u, \d_v, \d_w in terms of x,y,z\d_x, \d_y, \d_z.
  2. (3pt) Write the 1-forms (co-vector fields) du,dv,dwdu,dv,dw in terms of dx,dy,dzdx, dy, dz.
  3. (4pt) Write down the standard metric of R3\R^3 in coordinates (u,v,w)(u,v,w).

Solution: We write x,y,zx,y,z in terms of u,v,wu,v,w x=u,y=v,z=w+u2v2 x = u, y = v, z = w + u^2 - v^2 then u=u(x)x+u(y)y+u(z)z=x+2uz=x+2xz\d_u = \d_u(x) \d_x + \d_u(y) \d_y + \d_u(z) \d_z = \d_x + 2u \d_z = \d_x + 2x \d_z The others are similar.

dw=dz2xdx+2ydydw = dz - 2x dx + 2y dy

Then finally g=(dx)2+(dy)2+(dz)2=(du)2+(dv)2+(dz2xdx+2ydy)2 g = (dx)^2 + (dy)^2 + (dz)^2 = (du)^2 + (dv)^2 + (dz - 2x dx + 2y dy)^2 open up the parenthesis if you wish.

2. Special Functions and Differential Equations (50 pts)

1. (10 pt) Orthogonal polynomials. Let I=[1,1]I = [-1,1] be a closed interval. w(x)=x2w(x) = x^2 a non-negative function on II. For functions f,gf,g on II, we define their inner products as $$ \la f, g \ra = \int_{-1}^1 f(x) g(x) w(x) dx $$ The normalized orthogonal polynomials P0,P1,P_0, P_1, \cdots are defined by

  1. Pn(x)P_n(x) is a degree nn polynomial.
  2. $\la P_n, P_n \ra = 1$
  3. $\la P_i, P_j \ra = 0$ if iji \neq j.

Find out P0,P1,P2P_0, P_1, P_2.

P0(x)=±3/2,P1(x)=±5/2x,P2(x)=±14/4(3+5x2).P_0(x) = \pm \sqrt{3/2}, \quad P_1(x) = \pm \sqrt{5/2} x, \quad P_2(x) = \pm \sqrt{14}/4(-3+5 x^2).

2. (10 pt) Find eigenvalues and eigenfunctions for the Laplacian on the unit sphere S2S^2, i.e., solve ΔF(θ,φ)=λF(θ,φ) \Delta F(\theta, \varphi) = \lambda F(\theta, \varphi) for appropriate λ\lambda and FF. The Laplacian on a sphere is Δf=1sinθθ(sinθθ(f))+1sin2θφ2f. \Delta f = \frac{1}{\sin \theta} \d_\theta(\sin \theta \d_\theta(f)) + \frac{1}{\sin^2 \theta} \d_\varphi^2 f.

Solution: Eigenvalue λ=l(l+1)\lambda = -l(l+1) and eigenfunctions F(\theta, \varphi) = P_l^m(\cos \theta) \cos (m \varphi) , P_l^m(\cos \theta) \sin (m \varphi)

3. (5 pt) Find eigenvalues and eigenfunctions for the Laplacian on the half unit sphere S2S^2 with Dirichelet boundary condition, i.e., solve ΔF(θ,φ)=λF(θ,φ),F(θ=π/2,φ)=0. \Delta F(\theta, \varphi) = \lambda F(\theta, \varphi), \quad F(\theta=\pi/2, \varphi)=0. for appropriate λ\lambda and FF.

Here the trick is that, any eigenfunction on the upper-semisphere, by reflection, can be extended to an eigenfunction on the whole sphere F(π/2θ,φ)=F(π/2+θ,φ) F(\pi/2 - \theta, \varphi) = - F(\pi/2 + \theta, \varphi) hence, we want those eigenfunction on the whole sphere that satisfies P_l^m(x) = - P_l^m(-x) this turns out to be satisfies if l+ml+m is odd.

Note that, even though the boundary condition is rotational symmetric, it does not mean the solution is rotational symmetry (after all, the S2S^2 itself is symmetric, but the eigenfunction can have fluctuations).

4. (15 pt) (Heat flow). Consider heat flow on the closed interval [0,1][0,1] tu(x,t)=x2u(x,t), \d_t u(x,t) = \d_x^2 u(x,t), where u(x,t)u(x,t) denote the temperature. \ Let u(0,t)=u(1,t)=0u(0, t) = u(1, t) = 0 for all tt. Let the initial condition be u(x,0)={2xx[0,1/2]2(1x)x[1/2,1] u(x, 0) = \begin{cases} 2x & x \in [0, 1/2] \cr 2(1-x) & x \in [1/2, 1] \end{cases}

  • (12pt) Solve the equation for t>0t > 0.
  • (3 pt) Does the solution make sense for any negative tt? Why or why not?

The problem is standard, I won't repeat the solution.

The solution will diverge for any negative tt, no matter how small t|t| is, since ncnen2t=ncnen2t \sum_n c_n e^{-n^2 t} = \sum_n c_n e^{n^2 |t|} will diverge quite fast due to en2e^{n^2}, and cnc_n is only decaying as 1/nc1/n^c for some constant cc.

Note that, how much you can go negative in time, depends on how smooth the initial condition is. Here the inital condition is already non-smooth, hence you cannot extend the solution to t(ϵ,+)t \in (-\epsilon, +\infty) from t(0,+infty)t \in (0, +infty).

5. (10 pt) (Steady Heat equation). Let DD be the unit disk. We consider the steady state heat equation on DD Δu(r,θ)=0 \Delta u(r, \theta) = 0

  • (3 pt) Write down the Laplacian Δ\Delta in polar coordinate
  • (5 pt) Show that, if the boundary value is u(r=1,θ)=0u(r=1, \theta) = 0, then u=0u=0 on the entire disk.
  • (2 pt) Is it possible to have a boundary condition u(r=1,θ)=f(θ)u(r=1, \theta) = f(\theta), such that there are two different solutions u1(r,θ)u_1(r,\theta) and u2(r,θ)u_2(r,\theta) to the problem?

It is impossible to have a boundary condition u(r=1,θ)=f(θ)u(r=1, \theta) = f(\theta), such that there are two different solutions u1(r,θ)u_1(r,\theta) and u2(r,θ)u_2(r,\theta) to the problem. Otherwise, let u1u_1 and u2u_2 be the two solutions, and we may take their difference u=u1u2 u = u_1 - u_2 then Δu=0,uD=0 \Delta u = 0, \quad u|_{\d D} = 0 by part (b), u=0u=0.

3. Probability and Statistics (20 pts)

1. (5 pt) Throw a die 100 times. Let XX be the random variable that denote the number of times that 44 appears. What distribution does XX follow? What is its mean and variance?

Binomial distribution, with n=100,p=1/6n=100, p=1/6. Mean is npnp, variance is np(1p)np(1-p).

2. (5 pt) Let XN(0,1)X \sim N(0,1) be a standard normal R.V . Compute its moment generating function E(etX). \E(e^{t X}). Use the moment generating function to find out E(X4)\E(X^4). Let Y=X2Y = X^2. What is the mean and variance of YY?

Do the integral, we get E(etX)=et2/2\E(e^tX) = e^{t^2/2}.

To compute E(X4)\E(X^4), we note that E(etX)=m0+m1t+t22!m2+t33!m3+t44!m4+ \E(e^tX) = m_0 + m_1 t + \frac{t^2}{2!} m_2 + \frac{t^3}{3!} m_3+ \frac{t^4}{4!} m_4 + \cdots where E(Xk)=mk\E(X^k) = m_k. Hence, we may get the coefficients of Taylor expansion et2/2=1+(t2/2)+(t2/2)22!+=1+t2/2+t4/8+ e^{t^2/2} = 1 + (t^2/2) + \frac{(t^2/2)^2}{2!} + \cdots = 1 + t^2/2 + t^4/8 + \cdots comparing the t4t^4 coefficients, we see m4/4!=1/8m4=3 m_4/4! = 1/8 \Rightarrow m_4 = 3

3. (5 pt) There are two bags of balls. Bag A contains 4 black balls and 6 white balls, Bag B contains 10 black balls and 10 white balls. Suppose we randomly pick a bag (with equal probability) and randomly pick a ball. Given that the ball is white, what is the probability that we picked bag A?

Note that P(white ball)=P(whitebag A)P(bag A)+P(whitebag B)P(bag B)=(6/10)(1/2)+(10/20)(1/2)P(\z{white ball}) = P(\z{white} | \z{bag A}) P( \z{bag A}) + P(\z{white} | \z{bag B}) P( \z{bag B}) = (6/10)(1/2) + (10/20)(1/2) instead of total number of white ball divided by total number of balls.

4. (5 pt) Consider a random walk on the real line: at t=0t=0, one start at x=0x=0. Let SnS_n denote the position at t=nt=n, then Sn=Sn1+XnS_n = S_{n-1} + X_n, where Xn=±1X_n = \pm 1 with equal probability.

  • (3pt) What is the variance of SnS_n?
  • (2pt) Use Markov inequality, prove that for any c>1c > 1, we have

P(Sn>cn)1/c2 \P(|S_n| > c \sqrt{n}) \leq 1/c^2

Since Sn=X1++XnS_n = X_1 + \cdots + X_n and XiX_i are independent, we have Var(Sn)=i=1nVar(Xi)=n(12(1/2)+(1)2(1/2)=nVar(S_n) = \sum_{i=1}^n Var(X_i) = n (1^2 (1/2) + (-1)^2 (1/2) = n

Then, by Markov inequality, we have P(X>cVar(X))1c2 \P(|X| > c \sqrt{Var(X)}) \leq \frac{1}{c^2} take X=SnX = S_n.

math121b/final-sol.txt · Last modified: 2020/05/17 18:55 by pzhou