Problem:
Solution:
This is a very long question, and we will tackle this part-by-part. This is for part (a)
Let's start with the formula:
xu=(β′(u)+vδ′(u))
xv=(δ(u))
xuu=(β″(u)+vδ″(u))
xuv=(δ′(u))
xvv=(0)
The observation is that n=U⋅xvv=U⋅0=0. So we have got the first equality.
Now we need to compute xu×xv
xu×xv=(β′(u)+vδ′(u))×δ(u)=β′(u)×δ(u)+vδ′(u)×δ(u)
So that explain the denominator, as per the hint. Finally we compute m
U⋅xuu=(β′(u)×δ(u)+vδ′(u)×δ(u))⋅δ′(u)=β′(u)×δ(u)⋅δ′(u). So we finally also explain the numerator!
Notice the numerator is not exactly the same form I had, but it is the same because it is a circular shift of the scalar triple product.
Now we moved on to part (b), we can parametrize (x,y,xy)=(u,0,0)+v(0,1,u) to make it a ruled surface.
The code for computing the Gaussian curvature is as follow:
syms u
syms v
beta = [u; 0; 0];
delta = [0; 1; u];
beta_u = diff(beta, u);
delta_u = diff(delta, u);
n = -(beta_u.' * cross(delta, delta_u))^2
D = cross(beta_u, delta) + v * cross(delta_u, delta);
d = (D.' * D)^2;
K = simplify(n/d)
We get the answer as −1(u2+v2+1)2.
For part (c) and (d), the Gaussian curvatures for both cone and cylinder are 0, this because either β′=0 or δ′=0.
For part (e), the helicoid has the parametrization as (vcosu,vsinu,bu)=(0,0,u)+v(cosu,sinu,0), so we use essentially the same code above except
syms b;
beta = [0;0;b*u];
delta= [cos(u);sin(u);0];
So we get the answer as −b2(b2+v2)2
For part (f), we will use the ruling patch we found in Exercise 2.1.22, so we simply put in yet another β and δ into the program and get the answer:
−a2b2c2(a2b2v2+a2c2(vcosu−sinu)2+b2c2(vsinu+cosu)2)2
For part (g), in some sense, we have already done with it, for the saddle z=xy is a hyperbolic paraboloid.
For a more general hyperbolic paraboloid, we consider
zc=y2b2−x2a2
Now let u=yb+xa and v=yb−xa
(a2(u−v),b2(u+v),cuv)
Now we obtain the ruled patch (a2u,b2u,0)+v(−a2,b2,cu)
Because we wanted the expression to show values in terms of x and y this time, so we modified the program a bit as follow:
% defining the ruled patch
syms a
syms b
syms c
beta = [ 0.5 * a * u; 0.5 * b * u; 0 ];
delta = [-0.5 * a ; 0.5 * b ; u/c];
% back substitute the x, y values
syms x
syms y
simplify(subs(subs(K, u, y/b + x/a), v, y/b - x/a))
So the answer −4a6b6c2(a4b4c2+4a4y2+4b4x2)2
Phew, finally!
Solution:
This is a very long question, and we will tackle this part-by-part. This is for part (a)
Let's start with the formula:
xu=(β′(u)+vδ′(u))
xv=(δ(u))
xuu=(β″(u)+vδ″(u))
xuv=(δ′(u))
xvv=(0)
The observation is that n=U⋅xvv=U⋅0=0. So we have got the first equality.
Now we need to compute xu×xv
xu×xv=(β′(u)+vδ′(u))×δ(u)=β′(u)×δ(u)+vδ′(u)×δ(u)
So that explain the denominator, as per the hint. Finally we compute m
U⋅xuu=(β′(u)×δ(u)+vδ′(u)×δ(u))⋅δ′(u)=β′(u)×δ(u)⋅δ′(u). So we finally also explain the numerator!
Notice the numerator is not exactly the same form I had, but it is the same because it is a circular shift of the scalar triple product.
Now we moved on to part (b), we can parametrize (x,y,xy)=(u,0,0)+v(0,1,u) to make it a ruled surface.
The code for computing the Gaussian curvature is as follow:
syms u
syms v
beta = [u; 0; 0];
delta = [0; 1; u];
beta_u = diff(beta, u);
delta_u = diff(delta, u);
n = -(beta_u.' * cross(delta, delta_u))^2
D = cross(beta_u, delta) + v * cross(delta_u, delta);
d = (D.' * D)^2;
K = simplify(n/d)
We get the answer as −1(u2+v2+1)2.
For part (c) and (d), the Gaussian curvatures for both cone and cylinder are 0, this because either β′=0 or δ′=0.
For part (e), the helicoid has the parametrization as (vcosu,vsinu,bu)=(0,0,u)+v(cosu,sinu,0), so we use essentially the same code above except
syms b;
beta = [0;0;b*u];
delta= [cos(u);sin(u);0];
So we get the answer as −b2(b2+v2)2
For part (f), we will use the ruling patch we found in Exercise 2.1.22, so we simply put in yet another β and δ into the program and get the answer:
−a2b2c2(a2b2v2+a2c2(vcosu−sinu)2+b2c2(vsinu+cosu)2)2
For part (g), in some sense, we have already done with it, for the saddle z=xy is a hyperbolic paraboloid.
For a more general hyperbolic paraboloid, we consider
zc=y2b2−x2a2
Now let u=yb+xa and v=yb−xa
(a2(u−v),b2(u+v),cuv)
Now we obtain the ruled patch (a2u,b2u,0)+v(−a2,b2,cu)
Because we wanted the expression to show values in terms of x and y this time, so we modified the program a bit as follow:
% defining the ruled patch
syms a
syms b
syms c
beta = [ 0.5 * a * u; 0.5 * b * u; 0 ];
delta = [-0.5 * a ; 0.5 * b ; u/c];
% back substitute the x, y values
syms x
syms y
simplify(subs(subs(K, u, y/b + x/a), v, y/b - x/a))
So the answer −4a6b6c2(a4b4c2+4a4y2+4b4x2)2
Phew, finally!
No comments:
Post a Comment