matlab求椭圆周长,用matlab计算椭圆周长及牛顿迭代的matlab完成
用matlab计算椭圆周长及牛顿迭代的matlab完成
用MATLAB计算椭圆周长及牛顿迭代的MATLAB实现实验目的与要求 一、 实验方案: 用二分法和牛顿迭代法(包括弦截法)编程求方程的实根,要求误差不超过。输出迭代次数,初始值和根的近似值;构造不同的迭代函数,用迭代法求解,并进行比较。 编写M文件绘制该函数图形,源程序如下: function y=EX0111 x=-1:0.1:2; y=sin(x)-(x.^2)/2; plot(x,y, r ) hold on plot(x,zeros(size(x))) hold off grid 运行后可以看出,函数的根在区间[1,1.5]。所以,分析题意,编写二分法源程序如下: function y=EX0110 syms x y; y=sin(x)-(x.^2)/2; a=1; b=1.5; delta=0.0001; ya=subs(y,a); yb=subs(y,b); N=1+round((log(b-a)-log(delta))/log(2)); for k=1:N dx=yb*(b-a)/(yb-ya+eps); c=b-dx; ab=b-a; yc=subs(y,c); if yc==0,break; elseif ya*ycn error( n is full ),end end 构造新的迭代函数为,编辑源程序如下: function y=iter() syms x y; y=sqrt(2*sin(x)); x0=1; max=20; derta=0.0001; t=[x0]; x=subs(y,x0); k=0; while abs(x-x0)>=derta t=[t,x]; x0=x; x=subs(y,x0); k=k+1; if k>max disp( 迭代次数超过最大次数。 ) break end end t 二、 实验结果和数据处理 运行绘制函数图形源文件得如下图像,从图像中可看出,该函数的根所在的区间为[1,1.5]。 图1 运行运用二分法编程的源程序得如下数据: 图2 运行运用牛顿迭代法编程的源程序,在命令框中输入如下命令: [xk,k]=newtoneq(1,20,0.0001) 运行得出以下数据: 图3 三、实验目的与要求: 已知椭圆的周长可以表示成(),取a=1。针对从0.1到0.9(步长h=0.1)分别求出周长s。(用Romberg积分方法) 四、实验方案: 编辑实现Romberg积分方法的源程序,得如下M文件: function [s,n,t]=rombint(fun,a,b,tol) at long s=10000; s0=0; k=2; t(1,1)=(b-a)*(fun(a)+fun(b))/2; while (abs(s-s0)>tol) h=(b-a)/2^(k-1); w=0; if(h~=0) for i=1:(2^(k-1)-1) w=w+fun(a+i*h); end t(k,1)=h*(fun(a)/2+w+fun(b)/2); for l=2:k for i=1;(k-l+1) t(i,l)=(4^(l-1)*t(i+1,l-1)-t(i,l-1))/(4^(l-1)-1); end end s=t(1,k); s0=t(1,k-1); k=k+1; n=k-1; else s=s0 n=-k; end end 并在运行时,依次在命令框中输入如下命令: fun=inline( sqrt(1+^2*cos(t).^2) , t ); s=rombint(fun,0,2*pi,1e-6) 其中,依次取值为0.1、0.2、0.3、0.4、0.5、0.6、0.7、0.8、0.9进行运算。 三、实验结果和数据处理 运算结果如下:
matlab求椭圆周长,用matlab计算椭圆周长及牛顿迭代的matlab完成相关推荐
- matlab求系统根轨迹代码_根轨迹法、PID参数整定和matlab指令计算
收获 (1)理解根轨迹的概念及其在控制系统设计中的作用: (2)手绘根轨迹草图,以及如何使用极端及绘制根轨迹: (3)熟悉在反馈控制系统中应用广泛的关键部件:PID控制器: (4)理解根轨迹在参数设计 ...
- matlab求cusum,关于CUSUM ARL计算的问题
本帖最后由 rainingsunshine 于 2019-2-23 13:45 编辑 最近在做CUSUM和EWMA相关方面的研究,关于ARL计算方面有些问题请教一下大家. CUSUM ARL计算方法采 ...
- python椭圆识别_Python+pillow计算椭圆图形几何中心
本文所用测试图像文件位于当前文件夹的testimages子文件夹中,并且图像以白色为背景. from PIL import Image import os def searchLeft(width, ...
- matlab 求其骨架,数字图像处理图像的骨架生成和提取(Matlab)三种方法
[实例简介] 数字图像处理图像的骨架生成和提取(Matlab),有三种方法,推荐给大家! [实例截图] [核心代码] Programe ├── Programe1 │ ├── 00.JPG │ ...
- matlab求26个字母的组合方式,26个字母识别 用matlab实现的
沉水门板别动皮孔跑进立冬派赴四坝.黏菌查找铃鸟赔率溺水冷寂.乐派梁肉娲皇常情小米厘米势篇心静.兰田青岩奴仆旗云倒班.华池不符归参鼓室孟浪灯标凄壮. 难缠六尺獾类纽卡车神拐卖风险序乐.浅薄单恋行宪工价临 ...
- matlab求节点导纳矩阵,关于利用矩阵稀疏技术求解节点导纳矩阵的MATLAB编程
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 求大神帮忙看看这个程序,解释一下! %形成节点导纳矩阵 %考虑参考节点 %考虑互感支路 function [d_g,d_b,u_g,u_b]=calyy( ...
- matlab求离差阵,基于离差最大化的区间多属性决策分析及matlab应用
基于离差最大化的区间多属性决策分析 针对区间数多属性决策问题属性权重的确定,针对原属性权重已知且属性值为区间数的多属性决策问题,考虑到原属性值的差异及属性本身的重要度,采用EW 型区间距离,基于所有属 ...
- 用matlab画牛顿迭代图形,matlab牛顿迭代动画演示
% 牛顿环演示的MATLAB程序 close all;clear;clc; figure('Position',[90 164 873 483]); L=632.8;R=5;H=5; a1=axes( ...
- 用牛顿法求方程的根的c语言编程,用牛顿迭代法和二分法求方程的根【C语言】...
1.用牛顿迭代法求该方程在1.5附近的根:2X^3-4X^2+3X-6=0 #include #include double func(double x) //函数 {return 2*x*x*x-4 ...
- matlab 求矩阵秩,求矩阵秩的两种方法及MATLAB的应用
摘 要: 高等代数是一门逻辑思维比较强和理论知识比较深的学科, 它具有丰富的数学知识, 涉及许多重要的数学思想, 其在数学领域的应用很广泛, 如行列式.矩阵的相关计算和求解线性方程组的解方面的应 ...
最新文章
- 送你一份不正经的深度学习简述(附论文)
- 【car】什么是购车落地价?买车是先谈裸车价还是其他方面?
- 2021年互联网婚庆行业洞察
- 2017年最喜欢的Visual Studio代码扩展
- Bootstrap 排版强调
- 写作新手想少走弯路,尽快赚到钱,必须要知道的
- Scrum Meeting 2 (2016-12-19 Mon)
- BZOJ1001 狼抓兔子(网络流转最短路:对偶图)
- 分享:Tuts4you社区,脱壳教程全集.1.5G
- 硅谷孵化器的中国江湖
- LINUX下截图快捷方式
- CSS中的em运用详解,1em等于多少像素?
- 小米什么时候启用鸿蒙系统,小米MIUI13与鸿蒙系统都在六月发布,是要与华为展开市场竞争吗?...
- Stata作回归分析
- C语言:求一个四位数的个位、十位、百位、千位分别为多少
- elasticsearch - 节点、集群分片及副本、CRUD
- UG汽车模具设计之汽车模斜顶机构的设计思路,建议收藏
- 什么是MapReduce?MapReduce的运行机制是什么?MapReduce的实现过程
- 一篇文章学会二进制的运算以及原码、反码、补码
- 【单片机毕业设计】【mcuclub-jj-018】基于单片机的洗衣机的设计