题目:用牛顿法求方程x-cos(x)=0的实根(精确到1E-6)。

(1)要求用函数调用。

(2)进一步研究和弦截法作比较。

算法分析:

(1)       此题是利用牛顿方法解一元非线性方程的根。(牛顿法是把非线性方程局部线性化的一种方法,它在单根附近具有较高的收敛速度。)所以首先我们应先给出估计的根,先对方程x-cos(x)=0变形,令y1=x,y2=cos(x),则两函数图象的交点,就是方程x-cos(x)=0的根,这里利用Matlab作图估计根的值。

在Matlab命令行中输入,并运行:

>>x=-2:0.01:2;

>>y1=cos(x);

>>y2=x;

>>plot(x,y1,x,y2);

>>grid on;

可得下图:

从图中可以很容易得到根x的初值可选0.6。 (2)利用牛顿法的迭代公式x1=x0-(f(x0)/f’(x0));这里首先取x0=0.6,代入迭代公式,然后判断x1与x0之差的绝对值是否小于精度,如果小于精度,则停止,即得出根x的值,如果不小于精度,则继续迭代,直到符合精度为止。

Matlab代码如下:

(1)funNewt

matlab 一元方程程序,用牛顿方法解一元非线性方程的根(Matlab实现)相关推荐

  1. 问题四十:对ray tracing圆环图形进行debug(2)——C++,用“笛卡尔”方法解一元四次方程

    第七步:用"笛卡尔"方法解一元四次方程 "笛卡尔"方法也就是"待定系数法" 通过"式子八"可以解出s(可能有0到3个实根 ...

  2. 非线性方程组求解Matlab实现 (多元牛顿方法、Broyden方法、Broyden方法2)

    三元牛顿方法(非线性方程中三个未知数) function jie = multvarnewton3(g1,g2,g3,c) e=[inf,inf,inf]; syms x y z k1=g1(x,y, ...

  3. 二分法求非线性方程组Java_用C#编写二分法解一元非线性方程

    在计算方法里有一种求一元非线性方程的解法,叫做二分法. 简单介绍如下: 函数F(x)在区间[a,b] 上连续,假定在区间内有唯一实根,计为x* 二分法的思想为:首先确定有根区间,将区间二等分,通过判断 ...

  4. MATLAB用SOR方法解方程组,SOR迭代法 求MATLAB程序

    满意答案 vwbog 2014.11.02 采纳率:41%    等级:6 已帮助:1507人 function [n,x]=sor22(A,b,X,nm,w,ww) %用超松弛迭代法求解方程组Ax= ...

  5. 二元牛顿迭代法matlab,牛顿迭代法解二元方程组以及误差分析 matlab实现

    电子科大数值分析第一次作业 类似一元函数的泰勒展开: f(x) f(x0) f'(x0)(x x0) 二元函数的展开为: f(x,y) f(x0,y0) [(x x0) (y y0)]f(x,y)x ...

  6. MATLAB递归程序的调试方法

    MATLAB具有和VC类似的调试功能,各功能的按钮图标和VC也很相似,基本可以使用VC的调试思想. MATLAB自己的特点: 一. 查看中间变量. 1. 在调试状态下,鼠标放在变量上就可以显示变量的基 ...

  7. 牛顿迭代法解一元三次方程

    题目描述 有形如:ax^3+bx^2+cx^1+dx^0=0ax3+bx2+cx1+dx0=0 这样的一个一元三次方程.给出该方程中各项的系数(a,b,c,da,b,c,d均为实数),并约定该方程存在 ...

  8. matlab 纽马克 激励,用Newmark方法计算系统的动力学响应的matlab程序

    请大家帮忙看看这个程序有什么问题?用Newmark方法计算系统的动力学响应,结果大的惊人. function[Q,V,AA]=newmarkb E=2.1e11;P=7850;D1=0.405;d1= ...

  9. python迭代法求解方程_第一部分:趣味算法入门;第六题牛顿迭代法求一元三次方程的根...

    100个不同类型的python语言趣味编程题 在求解的过程中培养编程兴趣,拓展编程思维,提高编程能力. 第一部分:趣味算法入门:第六题SRE实战 互联网时代守护先锋,助力企业售后服务体系运筹帷幄!一键 ...

最新文章

  1. java实现qq_java实现的类似qq聊天系统
  2. Windows2003如何安装IIS
  3. 在Visual Studio 2010 中使用菱形向导对窗口进行布局
  4. OpenStack绝非昨日黄花 | 人物志
  5. c# 操作word中在右下角插入图片
  6. 浏览器升级怎么升_测评:iOS 13 升不?新老机型跑分对比!
  7. 【路径规划】基于matlab蚁群算法无人机巡检路径规划【含Matlab源码 138期】
  8. c语言学习宝典怎么样,C语言学习宝典
  9. Fiddler(FD)抓包工具汉化版及使用方法
  10. (附源码)PHP疫情上报管理系统 毕业设计 170948
  11. QQ VS MSN——QQ “只爱陌生人”,MSN “不和陌生人说话”
  12. spoj2142 Arranging Flowers
  13. python分析txt文件_python解析文本文件
  14. RabbitMQ第二话 -- Springboot基于四种Exchange(Direct、Fanout、Topic、Heders、延时队列)的实现和多虚拟主机下的生产消费者实现
  15. 面向物流行业的文档管理系统
  16. RAM、ROM、emmc、iNand、SD卡、mmc 与 Nandflash 的区别
  17. 计算机网络第七版(谢希仁) 第一章 概述 1-10,1-17作业答案
  18. 《网络传播概论》第二章 网络的属性与传播形态
  19. PyQt5在QGridLayout布局下动态删除增加控件
  20. Linux进程管理概述

热门文章

  1. 付出与回报——伍迷“我的2007”随想 (转载 注 这个文章再次告诉我,做人要厚道,只要付出了,总有...
  2. [Codeforces 940E]Cashback
  3. FastDFS,Redis,Solr,ActiveMQ核心技术整合二(1)
  4. AS 项目目录组成架构
  5. torch.zeros_like() 和 torch.zeros()的区别
  6. 加载组件Already included file name ‘e:/Vuejs/组件/单文件组件/vue_test/src/components/Message.vue‘ differs from
  7. nc文件分解为tiff文件
  8. 一日不刷算法题,明日编码临大敌
  9. linux ntp攻击python脚本,(转)如何发现 NTP 放大攻击漏洞
  10. 3个视频提取字幕的方法,简单又高效,建议收藏