您现在的位置:首页 >> ★免费资源 >> 各类算法程序 >> 内容

自适应Runge-Kutta-Fehlberg法求解微分方程组

时间:2019-7-15 13:30:25 点击:

  核心提示:MATLAB代做,FPGA代做,Coursework代做、Assignment代做,淘宝交易。QQ:1224848052 ...

function varargout=saxplaxltest31(varargin)
clear,clc
a=0;b=0.2;y0=[6 4];hmax=0.02;hmin=0.0001;e=10^(-4);
[x,y]=RKF1(a,b,y0,e,hmax,hmin);
n1=length(x);
fprintf(' xn        y1n         y2n\n');
for i=1:n1
    fprintf('%f   %f     %f \n',x(i),y(i,1),y(i,2));
end
function u=f(x,y)
         n=length(y);
         u=zeros(1,n);
         u(1)=y(1)+2*y(2);
         u(2)=3*y(1)+2*y(2);
function [x,y]=RKF1(a,b,y0,e,hmax,hmin)
t=a;
x=a;
h=hmax;
i=1;
y=y0;
while (t<b)
        k1=h*f(t,y0);   
        k2=h*f(t+h/4,y0+1/4*k1);
        k3=h*f(t+3/8*h,y0+3/32*k1+9/32*k2);   
        k4=h*f(t+12/13*h,y0+1932/2197*k1-7200/2197*k2+7296/2197*k3);
        k5=h*f(t+h,y0+439/216*k1-8*k2+3680/513*k3-854/4104*k4);   
        k6=h*f(t+h/2,y0-8/27*k1+2*k2-3544/2565*k3+1859/4104*k4-11/40*k5);
        R=abs(1/360*k1-128/4275*k3-2197/75240*k4+1/50*k5+2/55*k6)/h;       
        RL=max(R);
        delta=0.84*(e/RL)^(1/4);
      if R<=e
        t=t+h;
        x=[x;t];
        y0=y0+25/216*k1+1408/2565*k3+2197/4104*k4-1/5*k5;
        y=[y;y0];
     end
    if delta<=0.1
         h=0.1*h;
    elseif delta>=4
          h=4*h;
        else
            h=delta*h;
    end
    if h>=hmax
        h=hmax;
    end
    if h<hmin
        error('Minimum h exceeded!!!');       
    end
end
% disp('The procedure is complete!');
-----------------------------------------------------------------------------
MATLAB代做,fpga代做,Simulink代做,创新,淘宝交易,MATLAB毕业设计,FPGA毕业设计,

Simulink毕业设计,MATLAB专业代做,fpga专业代做,Simulink专业代做,博士,硕士,本科,

毕业设计,硕士毕业论文,Coursework,Assignment,Verilog

QQ:1224848052

作者:网络 来源:网络
  • 您是如何找到本站的?
  • 百度搜索
  • Google搜索
  • 查阅资料过程中
  • 论坛发现
  • 百度贴吧发现
  • 朋友介绍
本站最新成功开发工程项目案例
相关文章
  • 没有相关文章
相关评论
发表我的评论
  • 大名:
  • 内容:
  • matlab代做|matlab专业代做|matlab淘宝代做|python人工智能代做|FPGA项目合作|C代做|C++代做(www.hslogic.com) © 2021 版权所有 All Rights Reserved.
  • Email:highspeed_logic@163.com 站长QQ: 1224848052

    专业代做/代写/承接、MATLAB、SIMULINK、FPGA项目、博士/硕士/本科毕业设计、课题设计、论文,毕业论文,Coursework、Eassy、Assignment