文档介绍：学生实验报告实验课程名称偏微分方程数值解开课实验室数统学院学院数统年级2013专业班信计2班学生姓名学号开课时间2015至2016学年第2学期总成绩教师签名数学与统计学院制开课学院、实验室:数统学院实验时间:,要求学生掌握求解两点问题的有限差分法,并能通过计算机语言编程实现。:(1)(2)其中,,,,是给定常数。将区间等分,设,网点。(1)和(2)的差分格式,并给出该格式的局部截断误差。,编写一个有限差分法程序。将所写程序放到第四部分。,,,,问题(1)的精确解,其中将及带入方程(1)可得。分别取,用所编写的程序计算问题(1)和(2)。将数值解记为,,网点处精确解记为,。然后计算相应的误差,及收敛阶,将计算结果填入第五部分的表格,并对表格中的结果进行解释?,每种网格上的解画在一图。、方法(算法)、步骤差分格式:Lhui=-1/h^2(pi+1/2ui+1-(pi+1/2+pi-1/2)ui+pi-1/2ui-1)+ri(ui+1-ui-1)/2h+qiui=fiu0=A,uN=:Ri(u)=O(h^2)=10;a=0;b=1;p=(x)1;r=(x)2;q=(x)3;alpha=0;beta=1;f=(x)(4*x^2-2)*exp(x-1);h=(b-a)/N;H=zeros(N-1,N-1);g=zeros(N-1,1);%fori=1H(i,i)=2*(p(a+(i+1/2)*h)+p(a+(i-1/2)*h))/h+2*h*q(a+i*h);H(i,i+1)=-(2*p(a+(i+1/2)*h)/h-r(a+i*h));g(i)=2*h*f(a+i*h)+(2*p(a+(i-1/2)*h)/h+r(a+i*h))*alpha;endfori=2:N-2H(i,i-1)=-(2*p(a+(i-1/2)*h)/h+r(a+i*h));H(i,i)=2*(p(a+(i+1/2)*h)+p(a+(i-1/2)*h))/h+2*h*q(a+i*h);H(i,i+1)=-(2*p(a+(i+1/2)*h)/h-r(a+i*h));g(i)=2*h*f(a+i*h);endfori=N-1H(i,i-1)=-(2*p(a+(i-1/2)*h)/h+r(a+i*h));H(i,i)=2