- 1、本文档共27页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
数值分析实验报告
第一题 实验题1.2
1、实验内容
实验1.2 体会稳定性在选择算法中的地位,误差扩张的算法不稳定,而误差衰竭的算法是稳定的。分别采用E.1.6(即E.1.4)E1=1e,En=1-nEn-1,n=2,
2、源程序
%function t_charpt1
%数值试验1.2:误差传播与算法稳定性
%输入:递推公式选择与递推步数
%输出:各步递推值及误差结果,以及递推值和误差与递推步数的关系图
clear;
clc;
promps = {请选择递推关系式,若选(1.4),请输人1,否则输入2:};
result = inputdlg(promps,charpt 1_2,1,{1});
Nb = str2num(char(result));
if((Nb~= 1)(Nb~= 2))
errordlg(请选择递推关系式,若选(1.4),请输人1,否则输人2!);
return;
end
result = inputdlg({请输人递推步数 n:},charpt 1_2,1,{ 10});
steps = str2num(char(result));
if(steps1)
errordlg(递推步数错误!);
return;
end
result = inputdlg({请输入计算中所采用的有效数字位数:},charpt 1_2,1,{5});
Sd = str2num(char(result));
format long %设置显示精度
result=zeros(1,steps) ; %存储计算结果
err=result; %存储计算的绝对误差
func=result; %存储用库函数quadl计算出的积分的近似值
%用库函数quadl计算积分的近似值
for n= 1:steps
fun=@(x) x.^n.* exp(x-1);
func(n) = quadl(fun,0,1);
end
if(Nb==1)
%用算法(1.4)计算
digits(Sd); %控制有效数字位数
result(1) = subs(vpa(1/exp(1)));
for n=2:1:steps
result(n)=subs(vpa(1-n * result(n-1)));
end
err=abs(result-func);
elseif(Nb==2)
%用算法(1.5)计算
digits(Sd); %控制有效数字位数
result(steps)=0;
for n=steps:-1:2
result(n-1) = subs(vpa((1-result(n))/n));
end
err=abs(result-func);
end
clf;%清除当前图像窗口
disp(递推值:);
disp(sprintf(%e , result));
disp(误差:);
disp(sprintf( %e ,err));
plot([1:steps],result,-,LineWidth,2);
set(gca,linewidth,0.5,fontsize,16);
grid on
hold on;
plot([1:steps],err,r--,LineWidth,2);
xlabel(steps n,FontSize,18);
ylabel(En-and ERR n--,FontSize,18);
legend(En,err(n));
title([Algorithm (1., num2str(Nb+3),) Significant Digits , num2str(Sd)],FontSize,18);
% text(2,err(2),\uparrow err(n));
% text(4,result(4),\downarrow En);
3、实验结果
(1)算法E1.6,有效数字5位
递推值:
3.678800e-01 2.642400e-01 2.072800e-01 1.708800e-01 1.456000e-01 1.264000e-01 1.152000e-01 7.840000e-02 2.944000e-01 -1.944000e+00
误差:
5.588280e-07 1.117662e-06 3.352927e-06 1.341222e-05 6.705713e-05 4.023702e-04 2.816427e-03 2.253226e-02 2.027877e-01 2.02
(2)算法E1.6,有效数字6位
递推值:
您可能关注的文档
- 人教版新教材地理1 海水的运动.pptx
- 《卫生监督学》模考试题和答案三.doc
- 2020中考(冀)道德与法治总复习练习题:小题狂练4 认识国情,爱我中华.doc
- 人教版必修2 地理:2.3 城市化 .pptx
- 外研版九年级上册 module 4综合能力测试.doc
- 2019年10月22日 文言文阅读之分析重点虚词的意义和用法——《每日一题》2019-2020年八年级语文人教(上).doc
- 四年级上册第一单元1.7数的产生与十进制记数法.ppt
- 湖南省攸县三中2020届高三上学期第二次月考数学(理)试卷和答案.doc
- 北师大版九年级数学上册 第二章 一元二次方程 单元评估测试试题(无答案).docx
- 小学生法制教育发言稿.docx
文档评论(0)