牛顿科茨matlab,牛顿-科茨(Newton-Cotes)公式算法.doc
牛顿-科茨(Newton-Cotes)公式算法.doc
牛顿-科茨(Newton-Cotes)公式算法:
Step 1:判断type类型,1转Step 2,2转Step 3,3转Step 4;否则输出值为0;
Step 2:计算科茨公式;
Step 3:计算牛顿-科茨六点公式;
Step 4:计算牛顿-科茨七点公式;
流程图:
开始type=2type=1计算科茨六点公式计算科茨公式type=3计算科茨七点公式输入type
开始
type=2
type=1
计算科茨六点公式
计算科茨公式
type=3
计算科茨七点公式
输入type
N N N
输出结果为0输出结果输出结果输出结果
输出结果为0
输出结果
输出结果
输出结果
定义function函数如下:
function I = NewtonCotes(f,a,b,type)
%type = 1 科茨公式
%type = 2 牛顿-科茨六点公式
%type = 3 牛顿-科茨七点公式
I=0;
switch type
case 1,
I=((b-a)/90)*(7*subs(sym(f),findsym(sym(f)),a)+...
32*subs(sym(f),findsym(sym(f)),(3*a+b)/4)+...
12*subs(sym(f),findsym(sym(f)),(a+b)/2)+...
32*subs(sym(f),findsym(sym(f)),(a+3*b)/4)+7*subs(sym(f),findsym(sym(f)),b));
case 2,
I=((b-a)/288)*(19*subs(sym(f),findsym(sym(f)),a)+...
75*subs(sym(f),findsym(sym(f)),(4*a+b)/5)+...
50*subs(sym(f),findsym(sym(f)),(3*a+2*b)/5)+...
50*subs(sym(f),findsym(sym(f)),(2*a+3*b)/5)+...
75*subs(sym(f),findsym(sym(f)),(a+4*b)/5)+19*subs(sym(f),findsym(sym(f)),b));
case 3,
I=((b-a)/840)*(41*subs(sym(f),findsym(sym(f)),a)+...
216*subs(sym(f),findsym(sym(f)),(5*a+b)/6)+...
27*subs(sym(f),findsym(sym(f)),(2*a+b)/3)+...
272*subs(sym(f),findsym(sym(f)),(a+b)/2)+...
27*subs(sym(f),findsym(sym(f)),(a+2*b)/3)+...
216*subs(sym(f),findsym(sym(f)),(a+5*b)/6)+41*subs(sym(f),findsym(sym(f)),b));
end
运算的输入方式及结果为:
>> y=NewtonCotes('sin(x)',0,10,1)
y =
3.7613
>> y=NewtonCotes('sin(x)',0,10,2)
y =
2.7865
>> y=NewtonCotes('sin(x)',0,10,3)
y =
1.5296
2、编程解决以下科学计算问题。
1.编程解决以下科学计算问题。
求编写求解方程函数文件,然后调用该函数文件求30000个点上的x,y,最后在所有的(x、y)坐标处标记一个点(不要连线)绘出图形。这种图形又称埃农(Henon)引力线图,他将迭代迟来的随机点吸引到一起,最后得出貌似连贯的引力线图。
流程图:
开始
开始
创建初始矩阵,并赋初值(0,0),固定点n=30000
创建初始矩阵,并赋初值(0,0),固定点n=30000
将第K行的值通过迭代关系赋值给第(K+1)行变量为k,k<3000
将第K行的值通过迭代关系赋值给第(K+1)行
变量为k,k<3000
N
牛顿科茨matlab,牛顿-科茨(Newton-Cotes)公式算法.doc相关推荐
- 牛顿迭代法例题 matlab,牛顿迭代法-matlab程序(解线性方程组)
牛顿迭代法 matlab程序(解线性方程组) 作者:佚名来源:转载发布时间:2009-3-7 16:55:53 减小字体 增大字体 1.功能 本程序采用牛顿法,求实系数高次代数方程 f(x)=a0xn ...
- 牛顿迭代法例题 matlab,牛顿迭代法求方程根的MATLAB程序
function [x_reality,n_reality] = Newt( f_name,x_start,tolerance,n_limit) %% %牛顿迭代法(切线法)求解方程f_name = ...
- 力科示波器 matlab,力科周末文章四周年180期合集目录
力科周末文章四周年180期合集 一.示波器基础系列 01 关于示波器的带宽-Frankie 01-1calibrating scopes 02 采样率和存储深度-Rick 03 DSO中的内插技术-P ...
- 用matlab画牛顿迭代图形,matlab牛顿迭代动画演示
% 牛顿环演示的MATLAB程序 close all;clear;clc; figure('Position',[90 164 873 483]); L=632.8;R=5;H=5; a1=axes( ...
- 力科示波器 matlab,力科示波器手册.pdf
LECROY LECROY X-STREAM™™ X-STREAM 系列示波器 系列示波器 用户手册 用户手册 2 00 3 年 12 月 2 00 3 年 12 月 LeCroy Corporati ...
- 双边滤波去雾matlab,基于双边滤波的去雾算法.doc
基于双边滤波的去雾算法 中图法分类号:TP391.41 文献标识码:A 文章编号:1006-8961(2013) - - 论文引用格式: 基于双边滤波的图像去雾 王一帆1 ,尹传历2,黄义明1 , 王 ...
- 模式识别感知器算法matlab,模式识别第三章-感知器算法.doc
模式识别第三章-感知器算法 模式识别第三章 感知器算法 一.用感知器算法求下列模式分类的解向量w: 将属于的训练样本乘以(-1),并写成增广向量的形式: ,,, ,,, 第一轮迭代:取, 因不大于0, ...
- 利用matlab实现复数域空间牛顿迭代法的分形图案展示(newton法)
利用matlab实现复数域空间牛顿迭代法的分形图案展示(newton法) 1 一维函数的牛顿迭代法 2 复平面的牛顿迭代法 2.1 简单方程结果 2.2 其它非线性方程结果 本文首发于 matlab爱 ...
- matlab差商chaS,牛顿插值MATLAB算法.doc
MATLAB程序设计期中作业 --编程实现牛顿插值 成员:刘 川(P091712797) 签名_____ 汤 意(P091712817) 签名_____ 王功贺(P091712799) 签名_____ ...
- matlab中牛顿下山法实例,非线性方程的数值解法牛顿下山法matlab
非线性方程的数值解法牛顿下山法matlab 1 非线性方程的数值解法 --计算物理实验作业九 陈万 物理学2013级 130******** ● 题目: 用下列方法求0133=--=x x f(x)在 ...
最新文章
- plsql视图添加表字段_教你不到两分钟,快速创建数据透视图
- ESL:我们如何使用首云混合云产品实现提效降本
- DevExpress中获取GridControl排序之后的List
- cpout引脚是干什么的_FPGA中差分信号的定义和使用(一)
- CodeVS 3027 线段覆盖2(DP)
- 【word使用技巧】删掉某一行参考
- 2020统计局的行政划分表_天津市第七次全国人口普查区域划分与地图绘制试点工作在西青区开展...
- [转载] Python数学实验与建模 课后习题第1章解析
- Android 通过StickyBroadcast获取电量
- Windows10 phpStudyV8.1的基本使用
- microsoft visual studio 2010 旗舰版的产品密
- 网页flash加载不出来
- 动手实验:使用 jstat 摸清线上系统的JVM运行状况
- 南京邮电大学java第三次实验报告
- java ico图片转png_Java 图片处理: ico 格式转 PNG/JPG 等格式
- 如何在vue项目中使用styl
- 中国大学慕课python答案 北京理工大学_中国大学慕课mooc_Python 语言程序设计_章节测验答案...
- 安科瑞无线测温方案part1
- element rules不生效
- excel2021 打印圆不圆