数值计算方法实验报告(西安工程大学附有程序).docxVIP

数 值 分 析上机实验报告 学院：理学院 专业：数学与应用数学 学号：41108040110 姓名：吕炳林实验报告1例1实验程序I0 = exp(-1)*pond(‘x.^0.*exp(x.^2)’,0,1) I0 =0.5381 vpa(I0,10) 实验结果 ：ans =0.5380795164例1.9实验程序：function p=p(n,x)a(1)=13;fori=1:n+1a(i+1)=2*a(i)+3;end S=a(n+1);for j=1:n S=x*S+a(n+1-j);end p=S;End在命令窗口进行计算：实验结果 p(100,0.5)ans = 600 p(150,13)ans = 1.0995e+213习题3实验程序Function y=y(n)Y=28fori=1:ny=y-1/100*sqrt(783);End在命令窗口进行计算实验结果y(100)y= 28ans=0.0179y(500)y= 28ans= -111.9107 实验报告2 习题3实验原理：一次插值：二次插值：实验运行环境： 本实验采用Matlab编写。实验程序： 线性插值： x=0.4:0.1:0.8; f=[-0.916291,-0.693147,-0.510826,-0.357765,-0.223144];format long interp1(x,f,0.54)实验结果：ans =-0.620218600000000二次插值：x=0.54;a=[0.4,0.5,0.6];b=[-0.916291,-0.693147,-0.510826];l=b(1)*(x-a(2))*(x-a(3))/((a(1)-a(2))*(a(1)-a(3)));m=b(2)*(x-a(1))*(x-a(3))/((a(2)-a(1))*(a(2)-a(3)));n=b(3)*(x-a(1))*(x-a(2))/((a(3)-a(1))*(a(3)-a(2)));y=l+m+n实验结果：y =-0.61531984000000 习题21实验运行环境： 本实验采用Matlab编写实验程序1：（输入函数y）for x=-4.5:4.5y=1/(x^2+1)end显示结果y=0.04705882352941；y=0.07547169811321；y=0.13793103448276； y=0.30769230769231；y=0.80000000000000；实验程序2（求的值）x=input(请输入x的值);a=[x-0.5,x+0.5];y=[1/(1+(x-0.5)^2),1/(1+(x+0.5)^2)];I=y(1)*(x-a(2))/(a(1)-a(2))+y(2)*(x-a(1))/(a(2)-a(1))显示结果当分别输入时I=0.0486 I=0.0794 I=0.1500I=0.3500 I=0.7500习题241(2)实验报告3习题22实验程序x=[19,25,31,38,44]y=[19.0,32.3,49.0,73.3 97.8][p,s]=polyfit(x,y,2)实验结果 p = 0.0497 0.0193 0.6882s = R: [3x3 double]df: 2normr: 0.1145习题23。实验程序t=[0,0.9,1.9,3.0,3.9,5.0];s=[0,10,30,50,80,110];[P,S]=polyfit(t,s,5)显示结果P = -0.5432 6.4647 -26.5609 46.1436 -13.2601 -0.0000S = R: [6x6 double]df: 0normr: 1.2579e-012实验结果：习题 24.实验程序： x=0:5:55;y=[0,1.27,2.16,2.86,3.44,3.87,4.15,4.37,4.51,4.58,4.62,4.64]; [P,S]=polyfit(x,y,5)显示结果P = 0.0000 -0.0000 0.0002 -0.0084 0.2851 0.0082S = R: [6x6 double]df: 6normr: 0.0487实验结果：实验报告4习题4实验程序先求精确值：x=sym(x);f=exp(-x);I=int(f,x,0,1)I = 1-exp(-1)利用simpson公式得；a=0;b=1;S=(b-a)/6*[exp(-a)+4*ex