0

多くのイプシロンのファンデルポール方程式の解をプロットしたい私のコードは次のとおりです。

tspan = [0, 10];
y0 = [0.5; 0];  % Initial location

for ep = 0.1:0.2:2.5  % Loop through a few epsilons
    ode = @(t,y) vanderpol(t,y,ep);   % Call vanderpol.m for the points (t,y)
    [t,y] = ode45(ode, tspan, y0);    % solve Van der Pol equation

    % Plot of the solution

    plot(t,y(:,1)); drawnow;
    %xlabel('Time');
    %ylabel('Van der Pol Solution');
    %title('Solutions to van der Pol equation for many \epsilon');
end

各プロットの各イプシロンのファン デル ポール方程式のプロットが必要です。これを行う方法が本当にわからない場合は、助けていただければ幸いです。

4

1 に答える 1

1

figureさまざまなプロットを作成するために使用します。

tspan = [0, 10];
y0 = [0.5; 0];  % Initial location

for ep = 0.1:0.2:2.5  % Loop through a few epsilons
    ode = @(t,y) vanderpol(t,y,ep);   % Call vanderpol.m for the points (t,y)
    [t,y] = ode45(ode, tspan, y0);    % solve Van der Pol equation

    % Plot of the solution

    figure,  

    plot(t,y(:,1)); drawnow;
    %xlabel('Time');
    %ylabel('Van der Pol Solution');
    %title('Solutions to van der Pol equation for many \epsilon');
end

単一プロットの使用hold on

tspan = [0, 10];
y0 = [0.5; 0];  % Initial location

for ep = 0.1:0.2:2.5  % Loop through a few epsilons
    ode = @(t,y) vanderpol(t,y,ep);   % Call vanderpol.m for the points (t,y)
    [t,y] = ode45(ode, tspan, y0);    % solve Van der Pol equation

    % Plot of the solution

    plot(t,y(:,1)); drawnow;
    %xlabel('Time');
    %ylabel('Van der Pol Solution');
    %title('Solutions to van der Pol equation for many \epsilon');
    hold on
end
于 2013-03-07T15:37:56.487 に答える