m0_58956357 2021-06-25 17:15 采纳率: 100%
浏览 100
已采纳

编制matlab程序,R-K公式求解初值问题,matlab编程

分别用2阶和3阶R-K公式求解初值问题(步长h=0.05)

  • 写回答

1条回答 默认 最新

  • CSDN专家-Matlab_Fans 2021-06-25 19:42
    关注
    dfun = @(x,y) cos(x)*sqrt(y);
    h = 0.05;
    
    h = 0.05;
    X = 1.2:h:2.05;
    %% RK2
    x = 1.2; % 初始值
    y = 3.2; 
    for ii = 1:length(X)
       x = X(ii); 
       k1 = dfun(x,y);
       k2 = dfun(x+h,y+h*k1);
       y = y + h*(k1+k2)/2;
       Y1(ii) = y;
    end
    
    %% RK3
    x = 1.2; % 初始值
    y = 3.2; 
    for ii = 1:length(X)
       x = X(ii); 
       k1 = dfun(x,y);
       k2 = dfun(x+h/2,y+h*k1/2);
       k3 = dfun(x+h,y+h*(2*k2-k1)/2);
       y = y + h*(k1+4*k2+k3)/6;
       Y2(ii) = y;
    end
    
    figure
    plot(X,Y1,X,Y2,'r--')
    legend('RK2','RK3')
    ylabel('y')
    xlabel('x')

    结果

     

    本回答被题主选为最佳回答 , 对您是否有帮助呢?
    评论

报告相同问题?

悬赏问题

  • ¥15 MATLAB yalmip 可转移负荷的简单建模出错,如何解决?
  • ¥15 数学的三元一次方程求解
  • ¥20 iqoo11 如何下载安装工程模式
  • ¥15 本题的答案是不是有问题
  • ¥15 关于#r语言#的问题:(svydesign)为什么在一个大的数据集中抽取了一个小数据集
  • ¥15 C++使用Gunplot
  • ¥15 这个电路是如何实现路灯控制器的,原理是什么,怎么求解灯亮起后熄灭的时间如图?
  • ¥15 matlab数字图像处理频率域滤波
  • ¥15 在abaqus做了二维正交切削模型,给刀具添加了超声振动条件后输出切削力为什么比普通切削增大这么多
  • ¥15 ELGamal和paillier计算效率谁快?