- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
第一章
一、题目
设,其精确值为。
1) 编制按从大到小的顺序,计算SN的通用程序。
2) 编制按从小到大的顺序,计算SN的通用程序。
3) 按两种顺序分别计算,并指出有效位数。(编制程序时用单精度)
4) 通过本次上机题,你明白了什么?
二、通用程序
N=input(Please Input an N (N1):);
AccurateValue=single((0-1/(N+1)-1/N+3/2)/2);
Sn1=single(0);
for a=2:N;
Sn1=Sn1+1/(a^2-1);
end
Sn2=single(0);
for a=2:N;
Sn2=Sn2+1/((N-a+2)^2-1);
end
fprintf(The value of Sn (N=%d)\n,N);
fprintf(Accurate Calculation %f\n,AccurateValue);
fprintf(Caculate from large to small %f\n,Sn1);
fprintf(Caculate from small to large %f\n,Sn2);
disp(____________________________________________________)
三、结果
从结果可以看出有效位数是6位。
感想:可以得出,算法对误差的传播有一定的影响,在计算时选一种好的算法可以使结果更为精确。从以上的结果可以看到从大到小的顺序导致大数吃小数的现象,容易产生较大的误差,求和运算从小数到大数所得到的结果才比较准确。
第二章
一、题目
(1)给定初值及容许误差,编制Newton法解方程f(x)=0的通用程序。
(2)给定方程,易知其有三个根
a) 由Newton方法的局部收敛性可知存在当时,Newton迭代序列收敛于根x2*。试确定尽可能大的。
b)试取若干初始值,观察当时Newton序列的收敛性以及收敛于哪一个根。
(3)通过本上机题,你明白了什么?
二、通用程序
1、定义函数和导函数
%% 定义函数f(x)
function Fx=fx(x)
Fx=x^3/3-x;
----------------------------------
%% 定义导函数df(x)
function Fx=dfx(x)
Fx=x^2-1;
2、寻找最大的
clear
flag=1;
k=1;
x0=0;
while flag==1
delta=k*10^-6;
x0=delta;
k=k+1;
m=0;
flag1=1;
while flag1==1 m=10^3
x1=x0-fx(x0)/dfx(x0);
if abs(x1-x0)10^-6 flag1=0;
end
m=m+1;
x0=x1;
end
if flag1==1||abs(x0)=10^-6 flag=0;
end
end
fprintf(The maximun delta is %f\n,delta);
3、Newton法求方程的根
clear
ef=10^-6; %% 给定容许误差10^-6
k=0;
x0=input(Please input initial value Xo:);
disp(k Xk);
fprintf(0 %f\n,x0);
flag=1;
while flag==1 k=10^3
x1=x0-fx(x0)/dfx(x0);
if abs(x1-x0)ef
flag=0;
end
k=k+1;
x0=x1;
fprintf(%d %f\n,k,x0);
end
三、结果
寻找最大的δ值:
在题目给出的不同区间内进行牛顿迭代:
综上所述:首先得到最大的δ值为0.774597。在各区间任意取初值进行迭代得到(-∞,-1)区间收敛于-1.73205,(-1,-δ)区间局部收敛于0, (-δ,δ)区间收敛于0,(δ,1)区间类似于(-1,δ)区间,收敛于0,(1,∞)收敛于1.73205。
感想:通过本上机题,了解到用Newton法求多根方程的根时,迭代序列收敛于某一个根有一定的区间限制。如果不清楚这个限制随意取值的话,会出现在一个区间上局部收敛于不同的根的情况。如下面的迭代:
第三章
一、
文档评论(0)