两点边值问题的有限差分法.doc

《两点边值问题的有限差分法.doc》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、.学生实验报告实验课程名称偏微分方程数值解开课实验室数统学院学院数统年级2013专业班信计2班学生姓名学号开课时间2015至2016学年第2学期总成绩教师签名word范文.数学与统计学院制开课学院、实验室：数统学院实验时间：2016年月日实验项目名称两点边值问题的有限差分法实验项目类型验证演示综合设计其他指导教师曾芳成绩是一．实验目的通过该实验，要求学生掌握求解两点问题的有限差分法，并能通过计算机语言编程实现。二．实验内容考虑如下的初值问题：(1)(2)其中，，，，是给定常数。将区间等分，设，网点。1．在第三部分写出问题（1）和（2）

2、的差分格式，并给出该格式的局部截断误差。2．根据你写出的差分格式，编写一个有限差分法程序。将所写程序放到第四部分。3．给定参数，，，，问题(1)的精确解，其中将及带入方程(1)可得。分别取，用所编写的程序计算问题(1)和(2)。将数值解记为，，网点处精确解记为，。然后计算相应的误差word范文.，及收敛阶，将计算结果填入第五部分的表格，并对表格中的结果进行解释？4.将数值解和精确解画图显示，每种网格上的解画在一张图。三．实验原理、方法（算法）、步骤1.差分格式：=-1/h^2(-()+)+()/2h+=A,2.局部阶段误差：(u)=O

3、(h^2)3.程序clearallN=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

4、)*h)/h+r(a+i*h))*alpha;endfori=2:N-2word范文.H(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*(p(a+(i+1/2)*

5、h)+p(a+(i-1/2)*h))/h+2*h*q(a+i*h);g(i)=2*h*f(a+i*h)+(2*p(a+(i+1/2)*h)/h-r(a+i*h))*beta;endu=Hg;u=[alpha;u;beta];x=a:h:b;y=(x.^2).*exp(x-1);plot(x,u);holdonplot(x,y);y=y'z=y-u四．实验环境（所用软件、硬件等）及实验数据文件Matlabword范文.五．实验结果及实例分析N收敛阶收敛阶100.00104256……0.00073524……200.000261681.9

6、3410.000183481.4530400.000065412.00010.000045852.0000800.000016361.99930.000011462.00001600.000004092.00000.000002872.0000N越大只会使绝对误差变小，方法没变，所以收敛阶一致。图示为：(绿线为解析解，蓝线为计算解)N=10N=20word范文.N=40N=80word范文.N=160word范文.教师签名年月日word范文.欢迎您的光临，word文档下载后可以修改编辑。双击可以删除页眉页脚。谢谢！单纯的课本内容，并不

7、能满足学生的需要，通过补充，达到内容的完善教育之通病是教用脑的人不用手，不教用手的人用脑，所以一无所能。教育革命的对策是手脑联盟，结果是手与脑的力量都可以大到不可思议。word范文