matlab四阶龙格库塔+打靶法

最新推荐文章于 2024-08-22 14:13:29 发布

AdaLeery

最新推荐文章于 2024-08-22 14:13:29 发布

阅读量8k

点赞数 3

分类专栏： matlab 文章标签：解微分方程

本文链接：https://blog.csdn.net/qq_37312720/article/details/93655535

版权

本文介绍了利用MATLAB的四阶龙格库塔算法解决初始值问题（IVP）解微分方程，并探讨了打靶法在边界值问题（BVP）中的应用，特别是Newton-Raphson射击法和放松法。此外，还详细阐述了MATLAB中floor、round、ceil和fix四个取整函数的区别。

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

function ys=dbf(f,a,b,a1fa,beta,h,eps)
ff=@(x,y)[y(2),f(y(1),y(2),x)];
    xvalue=a:h:b;
    n=length(xvalue)
    s0=a-0.01;
    x0=[a1fa,s0];
    flag=0;
    y0=rk4(ff,a,x0,h,a,b);
    if abs(y0(1,n)-beta)<=eps
    flag=1;
    y1=y0;
    else
        s1=s0+1;
        x0=[a1fa,s1];
        y1=rk4(ff,a,x0,h,a,b);
         if abs(y1(1,n)-beta)<=eps
             flag=1;
    end
    end
if flag~=1
    while abs(y1(1,n)-beta)>eps
        s2=s1-(y1(1,n)-beta)*(s1-s0)/(y1(1,n)-y0(1,n));
        x0=[a1fa,s2];
        y2=rk4(ff,a,x0,h,a,b);
        s0=s1;
        s1=s2;
        y0=y1;
        y1=y2;
    end
end
xvalue=a:h:b;
yvalue=y1(1,:);
ys=[xvalue',yvalue'];
function x=rk4(f,t0,x0,h,a,b)
t=a:h:b;
m=length(t);
t(1)=t0;
x(:,1)=x0;
for i=1:m-1
    L1=f(t(i),x(:,i)); 
    L2=f(t(i)+h/2,x(:,i)'+(h/2