微分方程数值解-打靶法求边值问题.docVIP

微分方程数值解 实验七 学号 班级 姓名 指导教师 实验题目 打靶法求边值问题 评 分 实验要求： 掌握打靶法求边值问题 问题： 用打靶法求如下边值问题 2、实验内容：（求解过程，包括简单的分析，求解程序代码，所得结果） 程序代码： （1）先猜测初值，调用T7函数龙格库塔算法计算出该初值下的边界值，与题中的边界值相减，计算差值，既F(s)=yb-y(b). function y=f(a) %利用龙格库塔算法计算F(s)=yb-y(b) qj=[1,2]; cz=[0,a]; [t,w]=T7(qj,cz,0.01); y=0-w(end,2); end （2）在ydot函数中输入题目中的微分方程，用龙格库塔算法计算预设初值下的边界值。 function [t,w]=T7(qj,cz,h) %求解实验7问题 %qj为所求区间 %cz=[y1,y2]为y1,y2的初值 %h为步长 %eg：[t,w]=T7([qj,cz,0.01); n=round((qj(2)-qj(1))/h); t=zeros(n,1); w=zeros(n,length(cz)); t(1)=qj(1); w(1,:)=cz; for i=1:n t(i+1)=t(i)+h; w(i+1,:)= RungeKutta(t(i),w(i,:),