本帖代码和教程有Matlab技术论坛原创,原帖参见http://www.matlabsky.com/viewthread.php?tid=3885一、数值积分基本公式数值求积基本通用公式如下Eqn1.gif(1.63KB)2009-11-2023:23xk:求积节点Ak:求积系数,与f(x)无关数值积分要做的就是确定上式中的节点xk和系数Ak。
可以证明当求积系数Ak全为正时,上述数值积分计算过程是稳定。
二、插值型数值积分公式对f(x)给定的n+1个节点进行Lagrange多项式插值,故Eqn2.gif(2.95KB)2009-11-2023:23即求积系数为Eqn3.gif(3.29KB)2009-11-2023:23三、牛顿-柯特斯数值积分公式当求积节点在[a,b]等间距分布时,插值型积分公式(先使用Lagrange对节点进行多项式插值,再计算求积系数,最后求积分值)称为Newton-Cotes积分公式。
由于Newton-Cotes积分是通过Lagrange多项式插值变化而来的,我们都知道高次多项式插值会出现Runge振荡现象,因此会导致高阶Newton-Cotes公式不稳定。
Newton-Cotes积分公式的求积系数为Eqn4.gif(3.38KB)2009-11-2023:28其中C(k,n)称为柯特斯系数。
(1)当n=1时,Newton-Cotes公式即为梯形公式Eqn5.gif(1.68KB)2009-11-2023:28容易证明上式具有一次代数精度(对于Newton-Cotes积分公式,n为奇数时有n次迭代精度,n为偶数时具有n+1次精度,精度越高积分越精确,同时计算量也越大)(2)当n=2时,Newton-Cotes公式即为辛普森(Simpson)公式或者抛物线公式Eqn6.gif(2.04KB)2009-11-2023:28上式具有3次迭代精度(3)当n=4时,Newton-Cotes公式称为科特斯(Cotes)公式Eqn7.gif(2.68KB)2009-11-2023:28上式具有5次迭代精度。
由于n=3和n=2时具有相同的迭代精度,但是n=2时计算量小,故n=3的Newton-Cotes积分公式用的很少(4)当≥8时,通过计算可以知道,在n=8时柯特斯系数出现负值由于数值积分稳定的条件是求积系数Ak必须为正,所以n>=8以上高阶Newton-Cotes公式,我们不能保证积分的稳定性(其根本原因是,Newton-Cotes公式是由Lagrange插值多项推导出来的,而高阶多项式会出现Rung现象)。
四、复化求解公式n阶Newton-Cotes公式只能有n+1个积分节点,但是高阶Newton-Cotes公式由不稳定。
为了提高大区间的数值积分精度,我们采用了分段积分的方法,即先将原区间划分成若干小区间,然后对每一个小区间使用Newton-Cotes积分公式,这就是复化Newton-Cotes求积公式。
(1)当n=1时,称为复化梯形公式。
将[a,b]等分为n份,子区间长度为h=(b-a)/n,则复化梯形公式为(注意:复化求解公式不需要求积子区间等间距,只是Newton-Cotes公式分段积分时自动对小区间进行等分,我们这里采用等分子区间是为了便于计算而已)Eqn8.gif(2.18KB)2009-11-2023:28(2)当n=2时,称为复化辛普森公式。
Eqn9.gif(2.96KB)2009-11-2023:28五、Newton-Cotes数值积分公式Matlab代码
2023/11/26 8:36:30 126KB Matlab 技术论坛 牛顿 科特斯
1
安川机器人梯形图教程,介绍的很详细。
只是现在DX200以后默认以语句形式显示了
2023/11/24 5:42:58 2.53MB 工业机器人
1
1、本程序运用C语言,根据三菱PLC_FX2N的通信协议和通信命令,基于主控芯片STM32F103XX(目前在STM32F103RC,STM32F103RD,STM32F103VC,STM32F103VD,STM32F103VE测试通过)上编写运行的程序,可以直接利用三菱编程软件编写梯形图下载运行,无需任何转换。
目前至少支持的指令有:(其他指令亲可以自己添加)RSTRSTSRSTTCOUTOUTSSETSETSADDSUBMULDIVLDLDILDPLDFANDANIORORIANDPANDFORPORFADDPSUBPMULPDIVPMOVMOVPENDFENDINCDECINCPDECPCJCALLRETINVLD=LD>LDAND<AND=编程语言梯形图程序容量8K步内部寄存器D8000个定时器T 256个记数器C256个输入点X256个输出点Y256个壮态继电器S600个辅助继电器M3071点M0-M3071特殊功能:M8000(运行监视触点)M8001(运行监视反触点).M8002(初始化脉冲触点)M8003(初始化脉冲反触点)M8004(错误指示触点)M8011(10毫秒时钟脉冲)M8012(100毫秒时钟脉冲)M8013(1秒时钟脉冲)M8014(1分时钟脉冲)M8020(零位标志)M8021(借位标志)M8022(进位标志)M8029(指令执行结束标志)M8033(内存保持触点)M8034(禁止输出触点).更多参考FX2N系列。
2023/11/10 0:46:10 6.92MB PLC源码
1
PLC水塔水位控制系统设计,包括梯形图,指令表,流程图,接线图,I/O地址分配,软硬件设计,组态仿真
1
此版本是机床厂新代系统编辑plc专用,新增加联机在线编辑功能,离线模拟功能,免安装,支持xp,win7及上版本。
2023/10/31 17:35:45 2.81MB LadEditor
1
西门子plc指令与梯形图速查手册
2023/10/27 6:10:44 6.65MB 西门子plc S7-200 PLC指令
1
西门子s7200变频一拖三恒压供水的工艺图包含梯形图分配图电路图
2023/10/6 2:24:27 146KB S7-200 变频一拖三
1
基于STM32F103ZET6的两路DAC模块输出两路相位相差120度的正弦波,方波,三角波以及梯形波,还有不对称的正弦波,叠加了高次谐波的正弦波等,利用按键切换,
2023/9/16 18:42:04 3.59MB STM32F103 正弦波 方波 三角波
1
里面有就是多种数据处理工具,非常便利。
1.文本文件操作1.1征地部标准坐标导出...21.2征地部标准坐标导入...61.3线封闭...61.4点集转面...72.    数据转换...82.1       SHP转数据库...82.2       批量数据库转数据库...92.3       栅格彩色转黑白...93.    数据检查和数据信息获得...103.1       锐角检查...103.2       获得内角并判断是否凸多边形...113.3       获得线(面)两个折点方向...123.4       四至和范围获得...143.4.1获得数据的XY范围...143.4.2获得数据的经纬度范围...153.4.3四邻信息获得...153.4.4地块四至点获得...173.4.5地块四至点坐标获得...183.4.6获得相对四至(适合大比例小地块)193.4.7获得绝对四至(根据四至点坐标)...203.5       属性赋值...223.5.1比例分析...223.5.2加权平均...233.6       椭球面积计算...233.7       计算点到线的距离...243.8       道路河流依次经过的地方...244.    裁剪和合并...254.1       按属性裁剪...254.2       矢量数据批量裁剪...254.3       矢量数据批量合库...264.4       影像批量裁剪...274.5       影像批量合并...285.    MXD文档处理...295.1       MXD批量裁剪...295.2       MXD批量导出图片...305.3       mxd压缩和版本另存...315.4       mxd文档相对路径和无效数据检查...316.    制图...326.1       梯形接幅表的创建...326.2       矩形接幅表的创建...336.3       公里网或方里网制作...346.4       经纬网制作...356.5       色带制作...357.    拓扑错误处理...367.1       删除线面直线上的点...367.2       点不在线面上处理...377.3       线部分或完全重叠处理...387.4       删除完全重复的点线面处理...387.5       面线边界不重合...397.6       面缝隙处理...397.7       面重叠处理...407.8       删除线面上重复点...407.9       检查多部件要素...417.10         删除伪节点...428.    编号工具...428.1       整库更新BSM..428.2       更新BSM..438.3       字符串前补零...438.4       按图形自动编号...449.    数据处理...459.1       两个图层按重叠度赋属性...459.2       分区域消除...469.3       批量压缩数据库...479.4       批量修复几何(修复前一定备份数据)479.5       按长度分割线...489.6       线分割面保留属性...489.7       融合时字段连接...499.8       要素移动...5110.           业务相关...5110.1         上级行政区和下级行政区图形不一致处理...5110.2         修改面左上角点为第一个点...5210.3         修改面左上角点为第一个点根据点层...5310.4
2023/9/6 6:03:24 1.66MB 真好
1
很实用的资源,可以把电梯编程主要思想很简单、容易理解,有梯形图
2023/8/8 6:05:56 440KB 电梯 实例 毕业 论文
1
共 90 条记录 首页 上一页 下一页 尾页
在日常工作中,钉钉打卡成了我生活中不可或缺的一部分。然而,有时候这个看似简单的任务却给我带来了不少烦恼。 每天早晚,我总是得牢记打开钉钉应用,点击"工作台",再找到"考勤打卡"进行签到。有时候因为工作忙碌,会忘记打卡,导致考勤异常,影响当月的工作评价。而且,由于我使用的是苹果手机,有时候系统更新后,钉钉的某些功能会出现异常,使得打卡变得更加麻烦。 另外,我的家人使用的是安卓手机,他们也经常抱怨钉钉打卡的繁琐。尤其是对于那些不太熟悉手机操作的长辈来说,每次打卡都是一次挑战。他们总是担心自己会操作失误,导致打卡失败。 为了解决这些烦恼,我开始思考是否可以通过编写一个全自动化脚本来实现钉钉打卡。经过一段时间的摸索和学习,我终于成功编写出了一个适用于苹果和安卓系统的钉钉打卡脚本。
2024-04-09 15:03 15KB 钉钉 钉钉打卡