21 Mart 2024 Perşembe

Konvekslik ve norm arasındaki bağ - Kreyszig'in örneği

Kreyszig'in fonksiyonel analiz kitabında Banach uzayları anlatılırken güzel bir önerme iki alıştırmaya paylaştırılarak veriliyor. Önerme, norm olmaya aday bir fonksiyonun gerçekten de norm olup olmadığının nasıl anlaşılacağına dair. Örneğin $x \in \mathbb{R}^{2}$ olmak üzere \begin{equation*} \varphi(x) := \left( \sqrt{|x_{1}|} + \sqrt{|x_{2}|} \right)^{2} \end{equation*} tanımlansın. Bu fonksiyonun $\mathbb{R}^{2}$ üzerinde bir norm tanımlayıp tanımlamadığını nasıl anlarız? Genel olarak bu tip sorgularda norm olmanın aksiyomlarından üçgen eşitsizliği hariç diğerleri kolayca ispatlanırken üçgen eşitsizliği ispatlanması en zor kademeyi teşkil eder. Kreyszig burada konveks analizden alternatif bir rota öneriyor. Ama öncelikle konvekslikten ne kastettiğimizi tanımlayalım.

Tanım: $X$ bir lineer uzay olsun. $A \subset X$ olacak şekilde \begin{equation*} x,y \in A \ \Rightarrow \ M := \left\{ z \in X \ | \ z = \alpha x + (1-\alpha)y, \ \alpha \in [0,1] \right\} \subset A \end{equation*} önermesi doğru ise $A$ konveks bir kümedir.

Bu temel tanımda $M$, başı $x$ ucu $y$ olan bir doğru parçası. Sözel olarak ifade edildiğinde verilen bir kümeden herhangi iki nokta alıp bunları bir doğru parçası oluşturacak şekilde birleştiriyoruz. Doğru parçasının tamamı bahsi geçen kümede kalıyorsa o zaman o kümenin konveks olduğunu söylüyoruz. Örneğin çember, elips, kare, küp ve küre konveks kümeler. Hilal, yıldız ve amip konveks kümeler değil. İkinci bir tanımımız daha var. O da kapalı birim topun tanımı. Bu da geometrik bir kavramı genelleştiriyor.

Tanım: Normlu bir lineer uzayda $\| x \| \leq 1$ eşitsizliğini sağlayan tüm $x$ vektörlerinin oluşturduğu kümeye kapalı birim top denir ve $B_{1}(0)$ ile gösterilir.

Geleneksel geometride "top" denilen objeler konveks. Kreyszig bu noktada basit bir önermeyle yol yapıyor. O da şu: Kapalı birim top gerçekten de konveks bir küme. İspatlamak için $x,y\in B_{1}(0)$ ve $\alpha \in [0,1]$ alalım. $z = \alpha x + (1-\alpha)y$ olsun. O zaman \begin{eqnarray}\nonumber \|z\| &=& \| \alpha x + (1-\alpha)y \| \\ \nonumber &\leq& \| \alpha x \| + \| (1-\alpha)y \| \\ \nonumber &=& |\alpha| \| x \| + |1-\alpha| \|y\| \\ \nonumber &=& \alpha \|x\| + (1-\alpha) \|y\| \\ \nonumber &\leq& \alpha + 1 - \alpha \\ \nonumber &=& 1 \end{eqnarray} olur. Son basamakta $x$ ve $y$ vektörlerinin kapalı birim top içinde olduklarını kullandık.

Gelelim norm olmaya aday $\varphi$ fonksiyonuna. $\varphi(x) \leq 1$ ile tanımlanan kapalı birim top içinde $x=(1,0)$ ve $y=(0,1)$ noktalarını alalım. $\alpha = 1/2$ olsun. O zaman $z = \alpha x + (1-\alpha)y = (1/2,1/2)$ olur. Ama $\varphi(z) = 2 > 1$ olduğundan bu nokta kapalı birim top içinde değildir. Demek ki $\varphi$ tarafından üretilen kapalı top konveks bir küme değilmiş. O zaman $\varphi$, $\mathbb{R}^{2}$ üzerinde bir norm tanımlamaz.

15 Mart 2024 Cuma

Makara sistemlerinde yol kaybı dinamik argumanlara girmeden de hesaplanabilir

Soru: Şekilde görülen ideal makara sistemi serbest bırakıldığında $M$ yukarı doğru ivmelenirken $m$ aşağı yönde hareket etmektedir. $M$ kütlesi $\xi$ kadar yükseldiğinde $m$ ne kadar aşağı doğru inmiştir?

Bu soruyu liselere geçiş sınavına hazırlanan kızım dün akşam getirdi ve ben de geometrik argumanlarla bir çözüm yaptım. Onu bu postada sizinle paylaşacağım. Yaptığım çözümü kızım çok uzun buldu. Onun yerine $M$'nin ağırlığı için $60$ N dedikten sonra ancak statik dengede geçerli olan iplerdeki gerilimleri heaplayarak kuvvet kazancının $4$ kat olduğunu hesapladı. Sonra da yol kaybından çözüme ulaştı. Halbuki bu pür geometrik bir problem. Dinamikle, kuvvetle, ip gerilimleriyle bir alakası yok!

Sağdaki şekli kitaptaki soruyu taklit ederek ben çizdim. İp parçalarının üzerindeki uzunluklar tabii ki kitapta yoktu. (Yeri gelmişken ilgili TikZ komutlarını web'de paylaşarak bizi günlerce debelenip durmaktan kurtaran Izaak Neutelings'e teşekkürler.) Bu problemi geometriye ve lineer cebirin kurallarına uyarak çözeceğiz.

Şekilde aslında iki parça ip var. İpler hareket esnasında esnemediklerinden bu ip parçalarının uzunlukları korunuyor. Özdeş makaraların çevresine $C$ diyelim. O zaman iplerin uzunlukları aşağıdaki denklemlerle verilir. \begin{eqnarray}\nonumber l & := & l_{1}(t) + l_{2}(t) + l_{3}(t) + C \\ \nonumber k & := & k_{1}(t) + k_{2}(t) + \tfrac{1}{2}C \end{eqnarray} Burada ipin bir makara etrafında yarım devir tur attığını varsaydık. Bu iki denklemi görmesi kolaydı. Şimdi, ortadaki makarayı $k$ ipine bağlayan aparatın sabit kalan uzunluğuna $A$ dersek, bütün hareket boyunca geçerli üçüncü bir denkleme daha ulaşırız. \begin{equation*} k_{1}(t) + l_{3}(t) + A = k_{2}(t) \end{equation*} Bu denklem $M$ kütlesinin bağlı olduğu makaranın tavana göre konumunu iki farklı yolla ölçüyor. Son olarak birinci makarayı tavana bağlayan aparatın uzunluğuna da $A$ dersek ortadaki makaranın tavana göre konumunu sağındaki ve solundaki ip uzunluklarıyla ölçen bir denklem daha bulunur. \begin{equation*} l_{2}(t)+A=l_{3}(t) \end{equation*} İşte bu dört geometrik kısıtlamayı kullanarak yol kaybını bulmaya çalışacağız.

$k_{2}$ sadece bir ip segmentinin uzunluğu değil aynı zamanda $M$ kütlesinin tavana göre konumunu da ölçüyor. Bu kütlenin $\xi(t)$ kadar yükselmesi demek $k_{2}(t)=k_{2}(0)-\xi(t)$ demektir. Ama $k$ ipinin toplam uzunluğu korunduğundan $k_{1}(0)+k_{2}(0)+\tfrac{1}{2}C=k_{1}(t)+k_{2}(0)-\xi(t)+\tfrac{1}{2}C$ denklemi geçerlidir ve buradan $k_{1}(t) = k_{1}(0)+\xi(t)$ ara sonucu bulunur.

$k$ ipindeki segmentlerin uzunluklarını bulduk. Şimdi aynı işlemleri $l$ ipi için de yapacağız. Şimdi, $l_{3}(t)=k_{2}(t)-k_{1}(t)-A$ ve $A=k_{2}(0)-k_{1}(0)-l_{3}(0)$ denklemleri beraber kullanıldığında $l_{3}(t) = l_{3}(0)-2\xi(t)$ eşitliği bulunur. Çorabın nasıl söküldüğünü bu aşamada artık okur kavramıştır. Bu yüzden sonuçları vererek çözümü bitirelim. $l_{2}(t) = l_{2}(0)+2\xi(t)$ ve $\boxed{l_{1}(t) = l_{1}(0)+4\xi(t)}$. Demek ki $m$ kütlesi $4\xi$ kadar aşağı inmiş.

İşaret: Öğrencimizin problemi hızlıca çözmek için neden büyük ağırlık yerine $60$ N koyduğunu ilk başta ben de merak ettim. Sonra Sümerliler'in neden aritmetiği 60 tabanında yaptığı aklıma geldi. Çünkü 60 sayısının bölünebilirliği 10'a kıyasla daha fazla! İp gerilimlerini hesaplarken bu ağırlığı böle böle ilerliyorsunuz.

10 Mart 2024 Pazar

Bir makara probleminde kütlelerin olası tüm hareket yönleri

Aşağıdaki ideal makara düzeneğinde genel olarak birbirinde farklı kütlelere sahip yükler esnemeyen ve gergin bir iple bağlanmıştır. Bu postada olası tüm kütle kombinasyonları için parametre uzayında kütlelerin ivmelenme yönlerini sorgulayacağız.

PDF formatında çözümü indirmek için tıkla.

3 Ocak 2024 Çarşamba

İki quadratiğin ortak kökü bulunması için gerek ve yeter şart

Soru: $q(X):=aX^{2}+bX+c$ ve $Q(X):=AX^{2}+BX+C$ ile iki kuadratik polinom veriliyor. Bu ikisinin ortak kökleri bulunması için katsayılar üzerindeki gerek ve yeter şartı bulunuz.

İkaz: Burada polinomların gerçel analitik düzlemde kesiştiği noktaları aramıyoruz. Örneğin $x^{2}+1$ ve $-x^{2}-1$ polinomları gerçel analitik düzlemde asla kesişmez ama tüm kökleri ortaktır.

Çözüm: Buna benzer bir problemi geçen yılın son blog postasında çözmüştük. Orada kullandığımız temel fikir aradığımız özellik sağlanınca sıfır olan ve kökler cinsinden yazılabilen simetrik bir indikatör fonksiyonu icat etmekti. Her ne kadar indikatör fonksiyon kökler cinsinden yazılsa da permütasyon simetrisinden dolayı kökleri hesaplamadan da değeri bulunabiliyordu. Bu problemi çözerken de böyle yapacağız ve adına resultant denilen bir fonksiyon tanımlayacağız. $q$ ve $Q$ polinomlarının kökleri sırasıyla $\{r_{1},r_{2}\}$ ve $\{s_{1},s_{2}\}$ olsun. O zaman ilk defa J. J. Sylvester tarafından tanıtılan resultant fonksiyonu \[ R(r_{1},r_{2};s_{1},s_{2}) := A^{2}a^{2}(r_{1}-s_{1})(r_{1}-s_{2})(r_{2}-s_{1})(r_{2}-s_{2}) \] ile verilir.

İşaret: Polinomların derecesi 2 olduğundan $a \ne 0$ ve $A \ne 0$ olmalıdır. O zaman $R$ fonksiyonunun sıfır olması için gerek ve yeter şart $r_{i}=s_{j}$ biçiminde bahsi geçen polinomların en az bir ortak köke sahip olmalarıdır.

İşaret: Resultant fonksiyonu $\{r_{1},r_{2}\}$ arasındaki permütasyonlara göre simetriktir. Benzer şekilde $\{s_{1},s_{2}\}$ değişkenleri arasındaki permütasyonlara göre de simetriktir. Ancak tüm değişkenleri alıp bir permütasyon yaparsanız simetri bozulabilir.

Resultantın simetrik olması onu Viete-Girard polinomlarını kullanarak katsayılar cinsinden ifade edebileceğimiz anlamına geliyor. Cebirin temel teoremine göre $Q(X)=A(X-s_{1})(X-s_{2})$ olduğundan \[ R = a^{2} Q(r_{1})Q(r_{2})\] yazabiliriz.

İşaret: Bize resultantın değerinden ziyade işareti lazım ve işaret $Q$ polinomunun $q$ polinomunun köklerindeki işaretine bağlı. Bir polinomun işaretinin başka bir polinomun köklerindeki değerinin sorgusuna literatürde Tarski sorgusu denir ve cebirsel geometeri literatüründe bu problemin kendine özgü bir yeri vardır.

Resultant, Sylvester matrisleri ve determinantları kullanılarak etkin bir şekilde hesaplanabilir. Biz burada problemin küçüklüğünden dolayı sadece simetrik fonksiyonları kullanacağız. Hatırlatalım: $\sigma_{1} := r_{1}+r_{2} = -b/a$ ve $\sigma_{2} := r_{1}r_{2} = c/a$ ile tanımlanmaktaydı. \begin{eqnarray}\nonumber R &=& a^{2}(Ar_{1}^{2}+Br_{1}+C)(Ar_{2}^{2}+Br_{2}+C) \\ \nonumber &=& a^{2} \left(A^{2}r_{1}^{2}r_{2}^{2} + B^{2}r_{1}r_{2} + C^{2} + ABr_{1}r_{2}(r_{1}+r_{2}) + AC(r_{1}^{2}+r_{2}^{2}) + BC(r_{1}+r_{2})\right) \\ \nonumber &=& a^{2}\left(A^{2} \frac{c^{2}}{a^{2}} + B^{2}\frac{c}{a} + C^{2} - AB \frac{bc}{a^{2}} + AC((r_{1}+r_{2})^{2}-2r_{1}r_{2}) - BC \frac{b}{a} \right) \\ \nonumber &=& A^{2}c^{2} + B^{2}ac + C^{2}a^{2} - ABbc + ACb^{2} - 2ACac - BCab \\ \nonumber &=& (Ac-Ca)^{2} - (Ab-Ba)(Bc-Cb) \end{eqnarray}