- 1、本文档共18页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
数值实验题new.doc
大学的 数值实验题1
实验1.1 病态问题
实验目的:
算法有“优”与“劣”之分,问题也有“好”与“坏”之别。对数值方法的研究而言,所谓坏问题是指问题本身对扰动敏感,反之属于好问题。本实验通过对一个高次多项式方程的求解,初步认识病态问题。
实验内容:
考虑一个高次的代数多项式
(E.1.1)
显然该多项式的全部根为1,2,…,20,共计20个,且每个根都是单重的(也称为简单的)。现考虑该多项式的一个扰动 , (E.1.2)
其中,ε是一个非常小的数。这相当于是对方程(E.1.1)中x19的系数作一个小的扰动。比较方程(E.1.1)和方程(E.1.2)根的差别,从而分析方程(E.1.1)的解对扰动的敏感性。
实验步骤与结果分析:
实验源程序
function t_charpt1_1
% 数值实验1.1病态问题
% 输入:[0 20]之间的扰动项及小的扰动常数
% 输出:加扰动后得到的全部根
clc
result=inputdlg({请输入扰动项:在[0 20]之间的整数:},charpt 1_1,1,{19});
Numb=str2num(char(result));
if((Numb20)|(Numb0))errordlg(请输入正确的扰动项:[0 20]之间的整数!);return;end
result=inputdlg({请输入(0 1)之间的扰动常数:},charpt 1_1,1,{0.00001});
ess=str2num(char(result));
ve=zeros(1,21);
ve(21-Numb)=ess;
root=roots(poly(1:20)+ve);x0=real(root); y0=imag(root);
plot(x0,y0, *);
disp([对扰动项 ,num2str(Numb),加扰动,num2str(ess),得到的全部根为:]);
disp(num2str(root));
对于x19项的扰动ess,不同的取值对应的结果如下所示。
对扰动项 19加扰动1e-010得到的全部根为:
19.9961,19.0257,17.9085,17.1508,15.7982,15.181,13.8995,13.0571,11.9753,11.0109,9.99608,9.00111,7.99978,7.00003,6,5,4,3,2,1。
对扰动项 19加扰动1e-009得到的全部根为:
19.952,19.2293,17.6573+0.692896i,17.6573-0.692896i,15.4524+0.875524i,15.4524-0.875524i,13.3527+0.486992i,13.3527-0.486992i,11.8578,11.0427,9.9916,9.00201,7.99952,7.00009,5.99999,5,4,3,2,1。
对扰动项 19加扰动1e-007得到的全部根为:
20.422+0.999203i,20.422-0.999203i,18.1572+2.4702i,18.1572-2.4702i,15.3149+2.69865i,15.3149-2.69865i,12.8466+2.06246i,12.8466-2.06246i,10.9216+1.10366i,10.9216-1.10366i,9.56629,9.11508,7.99387,7.00027,6,5,4,3,2,1。
对扰动项 19加扰动1e-005得到的全部根为:
22.5961+2.3083i,22.5961-2.3083i,18.8972+5.00563i,18.8972-5.00563i,14.9123+4.95848i,14.9123-4.95848i,12.0289+3.73551i,12.0289-3.73551i,10.059+2.33021i,10.059-2.33021i,8.63828+1.0564i,8.63828-1.0564i,7.70896,7.028,5.99942,5.00001,4,3,2,1。
根在复平面上的位置如图所示:
图 ess=1e-010 图 ess=1e-009
图 ess= 1e-007 图 ess=1e-005
从实验的图形中可以看出,当ess充分小时,方程E.1.1和方程E.1.2的解相差很小,当ess逐渐增大时,方程的解就出现了病态解,这些解都呈现复共轭性质。并且,病态解首先出现在x=16这个解附近,如ess=1e-009时,x=20,19,12,11,…,2,1的解基本误差不大。在x
您可能关注的文档
- 实用英语作文模板2.doc
- 《计算机应用基础 》课程教学大纲new.doc
- 如何编写牛顿盘项目可行性研究报告方案(可用于发改委立项及银行贷款+2013详细案例范文).docx
- 开题报告49982.doc
- 滨州医学院大学生科技创新活动基金项目申请书1.doc
- 中国城镇化水平和速度的实证分析与前景预测new.doc
- 《数值计算方法》试验报告册new.doc
- 数字化定量分析new.doc
- 张坤鹏 -外文翻译.doc
- 上海体育学.doc
- 四川省德阳市罗江中学2025届高三考前热身化学试卷含解析.doc
- 山东省枣庄现代实验学校2025届高三下学期第五次调研考试化学试题含解析.doc
- 吉林省长春市十一高中等九校教育联盟2025届高三一诊考试生物试卷含解析.doc
- 2025届江苏省盐城市伍佑中学高考仿真模拟化学试卷含解析.doc
- 2025届广西贺州中学高考冲刺押题(最后一卷)生物试卷含解析.doc
- 安徽省池州市贵池区2025届高三第一次模拟考试生物试卷含解析.doc
- 宁夏银川一中2025届高三(最后冲刺)化学试卷含解析.doc
- 广东省广州市增城区四校联考2025届高考压轴卷化学试卷含解析.doc
- 2025届邯郸市第一中学高考生物必刷试卷含解析.doc
- 2025届安徽省安庆市石化第一中学高考仿真卷化学试卷含解析.doc
文档评论(0)