Diferansiyel denklem fikri

Açılış sahnesi

Bir fizyoloji laboratuvarında öğrenci, buzdolabından çıkan 10 derecelik bir tüpü 50 derecelik su banyosuna koyuyor ve her dakika termometreyi okuyor. Birinci dakikanın sonunda tüp 18 derece, ikincide 24,4, üçüncüde 29,5. Asistan soruyor: "On dakika sonra kaç derece olacak?" Öğrenci önce ilk dakikadaki artışı, sekiz dereceyi, on ile çarpmayı düşünüyor ve 90 dereceyi buluyor; sonra banyonun kendisinin 50 derece olduğunu hatırlayıp duruyor.

Öğrencinin elinde tüpün on dakika sonraki sıcaklığını veren bir formül yok. Ama tüpün şu anda hangi kurala göre ısındığını biliyor: tüp ile banyo arasındaki fark ne kadar büyükse, tüp o kadar hızlı ısınıyor. Bu kural yeter; seyri dakika dakika kurmak mümkün. Diferansiyel denklem fikrinin tamamı bu: durumu tarif etmek yerine değişimi tarif etmek. Bu ders o bakış değişikliğini elle yürütüyor ve ikinci yarısında o bakışın nerede kırıldığını gösteriyor.

Bu derste ne olacak

Önce fikri kuracağız ve birikim ve alan dersindeki denge denkleminin aslında bir diferansiyel denklem olduğunu göreceğiz. Sonra tüpün kuralını küçük adımlarla elle yürütecek, aynı kuralı iki farklı başlangıçtan çalıştıracağız. Ardından üç temel davranışı tanıyacak ve ilk ikisini iki kısa seride yan yana koyacağız. Son bölümde modelin sınırlarına geleceğiz: varsayımlardaki küçük farkların, verinin bittiği yerden sonra büyük farklar üretmesine ve bir modelin ne zaman sınanmış sayılacağına.

Sıfır noktası

Fikrin gücü, değişim kuralının çoğu zaman durumdan daha kolay bilinmesinde. Bir havuzun on gün sonra ne kadar dolu olacağını bilmiyor olabilirsin, ama giren ve çıkan akışların kurallarını biliyorsan seyri kurabilirsin. Birikim ve alan dersindeki denge denklemi tam olarak buydu: havuzdaki değişim hızı eşittir giren eksi çıkan. Bu modülde diferansiyel denklemleri zaten kullandın; burada yalnız adı konuyor.

Tüpün kuralını da aynı biçimde yazabilirsin: değişim hızı = k × (hedef − şimdiki). Burada hedef banyonun 50 derecesi, k ise dakikada 0,2. Bu denklem tüpün hangi sıcaklıkta olduğunu söylemiyor; yalnız her an nasıl değiştiğini söylüyor. Seyri kurmak için bir bilgi daha gerekiyor: nereden başladığı.

Mekanizma

Kuralı adım adım yürütmek

Sayısal çözümün mantığı basit: şimdiki değerden değişim hızını hesapla, küçük bir adım ilerle, yeni değerle tekrarla. Birer dakikalık adımlarla tüpü yürüt. Başta fark 50 − 10 = 40, bir dakikalık değişim 0,2 × 40 = 8, yeni sıcaklık 18. İkinci dakikada fark 32, değişim 0,2 × 32 = 6,4, sıcaklık 24,4. Üçüncüde fark 25,6, değişim 5,12, sıcaklık 29,52. Dördüncüde fark 20,48, değişim yaklaşık 4,1, sıcaklık yaklaşık 33,6.

Artışlar 8, 6,4, 5,12 diye küçülüyor, çünkü her adımda fark bir öncekinin 0,8'ine iniyor. Tüp hedefe yaklaştıkça yavaşlıyor ve 50'yi hiçbir adımda geçmiyor. Öğrencinin 90 derecelik tahmini, ilk dakikanın hızını sabit sanmaktan geliyordu; kural ise hızın her dakika yeniden hesaplanmasını istiyor.

Şimdi aynı kuralı başka bir başlangıçtan çalıştır. Laboratuvarın öbür masasında 80 derecelik bir tüp aynı banyoya konuyor. Fark 50 − 80 = −30, değişim 0,2 × (−30) = −6, sıcaklık 74. Sonra fark −24, değişim −4,8, sıcaklık 69,2; ardından 65,36. Kural aynı, banyo aynı, k aynı; biri yukarıdan öteki aşağıdan aynı 50'ye yaklaşıyor. Kural ile başlangıç ayrı bilgilerdir: iki kişi aynı kurala göre değişiyor olabilir ve yine de çok farklı yerlerde durabilir.

Adım büyüklüğü de bir karar. Bu kuralın tam çözümü ikinci dakikada tüpü yaklaşık 23,2 derecede buluyor; birer dakikalık adımlar 24,4 verdi. Adımı yarım dakikaya indirirsen her adımda fark 0,9'una iner ve dört adımda 50 − 40 × 0,9⁴ ≈ 23,8 bulursun. Hata 1,2 dereceden 0,6 dereceye düştü. Adım küçüldükçe çözüm doğrulaşıyor ama hesap artıyor; sayısal bir çözümün doğruluğu bir ayar meselesidir ve bildirilmesi gerekir. İntegral sezgisinde dilimleri incelterek bir alanı toplamıştın; burada aynı inceltme ters yönde, bir hızdan bir seyir kurmak için işliyor.

Üç temel davranış

Karmaşık sistemlerin çoğu birkaç temel davranışın birleşimidir ve üçü tanınırsa çoğu seyir okunur. Birincisi kendi büyüklüğüyle orantılı değişim: değişim hızı mevcut miktarla orantılıysa ortaya üstel davranış çıkar, büyüme ya da azalma. Üstel fonksiyonu işleyen L02-P01-M02-D03 bunu kurmuştu ve ayırt edici işareti sabit bir iki katına çıkma ya da yarılanma süresidir.

İkincisi hedefe yaklaşma: değişim hızı mevcut değer ile bir hedef arasındaki farkla orantılıysa sistem hedefe yaklaşır ve giderek yavaşlar. Tüp tam olarak bunu yaptı. Ayırt edici işareti başta hızlı, sonra yavaşlayan ve bir düzeye oturan seyirdir; birikim ve alan dersindeki kararlı durum da bu davranışın sonucudur.

İki davranışı yan yana koymak farkı hemen gösteriyor. Aynı 10'dan ve aynı 0,2 katsayısıyla başla. Hedefe yaklaşmada seri 10, 18, 24,4, 29,52; adımlar 8, 6,4, 5,12 diye küçülüyor. Kendi büyüklüğüyle orantılı değişimde kural değişim = 0,2 × şimdiki olur ve seri 10, 12, 14,4, 17,28; adımlar 2, 2,4, 2,88 diye büyüyor. Birincisi bir tavana yavaşlayarak yaklaşıyor, ikincisinin tavanı yok ve hızlanıyor. Denklemde tek fark, parantezin içinde "hedef − şimdiki" mi yoksa yalnız "şimdiki" mi yazdığı.

Üçüncüsü geri bildirimli değişim: sistemin çıktısı kendi girdisini etkiler. Negatif geri bildirim sapmayı geri çeker ve sistem kararlı kalır; pozitif geri bildirim sapmayı büyütür ve sistem uçlara gider. Tüpün kuralı aslında en basit negatif geri bildirimdir, çünkü sıcaklık hedefi geçerse fark eksiye döner ve tüpü geri çeker. Bu üçüncü davranışta bir ek var ve dersin en önemli teknik noktası o.

Tüpü bu kez bir banyo yerine, termometreye göre açılıp kapanan küçük bir ısıtıcıyla ısıttığını düşün; termometre ekrana her değeri iki dakika geç yansıtsın. Ekran hâlâ düşük gösterirken tüp çoktan hedefe varmıştır ve ısıtıcı açık kalır; tüp hedefi aşar. Geri dönüşte aynı gecikme onu bu kez aşağıda aşırır. Sonuç sabit bir değer yerine hedef etrafında bir dalgalanmadır. Glikoz ve insülin dinamiklerini işleyen ders bu dalgalanmayı sayılarla kuruyor ve biyolojik denetim sistemlerinin neden dalgalandığını oradan göreceksin.

Modelin sınırı

Bir model bir varsayımlar kümesinin sonucudur ve varsayımlar küçük ölçüde değiştiğinde sonuç büyük ölçüde değişebilir. Bu bölüm dersin en önemlisi, çünkü modeller ikna edici görünür.

Öğrencinin ilk üç okumasına, başlangıçtaki 10'a ve ilk iki dakikanın 18 ile 24,4'üne, iki model uydur. Birincisi bildiğin model: hedef 50, k = 0,2. İkincisinde hedef 70 ve k yaklaşık 0,13 olsun; ilk adımda o da 0,13 × 60 ≈ 8 derece ekler ve 18'i verir, ikinci dakikada 24,9 der. Termometrenin okuma hatası yarım derece kadarsa iki model bu üç okumaya eşit iyi uyuyor. Onuncu dakikada ise birincisi yaklaşık 45,7, ikincisi yaklaşık 55,7 diyor. On derecelik fark, modellerin kurulduğu verinin bittiği ikinci dakikadan sonra açılıyor.

Kartın son cümlesi bir okuma kuralı veriyor: tek bir sayı olarak sunulan bir model öngörüsü eksik sunulmuştur. Öğrencinin dürüst cevabı "onuncu dakikada 45 ile 56 derece arasında, hangi modelin doğru olduğuna bağlı" olurdu. Bu belirsizliği gidermenin yolu da modeli biraz daha uzun gözlemek.

İkinci sınır, bir modelin makul görünmesinin doğru olduğunu göstermemesi. Model bir mekanizmayı tarif eder ve mekanizmacı akıl yürütmenin sınırlarını işleyen L01-P03-M01-D09'daki kural burada da geçerli: bir mekanizmanın tutarlı olması, sonucun doğru olduğunu göstermez. Bir model ancak öngörüsü sınandığında değer kazanır.

Tüp örneğinde sınama elinin altında: iki model de ilk üç okumadan kuruldu, üçüncü dakikanın okuması ise ikisinin de kullanmadığı yeni bir veri. Birincisi 29,5, ikincisi yaklaşık 30,9 diyor ve fark okuma hatasının üç katına yakın. Termometre 29,5 okumuştu; birinci model sınavı geçiyor, ikincisi düşüyor. Bu, rakip hipotezleri işleyen L01-P06-M03-D03'teki ayırt edicilik kuralının model hâli: uyum yeterli değildir, iki modelin farklı şey öngördüğü bir gözlem gerekir.

Klinik gerçeklik

Bir model öngörüsüyle karşılaştığında önce değişim kuralını sor: bu seyir kendi büyüklüğüyle orantılı değişim mi, hedefe yaklaşma mı, geri bildirimli değişim mi, yoksa bunların bir birleşimi mi? Sonra başlangıç koşulunu sor, çünkü aynı kural farklı başlangıçtan farklı seyir üretir ve iki kişinin aynı kurala göre değişmesi aynı yerde olduklarını göstermez. Gecikme varsa, dalgalanan bir ölçümün kendiliğinden bir bozukluk işareti olmayabileceğini aklında tut.

Ardından öngörünün bir aralık olarak mı yoksa tek bir sayı olarak mı verildiğine bak; tek sayıysa eksiktir. Son olarak modelin, kurulurken kullanılmayan bağımsız bir veride sınanıp sınanmadığını sor. Bu dersin ölçütü bu son sorudur: uydurulmuş bir model ile sınanmış bir model aynı şey değildir, ve aradaki fark çoğu zaman yalnız yöntem bölümünün bir cümlesinde yazılıdır.

Kendini sına

1. Diferansiyel denklem fikrini bir cümlede yaz. Çözümü neden bir sayı değil de bir seyir?
2. Birikim ve alan dersindeki denge denklemi neden bir diferansiyel denklem?
3. Aynı kuralı (hedef 50, k = 0,2) 30 derecelik bir tüpten başlatarak üç adım yürüt. Adımlar nasıl değişiyor?
4. Üç temel davranışı ve ayırt edici işaretlerini yaz. Birinci ile ikinciyi denklemde ne ayırıyor?
5. Gecikme geri bildirimli bir sistemde ne üretiyor? Termometre örneğiyle açıkla.
6. Adımı bir dakikadan yarım dakikaya indirmek sonucu nasıl değiştirdi? Bu hangi uyarıyı getiriyor?
7. Modelin sınanması nedir? Tüp örneğinde hangi ölçüm iki modeli ayırıyor?
8. Vakadaki ayrışma nerede başlıyor ve neden tek sayıyla bildirilince görünmüyor?

Bir sonraki derse köprü

Bu derste değişim kurallarını tek tek kurdun ve iki davranışı ayrı serilerde gördün. Sıradaki ders ikisini tek bir mekanizmada birleştiriyor: vücuda giren bir maddenin miktarının zamanla nasıl değiştiği. Tek bir girişten sonra birinci davranış, yani kendi büyüklüğüyle orantılı azalma ve sabit bir yarılanma süresi çıkıyor. Giriş tekrarlandığında ikinci davranış, yani bir düzeye yaklaşma devreye giriyor. Orada, o düzeye ulaşmanın neden verilen miktara değil yalnız yarılanma süresine bağlı olduğunu ve bir mekanizma modelinin neyi göstermediğini göreceksin.