计算方法第七章上机报告..docVIP

实验报告名称 求解常微分方程 班级：020991 学号 姓名：杜凡 成绩： 1实验目的 2 实验内容 实习题7.1用改进欧拉法与四阶龙格-库塔公式求解所给微分方程初值问题； 7.2 用四阶龙格-库塔公式解下列微分方程初值问题； 7.3用阿当姆斯方法解微分方程初值问题； 3实验步骤 // 实验环境的配置，例如添加什么函数，库，头文件等，以及你的思路都可以写。 3 // 程序流程图、代码function Heun2(f,a,b,y0,n) h=(b-a)/n; x=a:h:b; %ytrue=f1(-1*x); y=y0*ones(1,n+1); for j=2:n+1 yp=y(j-1)+h*f(x(j-1),y(j-1)); yc=y(j-1)+h*f(x(j),yp); y(j)=(yp+yc)/2; end for i=1:n+1 fprintf(x[%d]=%f\t y[%d]=%f\n,i-1,x(i),i-1,y(i)); %fprintf(x[%d]=%f\t y[%d]=%f\t ytrue[%d]=%f\n,i-1,x(i),i-1,y(i),i-1,%ytrue(i)); end 4实验结果 // 程序运行的结果，可以添加截图以说明问题。//实验结果分析，包括误差分析和结论。 2）实验结果分析 改进欧拉公式的局部截断误差O（h^3），h=0.1,则绝对误差e1.0*10^2. 四阶龙格库塔方法的局部截断误差为O(h^5),h=0.1,则绝对误差e1.0*10^4. 阿当姆斯方法的局部截断误差为O(h^5),h=0.1,则绝对误差e1.0*10^4. 5总结 // 通过掌握的内容，以及在实验中遇到的问题及解决方法。 6参考资料 // // 英文参考文献格式 // 期刊 // [序号] 主要责任者. 文献题名[J]. 刊名, 年, 卷(期): 起止页码. // 专著、论文集、学位论文、报告 // [序号] 主要责任者. 文献题名[文献类型标识]. 出版地: 出版者, 出版年. 起止页码.  计算方法与实习实验报告