计算方法综合实践.docVIP

应用自己熟悉的算法语言编写程序，使之尽可能具有通用性。 上机前充分准备，复习有关算法，写出计算步骤，反复检查，调试程序。（注：在练习本上写，不上交） 完成计算后写出实验报告，内容包括:算法步骤叙述，变量说明，程序清单，输出计算结果，结构分析和小结等。（注：具体题目具体分析，并不是所有的题目的实验报告都包含上述内容！） 独立完成，如有雷同，一律判为零分！ 上机期间不允许做其他任何与课程设计无关的事情，否则被发现一次扣10分，被发现三次判为不及格！非特殊情况，不能请假。旷课3个半天及以上者，直接判为不及格。 目 录 一、基本技能训练 4 1、误差分析 4 2、求解非线性方程 4 3、插值 4 4、数值积分 4 二、提高技能训练 4 1、 4 2、 4 三、本课程设计的心得体会（500字左右） 4 一、基本技能训练 1、误差分析 实验 1.2 误差传播与算法稳定性 实验目的：体会稳定性在选择算法中的地位。误差扩张的算法是不稳定的，是我们所不期望的；误差衰减的算法是稳定的，是我们努力寻求的，这是贯穿本课程的目标。 问题提出： 4、数值积分 实验目的： 掌握Romberg积分法 实验内容： 用Romberg积分法计算下列定积分： function [R,y]=Romberg(a,b,n) k=1; while 1 T=zeros(k+1,k+1); T(1,1)=1/2*(b-a)*(f(a)+f(b)); for i=1:k h=(b-a)/2^i; s=0; for j=1:2^(i-1) s=s+f(a+(2*j-1)*h); end T(i+1,1)=T(i,1)/2+h*s; end for j=1:k c=1/(4^j-1); for m=j:k T(m+1,j+1)=T(m+1,j)+c*(T(m+1,j)-T(m,j)); end end if abs(T(k,k)-T(k+1,k+1))n break; end k=k+1; end y=T(k,k); R=T; 输入积分公式 （1） function y=f(x) y=1/(1+x^4); [R,y]=Romberg(0,1,10) R = 0.7500 0 0.8456 0.8775 y = 0.7500 （2） function y=f(x) y=cos(sin(x)); [R,y]=Romberg(0,pi/4,10) R = 0.6912 0 0.7099 0.7161 y = 0.6912 （3） function y=f(x) y=exp(x)/x; [R,y]=Romberg(1,2,10) R = 3.2064 0 3.0971 3.0607 y = 3.2064 （4） function y=f(x) y=x.^x*log(3+x); [R,y]=Romberg(2,3,10) R = 27.4076 0 22.1271 20.3669 y = 27.4076 龙贝格计算量较小 二、提高技能训练 1、（2）实验目的与要求 （1）了解 Kepler 方程； （2）能够用牛顿迭代方法计算 Kepler 方程； （3）通过调节轨道偏心率 e 查看运算收敛情况。 实验内容及数据来源 编写解 kepler 方程的方法，给定偏心率为 0.01、平近点角为 32 度时计算出偏近点角。 clear; x0=32; syms E; f=inline(E-0.01*sin(E)-32); g=inline(1-0.01*cos(E)); x1=x0-f(x0)/g(x0); y=abs(x1-x0); x0=x1;i=i+1; while y10^(-5) x1=x0-f(x0)/g(x0); y=abs(x1-x0); x0=x1; i=i+1; end E=vpa(x1) 结果：E =32.005560568887858607922680675983 通过调节轨道偏心率 e 查看运算收敛情况； e 0.001