最近在学电动力学,要手算各种矢量微分,耗时又易错。为此,我编写了MATLAB符号运算矢量微分的函数,用MATLAB算微分,省了很多精力。这里分享给大家。

函数名:vd
作者:Freeman
日期:2020.5
应用:在笛卡尔坐标系、柱坐标系和球坐标系下计算梯度、散度、旋度与拉普拉斯量
输入:
output=vd(u,v,difftype,cood)
u - 被求导函数
v - 函数自变量,笛卡尔坐标系下为[x y z],柱坐标系下为[s phi z],球坐标系下为[r theta phi]
difftype - 求导类型,梯度为’g’,散度为’d’,旋度为’c’,拉普拉斯量为’l’
cood - 坐标系,笛卡尔坐标系为’car’,柱坐标系为’cyl’,球坐标系为’sph’
输出:散度和拉普拉斯量为符号标量,梯度和旋度为三维符号向量

举例:
1、在笛卡尔坐标系下求 u(x,y,z)=x+2y+3zu(x,y,z) =x+2y+3zu(x,y,z)=x+2y+3z 的梯度
输入:

clc;clear;
syms x y z
vd(1*x+2*y+3*z,[x y z],'g','car')

输出:

ans =
[ 1, 2, 3]

2、在球坐标系下求 u⃗(r,θ,ϕ)=r2r^\vec{u}(r,\theta,\phi)=r^2\hat{r}u(r,θ,ϕ)=r2r^ 的散度
输入:

clc;clear;
syms r the phi
u=[r^2,0,0]; % u是矢量函数,写成三个分量的形式。注意这三个分量必须对应于该坐标系,不一定是x,y,z分量。
vd(u,[r,the,phi],'d','sph') % 即使表达式中没有出现the和phi,也要补全在vd函数里。

输出:

ans =
4*r

(m文件中有英文说明)

下载链接:
https://download.csdn.net/download/FreemanThe/12426455

MATLAB自定义函数 计算三元函数矢量符号微分相关推荐

  1. matlab solve函数计算三元一次方程组

    Matlab solve函数计算三元一次方程组 %matlab代码 clc; clear; Lambda = 0.0001; %失效率 Mu = 2; %修复率 syms p0 p1 p2; %定义三 ...

  2. 对ndarray的操作:一元函数、二元函数、聚合函数、三元函数

    python高级应用与数据分析学习笔记 10 1.概念笔记 # ndarray的一元函数 # abs fabs 计算整数.浮点数.复数的绝对值,对于非复数,fabs更快 # sqrt 计算各个元素的平 ...

  3. matlab:变计算精度函数vpa,digits

    vpa vpa(x)使用可变精度浮点运算(vpa)对符号输入x中的每个元素求值到至少d位有效数字,其中d是digits函数的值. 数字的默认值为32. vpa (x,d)使用至少d位有效数字,而不是数 ...

  4. Python基础知识3:函数、三元函数、lambda函数

    函数:包含自定义函数.内置函数 1)自定义函数结构:包含五个部分:

  5. 使用函数计算分段函数的值

    #include<stdio.h> int sign(int x) {int z;if(x>0)z=1;if(x==0)z=0;if(x<0)z=-1;return z; } ...

  6. matlab学习—分段函数计算

    分段函数计算 分段函数用if函数,代码如下(示例): function f=fun(x) if x==10f=cos(x+1)+sqrt(x^2+1); else f=x*sqrt(x+sqrt(x) ...

  7. 【函数计算月报】2018年10月刊

    重点新闻 阿里新任 CEO 张勇首次电视采访:云计算将是阿里未来主要业务 阿里巴巴新任 CEO 张勇在接受 CNBC 的电视采访时表示坚信未来的每项业务都将由云计算提供支持,非常乐意在新的数字时代建立 ...

  8. 函数计算的开发与配置

    作者 | 夏莞 阿里云函数计算开发工程师 导读: 在本篇文章中"基本概念"部分主要对函数计算最核心的概念进行详细介绍,包括服务.函数.触发器.版本.别名以及相关的配置:" ...

  9. 开发函数计算的正确姿势 —— 移植 next.js 服务端渲染框架

    为什么80%的码农都做不了架构师?>>>    首先介绍下在本文出现的几个比较重要的概念: 函数计算(Function Compute): 函数计算是一个事件驱动的服务,通过函数计算 ...

最新文章

  1. Bootstrap响应式与自适应区别
  2. 2017-2018-1 20155338 《信息安全系统设计基础》 第三周学习总结
  3. Windows下更改pip镜像源
  4. boost::mpl模块实现empty_sequence相关的测试程序
  5. CPU资源和可用内存大小
  6. python画两条曲线_查找在matplotlib中绘制的两条曲线之间的区域(在区域之间填充)...
  7. java sql2005驱动_java.sql.SQLException:找不到适用于jdbc:microsoft:sqlserver的驱动程序...
  8. 成毅高清电脑桌面壁纸
  9. 备份自己的myeclipse快捷键到需要的环境
  10. 大话云时代rac_网易Q1财报中的增长信号:有道和云音乐如何打通“任督二脉”?...
  11. linux系统 mysql日志文件太大。造成数据库无法正常启动怎么解决
  12. 小数点后最多几位_一个“90后”女生的月账单,花得最多的果然是这个
  13. 买了小区一楼的感受是怎样的?
  14. 用PHP生成随机数的函数
  15. BZOJ5221[Lydsy2017省队十连测] 偏题
  16. 最新QQ搭建代挂系统2.0全新打包内带i优支付在线充值系统
  17. web前端学习13-19(HTML常用标签)
  18. word中表格出现无法自动换页问题
  19. 51单片机物联网开发
  20. 在疫情期间微分销商城如何运营

热门文章

  1. 云计算、大数据、物联网、人工智能都是什么关系?
  2. 激光雷达与视觉联合标定综述!(系统介绍/标定板选择/在线离线标定等)
  3. 浏览器缓存机制详细图解
  4. 电网窃电用户数据检测与识别数据集
  5. Sublime使用详细总结
  6. 点到直线垂足坐标的计算
  7. 学习英语01——听力——剑雅12 T7S3
  8. Oracle数据库中查询结果集中保留两位小数的问题
  9. bat脚本 获取当前所在文件夹的路径
  10. for循环多重嵌套的流程图