Ana Sayfa
Matematik

Diferansiyel Geometri Serisi 9. Bölüm

Bu bölümde birinci ve ikinci temel formdan ve eğrilikten bahsettik.

Gauss16 Mart 202610 dk okuma süresi
Diferansiyel Geometri Serisi 9. Bölüm

Yüzeylerin Sırları: İçsel Metrikler ve Dışsal Bükülmeler

Giriş: Karıncanın ve Dev'in Bakış Açısı

Herkese tekrar merhaba! Kahveleriniz hazırsa, diferansiyel geometrinin en heyecanlı virajlarından birini almaya hazırlanıyoruz. Geçen bölümde yüzeylere "yamalar" yapıştırıp onları tanımlamıştık. Ama bir sorunumuz var: Bu yüzeyin üzerinde yaşayan, dış dünyayı göremeyen bir "karınca" ne kadar yol gittiğini, üçgenlerinin iç açılarının toplamını veya bulunduğu alanın büyüklüğünü nasıl bilecek?

Biz dışarıdan bakan "devler" (gözlemciler) olarak 3 boyutlu uzayda (E3\mathbb{E}^3) cetvelimizi kullanıp ölçüm yapabiliriz. Ama diferansiyel geometri, sadece dışarıdan bakmak değil, yüzeyin içsel karakterini de anlamaktır. Bu bölümde hem karıncanın cetvelini (Birinci Temel Form) hem de dışarıdan bakan devin gözlemlerini (İkinci Temel Form) birleştireceğiz.

Birinci Temel Form: Yüzeyde Mesafe Ölçmek

Hatırlayın, bir SS yüzeyimiz ve üzerindeki bir pp noktasında teğet düzlemimiz TpST_pS vardı. Bu düzlemi xu\mathbf{x}_u ve xv\mathbf{x}_v baz vektörleri geriyordu. Teğet düzlemdeki herhangi bir w=axu+bxvw = a \mathbf{x}_u + b \mathbf{x}_v vektörünün boyunu hesaplamak istersek, kendisiyle iç çarpımını alırız:

w2=axu+bxv,axu+bxv|w|^2 = \langle a \mathbf{x}_u + b \mathbf{x}_v, a \mathbf{x}_u + b \mathbf{x}_v \rangle

Parantezleri açtığımızda karşımıza Gauss'un o meşhur üç katsayısı çıkar:

Tanım (Birinci Temel Form Katsayıları)

E=xu,xu,F=xu,xv,G=xv,xvE = \langle \mathbf{x}_u, \mathbf{x}_u \rangle, \quad F = \langle \mathbf{x}_u, \mathbf{x}_v \rangle, \quad G = \langle \mathbf{x}_v, \mathbf{x}_v \rangle

Bu katsayılar yüzeyin o noktadaki geometrisini (metriğini) belirler. Yüzey üzerindeki diferansiyel uzunluk karesi (ds2ds^2) şu formülle ifade edilir:

Ip(w)=Edu2+2Fdudv+Gdv2I_p(w) = E du^2 + 2F dudv + G dv^2

Bu formül sayesinde, x,y,zx, y, z koordinatlarına hiç bulaşmadan, sadece parametreler (u,vu, v) ve bu katsayılar ile yüzey üzerinde yay uzunluğu ve alan hesaplayabiliriz. Buna İçsel Geometri (Intrinsic Geometry) denir.

Şaşırtıcı Bir Sonuç: Silindir Aslında Düzdür!

Şimdi beyin yakan kısma, yani "karşıt örneklere" geliyoruz. İki farklı yüzeyi ele alacağız ve Birinci Temel Formlarının bize ne fısıldadığına bakacağız.

Düzlem (Referans Noktamız)

Düzlemi x(u,v)=(u,v,0)\mathbf{x}(u, v) = (u, v, 0) olarak parametrize edelim.

  • xu=(1,0,0)    E=1\mathbf{x}_u = (1, 0, 0) \implies E = 1
  • xv=(0,1,0)    G=1\mathbf{x}_v = (0, 1, 0) \implies G = 1
  • xuxv=0    F=0\mathbf{x}_u \cdot \mathbf{x}_v = 0 \implies F = 0

Metrik: ds2=du2+dv2ds^2 = du^2 + dv^2 (Bildiğimiz Pisagor!)

Silindir (Counter-Intuitive Example)

Yarıçapı 1 olan bir silindiri düşünün: x(u,v)=(cosu,sinu,v)\mathbf{x}(u, v) = (\cos u, \sin u, v).

  • xu=(sinu,cosu,0)    E=sin2u+cos2u=1\mathbf{x}_u = (-\sin u, \cos u, 0) \implies E = \sin^2 u + \cos^2 u = 1
  • xv=(0,0,1)    G=1\mathbf{x}_v = (0, 0, 1) \implies G = 1
  • F=0F = 0

Metrik: ds2=du2+dv2ds^2 = du^2 + dv^2

BİR DAKİKA! Düzlem ile Silindirin Birinci Temel Formları AYNI çıktı! 
Bu şu demek: İçsel geometri açısından Silindir ile Düzlem arasında HİÇBİR fark yoktur. Bir karıncayı silindirin üzerine koyun, bir diğerini de masanın üzerine. İkisi de üçgenler çizerse iç açıları toplamını 180 derece bulur. Silindir, bükülmüş olabilir ama eğrilmiş (curved) değildir. Bir kağıdı büküp rulo (silindir) yapabilirsiniz ama yırtmadan küre yapamazsınız.

Küre (Gerçek Eğrilik ve Haritacıların Çilesi)

Yarıçapı RR olan kürede işler değişir: x(u,v)=(Rsinucosv,Rsinusinv,Rcosu)\mathbf{x}(u, v) = (R \sin u \cos v, R \sin u \sin v, R \cos u)
Hesaplarsanız E=R2,F=0E=R^2, F=0 ama G=R2sin2uG=R^2 \sin^2 u bulursunuz. 
Katsayılar konuma (uu) göre değişiyor! Bu yüzden düzlemdeki geometri burada işlemez. Küreyi düzleme, mesafeleri bozmadan açamazsınız.

İkinci Temel Form: Dışarıdan Bakış ve Bükülme

Silindir ve düzlemin içsel olarak aynı olduğunu gördük. Ama göz var nizam var; silindir uzayda kıvrılıyor! İşte bu "bükülme" farkını anlamak için yüzeyin dışına çıkıp, yüzeyin normal vektörünün (NN) değişimine bakmamız gerekir.

Yüzeyin ikinci mertebeden türevlerini (xuu,xuv,xvv\mathbf{x}_{uu}, \mathbf{x}_{uv}, \mathbf{x}_{vv}) alıp, bunları yüzeyin birim normali (NN) üzerine izdüşürürsek, yüzeyin o yönde "teğet düzleminden ne kadar saptığını" buluruz.

Tanım (İkinci Temel Form Katsayıları)

L=xuu,N,M=xuvN,L = \langle \mathbf{x}_{uu}, N \rangle, \quad M = \langle \mathbf{x}_{uv} N \rangle, Nc=xvv,N N_c = \langle \mathbf{x}_{vv}, N \rangle

Şimdi silindir ve düzleme tekrar bakalım:

  • Düzlem: İkinci türevlerin hepsi 0'dır. Dolayısıyla L=M=Nc=0L=M=N_c=0. Düzlem bükülmez.
  • Silindir: xuu=(cosu,sinu,0)\mathbf{x}_{uu} = (-\cos u, -\sin u, 0) vektörü sıfır değildir ve normal vektörle aynı doğrultudadır. Bu yüzden L0L \neq 0 çıkar.

İşte fark burada! İkinci Temel Form, yüzeyin uzay ( E3\mathbb{E}^3 ) içindeki şeklini ve bükülmesini algılar. Gauss eğriliğine değinmeye ufak ufak başlayacağız, ancak öncesinde çok meşhur bir şekil olan hiperbolik paraboloid ya da nam-ı diğer Pringles cipsi'ne görsel olarak bir bakalım. Çünkü bu şekil normal vektörlerin oryantasyonunu gösterecek ve dolayısıyla eğrilik kavramına güzel bir giriş verecek.

Bir Örnek: Hiperbolik Paraboloid (Semer Yüzeyi)

Bükülmenin en güzel örneklerinden biri "semer" veya "Pringles cipsi" şeklinde olan Hiperbolik Paraboloiddir. Denklemi z=x2y2z = x^2 - y^2 (veya parametrik olarak x(u,v)=(u,v,u2v2)\mathbf{x}(u,v) = (u, v, u^2-v^2)) şeklindedir.

Bu yüzey üzerinde hareket ettiğinizi düşünün. Bir yönde yukarı doğru kıvrılırken, diğer yönde aşağı doğru kıvrılırsınız. İşte bu zıtlık, normal vektörlerin davranışında kendini gösterir.

Aşağıdaki Python kodunda bu yüzeyi ve üzerindeki normal vektörleri çizdiriyoruz. Normallerin nasıl farklı yönlere dağıldığına dikkat edin; bu bize yüzeyin "negatif eğriliğe" sahip olduğunu haykırır.

import numpy as np  
import matplotlib.pyplot as plt  
from mpl_toolkits.mplot3d import Axes3D  
  
fig = plt.figure(figsize=(10, 8))  
ax = fig.add_subplot(111, projection='3d')  
  
# Semer Yuzeyi: z = x^2 - y^2  
u = np.linspace(-1, 1, 15)  
v = np.linspace(-1, 1, 15)  
U, V = np.meshgrid(u, v)  
X, Y = U, V  
Z = U**2 - V**2  
  
ax.plot_surface(X, Y, Z, cmap='viridis', alpha=0.7)  
  
# Normallerin Hesabi ve Cizimi  
# xu = (1,0,2u), xv = (0,1,-2v) -> Normal ~ (-2u, 2v, 1)  
u_sample = np.linspace(-0.8, 0.8, 6)  
v_sample = np.linspace(-0.8, 0.8, 6)  
Us, Vs = np.meshgrid(u_sample, v_sample)  
Zs = Us**2 - Vs**2  
  
Nx = -2 * Us  
Ny = 2 * Vs  
Nz = np.ones_like(Us)  
Norm = np.sqrt(Nx**2 + Ny**2 + Nz**2) # Birimlendirme  
  
ax.quiver(Us, Vs, Zs, Nx/Norm, Ny/Norm, Nz/Norm, length=0.3, color='red')  
ax.set_title("Hiperbolik Paraboloid ve Normal Vektorleri")  
plt.show()

Hiperbolik Paraboloid ve Normal Vektörleri.

Büyük Finalin Fragmanı: Eğrilik

Elimizde iki tane form var:

  • II: İçsel ölçüm (Metrik)
  • IIII: Dışsal bükülme (Şekil)

Bu ikisinin oranı bize diferansiyel geometrinin en kutsal kavramını, Gauss Eğriliği'ni (KK) verir.

K=det(II)det(I)=LNcM2EGF2K = \frac{\det(II)}{\det(I)} = \frac{LN_c - M^2}{EG - F^2}

Bu formülde büyüleyici olan şudur: L,M,NcL, M, N_c dışarıdan bakışı gerektirse de, Gauss'un Theorema Egregium'u (Muhteşem Teorem) sayesinde KK'nın sadece E,F,GE, F, G cinsinden de hesaplanabileceği kanıtlanmıştır. Yani bir yüzeyin eğriliği, onun kaderidir; bükmekle değişmez!

  • K>0K > 0: Küre (Kapalı, bombeli)
  • K<0K < 0: Semer / Hiperbolik (Eyer şeklinde)
  • K=0K = 0: Düzlem veya Silindir (Düzlemsel geometri geçerli)

Python ile Metrik ve Bükülme Dedektifliği

Aşağıdaki kod, verdiğiniz herhangi bir yüzeyin hem içsel (E,F,GE,F,G) hem dışsal (L,M,NcL,M,N_c) analizini yapar.

import sympy as sp

u, v = sp.symbols('u v', real=True)
R = sp.symbols('R', real=True)

# Ornek: Kure Yuzeyi
surf = [R * sp.sin(u) * sp.cos(v), 
        R * sp.sin(u) * sp.sin(v), 
        R * sp.cos(u)]

def analyze_surface(pos, u, v):
    # 1. Kisim: Tegetler ve Birinci Form
    xu = [sp.diff(i, u) for i in pos]
    xv = [sp.diff(i, v) for i in pos]
    
    E = sp.simplify(sum(i*i for i in xu))
    F = sp.simplify(sum(i*j for i,j in zip(xu, xv)))
    G = sp.simplify(sum(j*j for j in xv))
    
    # 2. Kisim: Normal ve Ikinci Form
    # Normal vektor (Cross Product)
    n_raw = [xu[1]*xv[2]-xu[2]*xv[1], xu[2]*xv[0]-xu[0]*xv[2], xu[0]*xv[1]-xu[1]*xv[0]]
    n_mag = sp.sqrt(sum(k*k for k in n_raw))
    N = [k/n_mag for k in n_raw] # Birim Normal
    
    xuu = [sp.diff(i, u) for i in xu]
    xuv = [sp.diff(i, v) for i in xu]
    xvv = [sp.diff(i, v) for i in xv]
    
    L = sp.simplify(sum(i*j for i,j in zip(xuu, N)))
    M = sp.simplify(sum(i*j for i,j in zip(xuv, N)))
    Nc = sp.simplify(sum(i*j for i,j in zip(xvv, N)))
    
    return (E, F, G), (L, M, Nc)

# Hesapla
(E, F, G), (L, M, Nc) = analyze_surface(surf, u, v)

print(f"I. Form: E={E}, F={F}, G={G}")
print(f"II. Form: L={L}, M={M}, N={Nc}")

Gelecek Bölümde...

Artık elimizde hem cetvelimiz hem de bükülme ölçerimiz var. Gelecek bölümde bu eğrilik kavramının (KK) derinliklerine inecek, neden haritaların dünyayı hep yanlış gösterdiğini ispatlayacak ve Gauss'un dehasına şapka çıkaracağız.

Matematikle ve eğriyle kalın!

G

Gauss

Yazar