MATLAB

Question

'da ode45 kullanılırken bir işlevdeki parametreler nasıl değiştirilir Bir işlev (bir SIR modeli için) ve sonra bu işlevi çözen ve bu modele uymaya çalıştığım verilerle karşılaştıran bir komut dosyası var. Böylece fit'i optimize etmek için fonksiyondaki bir parametreyi değiştirmek için bir for döngüsü çalıştırmaya çalışıyorum. Ben onları eliyle değiştirmek zorunda kalmadan döngüsü için benim (r) ve (a) parametreleri değiştirmek için merak ediyorum:

function ydot=epidemic(t,y) 
r=0.000001; 
a=1/3; 
ydot=zeros(3,1); 
ydot(1)=-r*y(1)*y(2); 
ydot(2)=r*y(1)*y(2)-a*y(2); 
ydot(3)=a*y(2); 
end 

ve

[t,y]=ode45('epidemic',[0:222], [70500,1,0]) 

Teşekkür

Answer 1

You aşağıdaki kullanabilirsiniz: Eğer fonksiyonu

function ydot=epidemic(t,y,r,a) 
    ydot=zeros(3,1); 
    ydot(1)=-r*y(1)*y(2); 
    ydot(2)=r*y(1)*y(2)-a*y(2); 
    ydot(3)=a*y(2); 
end 

için r ve bir parametre ekleyin ve sonra gibi ode45 fonksiyonu geçmesi bu

r = 0.000001 ; 
a = 1/3 ; 
[t,y]=ode45(@(t,y)epidemic(t,y,r,a),[0:222], [70500,1,0]) 

Temel olarak, @ (t, y), epidemik (t, Y, R, a) bağımsız değişken (t, y) ve burada r ve hemen yukarıda tanımlandığı değerleri kullanan yeni bir fonksiyon tanımlandığı gibidir.

Tüm bunları bir for döngüsüne yerleştirebilirsiniz.