本帖代码和教程有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
matlab分段灰度线性变换代码
2023/11/18 14:36:30 543B matlab 分段灰度变换
1
本程序利用MFC实现了B样条曲线曲面的绘制:均匀B样条曲线准均匀B样条曲线分段Bezier曲线非均匀B样条曲线曲面
2023/10/28 16:41:29 1.98MB B样条 B样条曲面 曲线 MFC
1
单片机最小系统,时钟芯片,根据用电的峰时、平时、谷时分段计价,存储芯片,用于存储充值余额,脉冲输入,用于模拟用电量
2023/10/26 10:27:19 858KB 单片机
1
该算法是对STOMP算法的完整实现,使用者可改变图像后直接使用
2023/10/10 2:21:22 4KB MATLAB CS
1
选取windows系统自带的ding.wav信号作为分析对象,在Matlab软件平台下,利用函数wavread对音频信号进行采样,记住采样频率和采样点数,听一下原始声音sound(y,fs,bits)。
(2)音频信号的频谱分析,先画出音频信号的时域波形;
然后对音频号进行快速傅里叶变换fft(y,N),N取32768,画出信号的频谱特性,加深对频谱特性的理解。
(3)根据频谱,反演时域特性,画出时域波形。
寻找幅值最大的两个频率,此频率除以fft点数在乘以采样频率就是信号的主频,即可合成信号的时域图形,听一下声音。
(4)对原音频信号进行1024点的分段付立业分析meshgrid(5)根据主要频线合成音频,并画出时域图形,试听合成效果。
(6)采用线性插值(linspace)和傅立业反变换(fliplr,ifft)分别合成音频,并画出时域图形,试听效果。
2023/10/9 9:49:25 48KB ding 音频信号的频谱分析
1
excel常用函数公式及技巧搜集从身份证号码中提取出生年月日从身份证号码中提取出性别从身份证号码中进行年龄判断按身份证号号码计算至今天年龄以2006年10月31日为基准日,按按身份证计算年龄(周岁)的公式按身份证号分男女年龄段【年龄和工龄计算】根据出生年月计算年龄根据出生年月推算生肖如何求出一个人到某指定日期的周岁?计算距离退休年龄的公式求工齡计算工龄年龄及工龄计算自动算出工龄日期格式为(yyyy.mm.dd)【时间和日期应用】自动显示当前日期公式如何在单元格中自动填入当前日期如何判断某日是否星期天某个日期是星期几什么函数可以显示当前星期求本月天数显示昨天的日期关于取日期如何对日期进行上、中、下旬区分如何获取一个月的最大天数日期格式转换公式【排名及排序筛选】一个具有11项汇总方式的函数SUBTOTAL自动排序按奇偶数排序自动生成序号如何自动标示A栏中的数字大小排序?如何设置自动排序重复数据得到唯一的排位序列按字符数量排序排序字母与数字的混合内容随机排序排序的问题怎样才能让数列自动加数一个排序问题数字的自动排序,插入后不变?根据规律的重复的姓名列产生自动序号排名的函数自动排名公式百分比排名的公式写法为:平均分及总分排名求名次排名排名次根据分数进行普通排名对于普通排名分数相同时,按顺序进行不重复排名依分数比高低名次成绩排名美国式排名中国式排名求最精简的自动排名公式排序后排名位次排名根据双列成绩进行共同排名在双列间排名等次排名不等次排名(行小排先)不等次排名(行大排先)顺次排名有并列排名无并列排名有并列分段排名无并列分段排名成绩排名如何排名数据排名(隔几行排名)根据分数进行倒排名倒数排名函数是什么如何实现每日各车间产量的排名分数相同时按照一科的分数进行排名筛选后自动产生序列号并汇总如何筛选奇数行函数筛选姓名名次筛选如何实现快速定位(筛选出不重复值)如何请在N列中列出A1:L9中每列都存在的数值自动为性别编号的问题【文本与页面设置】EXCEL中如何删除*号将字符串中的星号“*”替换为其它字符去空格函数如何去掉字符和单元格里的空格怎样快速去除表中不同行和列的空格如何禁止输入空格代替单元格中字符串把单元格中的数字转变成为特定的字符格式把有六百多个单元格的一列,变成一页的多列将N列变M列公式归纳为一列变四列四列变一列重复四次填充多行数据排成一列将单元格一列分为多列首写字母大写把单元格编号中的小写字母变成大写字母让姓名左右对齐数字居中而小数点又对齐计算指定单元格编号组中非空单元格的数量比较两个单元格内容是否一致怎么样设置才能让这一列的每个单元格只能输入12位如何让工作表奇数行背景是红色偶数行背景是蓝色计算特定的一组单元格中,满足条件的单元格的个数把文本格式的数字转换成真正的数字设置页码Excel表格里如何插入页码的?如何设置页脚首页为第5页表格的页脚问题无拘无束的页眉打印表头Excel打印中如何不显示错误值符号对于一些不可打印的字符的处理用那个函数可将个位数前面的零值显示出来?如果你要在A3的前面插入100行请问如何每隔30行粘贴一新行在工作表里有连续10行数据,现在要每行间格2行一个大表每一行下面需要加一行空行,怎么加最方便Excel中插入空白行快速删除工作表中的空行快速删除空行一次删完Excel里面多出很多的空白行每30行为一页并加上一个标题如何实现如何实现隔行都加上标题项如何把标签页去掉的?去掉默认的表格线(网线)表格的框线列标的标识变了符号的意义双击格式刷竟也能COPY文本(不是文本格式)查找+格式刷的妙用光标移动最后一行为文本最后一行为数字如何在EXCEL中快速定位最后一行数据用SUN函数快速求和在Excel中快速查看所有工作表公式在Excel中设置行间距怎样同时改变多行行高快速换行让文本换行在Excel中行列快速转换将原有列中的内容倒置过来快速回到A1单元格复制粘贴中回车键的妙用一次选中批注单元格一次在所有单位格中插入批注在公式中插入批注不连续单元格填充同一数据空白行的填充怎样用函数向下实现自动填充怎么设置自动保存避免输入网址和电子邮件地址时的超链接单元格前面自动加了等号加盖公章查找+格式刷的妙用Excel中鼠标双击妙用Excel中快速定位的技巧实例在Excel中插入Flash时钟的步骤小写数字转换成人民币大写轻轻松松制作超复杂Excel表头【字符截取与增减】截取单元格里某个字符后的字符截取字符的公式如何确定*号
2023/10/8 7:32:08 965KB Excel 公式 实例
1
变电站电压/无功控制一直以来都是众多学者和供电部门关注的问题。
目前广泛使用的基于“九区图”电压无功控制策略,因为电压无功控制装置和控制策略的矛盾,在实际运行中暴露出严重的问题。
从根本上改善电压无功调节特性,不仅要从控制策略上着手,还要从电压无功调节装置上着手。
寻找一种满足工程实际需要的变电站电压无功控制方式具有重要的理论意义和工程实用价值。
静止同步补偿器(STATCOM)具有响应时间短、产生谐波含量少,在系统电压下降时,输出无功的能力不受母线电压的影响。
在系统故障或负荷突增时,能够快速的交换无功,动态提供无功支撑,从而较好地改善电压水平,抑制冲击负荷造成的电压波动。
因此本文将STATCOM用于变电站电压无功控制中,改善变电站调压装置与控制策略“九区图”之间存在的不足,提出了一种在线灵敏度计算的基于“离散设备优先动作,连续设备精细调节”原则的电压无功控制策略,协调有载调压变压器、电容器以及STATCOM之间的运行。
一方面,变电站侧的电容器作为主要的无功输出更加接近无功负荷端,对整个电网起到一个基础性的无功支撑作用,而STATCOM的无功出力保持一定裕量的状态,使之有足够的可调无功储备,以应对变电站紧急情况,提高运行的安全性,同时有效的减少变压器分接头和电容器组的动作次数;另一方面,离散控制的有载调压变压器和电容器只能实现阶跃、分段的控制,而且其调节容量是一个相对较大的数值,通过sTATc0M快速的无功交换能力,降低电容器投切动作对系统造成的冲击。
通过电磁暂态仿真软件EMTDC/PSCAD对EPR工36节点系统进行仿真分析,验证本文所提出的电压无功控制策略的可行性和正确性。
2023/9/29 7:05:29 3.95MB 无功补偿
1
一个基于android的秒表计时器,需实现以下功能:能够显示时间并且要精确到毫秒。
能够实现人机的简单交互功能,清零,暂停等功能。
能够实现多个时间的计时功能即分段次数的计数。
能够实现最大显示时间到小时。
2023/9/29 4:02:06 1.42MB android 秒表
1
设计程序sy2-2.c,要求从键盘输入x的值,计算以下分段函数的值y。
2023/9/24 0:36:44 173B c语言
1
共 119 条记录 首页 上一页 下一页 尾页
在日常工作中,钉钉打卡成了我生活中不可或缺的一部分。然而,有时候这个看似简单的任务却给我带来了不少烦恼。 每天早晚,我总是得牢记打开钉钉应用,点击"工作台",再找到"考勤打卡"进行签到。有时候因为工作忙碌,会忘记打卡,导致考勤异常,影响当月的工作评价。而且,由于我使用的是苹果手机,有时候系统更新后,钉钉的某些功能会出现异常,使得打卡变得更加麻烦。 另外,我的家人使用的是安卓手机,他们也经常抱怨钉钉打卡的繁琐。尤其是对于那些不太熟悉手机操作的长辈来说,每次打卡都是一次挑战。他们总是担心自己会操作失误,导致打卡失败。 为了解决这些烦恼,我开始思考是否可以通过编写一个全自动化脚本来实现钉钉打卡。经过一段时间的摸索和学习,我终于成功编写出了一个适用于苹果和安卓系统的钉钉打卡脚本。
2024-04-09 15:03 15KB 钉钉 钉钉打卡