这是一个由17个MATLAB函数组成的软件包,用于通过谱配置(伪谱)方法在有界、周期和无限区间上求解微分方程。附加的函数包括纳入边界条件,使用重心公式执行插值,和计算正交多项式的根。另外的程序演示了如何使用包解决特征值,边值和初值问题产生在特殊函数,量子力学,非线性波和水动力稳定性领域。
详情见发表于
JAC Weideman, SC Reddy, MATLAB微分
矩阵套件,ACM数学软件汇刊,
26卷,465-519页(2000)。
网站:
http://dip.sun.ac.za/~weideman/research/differ.html
在MATLAB R13下测试。
J.A.C Weideman(2021)。DMSUITE(//www.tatmou.com/matlabcentral/fileexchange/29-dmsuite), MATLAB中央文件交换。检索.
一个称为MaxPol的通用框架最近已经发布并在这里提供
//www.tatmou.com/matlabcentral/fileexchange/63294-maxpol-smoothing-and-differentiation-package
MaxPol提供了一个框架来设计具有如下性质的各种数值微分核:
(1)无旁瓣伪影的截止(低通)设计(用于噪声鲁棒的情况)
(2)任意阶微分
(3)任意多项式精度
(4)导数矩阵设计
(5)带有转向矩的2D导数核
(6)信号和图像处理中的直观例子
这是一个有用的包,主要是在使用ODE求解器进行及时演化时。
祝贺你
这是poldif函数背后的算法。M在数值上不稳定。考虑以下代码
N = 30;
x0 = linspace(0,1, N)';
x1 = x0。^ 2;
x2 = x1。* x0;
%
D0 = poldif(x0, 1);
D1 = poldif(x1, 1);
drawtext (D2 = poldif(x2, 1);
这是节点0 <= x <= 1的三种不同分布。
给定N个点,每个微分矩阵的秩为N - 1,微分矩阵的零空间是常数向量。
计算
%
rD0 = rank(D0)
rD1 = rank(D1)
rD2 = rank(D2)
%
分别排名29,6和3 !首先,这表明算法失败了,即使是适当数量的30个节点。现在如果我们计算D乘以一个常数向量它应该得到0向量。
%
vnD0 = D0 * ones(N, 1);
vnD1 = D1 * ones(N, 1);
vnD2 = D2 * ones(N, 1);
%
它不。这是一个更严重的错误。由于误差集中在支撑的末端。金宝app如果要精确地解决初值和边值问题,我们需要精确。
尽管存在这个问题,但一揽子计划作为一个整体是一个重要的贡献。
直观,良好的文档和非常有用!
优秀的! !
我得到你的网站在我的工作上FORTRAN因为我是学生
很好地放在一起包装欢呼!
非常感谢!!
好工作;
非常感谢。
一个伟大的深思熟虑的图书馆
非常好的包!然而,函数poldif。当节点数(根)增加时,m失效(例如,埃尔米特多项式为500,拉盖尔多项式为130)。这是由于计算条目c(k)/c(j):首先取乘积,然后取除法。产品生产过剩。
简单的和强大的。
对于执行动力系统的输入/输出(瞬态放大)分析非常有用。易于使用,文档齐全
对我的研究不可或缺!
优秀的工具!
我希望是否可以帮助我解决生物技术中的非线性BVP
很好的工具
伟大的工具箱!,我使用非线性Schrödinge方程。
伟大的工具箱,
非常有用的。
伟大的工作! !
理解简单概念非常有用
理论背后&形象化它们。
推荐给任何一个大学生。
太好了! ! !
解决了很多关于分化的问题!!
这东西太棒了!
用谱和伪谱方法解决微分方程的有用和漂亮的代码。
我对这个方法很感兴趣,谢谢
的帮助