本帖代码和教程有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代码

歌华链路由器专用Breed，代替不死Uboot，该版本支持固件稳定不掉无线信号，而且比M1版本更完善，不存在除电源灯以外其它四个LED微亮情况，出于商业目的该版本H大已删除。

星云稳定器4100参数，调乱参数的朋友肯定需要的，直接在软件重置就可以了

设计出了一种用于光强检测的前置放大及量程自动转换电路。
许多光强信号放大电路仅追求高增益，忽略了对测量范围的考虑。
本文采用同轴尾纤型光电探测器把光强信号转换成光电流信号，精密截波稳定型运算放大器ICL7652把光电流信号转化为电压信号，量程转换电路74HC4052受单片机控制可在4个量程之间自动转换，通过调节暗电流补偿电路减小光电二极管暗电流所产生的影响。
仿真测试结果表明，电路参数选择合理、电路模块性能稳定，并且很好地降低了噪声的影响，设计的电路具有低噪声、高增益、高共模抑制比、失调小等优点，探测光强动态范围可达76dB。

一款TCP转串口的转发工具，本人亲自开发，提学习和交流。
程序将来自TCP服务端发送过来的socket数据转到串口发送出去，架起了TCP与串口之前的通讯，程序测试稳定可靠。
使用vs2017开发。

用matlab实现扩展有限元计算，有限元法（FEM，FiniteElementMethod）是一种为求得偏微分方程边值问题近似解的数值技术。
它通过变分方法，使得误差函数达到最小值并产生稳定解。
类比于连接多段微小直线逼近圆的思想，有限元法包含了一切可能的方法，这些方法将许多被称为有限元的小区域上的简单方程联系起来，并用其去估计更大区域上的复杂方程。
它将求解域看成是由许多称为有限元的小的互连子域组成，对每一单元假定一个合适的(较简单的）近似解，然后推导求解这个域总的满足条件(如结构的平衡条件），从而得到问题的解。
这个解不是准确解，而是近似解，因为实际问题被较简单的问题所代替。
由于大多数实际问题难以得到准确解，而有限元不仅计算精度高，而且能适应各种复杂形状，因而成为行之有效的工程分析手段

链家网和贝壳网房价爬虫，采集北京上海广州深圳等21个中国主要城市的房价数据（小区，二手房，出租房，新房），稳定可靠快速！支持csv,MySQL,MongoDB,Excel,json存储，支持Python2和3，图表展示数据，注释丰富

在长流程浮选过程中,生产指标难以在线检测,造成操作不及时,影响系统的稳定运行.本文提出了一种基于多源数据的铝土矿浮选过程生产指标集成建模方法.首先结合浮选机理和现场工人经验,分析影响和反映生产指标的多源数据(生产数据和泡沫图像特征数据);然后分别建立各生产指标预测子模型和同步误差补偿子模型;最后采用信息熵和智能协调策略分别构建精矿品位和尾矿品位的集成预测模型.工业验证和工况分析表明,本文集成建模方法具有良好的预测性能和较强的泛化性,为基于生产指标的浮选过程操作参数控制和全流程优化奠定基础.

研究了一种对数螺线柱面晶体配接针孔对Z箍缩铝等离子体进行单色谱成像的摄谱成像仪，摄谱仪具有结构简单、外形尺寸紧凑的特点。
由于对数螺线晶体的保角特性，摄谱成像仪可在较大视场范围内对Z箍缩内爆等离子体进行单色谱成像。
在“阳”加速器上，针对Z箍缩铝等离子体K壳层的X射线辐射进行了成像实验，得到了铝丝阵内爆等离子体的类氢(1727.7eV)和类氦线(1588.3eV)单色图像。
在箍缩单色图像上观察到了磁瑞利泰勒不稳定性引起的“热点”及内爆不稳定性造成的螺旋形结构，反映了等离子体的内爆形态，为进一步理解Z箍缩物理过程和确定等离子体的辐射特性提供了参考。

具有语音播报的倒车雷达本系统以AT89S51控制器为核心、40KHz方波信号升压电路、超声波发射换能器、超声波接收换能器、信号放大与整形电路、液晶显示模块、键盘控制电路、语音播报电路、测温电路。
由本系统构建的超声波测距仪具有测量准确，显示便捷，操作灵活，反应迅速，使用方便，系统工作稳定，耗电量低等许多特点。
一、倒车雷达的主要功能概述　●实时显示测量距离；
●实时显示当前环境温度；
●距离1.5m开始语音播报测量距离,尔后每变化0.2m均报告之●当距离小于1.5m蜂鸣器报警●当距离小于0.5m输出紧急停止的灯光警示●具有开机音乐及语音提示●倒车语音注意提示,0.3m时则紧急告示●具有语音播报使能控制；
●使用4节AA干电池供电；

在日常工作中，钉钉打卡成了我生活中不可或缺的一部分。然而，有时候这个看似简单的任务却给我带来了不少烦恼。 每天早晚，我总是得牢记打开钉钉应用，点击"工作台"，再找到"考勤打卡"进行签到。有时候因为工作忙碌，会忘记打卡，导致考勤异常，影响当月的工作评价。而且，由于我使用的是苹果手机，有时候系统更新后，钉钉的某些功能会出现异常，使得打卡变得更加麻烦。 另外，我的家人使用的是安卓手机，他们也经常抱怨钉钉打卡的繁琐。尤其是对于那些不太熟悉手机操作的长辈来说，每次打卡都是一次挑战。他们总是担心自己会操作失误，导致打卡失败。 为了解决这些烦恼，我开始思考是否可以通过编写一个全自动化脚本来实现钉钉打卡。经过一段时间的摸索和学习，我终于成功编写出了一个适用于苹果和安卓系统的钉钉打卡脚本。