matlab:打靶法求解边值问题

函数1:龙格库塔方法
%书籍:常用数值算法及其matlab实现
%第10章 常微分方程边值问题的数值解法,例11.1使用
%打靶法+四阶龙格库塔方法
%调用:function [t,z] = rk4symeq(fun, t0, tf, Za, h)
function [t,w] = shootin(fun1, fun2,t0, tf, alpha, beta, h)
%fun1, fun2:边值问题转化为初值问题的函数
%t0, tf:区间
%alpha, beta,初值和终值,是一个数
%h:步长
M = floor(tf-t0)/h ;      %离散点的个数M+1
if t0 >= tf
    printf('左端点必须小于右端点');
    return;
end
N = 2;           %获得变量个数,N =2 ,二阶方程
z1 = zeros(M+1, N);
z2 = zeros(M+1, N);
t = zeros(M +1, 1);
t =[t0 : h :tf]';
z1(1,:) = [alpha; 0] ;            %假设alpha为数字,补充导数的初值,变成列向量,符合rk4symeq的要求
z2(1,:) = [0; 1 ];            %假设alpha为数字,补充导数的初值,变成列向量,符合rk4symeq的要求
[t,z1] = rk4symeq(fun1, t0, tf, [alpha; 0], h);  %第1个微分方程
[t,z2] = rk4symeq(fun2, t0, tf, [0; 1], h);   %第2个微分方程

U = z1(:,1);        %取出y1,是一个列向量
V = z2(:,1);        %取出y2,是一个列向量
w = U + (
  • 0
    点赞
  • 10
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值