数学实验“微分方程组边值问题数值算法(打靶法,有限差分法)”实验报告(内含matlab程序).doc

西京学数学软件实验任务书课程名称数学软件实验班级数 0901 学号 0912020107 姓名李亚强实验课题微分方程组边值问题数值算法( 打靶法, 有限差分法) 实验目的熟悉微分方程组边值问题数值算法(打靶法,有限差分法) 实验要求运用 Matlab/C/C++/Java/Maple/Mathematica 等其中一种语言完成实验内容微分方程组边值问题数值算法( 打靶法, 有限差分法) 成绩教师-1- 实验二十七实验报告一、实验名称: 微分方程组边值问题数值算法(打靶法,有限差分法)。二、实验目的: 进一步熟悉微分方程组边值问题数值算法( 打靶法, 有限差分法)。三、实验要求: 运用 Matlab/C/C++/Java/Maple/Mathematica 等其中一种语言完成程序设计。四、实验原理: 1. 打靶法: 对于线性边值问题??????????????)(,)( ],[)()()(byay baxxfyxqyxpy (1) 假设 L 是一个微分算子使: ( ) ( ) Ly y p x y q x y ?? ?? ? ?则可得到两个微分方程: )( 1xf Ly?,??)( 1ay ,0)( 1??ay?)()()( 111xfyxqyxpy??????,??)( 1ay ,0)( 1??ay (2) 0 2? Ly ,0)( 2?ay ,1)( 2??ay?0)()( 222??????yxqyxpy ,0)( 2?ay ,1)( 2??ay (3) 方程( 2),(3 )是两个二阶初值问题. 假设 1y 是问题( 2) -2- 的解, 2y 是问题(3) 的解,且 2 ( ) 0 y b ?, 则线性边值问题(1) 的解为: 1 1 2 2 ( ) ( ) ( ) ( ) ( ) y b y x y x y x y b ??? ?。 2. 有限差分法: 基本思想是把连续的定解区域用有限个离散点构成的网格来代替,这些离散点称作网格的节点;把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似; 把原方程和定解条件中的微商用差商来近似,积分用积分和来近似,于是原微分方程和定解条件就近似地代之以代数方程组,即有限差分方程组, 解此方程组就可以得到原问题在离散点上的近似解。然后再利用插值方法便可以从离散解得到定解问题在整个区域上的近似解。五、实验内容: % 线性打靶法 function [k,X,Y,wucha,P]=xxdb(dydx1,dydx2,a,b,alpha,beta,h) n=fix((b-a)/h); X=zeros(n+1,1); CT1=[alpha,0]; Y=zeros(n+1,length(CT1)); Y1=zeros(n+1,length(CT1)); Y2=zeros(n+1,length(CT1)); X=a:h:b; Y1(1,:)= CT1; CT2=[0,1];Y2(1,:)= CT2; for k=1:n k1=feval(dydx1,X(k),Y1(k,:)) x2=X(k)+h/2;y2=Y1(k,:)'+k1*h/2; -3- k2=feval(dydx1,x2,y2); k3=feval(dy