Özdeğer ve özvektör sezgisi

Açılış sahnesi

Yıl sonunda kliniğin veri sorumlusu, aynı aylık geçiş oranlarıyla izlenen iki takip listesini yan yana koyuyor. İkisinde de 150 kişi var. Birinci liste yıla 140 belirtisiz ve 10 belirtili kişiyle başlamış, ikincisi tam tersine yakın, 20 belirtisiz ve 130 belirtiliyle. Aralık ayında birinci listede yaklaşık 51 belirtisiz, ikincide yaklaşık 50 belirtisiz var. İki liste on iki ay önce birbirinden bu kadar uzakken, şimdi neredeyse aynı yerde duruyor.

Veri sorumlusunun sorusu basit: başlangıç bu kadar farklıyken sonuç neden bu kadar benzer? Cevap, geçiş matrisinin içinde saklı iki sayıda ve iki doğrultuda. Biri listenin varacağı yeri, öteki oraya ne hızla varacağını söylüyor. Bu ders o sayılara ve doğrultulara ad veriyor ve onları elle bulmayı gösteriyor.

Bu derste ne olacak

Önce bir dönüşümün yönünü değiştirmediği özel doğrultuları tanımlayacak ve takip listesinin geçiş matrisinde bu doğrultuları elle çarparak bulacağız. Sonra her başlangıç durumunu bu iki doğrultunun toplamı olarak yazacak ve açılıştaki iki listenin neden aynı yere vardığını hesaplayacağız. Ardından özdeğerin büyüklüğünden kararlılığı okumaya ve salınımın nereden çıktığına geçeceğiz. Son bölümde bu güçlü sezginin üç yerde nasıl kırıldığını göreceksin.

Sıfır noktası

Matris çarpımı sezgisini işleyen L02-P01-M04-D02'de 150 kişilik listenin durumunu (belirtisiz, belirtili) diye iki sayılık bir vektörle yazmış ve aylık geçişi şu matrisle tarif etmiştik: birinci satır (0,8, 0,1), ikinci satır (0,2, 0,9). Bu matris çoğu vektörü hem kısaltır ya da uzatır hem de döndürür, yani iki sayının birbirine oranını değiştirir. (100, 50) durumu bir ay sonra (85, 65) olur; oran 2'ye 1'den yaklaşık 1,3'e 1'e iner.

Tanımı sınamanın yolu çarpmaktır. (1, 2) vektörünü, yani bir belirtisize karşı iki belirtiliyi al. Birinci satır 0,8 × 1 + 0,1 × 2 = 0,8 + 0,2 = 1, ikinci satır 0,2 × 1 + 0,9 × 2 = 0,2 + 1,8 = 2 verir. Sonuç yine (1, 2): vektör ne döndü ne değişti, özdeğeri 1. Liste bu oranda, örneğin 50 belirtisiz ve 100 belirtiliyle başlasaydı, hiçbir ay yerinden kıpırdamazdı. Dengenin adı budur.

İkinci özel doğrultu (1, −1) vektörüdür. Bir kişinin belirtili kutusundan belirtisiz kutusuna geçtiği bir değişikliği temsil eder, toplamı sıfırdır. Çarp: 0,8 × 1 + 0,1 × (−1) = 0,7 ve 0,2 × 1 + 0,9 × (−1) = −0,7. Sonuç (0,7, −0,7) = 0,7 × (1, −1). Yön aynı kaldı, uzunluk 0,7 katına indi; özdeğeri 0,7. Özvektörler sistemin doğal yönleridir, özdeğerler o yönlerdeki büyüme ya da sönme oranlarıdır.

Mekanizma

Her durum iki doğrultunun toplamı

Bu iki doğrultunun gücü, her başlangıç durumunun onların toplamı olarak yazılabilmesinden gelir. Matris çarpımı dersindeki başlangıç (100, 50) şöyle ayrılır: (50, 100) + 50 × (1, −1) = (50 + 50, 100 − 50) = (100, 50). Birinci parça denge doğrultusunda, ikinci parça sapma doğrultusunda. Toplamı 150 olan her liste için denge parçası hep (50, 100) çıkar, çünkü sapma vektörünün toplamı sıfırdır ve kişi sayısını değiştirmez.

Şimdi dönüşümü parça parça uygula. Denge parçasının özdeğeri 1 olduğu için yerinde kalır. Sapma parçası her ay 0,7 ile çarpılır: 50, 35, 24,5. Bir ay sonra liste (50, 100) + 35 × (1, −1) = (85, 65), iki ay sonra (74,5, 75,5). Bunlar matris çarpımı dersinde satır satır bulduğun sayıların aynısı, ama bu kez tek bir çarpmayla geliyor. Her adımda iki satırı yeniden hesaplamak yerine, sapmanın 0,7'nin kuvvetleriyle küçüldüğünü bilmek yetiyor.

Açılıştaki iki liste de böyle çözülür. Birinci liste (140, 10) = (50, 100) + 90 × (1, −1), ikinci liste (20, 130) = (50, 100) − 30 × (1, −1). On iki ay sonra sapmalar 90 × 0,7¹² ≈ 1,2 ve −30 × 0,7¹² ≈ −0,4 olur. Listeler yaklaşık 51,2 ve 49,6 belirtisize varır. Başlangıç farkının tamamı sapma doğrultusundaydı ve o doğrultu her ay sönüyordu.

Takip listesinde baskın yön (1, 2) doğrultusudur, çünkü özdeğeri 1, ötekinin 0,7'den büyük. Diferansiyel denklem fikrini işleyen L02-P01-M03-D06'da aynı kuralın farklı başlangıç koşullarından farklı seyirler ürettiğini görmüştün. Burada o gözleme bir sınır ekleniyor. Başlangıç koşulu kısa vadeyi belirliyor, yani ilk aylarda listelerin nerede olduğunu; baskın yön uzun vadeyi, yani hepsinin nereye vardığını belirliyor.

Özdeğer kararlılığı söylüyor

Sapmanın özdeğeri 1'den küçük olduğu için sapma söner ve sistem bozulunca kendine gelir. Özdeğer 1'den büyük olsaydı ne olurdu, bir satırla gör: 1,1'lik bir özdeğerle 50 kişilik sapma 50 × 1,1 = 55, sonra 60,5, sonra 66,55 olur ve on iki ay sonra yaklaşık 157'ye çıkar. Yani 150 kişilik listede taşacak kadar büyür. Böyle bir sistemde en küçük bozulma bile büyür; sistem kararsızdır. Özdeğer tam 1'e yakınsa sapma ne büyür ne belirgin biçimde söner ve sistem küçük değişikliklere duyarlı, sınırda bir yerde durur.

Glikoz ve insülin dinamiklerini işleyen dersteki negatif geri bildirim bu tablonun fizyolojik karşılığıdır. Döngü bir sapmayı her adımda belirli bir oranda geri çekiyorsa, o oran sapma doğrultusunun 1'den küçük özdeğeridir. Negatif geri bildirimin sapmayı söndürmesi, matris dilinde sapma özdeğerinin 1'den küçük olması demektir.

Sönme hızını tek bir sayıyla da okuyabilirsin: sapmanın yarıya inmesi için kaç adım gerektiği. 0,7'lik özdeğerde 0,7 × 0,7 = 0,49, yani yaklaşık iki ayda sapma yarıya iner. 0,95'lik bir özdeğerde bu yaklaşık on üç buçuk adım sürer. İkinci sistem de kararlıdır ama bozulduğunda çok uzun süre bozuk görünür.

Kartın notundaki uyarı bilinçli. Takip listesinde baskın yön olan 1'e 2 oranı, belirtisizliğin ve belirtililiğin iki ayrı mekanizmayla sürdüğünü göstermez; yalnız iki geçiş oranının birbirini dengelediği noktayı gösterir. İstatistiksel bir baskın yön bulmak, bir mekanizma bulmak değildir.

Salınım nereden çıkıyor

Bazı sistemlerde uzun vadeli davranış bir doğrultuya oturmak değil, bir çevrimdir. Bunun matristeki izi, yönü değişmeyen hiçbir gerçek doğrultunun bulunmamasıdır. En yalın örnek birinci satırı (0, −1), ikinci satırı (1, 0) olan matristir. (1, 0) vektörü bir adımda (0, 1)'e, sonra (−1, 0)'a, sonra (0, −1)'e döner ve dördüncü adımda başladığı yere gelir. Her vektör çeyrek tur döndüğü için hiçbiri yönünü korumaz.

Bu durumda özdeğerler karmaşık sayılar olarak çift hâlinde çıkar ve sistem bir denge noktasına gitmek yerine onun çevresinde dolaşır. Dönmeyle birlikte her adımda bir küçülme de varsa çevrim içe doğru bir sarmala dönüşür: sistem salınarak dengeye yaklaşır. Aynı diferansiyel denklem dersinde gecikmenin geri bildirimli bir sistemde salınım ürettiğini görmüştün. Matris tarafında gecikme, sistemin durumuna "iki adım önceki değer" gibi ek boyutlar ekler ve dönmeyi üreten de o boyutlardır.

Sezginin sınırı

Bu araç güçlüdür ama üç yerde kırılır. Birincisi doğrusallık varsayar. Özdeğer çözümlemesi dönüşümün doğrusal olduğunu kabul eder, oysa biyolojik süreçler çoğu zaman eşikli ve doygunluk gösteren, parçalı fonksiyonlarla anlatılan ilişkilerdir. Araç bu yüzden yalnız dengenin yakınındaki dar bir bölgede geçerlidir; doğrusal fonksiyon dersinde gördüğün yerel doğrusallığın sistem düzeyindeki hâli.

İkincisi, "uzun vade" ne kadar uzun olduğunu söylemez; bunu özdeğerlerin birbirinden ne kadar ayrıldığı belirler. Takip listesinde 1 ile 0,7 arasındaki fark geniş, liste bir yılda dengeye oturuyor. Sapma özdeğeri 0,98 olsaydı on iki ay sonra sapmanın yaklaşık dörtte üçü hâlâ yerinde dururdu ve denge gözlem penceresi içinde hiç görünmeyebilirdi. Üçüncüsü, dönüşümün kendisi zamanla değişiyorsa, ki biyolojide sık görülür, sabit bir baskın yön yoktur.

Klinik gerçeklik

Bir sistemin davranışıyla karşılaştığında önce uzun vadede oturduğu örüntüyü sor: bir dengeye mi yaklaşıyor, büyüyor mu, salınıyor mu? Sonra bir sapmanın sönüp sönmediğine bak; sönüyorsa sistem kararlıdır, büyüyorsa değildir. Ardından sönmenin ne kadar hızlı olduğunu sor. Yarılanma süresi uzunsa sistem kararlı ama tepkisizdir ve dinlenme ölçümleri bunu göstermez.

Son olarak çözümlemenin hangi aralıkta geçerli olduğunu sor, çünkü özdeğer hesabı dengenin yakınında yapılmış bir doğrusal yaklaşımdır. Büyük bir bozulmada sistemin aynı biçimde davranacağını söylemez. Bu dört sorunun en az sorulanı ve en çok bilgi vereni üçüncüsüdür: toparlanma hızı, dinlenme düzeyinden farklı bir şey söyler.

Kendini sına

1. Özvektör ve özdeğeri tanımla. (1, 2) ve (1, −1) vektörlerini geçiş matrisiyle çarparak özdeğerlerini bul.
2. (30, 120) başlangıcını denge ve sapma parçalarına ayır. Bir ay sonra kaç belirtisiz olur? Sonucu doğrudan matris çarpımıyla kontrol et.
3. Sapma vektörünün toplamı neden sıfırdır? Bu, denge parçası için ne anlama gelir?
4. Başlangıç koşulu neyi belirliyor, baskın yön neyi? Açılıştaki iki liste üzerinden açıkla.
5. Sapma özdeğeri 1,1 olsaydı sistem nasıl davranırdı? 0,98 olsaydı ne değişirdi?
6. Sönme oranı nedir? 0,7 ile 0,95 özdeğerlerini yarılanma süresiyle karşılaştır.
7. Dönme matrisinde neden yönü değişmeyen gerçek bir doğrultu yok? Bu hangi davranışı üretir?
8. Vakadaki iki sistem nasıl ayırt ediliyor? Dinlenme ölçümü neden yetmiyor?

Bir sonraki derse köprü

Bu derste sistemi bir matris olarak yazdın ve her sayının neyi neye bağladığını bildiğin için hesap yapabildin. Ama matrisi kuran biri önce hangi niceliklerin ayrı birer bileşen sayılacağına ve hangi ikisinin arasında bir bağlantı bulunduğuna karar verdi. Sıradaki derste bu kararların kendisine geçiyoruz: sistemi noktalar ve çizgilerle, yani düğümler ve kenarlarla çiziyoruz. O çizimin kenarını hangi eşikte, hangi anlamla çizdiğin, sonradan hesaplanacak her ölçüyü belirliyor.