Problem:
Show that the Gaussian curvature of the hyperboloid of one sheet x(u,v)=(acoshucosv,bcoshusinv,csinhu) maybe written in Cartesian coordinate as
K=−1a2b2c2[x2a4+y2b4+z2c4]2
Solution:
This code should verify the correctness of the identity. It would be too tedious for manual computation.
clear;
clc;
syms a;
syms b;
syms c;
syms x;
syms y;
z = sqrt(c^2 * (x^2/a^2 + y^2/b^2 - 1));
s = [x; y; z];
sx = diff(s, x);
sy = diff(s, y);
E = simplify(sx.' * sx);
F = simplify(sx.' * sy);
G = simplify(sy.' * sy);
normal = cross(sx, sy);
normal_norm = simplify(sqrt(normal.' * normal));
u = simplify(normal/normal_norm(1));
sxx = simplify(diff(sx, x));
sxy = simplify(diff(sx, y));
syy = simplify(diff(sy, y));
l = simplify(u.' * sxx);
m = simplify(u.' * sxy);
n = simplify(u.' * syy);
K = simplify((l * n - m * m)/(E * G - F * F))
Now we get from the program that:
K=−c2a6b6(c2x2b4+c2y2a4+a2b4x2+a4b2y2−a4b4)2
PK represent the K value we need to prove, and the program output 0 meaning the K we found matches the PK we need to reach!
Show that the Gaussian curvature of the hyperboloid of one sheet x(u,v)=(acoshucosv,bcoshusinv,csinhu) maybe written in Cartesian coordinate as
K=−1a2b2c2[x2a4+y2b4+z2c4]2
Solution:
This code should verify the correctness of the identity. It would be too tedious for manual computation.
clear;
clc;
syms a;
syms b;
syms c;
syms x;
syms y;
z = sqrt(c^2 * (x^2/a^2 + y^2/b^2 - 1));
s = [x; y; z];
sx = diff(s, x);
sy = diff(s, y);
E = simplify(sx.' * sx);
F = simplify(sx.' * sy);
G = simplify(sy.' * sy);
normal = cross(sx, sy);
normal_norm = simplify(sqrt(normal.' * normal));
u = simplify(normal/normal_norm(1));
sxx = simplify(diff(sx, x));
sxy = simplify(diff(sx, y));
syy = simplify(diff(sy, y));
l = simplify(u.' * sxx);
m = simplify(u.' * sxy);
n = simplify(u.' * syy);
K = simplify((l * n - m * m)/(E * G - F * F))
PK = -1/(a^2 * b^2 * c^2 * (x^2 / a^4 + y^2 / b^4 + z^2/c^4)^2);
simplify(K - PK)
K=−c2a6b6(c2x2b4+c2y2a4+a2b4x2+a4b2y2−a4b4)2
PK represent the K value we need to prove, and the program output 0 meaning the K we found matches the PK we need to reach!
No comments:
Post a Comment