河南大学数值分析作业.docVIP

微分方程数值解法实验报告 主讲教师 葛志昊 指导教师 葛志昊 实验日期 2017.3.16 专业 研究生 课程名称 微分方程数值解法 同组实验者 一、实验名称： 实验一、两点边值问题有限差分法 实验目的： 掌握两点边值问题的有限差分法的程序设计 求解初值问题 其中，q(x)取 要求 给出该问题的有限差分格式，并给出该格式的局部截断误差 解：将求解区间进行等分，设步长,节点， 该问题在点的有限差分格式为： 局部截断误差为：，其中表示在处的函数值。 （2）取 ，j=2, 3, 4, 5, 6, 7, 8,写程序计算该问题 %两点边值问题的有限差分方法 % clear all;clc; %% x0=0;xN =pi; %求解区间 time=1; %43行会用到 for n1=2:8; hx=pi/2^n1; x=x0:hx:xN; n2=length(x)-2; %% a1=1+2/hx^2; a2=-1/hx^2; v1=sparse(1:n2,1:n2,a1*ones(1,n2),n2,n2)... +sparse(1:n2-1,2:n2,a2*ones(1,n2-1),n2,n2)... +sparse(2:n2,1:n2-1,a2*ones(1,n2-1),n2,n2); %系数矩阵 u(1)=0; u(n2