matlab上机实验报告.docVIP

上绛电力嗲紇 数值计算方法上机实习 报告 院 系 专业年级 学生姓名学 号指导教师能源与机械丁.程学院动力机械及丁?程2012级 院 系 专业年级 学生姓名 学 号 指导教师 能源与机械丁.程学院 动力机械及丁?程2012级 张亚杰 VS1210124014 黄建雄 2012年12月26日 数值计算方法上机实习题 由递推公式八=-5/^+1，从/0的几个近似值出发，计算/2(); n r\ 1 解：1{)= =0.1823 Jo 5 + x 计算120编輯matlab命令如下： 1=0.1823 for n=l:1:20, I=-5*I+l/n； fprintf(%.Id %.4f\n,n ,I); end 0.0885 0.0575 0.0458 0.0208 0.0958 -0.3125 1.7054 -8.4018 42. 1200 -210.5002 1052.5919 -5262.8762 26314.4578 -131572.2175 657861. 1542 -3289305.70836002 9455 411163214.7803 -2055816073.8513 粗糙估计/，。，用—丄八+丄，计算/。； 5 5/7 解：l20= — 5 + x 使川复合中点公式进行积分，相应的matlab程序如下: 1 = 0; for h=0:0.001:1, m=h+0.0005; 1=1+0.001*m^20/(5+m)； fprintf(1 %.Id %.4f\n1,m ,I)； end disp(I)； for k=l:20, n=21-k； 1=0.2*(1/n-I)； fprintf(%.Id %.4f\nz n ,I)； end disp(I) 结果： ae-ui U. UUl 8e-0l 0.0052 8e-0l 0.0054 8e-0l 0.0055 20 0 .0084 8e-0l 0.0056 19 0 .0089 8e-0l 0.0057 18 0 ,0093 8e-0l 0.0058 17 0 .0099 9e-0l 0.0060 16 0 .0105 9e-0l 0.0061 15 0 .0112 9e-0l 0.0062 14 0 .0120 9e-0l 0.0064 13 0 ? 0130 9e-0l 0.0065 12 0 .0141 9e-0l 0.0066 11 0 .0154 9e-0l 0.0068 10 0 .0169 9e-0l 0.0069 9 0. 0188 9e-0l 0. 0071 8 0. 0212 9e-0l 0. 0072 — 0. 0243 0e+00 0. 0074 6 0. 0285 0e+00 0.0075 5 0. 0343 0e+00 0. 0077 4 0. 0431 0e+00 0. 0078 3 0. 0580 0e+00 0.0080 2 0. 0884 0e+00 0.0082 1 0. 1823 0.0082 0. 1823 程序结朿时输出两个I值，第一个表示120,笫二个表示 分别为 12()=0.0082 1()=0.1823 分析结果的可靠性及产生此现象的原因(重点分析原因)。 从上述计算屮分析得到如果先得到10，再从IQ由递推公式得到I2o, I2e结果跟精 确值相比误差很人；如果先估算I2C,在从ho冇递推公式彳?到Ie，Ie的结果跟精确值 相比近似相等。 原因分析： 如果从10推I2Q的近似值，需要用到递推公式In=-5Ind+l/n, 10本身结果是有误 差的：经过递推公式计算20次，就等于误差被认为的放大5的20次方倍，所以得 到的I2Q与其精确值相差其远。 如果从I2o推Io的近似值,需要用到Ind=0.2 (l/n-In),尽管I2Q本身有误差，但 是经过20次运算，其误差缩小到原來的0.2的20次方估，所以得到的Io与其精确 值比较相近。 2.求方程eA +10x —2 = 0的近似根，要求|人+1-a|5x10并比较计算量。 (1) 在［0, 1］上用二分法； Matlab程序如下: a=0; b=l ； c=b-a； n=0 while c0.0005, x=(a+b)/2； f=exp(x)+10*x-2； if f0, b=x； c=b-a; elseif f 0, a=x； c=b-a； else x=x； c=0 ； end n=n+l； fprintf(1 %. Id % . 4f end %.4f\n1,n ,x,c); 结果如下： 1 0. 5000 0. 5000 2 0. 2500 0. 2500 3 0. 1250 0. 1250 4 0. 0625