偏微分方程实验1-张伟-20122058范文.doc

重 庆 大 学 学 生 实 验 报 告 实验课程名称 偏微分方程 数统学院 学 院 数 统 年级 2012 专业班 信计1班 学 生 姓 名 张 伟 学 号 开 课 时 间 2014 至 2015 学年第 2 学期 总 成 绩 教师签名 数学与统计学院制 开课学院、实验室： 数统学院 2015年 5月6日 实验项目 名 称 初值问题的Euler方法和梯形法 实验项目类型 验证 演示 综合 设计 其他 指导教师 李茂军 成 绩 是 一．实验目的 通过该实验，要求学生掌握初值问题的欧拉法和梯形法并能通过计算机语言编程实现两种算法。如下的初值问题： 该问题有解析解欧拉法求解，取步长将四种步长的结果时刻的解析结果的误差列表显示 用梯形法求解该问题，取步长将四种步长的结果时刻的解析结果的误差列表显示 在同一种方法下，请说明种网格大小的计算结果更加精确，并说明理由。相同的网格大小下，比较上述两种算法的计算结果，那种算法的结果要好一些，并说明理由。 ：放到部分将电子版的实验报告发送邮箱limj@cqu.edu.cn,标题为：方程实验报告-姓名学号纸质的实验报告在下次上机提交欧拉法迭代格式误差 ，。 2. 梯形法的迭代格式及其误差估计。，。 四．实验环境（所用软件、硬件等）及实验数据文件 Matlab 或者Fortran 1、用欧拉法求解f=exp(x); plot(x,f,g-*) box on grid on hold on [s,t]=euler(fun1,0,1,1,5) %这里5=1/0.2 plot(s,t,r-o) xlabel(x) ylabel(y) title(精确解与欧拉法数值解) legend(exp(x),欧拉法数值解) disp( 步长h 精确解u 数值解u 误差估计O(h)); disp([ x (exp(x)) t abs((t-exp(x)))]) 调用子程序一 %欧拉方法求解常微分方程 function [x,y]=euler(fun,x0,xfinal,y0,N) if nargin5,N=50;%区间长度 end h=(xfinal-x0)/N;%步长 x(1)=x0; y(1)=y0; for i=1:N x(i+1)=x(i)+h; y(i+1)=y(i)+h*feval(fun,x(i),y(i)); end 调用子程序二 function f=fun1(x,y) %函数定义,可以据情况更改 f=y; 2、用梯形法求解该问题f=exp(x); plot(x,f,g-*) box on grid on hold on [s,t]=tixingfa(fun1,0,1,1,5) plot(s,t,r-o) xlabel(x) ylabel(y) title(精确解与梯形法数值解) legend(exp(x),梯形法数值解) disp( 步长h 精确解u 数值解u 误差估计O(h)); disp([ x (exp(x)) t abs((t-exp(x)))]) 调用子程序一(子程序二同上) %梯形法求解常微分方程 function [x,y]=tixingfa(fun,x0,xfinal,y0,N) if nargin5,N=50;%区间长度 end h=(xfinal-x0)/N;%步长 x(1)=x0;y(1)=y0; for i=1:N x(i+1)=x(i)+h; y1=y(i)+h*feval(fun,x(i),y(i)); y2=y(i)+h*feval(fun,x(i+1),y1); y(i+1)=(y1+y2)/2; end 五．实验结果及实例分析 1、用欧拉法求解 0 1.0000 1.0000 0 0.2000 1.2214 1.2000 0.0214 0.4000 1.4918 1.4400 0.0518 0.6000 1.8221 1.7280 0.0941 0.8000 2.2255 2.0736 0.1519 1.0000 2.7183 2.4883