微分方程组：向后（隐）欧拉法（一、两步）matlab代码

weixin_66119178

已于 2022-02-10 19:12:34 修改

阅读量1.9k

点赞数 2

文章标签： matlab 矩阵开发语言

于 2022-02-10 19:08:24 首次发布

本文链接：https://blog.csdn.net/weixin_66119178/article/details/122852707

版权

该博客展示了如何利用后欧拉一阶法（BDF1）结合fsolve函数来求解两个耦合的一阶常微分方程。代码详细解释了时间步长设置、初始条件以及fsolve在反解过程中的应用。同时，提到了与BDF2方法的关联。

摘要由CSDN通过智能技术生成

function main() % 后欧拉一阶法
clear all
clc
hold off
T = 1
dt = 10e-5;
y0=[2,1]'; % y[0] X, y[1] Y
minx=0.2;
n = (T-minx)/dt;   %the number of required time steps
tspan = [minx,T]
a = -20201;
b = 20402;
c = 20000;
  options = optimoptions ( 'fsolve', 'Display', 'off' );
  m = length ( y0 );
  t = zeros ( n + 1, 1 );
  y = zeros ( n + 1, m );
  dt = ( tspan(2) - tspan(1) ) / n;
  t(1,1) = tspan(1);
  y(1,:) = y0(:);
  for i = 1 : n
    to = t(i,1);
    yo = y(i,:);
    tp = to + dt; 
    %size(yo) % 初始的yo
    yp = yo + dt *  ( f (yo ) );
    rhs= f(yp) ;
 [yp, fval] = fsolve(@(yp)fsolve_test_1_function_1(yp,yo,dt,rhs), yp, options);
   t(i+1,1) = tp;
    y(i+1,:) = yp;
%     return 
  end
  plot(t,y(:,1),'r-')
  hold on
  plot(t,y(:,2),'b-')
  title ("BFD1")
  
function  res= f( yp)
a = -20201;
b = 20402;
c = 20000;
res (1,1)= a*yp(1,1) + b *yp(1,2);
res (1,2)= c*yp(1,1) + a *yp(1,2) ;

function F = fsolve_test_1_function_1(yp,yo,dt, rhs)
F(1) = yp(1)-yo(1)- dt* rhs(1);
F(2)= yp(1)-yo(1)- dt* rhs(2);

技术难题：利用fsolve 反解y(n+1)

题目是：两个一阶微分方程

上面的是BDF1, 至于BDF2可联系博主。