原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
第一题
问题叙述：
构造算法并编程序精确计算二次方程的根
●设a≠0,b2
●包括b2≈b
问题分析：
对于一般情况，可以直接使用求根公式x=-b
clear
clear
a=input(a=);b=input(b=);c=input(c=);
x1=-b/2*a+1/2*a*sqrt(b*b-4*a*c)
x2=-b/2*a+1/2*a*sqrt(b*b-4*a*c)
但是当b2?4ac时，b2≈
解决方法：
（1）方法1：使用双精度计算
（2）方法2：变换求根公式
当b0时，
x
当b0时，
x
这样可以很好地避免出现减性抵销现象
实例分析：
以a=c=1,b=±1000000.000001为例
（1）方法1：使用双精度计算（与单精度作对比）
%calculateusingdoubleprecision
%calculateusingdoubleprecision
%forthepurposesofcomparison,alsousesingletocalculate
clear
a=input(a=);b=input(b=);c=input(c=);
%a,b,cisthecoefficientofequationax^2+bx+c=0
delta=b*b-4*a*c;
%first,usingdouble
d1=(-b+sqrt(delta))/2/a;
d2=(-b-sqrt(delta))/2/a;
d1=vpa(d1,12),d2=vpa(d2,12)
%than,usingsingle
delta=single(delta);
a=single(a);b=single(b);c=single(c);
s1=single((-b+sqrt(delta))/2/a);
s2=single((-b-sqrt(delta))/2/a);
s1=vpa(s1,12),s2=vpa(s2,12)
对于不同的系数，计算结果如下：
方程系数
单精度
双精度
a
b
c
x
x
x
x
1
7
5
-0.807417631
-6.192582130
-0.807417596
-6.192582403
1
10000
1
0.0
-10000.0
-0-9999.9999
1
1000000.000001
1
0.0
-1000000.0
1.0000076e-6
-1000000.0
1
-1000000.000001
1
0.0
1000000.0
1.0000076e-6
1000000.0
结论：
显然，在一般情况下，单精度还能保持一定准确度；但是当a=1，b=10000，c=1时，减性抵销现象很明显，单精度下的x1计算结果为0，在b=1000000.000001
（2）方法2：变换求根公式
当b0时，
x
当b0时，
x
%calculateusingdouble
%calculateusingdouble
%forthepurposesofcomparison,alsousesingletocalculate
clear
a=input(a=);b=input(b=);c=input(c=);
%a,b,cisthecoefficientofequationax^2+bx+c=0
delta=b*b-4*a*c;
%first,usingdouble
ifb0
d1=-2*c/(b+sqrt(delta));%usingnewformula
d2=(-b-sqrt(delta))/2/a;
else
d1=-2*c/(b-sqrt(delta));%usingnewformula
d2=(-b+sqrt(delta))/2/a;
end
d1=vpa(d1,12),d2=vpa(d2,12)
%than,usingsingle
delta=single(delta);
a=single(a);b=single(b);c=single(c);
ifb0
s1=single(-2*c/(b+sqrt(delta)));%usingnewformula
s2=single((-b-sqrt(delta))/2/a);
else
s1=single(-2*c/(b-sqrt(delta)));%usingnewformula
s2=single((-b+sqrt(delta))/2/a);
end
s1=vpa(s1,12),s2=vpa(s2,12)
计算结果如下：
方程系数
单精度
双精度
文档评论（0）