caputo分数阶导数程序_caputo分数阶导数

最新推荐文章于 2022-06-20 09:47:59 发布
最新推荐文章于 2022-06-20 09:47:59 发布
阅读量4.2k 收藏 3
点赞数
文章标签： caputo分数阶导数程序
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_39641463/article/details/112922112
版权
本文探讨了Caputo分数阶导数的理论与应用，包括其在微分方程求解、气候研究和工程实践中的重要性。通过分析，展示了Caputo分数阶导数与Riemann-Liouville导数的关系，并提出了数值求解方法。此外，还讨论了分数阶导数在处理具有记忆效应系统中的独特优势。
摘要由CSDN通过智能技术生成

Caputo 分数阶微分算子合成性质的推广① 秦彤晖, 张笛 【摘要】摘要: 运用微积分学的基本方法,得到了 Riemann-Liouville 分数 阶积分算子和 Caputo 分数阶微分算子......

而且对于分数阶序列导数的微分方程的求解 是比较实用的.本文的目的在于展示微分变换方法不仅适用于Caputo分数阶导数, 而且同样适用于Riemann—Liouville导数. 参考文献[......

定义 2.4 [26]:对 Caputo 分数阶导数的 Sumudu 变换 ? ...

4 现在,基础数学研究和工程应用研究中最 常用的有以下四种分数阶微积分的定义: Grunwald-Letnikov分数阶微积分, RiemannLiouville分数阶微积分,Caputo型分数阶导数和 ......

分数 阶导 数 袁 占斌 , 王俊 刚,聂 玉峰 ,孙浩 ( 西北工业大学理学院...

北京大学大气物理系刘式达教授将分数阶导数描述为 “天气与气候之间的桥 梁” , 气候的分数阶导数是天气,正是由于分数阶导数的存在,使得气候较天气 的记忆性好......

同时, 变导数作为分数阶导数的发展形式已经越来越受到研究人员的重视 . Caputo 形式的变导数定义如下: t 1 q( t) D0 f ( t ) = ( t -τ)-q( t) ......

当α= 2 时 , 由 ( 20) 式和 Caputo 分数阶导数的 定义 : ...

1, y (0) = y , (1) 这里 0c D αt y (t ) 表示 Caputo 分数阶导数,其定义形式为 c α 0 Dt t y ( l ) (τ ) dτ y ( t ) = Γ......

?m?1 ? 一.分数阶微积分定义和数值求解方法 1.1.3 Caputo 定义 进一步证明发现,在t>0时,如果考虑一类函数 ,它具有m+1阶连续的 导数,那么, 分数阶微积分......

该方程如下: 并服从下述边界条件 其中 cD

最低0.47元/天 解锁文章
weixin_39641463
关注
caputo分数阶导数程序_基于Caputo导数的分数阶非线性振动系统响应计算
weixin_30668061的博客
01-14 1010
基于Caputo导数的分数阶非线性振动系统响应计算李亚杰1,2，吴志强*1,2，章国齐1,2【摘要】摘要:研究了含分数阶Caputo导数的非线性振动系统响应的数值计算方法。首先，由Caputo分数阶导数算子的叠加关系，得到含分数阶导数项非线性振动系统状态方程的标准形式。其次，基于Caputo导数与Riemann-Liouville导数和Grunwald-Letnikov导数间的关系，推导计算了Ca...
Caputo分数阶快速算法——基于分数阶导数核函数的指数和逼近的Matlab程序
jstboyfriend的博客
02-20 1275
博主自己敲的分数阶微分方程的有限差分方法（第二版）附录部分的代码，该代码给出了caputo分数阶导数的核函数的指数和逼近， 可以根据公式进一步加工得到给出快速L1逼近，快速L2-1逼近和快速H2N2逼近。
caputo分数阶导数程序_第二章：黏弹性与波传播——2.5.2 微分形式的波动方程：分数阶导数...
weixin_28707365的博客
12-31 1028
Wave Fields in Real Media:Wave Propagation in Anisotropic,Anelastic,Porous and Electromagnetic MediaJose M. Carcione 著出版时间：2015年蜜酒厅通讯社 固体地球物理学部 译数值实现：Matlab 2019a仅供学习和参考，禁止转载封面及文中照片感谢 @CycleUser 友情提供前...
Caputo 分数阶导数的 H2N2 插值逼近 (附Matlab程序)
图灵的猫的博客
06-19 1039
本文介绍 H2N2 插值逼近方法, 结合数值算例验证了收敛情况,并给出了Matlab程序源代码.
Caputo 分数阶导数快速的 H2N2 插值逼近(附Matlab代码)
图灵的猫的博客
06-19 1577
众所周知, 分数阶微分方程由于其存在全局依赖(非局部性)加大了求解困难, 于是有学者给出了一系列的快速算法,本文我们讨论 Caputo 分数阶导数基于 H2N2 逼近的快速算法,并结合数值算例验证了其求解高效性,同时给出了相应的Matlab源代码.......
caputo分数阶导数程序_Caputo分数阶导数的一种二阶逼近方法
weixin_39885683的博客
12-24 826
A second-order accurate numerical approximation for Caputo fractional derivativesHuang Fenghui1黄凤辉,(1974-)，女，副教授，主要研究方向：微分方程数值解1、Department of Mathematics, School of Sciences, South China University o...
latex 分数_【重点】分数阶微积分和分数阶微分方程数值实验（8）——G-L分数阶微积分(2)——短时记忆效应...
weixin_39806818的博客
11-30 748
参考文献薛定宇《分数阶微积分学与分数阶控制》数值实现Matlab 2019a 主要基于薛定宇开发的FOTF工具箱蜜酒厅通讯社 固体地球物理学部前情回顾形式主义的居士：【重点】分数阶微积分和分数阶微分方程数值实验（8）——G-L分数阶微积分（1）​zhuanlan.zhihu.com短时记忆效应定义：一般情况下，如果计算步长 选择过小，或者 的值过大，则 参与求和的点数会特别庞大，可能导致计算量...
matlab实现分数阶
02-28
matlab实现分数阶 非常好的分数阶资源，能够帮助理解，加深印象，能够把分数阶融会贯通。
Caputo分数阶导数的稳定数值逼近 (2008年)
04-22
Caputo导数能够对初始值进行调整，以便在考虑分数阶导数的同时，保持系统的初始条件不变。这种导数在处理具有明显初始状态特征的物理问题时尤其重要。 文章提出的方法，是基于对函数进行傅里叶变换，并采用正则化...
分数阶matlab
04-27
非常好的分数阶资源，能够帮助理解，加深印象，能够把分数阶融会贯通
论文研究-含Caputo型分数阶导数的灰色预测模型.pdf
09-20
论文研究-含Caputo型分数阶导数的灰色预测模型.pdf,  针对缺乏统计规律的小样本预测系统, 如何挖掘其发展规律, 一直是学术界的难点. 本文依据分数阶微积分理论, 将整数阶导数灰色模型推广到分数阶导数灰色模型, 并从是否满足新信息优先原理、初值利用情况、还原误差大小和稳定性等方面说明了新模型的优势, 以期用Caputo型分数阶导数的记忆性描述小样本预测系统. 实例表明含有Caputo型分数阶导数的灰色预测模型的有效性与实用性.
分数阶方程的MATLAB代码
04-29
计算分数阶微分方程非常有效且准确的计算程序，并包含示例使用说明。
分数阶混沌程序
12-14
混沌分数阶仿真程序，matlab进行编程仿真，得出图形
关于分数阶的matlab程序
05-20
这是一个分数阶的程序，可以通过此函数来辨识出传递函数的系数，欢迎大家来下载~
caputo分数阶导数程序_Caputo分数阶导数,Caputo fractional derivative,音标,读音,翻译,英文例句,英语词典...
weixin_29492649的博客
12-24 943
补充资料：分数阶积分与微分分数阶积分与微分og fractional integration and differentia-分数阶积分的逆运算称为分数阶微分:若几介F，则f为F的:阶分数阶导数(na ctional deriVative).若00: ;、一上一f一工鱼一一添 r回几恤一t)’-(对f给予适当的限制;见!IL那里还包含算子人关于乌的估...
Caputo 分数阶一维问题基于 L1 逼近的空间二阶方法(附Matlab代码)
图灵的猫的博客
06-20 3371
Caputo 分数阶一维问题基于 L1 逼近的快速差分方法(附Matlab程序),考虑时间分数阶慢扩散方程初边值问题并结合数值算例进行分析
find函数matlab_【重点】分数阶微积分和分数阶微分方程数值实验（6）——Mittag-Leffler函数...
weixin_39640646的博客
11-28 1408
参考文献薛定宇《分数阶微积分学与分数阶控制》数值实现Matlab 2019a 主要基于薛定宇开发的FOTF工具箱蜜酒厅通讯社 固体地球物理学部封面及文中照片感谢 @CycleUser 友情提供前情回顾形式主义的居士：分数阶微积分和分数阶微分方程数值实验（5）——超几何函数​zhuanlan.zhihu.comMittag-Leffler函数一、单参数Mittag-Leffler函数1.定义： 其...
时间域的品质因子，在频率域可以通过引入复速度实现，请用c++程序体现黏声波动方程的复速度的实部与虚部
最新发布
04-07
黏声波动方程可以表示为： ∂u/∂t - β∇²u - α∇(∇·u) = iωρu 其中，u是声波的速度，β是介质的弹性模量，α是介质的黏性系数，ω是角频率，ρ是介质的密度。 引入复速度v，则有： v = u + iωρ/β u 实部为： Re(v) = u 虚部为： Im(v) = iωρ/β u 下面是用c程序实现黏声波动方程的复速度的实部与虚部的代码： #include <stdio.h> #include <stdlib.h> #include <math.h> #define N 100 // 网格数 #define L 1.0 // 区域长度 #define dx (L/N) // 网格间距 #define dt (0.01*dx*dx) // 时间步长 int main() { double u[N][N], v[N][N]; // 实部和虚部 double beta = 1.0, alpha = 0.1, rho = 1.0, omega = 2*M_PI; // 参数 int i, j, k; // 初始化 for(i=0; i<N; i++) { for(j=0; j<N; j++) { u[i][j] = exp(-(i*dx-L/2)*(i*dx-L/2)/0.01)*exp(-(j*dx-L/2)*(j*dx-L/2)/0.01); v[i][j] = 0.0; } } // 时间循环 for(k=0; k<1000; k++) { // 计算实部和虚部 for(i=1; i<N-1; i++) { for(j=1; j<N-1; j++) { v[i][j] = u[i][j] + I*omega*rho/beta*u[i][j]; } } // 计算黏声波动方程 for(i=1; i<N-1; i++) { for(j=1; j<N-1; j++) { u[i][j] += dt*(beta*(u[i+1][j]+u[i-1][j]+u[i][j+1]+u[i][j-1]-4*u[i][j])/dx/dx - alpha*(u[i+1][j]-u[i-1][j])*(u[i][j+1]-u[i][j-1])/dx/dx + I*omega*rho/beta*v[i][j]); } } // 输出结果 printf("%d %f %f\n", k, creal(v[N/2][N/2]), cimag(v[N/2][N/2])); } return 0; }
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值