MAT 3310 (2009-2010, Second Term) Suggested Solution For Assignment 2
MAT 3310 (2009-2010, Second Term) Suggested Solution For Assignment 2
MAT 3310 (2009-2010, Second Term) Suggested Solution For Assignment 2
(1)
y1 0,
(2)
y2 0.
(3)
y1 = 4 + 2y1 = 0
L
y = 4 + 2y1 = 0
L1
2
2
= y1 + 4y2 1 = 0
Solve the system, y1 =
2 5
5 ,
y2 =
5
10 ,
5.
2 5
5
(4)
+4
5
10
= 2 5.
2. For point A (y1 , y2 , y3 ) on the plane, if the distance from the point P (1, 1, 1) to A is
3. (a) Let A = (1, 2, 3)T , Ax and ~y are projections of vector ~a = (1, 1, 1)T on the line
y1 = y22 = y33 and its vertical line respectively, then
~a = Ax + ~y
(6)
(Ax) ~y = 0
(7)
I A
~y
~a
=
AT 0
x
0
solve the system, we get ~y = ( 47 , 17 , 27 )T , x = 73 .
the minimum distance from point (1, 1, 1) to the line is k~y k =
21
7 .
1 2 3
1
2 4 6
P = A(AT A)1 AT =
14
3 6 9
(c) By the definition of P , P T = A(AT A)1 AT = P ;
and P 2 = P P = A(AT A)1 AT A(AT A)1 AT = A(AT A)1 (AT A)(AT A)1 AT =
A(AT A)1 AT = P .
For any vector ~b,
kP bk2 = (P b)T (P b) = bT P T P b = bT P P b = bT P b
2
k(I P )bk = kb P bk = (b P b) (b P b)
= bT b 2bT P b + (P b)T (P b)
2
(8)
(9)
(10)
= kbk 2b P b + kP bk
(11)
(12)
= kbk kP bk
(13)
1
U = (K~u, ~u) = Pmin
2
2
(14)
(15)
(16)
By w = CA~u, K = AT CA,
1
Qmin = (C 1 w, w),
2
1
= (A~u, CA~u)
2
1
= (AT CA~u, ~u)
2
1
= (K~u, ~u) = U
2
(17)
(18)
(19)
(20)
(b) At equilibrium,
1
1
U = (K~u, ~u) = (f~, ~u)
2
2
U is NOT equal to the external potential energy (f~, ~u)!
If they are equal, then U = (f~, ~u) = 0, which means U = 21 (K~u, ~u) = 0 by (a).
Since K is positive symmetric matrix, ~u = ~0 and f~ = K~u = ~0, which means the bar
is not stretched, and there are no external force, it is impossible when the f~ is the
gravitational forces on the masses.
5. (a)
"
#
d(u + tv) 2
dx
dx
0
2 #
2
Z "
1 1 d2 u
du dv
dv
d2 v
du
dx
=
+ 2t
+ t2
+t 2 +
2
2 0 dx
dx
dx
dx dx
dx
Z
Z
t2 1 dv 2
t 1 d2 v
du dv
= P (u) +
dx +
+
2
dx
2 0
dx
2 0 dx2
dx dx
1
P (u + tv) =
2
d2 (u + tv)
+
dx2
(21)
(22)
(23)
For fixed u and v, we can treat P (u + tv) as a quadratic function on t. The first
variation P
u is achieved at
P (u + tv) P (u)
1
P (u) v = lim
=
t0
t
2
1 2
d v
du dv
+2
dx
dx2
dx dx
u00 (x) = 3
u(0) = 1
Z 1 2
du dv
d v
0
+2
dx = 0
2P (u) v =
dx2
dx dx
0
Z 1
(v 00 + 2u0 v 0 )dx = 0
0
(24)
(25)
(26)
(27)
Integration by part,
Z
Z
(28)
(29)
0
1
for any candidate (u + tv), from (u + tv)00 (x) = 3, and u00 (x) = 3, we know v 00 (x) = 0,
then v 0 (x) C, C is some constant. The above equation now become:
0 + 2C(u(1) u(0)) = 0
(30)
u(1) = u(0) = 1
(31)