다양한 파라미터

Question

와 매트랩과 미분 방정식을 해결 제가 단순 로지스틱 방정식을 말하게

DX/DT = 2ax (1 - X/N)N가 수용력이다 , 어떤 성장률이며, 내가 변화하고 싶습니다 매개 변수입니다 모두 및 N.

그래서 내가하고 싶은 것은 고정 점과 두 매개 변수의 3D 그래프를 그려 보는 것입니다.

단일 매개 변수의 고정 점을 찾는 방법을 이해합니다. 여기

날 다를 수 있습니다 내 샘플 코드 나는 두 가지 기능을 통해 루프에 대한 을 넣어 어떻게

% Initial Number 
x0 = 0.4; 

%Time Window 
tspan=[0 100]; 

[t,x]=ode45(@MyLogisticFunction,tspan,x0); 

clf 

x(end,1) % This gives me the fixed point for the parameters above. 

그래서 내 진짜 질문은

function xprime = MyLogisticFunction(t,X) %% The ODE 

% Parameters 
N = 10 % Carrying Capacity 
a = 0.5 % Growth Rate 

x1prime = 2*a*X(1)*(1 - X(1)/N); 

xprime = [x1prime ]'; 

end 

다음 내 해결사,이다 a 및 n이고, 따라서 및 N 및 고정 소수점 x *의 3D 그래프를 그릴 수 있습니다.

내가 일 개하는 .m 파일에 두 가지 기능을 결합하려고했습니다하지만 당신은 함수에 매개 변수를 전달해야

Answer 1

를 작동하지 않는 것 : 당신이 전화를 한 후

function xprime = MyLogisticFunction(t,X,a,N) %% The ODE 

% Parameters (passed as function arguments) 
% N = 10 % Carrying Capacity 
% a = 0.5 % Growth Rate 

x1prime = 2*a*X(1)*(1 - X(1)/N); 

xprime = [x1prime ]'; 

end 

과 ode 해결 자 :

% Initial Number 
x0 = 0.4; 

%Time Window 
tspan=[0 100]; 

a = 0.1:0.1:1; % or whatever 
N = 1:10; % or whatever 
x_end = zeros(length(a),length(N)); 

for ii = 1:length(a) 
    for jj = 1:length(N) 
     [t,x]=ode45(@(t,X)MyLogisticFunction(t,X,a(ii),N(jj)),tspan,x0); 
     x_end(ii,jj) = x(end,1); 
    end 
end