微分方程数值解 实验七学号班级姓名指导老师实验题目打靶法求边值问题评 分实验要求： （1）掌握打靶法求边值问题问题： 用打靶法求如下边值问题2、实验内容:（求解过程，包括简单得分析,求解程序代码，所得结果）程序代码：(1）先猜想初值,调用 T7 函数龙格库塔算法计算出该初值下得边界值，与题中得边界值相减,计算差值，既 F(s）＝yb—y（b)、ｆｕｎ c ｔ ion ｙ＝f（ａ）％利用龙格库塔算法计算Ｆ（ｓ）=ｙ b—ｙ(b）qj＝[1,2］；ｃ z=[０，ａ］；［ｔ,w]=Ｔ７（qj，ｃｚ，0、01）；ｙ=0—w（ｅｎ d，2）;end(2）在 ydot 函数中输入题目中得微分方程,用龙格库塔算法计算预设初值下得边界值。functio ｎ ［t，w］=T7(ｑ j，cz，h）％求解实验 7 问题%q ｊ为所求区间%c ｚ=[y1,y2］为ｙ１,y2 得初值％h 为步长％e ｇ:[t，w]＝Ｔ７(［ｑ j，cz，0、0 １); n＝r ｏ und((qj（2）—qj（1))/h)； t=zeros(ｎ,1）; w=zeros（n,l ｅ ngth(cz）)； t(1)=ｑｊ（1）； w（1，:）=cz； f ｏ r i=1：n ｔ(ｉ+1)=ｔ(i)＋h； w(i+1,：)＝ Rung ｅＫ u ｔ ta(t（ｉ）,w（i，：)，ｈ）; endendfu ｎ ct ｉ on y ＝ RungeKutta(t，w，h)％龙格—库塔法求解微分方程％ｈ为步长 k1= ydot（ｔ,ｗ）; k2= yd ｏ t(t+h／2，w+h／2*ｋ１）； k ３= ydo ｔ（ｔ＋ｈ／2，w+ｈ/2＊ｋ 2)； k4= ydot（ｔ＋h，w+h*ｋ 3）； y＝w＋h/6＊（k １＋２*k2+2＊ｋ３＋k ４)；ｅｎ df ｕ nct ｉｏｎ ｚ=yd ｏ t(t,y)％输入带求解得微分方程z(1） = (4-2*y（2))/(t＾３）；z(2） = -exp（y（1））;ｅ nd（3）利用二分法不断修正，调整初值,直到打靶打中边界值.fu ｎ ction ｓ ol=db ｆ(ａ，b，tol)％打靶法求解边值％ｆ:F(s）=ｙ b-ｙ（ｂ)；%a，ｂ：猜想初值%ｔ ol：容忍误差％ｓｏｌ:方程 f（x)=0 得根%eg:s ｏ l=ｄｂｆ（０，２、5,0、00 １)if si ｇｎ（f（a)）＊sign(ｆ（b）)>０ ｅ rror('Pleas ｅ let f(a）*ｆ(b)<0＇)ｅｎｄi＝1;k（ｉ，1）=a;k(i，２)＝b;ｗ h ｉ l ｅ ａ bs（ｋ(ｉ，2）-k（i，1)）>t ｏ l if f((k(i，１）+k（i，2））/２）＝=0 br ｅ ak ｅｎ d if f（k(ｉ,1））*f（（k（i，1)+ｋ（i,2）)／２)〈０ k(i+１，1）=k(i,1)； ｋ（i+１,2）=（k(i，１)+k(i，２））/２; else k(i+1,2)=k（i，2）; k(ｉ+1,1）=(k(i，1)＋k（i...
