用于二维域的快速、完全矢量化的Simpson方法版本。该代码避免使用任何for循环等,可在给定精度级别下比dblquad快一个数量级或更快。代码示例提供了如何使用。功能: ans = simp2D('func',xs,xe,ys,ye,NX,NY)输入参数: func - 接受向量输入的二维函数(否则可能导致错误结果) xs, xe - x积分极限 ys, ye - y积分极限 NX - x方向的积分区间数(应为偶数) NY - y方向的积分区间数(应为偶数)
2D Simpson方法的完全矢量化积分器基于Simpson方法的二维域快速积分器-matlab开发
相关推荐
使用积分器模拟物体运动
使用 Matlab ODE45 积分器和标准的 Runge-Kutta 4 积分器模拟物体运动。更多详情请查阅博客文章。代码摘要:https://mikescodeprojects.files.wordpress.com/2019/12/matlab_run-2.mp4?=1;代码演练:https://mikescodeprojects.files.wordpress.com/2019/12/matlab_code_walkthrough-1.mp4?=3。
Matlab
13
2024-05-12
IRIDirid_fcoi函数实现复分数阶积分器的新途径 - MATLAB开发
irid_fcoi函数设计用于离散时间或连续时间传递函数的计算,以近似形式表达((wgc/s)^lamtacos(mulog(wgc/s))的连续时间分数复阶积分器)^(-sign(mu))。这里,“s”代表拉普拉斯变换变量,“lamta”为分数复阶的实部,取值范围在0到2之间,“mu”为虚部,取值范围在-1到0之间,表示分数复数阶的特性,“wgc”则是增益交叉频率。
Matlab
7
2024-08-19
保持矢量化优化功能的矢量化版本开发 - MATLAB应用
VHOLD(multiax, onoff)用于设置多轴保持状态。 VHOLD(multiax, onoff)是函数hold的优化版本,利用句柄在矩阵中设置多个轴对象的状态multiax,并根据提供的onoff状态。参数onoff可以是字符串'on'或'off',将所有轴设置为相同的保持状态,或者是单元矩阵,以便各个轴可以设置为不同的状态。请注意,当onoff为单元矩阵时,矩阵multiax和单元矩阵onoff应具有相同的大小,即size(multiax)应等于size(onoff)。使用示例:VHOLD(多轴,开关)输入multiax =轴对象的句柄矩阵= [ax11,ax12,...,ax1
Matlab
6
2024-09-14
通用IIR数字微分/积分器新型加权算子基础下的连续分数展开
这种创新型通用IIR数字微分/积分器采用了加权算子,结合了梯形(Tustin)规则和后向差分(Euler)规则,实现了连续分数展开。该技术由Ivo Petras在他的著作《分数阶非线性系统:建模、分析和仿真》中详细描述。
Matlab
5
2024-07-16
MATLAB图像拼接代码-2D匹配二维匹配
MATLAB影像拼接代码图像马赛克和拼接-Yiren Lu (luyiren [at] seas [dot] upenn [dot] edu)图像拼接和拼接的MATLAB实现:哈里斯角检测器见corner_detector_impl.m哈里斯、克里斯和迈克·斯蒂芬斯。“组合角和边缘检测器。”阿尔维视觉会议。卷。15. 1988年。自适应非极大值抑制(ANMS)见anms.m布朗、马修、理查德·塞利斯基和西蒙·温德。“使用多尺度面向补丁的多图像匹配。”2005年IEEE计算机协会计算机视觉和模式识别会议(CVPR'05)。卷。1. IEEE,2005。几何模糊见geo_blur.m Berg、A
Matlab
9
2024-07-15
基于矢量化的压缩感知OMP算法
OMP算法的基本思想是从字典矩阵D(也称为过完备原子库)中选择与信号y最匹配的原子(即某列),构建稀疏逼近。然后将剩余残差减去所有已选择的原子组成的矩阵在空间上的正交投影,得到下一步的信号残差。随后,继续选择与信号残差最匹配的原子,反复迭代。信号y可以由这些原子的线性和,加上最后的残差值来表示。如果残差值在可以忽略的范围内,则信号y即为这些原子的线性组合。OMP分解过程实际上是依次对所选原子进行Schmidt正交化,然后将待分解信号减去在正交化后的原子上的各自分量,即可得到残差。
Matlab
5
2024-07-14
无网格指数积分器的实验实现基于紧支持径向基函数的瞬态偏微分方程数值解
针对瞬态偏微分方程的数值解问题,提出了一种无网格指数积分器的实验实现方法。该方法特别适用于微分方程解在计算域的局部区域随时间变化的情况。空间离散化阶段采用了具有紧支持的径向基函数。时间积分使用了exprb32指数Rosenbrock方法。所需的矩阵函数通过Leja点结合牛顿插值计算得出。这种积分器在空间和时间上表现出完全的自适应性。
Matlab
6
2024-09-01
完全加载的2D绘图仪定制数据集或函数的二维像素图生成工具 - MATLAB开发
此工具可根据输入类型绘制热图:1. 对于@FunctionHandle图示例: heatPlot([-100 100],[-250 250],@(x,y)(sqrt(x.*y)),'TITLE_','$$\sqrt{\chi\gamma}$$' ,'X_LABEL','\chi','Y_LABEL','\gamma','xTickRes',20,'yTickRes',25);2. 对于矩阵图示例: heatPlot([-100 100],[-250 250],randi(1e2,[2e2,5e2]),'TITLE_','$$Z$$','X_LABEL','x','Y_LABEL','y','x
Matlab
6
2024-08-12
Matlab开发解析地形积分方程求解器
在Matlab开发中,设计了一个用于解析沃尔泰拉积分方程的地形积分方程求解器。
Matlab
7
2024-09-30