3、80040001430145014701320128012001080940145014801500155015101430130012001460150015501600155016001600160013701500120011001550160015501380127015001200110013501450120011501230139015001500140090011001060118013201450142014001300700900Y/X040080012001600200024002800三、实验步骤1．开启软件平台——MATLAB，开启MATLAB编辑窗口
4、；2．根据各种数值解法步骤编写M文件3．保存文件并运行；4．观察运行结果(数值或图形)；5．写出实验报告，并浅谈学习心得体会。实验题目：1.在函数y=cos10x,-2£x£2上取n个结点（n不要太大，如5~11），用拉格朗日，分段线性，三次样条三种插值方法，计算m个插值点的函数值（m要适中，如50~100），通过数值和图形输出，将三种插值结果与精确值比较。适当增加n,再作比较。拉格朗日插值方法的函数m文件：functiony=lagr(x0,y0,x)n=length(x0);y=0.0;fork=1:np=1.0;forj=1:nifj~=kp=p.*(x-x0(j))/(
5、x0(k)-x0(j));endendy=p*y0(k)+y;end比较三种插值方法的程序m=99;x=-2:4/(m-1):2;y=cos(x).^10;z=0*x;k=1;forn=[5,7,10,11]figure(k)plot(x,z,'c',x,y,'k','Linewidth',2),x0=-2:4/(n-1):2;y0=cos(x0).^10;y1=interp1(x0,y0,x);y2=lagr(x0,y0,x);y3=interp1(x0,y0,x,'spline');holdon,plot(x,[y1;y2;y3],'Linewidth',2),pause,
6、holdoffk=k+1;end图形为：2.确定地球与金星之间的距离天文学家在1914年8月份的7次观测中，测得地球与金星之间距离（单位：米），并取其常用对数值，与日期的一组历史数据如表7.2。表7.2日期（号）18202224262830距离对数9.96177249.95436459.94680699.93909509.93122459.92319159.9149925由此推断何时金星与地球的距离（米）的对数值为9.9351799？算法程序：day0=18:2:30;s0=[9.96177249.95436459.94680699.93909509.93122459.9231
7、9159.9149925];day=18:0.001:30;sequ=9.9351799;s=interp1(day0,s0,day);plot(day0,s0,'r',day,s),holdonplot(day,sequ,'g'),title('地球与金星距离随时间变化图');xlabel('时间day'),ylabel('距离s');grid运行结果：图3：地球与金星距离随时间变化图3．山区地貌图在某山区（平面区域（0，2800）´（0，2400）内，单位：米）测得一些地点的高程（单位：
