Plot 人口增长 function 在 Matlab 中使用不同的 b 值使用 for 循环

Question

我需要 plot 人口增长 model 具有不同的 b 值，其中 b>d 和 b=d 使用 for 循环。 两个结果应该在具有不同 colors 的同一图表上可见。这是我的 function 和初始值，但我不知道如何获得 plot。

    [T,Y]=ode45(H,[0:1:time],[N0]);

    Intitial value 
   b= 0.1 0.001
   d=0.001
   N0=400````

I need two line for each b on the same plot.
```***
    Intitial value 
    b= 0.1 0.001;
    d=0.001;
    N0=400;

    H=@(t,N)[b*N-d*N];
    [T,Y]=ode45(H,[0:1:time],[N0]);

    plot(T,Y)
***```

Answer 1

首先，您需要以正确的格式放置代码，现在您复制了两次初始值块（我猜您在放置代码格式时遇到了一些问题。我希望您在提出新问题时可以关注这一点。

当前代码在 b 的定义上出错，因为您没有将其定义为向量。 要有两条线，最好将 output Y 存储为两列矩阵（每列是一个向量）。 MATLAB 然后将每一列输出为 plot 中的一行。

以下代码解决了您的问题：

% Initial value
b = [0.1, 0.001];
d = 0.001;
N0= 400;
simTime = 0:1:300;

% Preallocate output for the loop
Y = zeros(length(simTime), length(b));

% Run for loop
for i = 1:length(b)
    b_i = b(i);
    H=@(t,N)[b_i*N-d*N];
    [T,y]=ode45(H, simTime, [N0]);
    Y(:,i) = y;
end

% Plot output
plot(T,Y)