Glikoz–insülin dinamikleri

Açılış sahnesi

Bir sabah iki komşu aynı laboratuvarda açlık kan şekeri ölçtürüyor ve öğleden sonra sonuçları yan yana koyuyorlar. İki sayı birebir aynı ve ikisi de raporun referans aralığının ortasında duruyor. Birinin sistemi hiç zorlanmadan o düzeyde duruyor olabilir; ötekinin sistemi aynı düzeyi tutmak için her gün çok daha fazla çalışıyor olabilir. Kâğıttaki sayı bu iki durumu birbirinden ayırmıyor.

Bu ders o ayırmama durumunun nereden geldiğini, geri bildirimli bir denetim döngüsünün matematiğiyle gösteriyor. Döngüyü küçük bir kurgu modelle adım adım hesaplayacak, gecikmenin nasıl salınım ürettiğini sayılarla görecek ve sonunda açılıştaki iki komşuya döneceğiz.

Bu derste ne olacak

Önce döngüyü kuracağız: iki nicelik ve aralarındaki iki yönlü bağ. Sonra bu döngünün bir sapmayı nasıl geri çektiğini kurgu bir modelde elle hesaplayacak, kazancın hızı nasıl değiştirdiğini göreceğiz. Ardından modele gecikme ekleyip hedefin neden aşıldığını adım adım izleyeceğiz. Son bölümde dersin asıl sorusuna geleceğiz: denetlenen bir değişkeni ölçmek neyi göstermiyor?

Sıfır noktası

Kandaki glikoz yükseldiğinde pankreastan insülin salınır. İnsülin glikozun dokulara alınmasını ve depolanmasını artırır, karaciğerdeki üretimini azaltır; glikoz düştükçe insülin salınımı da azalır. Düzey düştüğünde ise karşıt yönde çalışan sinyaller devreye girer ve glikozu yükseltir. Yani sistem iki yönden, hem yükselmeye hem düşmeye karşı denetleniyor.

Birikim ve alan dersindeki havuz dili burada da işe yarıyor. Kandaki glikoz bir havuz; giren akışlar besinden gelen alım ve karaciğerin üretimi, çıkan akışlar dokulara alım ve depolama. İnsülin bu musluklardan bazılarını açıp bazılarını kısan denetleyici.

Mekanizma

Döngü sapmayı nasıl geri çekiyor

Şimdi bu döngüyü bir kurgu modele indirelim. Modeldeki sayılar birimsizdir ve hiçbir gerçek fizyolojik değeri temsil etmez; yalnız davranışı görmek için seçildi. Hedef düzey 100, başlangıç düzeyi 140. Kural tek: her adımda düzey, hedefle arasındaki farkın 0,3'ü kadar düşüyor.

İlk adımda fark 40, düşüş 40 × 0,3 = 12 ve düzey 140 − 12 = 128. İkinci adımda fark 28, düşüş 28 × 0,3 = 8,4, düzey 119,6. Üçüncü adımda fark 19,6, düşüş 5,88, düzey 113,72. Farklar 40, 28, 19,6, 13,72 diye gidiyor: her adımda bir öncekinin 0,7'si kalıyor. Diferansiyel denklem fikri dersinde "hedefe yaklaşma" diye adlandırılan davranış bu: başta hızlı, sonra giderek yavaşlayan ve hedefe oturan bir seyir.

Kazancı 0,6'ya çıkaralım. İlk adımda düşüş 40 × 0,6 = 24 ve düzey 116; sonra fark 16, düşüş 9,6, düzey 106,4; sonra 102,56. Üç adımda kalan fark 2,56, oysa 0,3 kazançla üç adım sonra hâlâ 13,72 vardı. Buraya kadar yüksek kazanç yalnız bir kazanç gibi görünüyor. Bunun bir bedeli olduğunu görmek için modele bir şey daha eklemek gerekiyor: zaman.

Gecikme neden dalgalanma üretiyor

Gerçek bir döngüde tepki anında oluşmaz. Sinyalin salınması, dolaşıma karışması, dokuya ulaşması ve dokunun tepki vermesi zaman alır. Modelde bunu en basit hâliyle yazalım: her adımdaki düşüş, şimdiki farka değil iki adım önceki farka göre hesaplanıyor. Sistemin başlangıçtan önce bir süre 140'ta beklediğini varsayalım.

Kazanç yine 0,3. İlk üç adımda denetleyici hâlâ 140'ı "görüyor", bu yüzden her adımda 12 düşüyor: 140, 128, 116, 104. Dördüncü adımda iki adım önceki 128'e tepki veriyor, düşüş 28 × 0,3 = 8,4 ve düzey 95,6. Hedef aşıldı, ama denetleyici bunu henüz görmüyor. Beşinci adımda 116'ya göre 4,8 düşüş, düzey 90,8; altıncıda 104'e göre 1,2 düşüş, düzey 89,6. Yedinci adımda nihayet hedefin altındaki 95,6'ya tepki geliyor: fark −4,4, düzeltme 4,4 × 0,3 = 1,32 yukarı ve düzey 90,92. Seyir önce hedefin on birim kadar altına iniyor, sonra geri dönerken hedefi yukarı doğru yeniden aşıyor.

Tanımdaki ikinci cümleyi sayılarla sına. Aynı iki adımlık gecikmeyle kazancı 0,6'ya çıkarınca seyir 140, 116, 92, 68 ve 58,4'e kadar iniyor, sonra 63,2, 82,4, 107,36, 129,44 diye tırmanıp dokuzuncu adımda yeniden 140'a varıyor. Gecikmesiz modelde en hızlı düzelten ayar, gecikmeli modelde sistemi başladığı yere geri götürdü. Maksimum ve minimum dersindeki çok amaçlı en iyileme burada somutlaşıyor: hızlı düzeltme ile kararlılık birbiriyle çatışıyor ve ikisini birden en iyileyen tek bir ayar yok.

Bu hesabın okuma sonucu şu: bir dalgalanma, sistemin bozuk olduğunu değil denetlendiğini gösterebilir. Denetimsiz bir sistem dalgalanmaz, sürüklenir. Dalgalanmanın varlığı bir denetim işareti; genliği ve sönme hızı ise kazanç ile gecikme hakkında bilgi taşıyor.

Kartın son cümlesi bir aşırı yorumu kesiyor. Kurgu modelin ürettiği salınım bir gözleme benzediğinde, bu modelin doğru olduğunu göstermez; aynı biçimi başka mekanizmalar da üretebilir.

Merkezî soru: denetlenen değişkeni ölçmek

Şimdi açılıştaki iki komşuya dönebiliriz. Bir denetim döngüsünde iki tür nicelik var ve ikisi çok farklı davranıyor.

Kurgu modeli biraz genişletelim. Havuza her adımda 20 birim giriyor ve düzeyin 100'de kalması için aynı miktarın çıkması gerekiyor. Birinci sistemde denetleyici sinyalin her birimi 4 birim uzaklaştırıyor, ikincide yalnız 1 birim. Birincinin dengede kalmak için ihtiyaç duyduğu sinyal 20 ÷ 4 = 5, ikincininki 20 ÷ 1 = 20. İki sistemde de denetlenen düzey 100, ama ikincinin denetleyicisi dört kat fazla çalışıyor.

Yedeğin neden önemli olduğunu görmek için denetleyiciye bir üst sınır koyalım: kurgu modelde sinyal en fazla 25 olabiliyor. Giriş 20'den 30'a çıktığında birinci sistemin ihtiyacı 30 ÷ 4 = 7,5 ve rahatça karşılanıyor. İkincinin ihtiyacı 30 ÷ 1 = 30, ama sinyal 25'te tavan yapıyor. Her adımda 30 − 25 = 5 birim havuzda kalıyor ve düzey artık adım adım yükseliyor. L02-P01-M02-D06'daki doygunluk ve eşik terimleri burada birleşiyor: denetleyici doygunluğa ulaşana kadar denetlenen değişken hiç kıpırdamıyor, eşik aşıldığında ise birden hareket etmeye başlıyor. O derste gizli değişim diye adlandırılan şey bu: ölçüm sabitken yedek tükeniyor.

Ayrım nasıl yapılıyor

İki sistemi ayırmanın ilk yolu, denetleyiciyi de ölçmek. Kurgu modelde düzeyin yanına sinyal düzeyi yazıldığında 5 ile 20 arasındaki fark hemen görünüyor; tek değişkene bakan ölçüm bu farkı hiçbir zaman gösteremez. İkinci yol sistemi zorlamak: dinlenmede ayırt edilemeyen iki sistem, girişin 20'den 30'a çıktığı bir yük altında tamamen farklı davrandı. Kapasite dinlenmede değil yük altında görünür ve dinamik değerlendirmelerin var olma nedeni budur.

Üçüncü yol zaman içinde izlemek. Anlık ve ortalama değişim dersinde iki noktanın aradaki yolu göstermediğini görmüştün; tek bir nokta hiç yol göstermiyor. İkinci sistemin tavana yaklaştığı dönemde bile düzey 100'de sabit görünecek, ama tekrarlanan ölçümlerde sinyalin yavaş yavaş tırmandığı okunabilir.

Modelin sınırı

Bu model dört yerde bilerek basitleştiriyor. Gerçek sistemde iki nicelik yetmez: çok sayıda sinyal, çok sayıda doku ve iç içe geçmiş birçok döngü birlikte çalışır, iki değişkenli model bunların yalnız iskeleti. Ayar noktası da sabit değildir; koşullara, günün saatine ve duruma göre kayar, oysa model 100'ü hiç oynatmadı.

Üçüncüsü, kazanç ve gecikme kişiden kişiye ve aynı kişide zaman içinde değişir; modeldeki 0,3, 0,6 ve iki adım yalnız davranışı göstermek için seçilmiş sayılar. Dördüncüsü, model bir davranışı tarif ediyor, nedenini değil. Bir salınımın modelin ürettiğine benzemesi, onu bu mekanizmanın ürettiğini göstermez; ilk kanıt kartının son cümlesi tam bunu söylüyor.

Bu ders neyi göstermiyor

Bu ders hiçbir eşik, aralık ya da ölçüm yorumu vermiyor ve bir durumu tanımlamıyor. Modeldeki 100, 140 ya da 25 hiçbir laboratuvar değerine karşılık gelmiyor; sistem dinamiği ile klinik değerlendirme ayrı işlerdir. Açılıştaki iki komşunun hangisinin hangi sistem olduğunu bu dersten çıkarmak mümkün değil.

Mekanizmanın tutarlı olması bir müdahaleyi de desteklemiyor. L01-P03-M01-D09'da kurulduğu gibi, bir döngüyü anlamak ona müdahalenin işe yarayacağını göstermez; döngü daha geniş bir ağın içindedir ve telafi yolları vardır. Son olarak model, kurulurken kullanılmayan bir veriyle sınanmadan bir öngörü üretmiyor; diferansiyel denklem fikri dersindeki model sınama ölçütü burada da geçerli.

Klinik gerçeklik

Denetlenen bir sistemle ilgili bir ölçüm okurken ilk soru, önündeki sayının denetlenen değişkene mi yoksa denetleyiciye mi ait olduğudur. Denetlenen değişkense, normal görünmesi denetimin şimdilik başarılı olduğunu söyler ve yük hakkında sessiz kalır. Ardından denetleyicinin de ölçülüp ölçülmediğini sor, çünkü telafiyi tek değişken göstermez.

Sonra ölçümün dinlenmede mi yük altında mı yapıldığına bak; kurgu modelde iki sistem ancak giriş arttığında birbirinden ayrıldı. Son olarak tek bir noktaya mı yoksa bir seyre mi baktığını düşün. Bu soruların hiçbiri bir sonucu yorumlamak için değil, bir sonucun neyi söyleyemeyeceğini görmek için; yorumun kendisi klinik görüşmenin işidir.

Kendini sına

1. Kurgu modelde hedef 100, başlangıç 150 ve kazanç 0,3 olsun. İlk üç adımın düzeylerini hesapla.
2. Aynı başlangıçla kazanç 0,5 olsa üç adım sonra fark kaç kalır? 0,3 kazançla karşılaştır.
3. İki adımlık gecikmeli modelde 0,3 kazançla düzey neden dördüncü adımda hedefin altına iniyor?
4. Kazanç 0,6'ya çıkınca gecikmeli modelin seyri nasıl değişti? Bu hangi çatışmayı gösteriyor?
5. Dalgalanma neden bir denetim işareti olabilir? Bir salınım gözlemek mekanizma hakkında neyi göstermez?
6. Denetlenen ve denetleyici değişkeni ayır. Hangisi en az değişiyor, neden?
7. Giriş 30 iken sinyal tavanı 25 olan ikinci sistemde düzey neden her adımda 5 birim artıyor? Bu hangi iki terimi birleştiriyor?
8. Telafi edilmiş durumu ayırt etmenin üç yolunu yaz ve her birinin kurgu modeldeki karşılığını göster.

Modülün kapanışı

Bu modül boyunca bir niceliği durduğu yerden değil değiştiği biçimden okumayı kurduk. Önce iki ölçüm arasındaki ortalama hızın aradaki yolu gizlediğini gördük, sonra türevle anlık değişimi ve integralle birikimi okuduk. Birikim ve alan dersinde sabit bir düzeyin arkasında büyük akışların olabileceğini hesapladık, uç noktalarda "en iyi"nin kimin için en iyi olduğunu sorduk. Diferansiyel denklem fikriyle değişimin kuralını yazdık, son iki derste bu kuralı bir maddenin kandaki seyrine ve bir denetim döngüsüne uyguladık.

Modülün ortak bulgusu tek bir cümleye iniyor: tek bir ölçüm bir sistemi tarif etmiyor. İki nokta yolu göstermedi, sabit depo göstergesi akışları göstermedi, bu dersteki normal düzey de denetleyicinin yükünü göstermedi. Hepsinin nedeni aynı, çünkü ölçtüğümüz şey bir durum, anlamak istediğimiz şey ise bir süreç.

Bu derste iki komşuyu ayırmanın ilk yolu, düzeyin yanına denetleyiciyi de yazmaktı: tek sayının yerine iki sayı. Gerçek bir değerlendirmede bu sayı onlarca olur ve hepsini birlikte tutmanın, birbirine karşı tartmanın ya da tek bir puana indirmenin kendi matematiği var. Sıradaki modül bu çok bileşenli bakışla açılıyor: farklı birimlerdeki ölçümleri tek bir nesnede tutmak ve onları bir puana indirirken hangi ağırlık kararının verildiği (L02-P01-M04-D01).