秦九韶算法matlab实验报告,数值分析上机实验报告.doc

实验报告一

题目: (绪论) 非线性方程求解及误差估计

摘要:非线性方程的解析解通常很难给出,因此线性方程的数值解法就尤为重要。本实验采用两种常见的求解方法二分法、Newton法和改进的Newton法。可以节省计算机的计算时间,还能减小不必要的误差。

前言:(目的和意义)

掌握二分法与Newton法的基本原理、应用以及熟练掌握用MATLAB求函数积分

数学原理:

函数的调用格式:quadl(filename,a,b,tol,trace)

其中filename是调用函数名,a和b分别为定积分的下限和上限。用来控制积分精度。

(2)秦九韶算法:

Sn=an

Sk=xSk+1+ak (k=n-1,n-2,...,0),

Pn(x)=S0

程序设计:

计算积分

利用MATLAB,下面给出主程序

>>g=inline('x.^10.*exp(x-1)'); %定义一个语句函数g(x)=exp(x^10*exp(x-1))

I=quadl(g,0,1)

I =

0.0098

例1.9 秦九韶算法

a0=3,

ak=2ak-1+3,

Pn(x)=anx^n+an-1x^(n-1)+...+a1x+a0

求I1=P100(0.5),I2=P150(13)

>>x=input('x=');

n=input('n=');

a=3;

for i=1:n

a=2*a+3;

end

s=z;

b=(a-3)/2;

for m=1:100

s=x*s+b;

b=(b-3)/2;

end

disp(s);

>>x=0.5

n=100

600.0000

>>x=3

n=100

4.7039e+078

结果分析和讨论:

结论:

对于二分法,只要能够保证在给定的区间内有根,使能够收敛的,当时收敛的速度和给定的区间有关,二且总体上来说速度比较慢。Newton法,收敛速度要比二分法快,但是最终其收敛的结果与初值的选取有关,初值不同,收敛的结果也可能不一样,也就是结果可能不时预期需要得结果。改进的Newton法求解重根问题时,如果初值不当,可能会不收敛,这一点非常重要,当然初值合适,相同情况下其速度要比Newton法快得多。

实验报告二

题目: (插值法)用各类插值方法法求解相关数值分析问题

摘要: 非线性方程的解析解通常很难给出,因此线性方程的数值解法就尤为重要。本实验采用两种常见的求解方法二分法和Newton法及改进的Newton法。

前言:(目的和意义)

掌握插值法的基本原理及用插值法求解相关数值分析问题的方法以及熟练掌握用MATLAB求函数积分。

数学原理: MATLAB提供了一、三、N维数据插值函数以及三次样条插值函数等。数据插值问题为一维插值,其采用的方法有线性方法、最近方法、三次样条和三次插值,在MATLAB中,实现这些插值的函数为interp1,其调用格式为inpert1(a,b,c,method);函数根据a、b的值,计算函数在c处的值,method为插值方法,'linear' 为线性插值,'cubic'为三次多项式插值,'spline' 为三次样条插值。

程序设计:

3.给出f(x)=lnx,用一、二、三次线形插值计算`ln0.54的近似值

x0.40.50.60.70.8-0.916291-0.693147-0.510826-0.357765-0.223144

>>x=0.4:0.1:0.8;

f=[-0.916291,-0.693147,-0.510826,-0.357765,-0.223144];

I1=interpl(x,f,0.54) %用一次线形插值计算f(x)

I1 =

-0.620218600000000

I3=interp1(x,f,0.54,'spline') %用3次样条插值计算f(x)

I3 =

-0.615977770000000

21.设f(x)=1/(1+x^2),在[-5,5]上取n=10,安等距节点求分段线性插值函数Ib(x),计算节点间中点处的 与f(x)的值,并估计误差。

>> x=linspace(-5,5,10);

>> y=1\(1+(x.^2));

>> p=polyfit(x,y,1)

p =

0.0000 11.1852

4.给定数据如下表,试求三次样条插值及Hermite插值S(x),并满足条

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/373250.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

Flex 布局教程:语法篇

网页布局(layout)是CSS的一个重点应用。 布局的传统解决方案,基于盒状模型,依赖 display属性 position属性 float属性。它对于那些特殊布局非常不方便,比如,垂直居中就不容易实现。 2009年,W3…

练习错误

form:阻止表单提交的方法一:在form标签中给出以下代码: 1 onsubmit "return False" 方法二:设置事件阻止 1 e.preventDefault() js中判断:只要非数字都应该表示为字符串 1 if(Email.indexOf("") -1){ 2 …

JavaFX 2中的PopupMenu

创建弹出菜单 要在JavaFX中创建Popupmenu,可以使用ContextMenu类。 您向其中添加MenuItems,也可以使用SeparatorMenuItem创建可视分隔符。 在下面的示例中,我选择子类ContextMenu并将MenuItems添加到其构造函数中。 public class Animatio…

matlab中CH指标聚类评价指标,MATLAB聚类有效性评价指标(外部)

MATLAB聚类有效性评价指标(外部)作者:凯鲁嘎吉 - 博客园 http://www.cnblogs.com/kailugaji/更多内容,请看标签:MATLAB、聚类前提:数据的真实标签已知!1. 归一化互信息(Normalized Mutual information)定义程序functio…

学习进度表

周数 专业学习目标 专业学习时/每分钟 新增代码量 知识技能总结 第六周 ps的图像处理 80 30 看书加以实践 第七周 数据结构的链式结构 100 50 多做习题加以巩固知识 第八周 网页设计 80 30 多多练习,学会用代码设计 第九周 图片美工 70 30 慢慢学会运用软…

Axis通过wsdd部署Web Service

axis网上的教程很多,不过搜来搜去,总是只有那么几篇。仔细看了一下那几篇文章,都感觉到不是自己想要的,所以自己整理了一篇分享一下。 本文介绍axis应用的一个小例子,没有麻烦的命令行操作,只需照下面的步骤…

弹簧特性

1.概述 本教程将展示如何通过XML或Java配置在Spring中设置和使用属性 。 在Spring 3.1之前 ,将新的属性文件添加到Spring并使用属性值并不像它那样灵活和健壮。 从Spring 3.1开始 ,新的Environment和PropertySource抽象大大简化了此过程。 2.通过XML名…

php-cgi cpu很高,php-cgi占用cpu资源过高的解决方法

转的网上的,不过对PHP-CGI菜鸟的人,还是有点帮助的。1. 一些php的扩展与php版本兼容存在问题,实践证明 eAccelerater与某些php版本兼容存在问题,具体表现时启动php-cgi进程后,运行10多分钟,奇慢无比&#x…

《做中学》读后有感

《做中学》读后有感 最近读了娄老师的“做中学”系列文章,有很大感触,今天想着重谈一谈我在学习方面收到的启发。 如何成功get一项技能 如果问到“如何开始get一项技能”,我想我们应该是最有发言权的一代。从小就被爸爸妈妈引导着参加各种课外…

多表之间关联查询

内连接 jion on 自连接 本表进行内连接的查询形式 外链接: 左链接 写法:select 字段 from 表1 t left join 表2 s on t.字段1 s.字段1 where 条件 或者 作用:保证左边的表的数据全部显示,包括空的 右链接 写法 :sele…

php文件夹0777,PHP代码mkdir(‘images’,’0777′)创建一个具有411权限的文件夹!为什么?...

我发誓这是昨天的工作.然而,现在下面的代码破坏文件夹没有问题,但创建一个具有411权限的新文件夹应该是777.我的代码昨天这样做.这样做的目的是压缩文件夹,传递文件夹,删除图像,然后为图像创建新目录.有人能告诉我我做错了什么或我应该做什么?谢谢function delete_d…

调查HashDoS问题

近一个月前,我就如何在不与供应商互动的情况下临时解决 28C3上出现的HashDoS问题或其他代码缺陷发表了一些想法。 现在是时候更深入地研究复杂性攻击并查看来源了。 我完全假设java.util.HashMap和java.util.Hashtable是受此攻击影响的最常用的Java数据结构&#xf…

Linq 和 EF Contains示例

List<int> unitIDListnew List<int>(); //此处添加int元素 var query DB.ElecConsumers.Where(c > unitIDList.Contains(c.ParentUnitID)); //EF方式 var query1 (from c in DB.ElecConsumers where unitIDList.Contains(c.ParentUnitID ) select c); //Linq方…

date 显示或设置系统时间和日期

显示或设置系统时间和日期 date [options] [format] date [options] [new date] date用来显示系统的时间和日期&#xff0c;超级用户可以使用date来更改系统时钟 选项 %H 小时&#xff0c;24小时制&#xff08;00~23&#xff09; %I 小时&#xff0c;12小时制&#xff…

Java 7:WatchService

在Java 7的所有新功能中&#xff0c;更有趣的是WatchService&#xff0c;它增加了监视目录更改的功能。 WatchService直接映射到本机文件事件通知机制&#xff08;如果有&#xff09;。 如果本机事件通知机制不可用&#xff0c;则默认实现将使用轮询。 结果&#xff0c;响应性&…

做一件事情的3个关键指标:兴趣、能力和回报

最近突然有了一点新的感悟&#xff0c;在原有的认识基础之上。关于找工作&#xff0c;大家说的最多的&#xff0c;根据自己的“兴趣”和“能力”。我觉得这是不够的&#xff0c;还应该加上一个“回报”。兴趣&#xff1a;对一件事有没有愿望去尝试&#xff0c;侧重“好奇心”。…

iOS应用内支付(IAP)详解

在iOS开发中如果涉及到虚拟物品的购买&#xff0c;就需要使用IAP服务&#xff0c;我们今天来看看如何实现。 在实现代码之前我们先做一些准备工作&#xff0c;一步步来看。 1、IAP流程 IAP流程分为两种&#xff0c;一种是直接使用Apple的服务器进行购买和验证&#xff0c;另一种…

vagrant box php,vagrant box php开发环境配置 -- nginx

centos7.3 直接用yum安装nginx的版本是1.10.2&#xff0c;当前的最新稳定版是1.10.3&#xff0c;暂时不更新&#xff0c;直接安装yum安装nginxsudo yum install -y nginx测试nginx -t启动sudo service nginx startps -ef|grep nginxcurl -i localhost在virtualbox设置网络的端口…

使用ASM 4处理Java类文件–第二部分:Tree API

什么是ASM树API&#xff1a; ASM树API是ASM的一部分&#xff0c;可让您创建/修改内存中的类。 该类被视为信息树。 像整个类一样&#xff0c;它是ClassNode的实例&#xff0c;其中包含FieldNode对象列表&#xff0c;MethodNode对象列表等。本文假设读者已经在这里阅读了第一部分…

php 去除 html 属性,用PHP 去掉所有html标签里的部分属性

用PHP 去掉所有html标签里的部分属性http://zhidao.baidu.com/question/418471924.html用PHP 去掉所有html标签里的部分属性 tppabsset_time_limit(0);function view_dir($dir){$dpopendir($dir); //打开目录句柄//echo "".$dir."";$path2;while ($file r…