logo

您所在位置网站首页 > 海量文档  > 教育文化 > 教育文化

数值计算方法上机实习题.docx 13页

本文档一共被下载: ,您可全文免费在线阅读后下载本文档。

  • 支付并下载
  • 收藏该文档
  • 百度一下本文档
  • 修改文档简介
全屏预览

下载提示

1.本站不保证该用户上传的文档完整性,不预览、不比对内容而直接下载产生的反悔问题本站不予受理。
2.该文档所得收入(下载+内容+预览三)归上传者、原创者。
3.登录后可充值,立即自动返金币,充值渠道很便利
下载地址
需要金币:300    文档纠错   收藏文档   下载帮助
支付并下载
特别说明: 下载前务必先预览,自己验证一下是不是你要下载的文档。
  • 上传作者 高江辉(上传创作收益人)
  • 发布时间:2017-04-21
  • 需要金币300(10金币=人民币1元)
  • 浏览人气
  • 下载次数
  • 收藏次数
  • 文件大小:620.38 KB
下载过该文档的会员
你可能关注的文档:
数值计算方法上机实习题设,由递推公式,从I0=0.1824, 出发,计算;,, 用,计算;(1)由I0计算I20递推子程序:function f=fib(n,i)if n>=1 f=fib(n-1,i)*(-5)+(1/(n));elseif n==0 f=i;end计算和显示程序:I=0.1824;I1=0.1823;fib1=fib(20,I);fib2=fib(20,I1);fprintf('I_0=0.1824时,I_20=%d\n',fib1);fprintf('I_0=0.1823时,I_20=%d\n',fib2);计算结果显示:I_0=0.1824时,I_20=7.480927e+09I_0=0.1823时,I_20=-2.055816e+09(2)由I20计算I0程序:n=21;i1=0;i2=10000;f1=i1;f2=i2;while n~=0; f1=f1*(-1/5)+(1/(5*n)); f2=f2*(-1/5)+(1/(5*n)); n=n-1;endfprintf('I_20=0 时,I_0=%4.8f\n',f1);fprintf('I_20=10000时,I_0=%4.8f\n',f2);计算结果显示:I_20=0 时,I_0=0.18232156I_20=10000时,I_0=0.18232156分析结果的可靠性及产生此现象的原因(重点分析原因)。答:第一个算法可得出 e0=I0-I0*en=In-In*=5ne0易知第一个算法每一步计算都把误差放大了5倍,n次计算后更是放大了5n倍,可靠性低。 第二个算法可得出en=In-In*e0=15nen可以看出第二个算法每一步计算就把误差缩小5倍,n次后缩小了5n倍,可靠性高。求方程的近似根,要求,并比较计算量。在[0,1]上用二分法;[0,1]上的二分法二分法子程序:function [root,n]=EFF3(f,x1,x2)%第二题(1)二分法f1=subs(f,symvar(f),x1);%函数在x=x1的值f2=subs(f,symvar(f),x2);%x=x2n=0;%步数er=5*10^-4;%误差if(f1==0) root=x1; return;elseif(f2==0) root=x2; return;elseif(f1*f2>0) disp('两端点函数值乘积大于0!'); return;else while(abs(x1-x2)>er)%循环 x3=(x1+x2)/2; f3=subs(f,symvar(f),x3); n=n+1; if(f3==0) root=x3; break; elseif(f1*f3>0) x1=x3; else x2=x3; end end root=(x1+x2)/2;%while循环少一步需加上end计算根与步数程序:fplot(@(x) exp(x)+10*x-2,[0,1]);grid on;syms x;f=exp(x)+10*x-2;[root,n]=EFF3(f,0,1);fprintf('root=%6.8f ,n=%d \n',root,n);计算结果显示:root=0.09057617 ,n=11 取初值,并用迭代;初值x0=0迭代迭代法子程序:function [root,n]=DDF(g,x0,err,max) (接下页)%root根,n+1步???,g函数,x0初值,err误差,max最大迭代次数X(1)=x0;for n=2:max X(n)=subs(g,symvar(g),X(n-1)); c=abs(X(n)-X(n-1)); root=X(n); if(c<err) break; end if n==max disp('超出预设迭代次数'); endendn=n-1;%步数减1计算根与步数程序:syms x;f=(2-exp(x))/10; (接下页)x0=0; err=5*10^(-4);max=100;[root,n]=DDF(f,x0,err,max);fprintf('root=%6.8f ,n=%d \n',root,n);计算结果显示:root=0.09051262 ,n=4加速迭代的结果;加速迭代加速迭代计算程序:x0=0;err=5*

请自觉遵守互联网相关的政策法规,严禁发布色情、暴力、反动的言论。
用户名: 验证码: 点击我更换图片

©2010-2013 max.book118.com在线文档投稿赚钱网. All Rights Reserved 蜀ICP备08101938号