回覆列表
  • 1 # 使用者5635123734932

    求解二階微分方程,初始條件還需要給出y1"(0)和y2"(0)。這裡暫時按照0處理。

    function zd530003514

    a=0.1;

    b=0.1;

    Y0 = [b-1; 0; b; 0];

    % 解方程

    [t,Y]= ode45(@ode,[0 10],Y0);

    y1=Y(:,1);

    y2=Y(:,3);

    % 繪圖

    subplot 211

    plot(t,y1);

    subplot 212

    plot(t,y2);

    % 微分方程定義

    function dY = ode(t, Y)

    L1=5;

    L2=0.01;

    a0=2;

    b0=2;

    c0=2;

    y1=Y(1);y2=Y(3);

    dY = [

    Y(2);

    -(a0*y2+b0*y2^2+c0*y2^3) - L1^2*L2*y1 - L1^2*y1;

    Y(4);

    -(a0*y2+b0*y2^2+c0*y2^3) - L1^2*L2*y1;

    ];

  • 中秋節和大豐收的關聯?
  • 咋解二元一次方程,要詳細點的?