matlab 微分符号,Matlab 符号微积分
6.3 符号微积分
6.3.1 符号序列的求和
【 * 例 6.3.1 -1 】求 ,
syms k t;f1=[t k^3];f2=[1/(2*k-1)^2,(-1)^k/k];
s1=simple(symsum(f1)) %f1 的自变量被确认为 t
s2=simple(symsum(f2,1,inf)) %f2 的自变量被确认为 k
s1 =
[ 1/2*t*(t-1), k^3*t]
s2 =
[ 1/8*pi^2, -log(2)]
6.3.2 符号微分和
矩阵
【 * 例 6.3.2 -1 】求
、
和
syms a t x;f=[a,t^3;t*cos(x), log(x)];
df=diff(f) % 求矩阵 f 对 x 的导数
dfdt2=diff(f,t,2) % 求矩阵 f 对 t 的二阶导数
dfdxdt=diff(diff(f,x),t) % 求二阶混合导数
df =
[ 0, 0]
[ -t*sin(x), 1/x]
dfdt2 =
[ 0, 6*t]
[ 0, 0]
dfdxdt =
[ 0, 0]
[ -sin(x), 0]
【 * 例 6.3.2 -2 】求
的
矩阵。
syms x1 x2 x3;f=[x1*exp(x2);x2;cos(x1)*sin(x2)];
v=[x1 x2];fjac=jacobian(f,v)
fjac =
[ exp(x2), x1*exp(x2)]
[ 0, 1]
[ -sin(x1)*sin(x2), cos(x1)*cos(x2)]
6.3.3 符号积分
6.3.3.3 符号积分示例
【 * 例 6.3.3 .3-1 】求
。演示:积分指令对符号函数矩阵的作用。
syms a b x;f=[a*x,b*x^2;1/x,sin(x)];
disp('The integral of f is');pretty(int(f))
The integral of f is
[ 2 3]
[1/ 2 a x 1/3 b x ]
[ ]
[ log(x) -cos(x) ]
【 * 例 6.3.3 .3-2 】求
。演示如何使用 mfun 指令获取一组积分值。
(1)求一般积分结果
F1=int('1/log(t)','t',0,'x')
F1 =
-Ei(1,-log(x))
(2)利用 mfun 指令求 x=0.5 , 0.6 , 0.7 , 0.8 , 0.9 时的定积分
x=0.5:0.1:0.9
F115=-mfun('Ei',1,-log(x))
x =
0.5000 0.6000 0.7000 0.8000 0.9000
F115 =
-0.3787 -0.5469 -0.7809 -1.1340 -1.7758
【 * 例 6.3.3 .3-3 】求积分
。注意:内积分上下限都是函数。
syms x y z
F2=int(int(int(x^2+y^2+z^2,z,sqrt(x*y),x^2*y),y,sqrt(x),x^2),x,1,2)
VF2=vpa(F2) % 积分结果用 32 位数字表示
F2 =
1610027357/6563700-6072064/348075*2^(1/2)+14912/4641*2^(1/4)+64/225*2^(3/4)
VF2 =
224.92153573331143159790710032805
【 * 例 6.3.3 .3-4 】利用 rsums 求
积分。(与例 6.3.3.3-2 结果比较)
syms x positive;px=0.5/log(0.5*x);rsums(px)
图 6.3.3 .3-4 交互式近似积分
6.3.4 符号卷积
【 * 例 6.3.4 -1 】本例演示卷积的时域积分法:已知系统冲激响应
,求
输入下的输出响应。
syms T t tao;ut=exp(-t); % 定义系统输入
ht=exp(-t/T)/T; % 定义系统冲激响应
uh_tao=subs(ut,t,tao)*subs(ht,t,t-tao); % 运用变量替换指令形成被积函数
yt=int(uh_tao,tao,0,t); % 实施卷积
yt=simple(yt)
yt =
-(exp(-t)-exp(-t/T))/(T-1)
【 * 例 6.3.4 -2 】本例演示通过变换和反变换求取卷积。系统冲激响应、输入同上例,求输出。
syms s;yt=ilaplace(laplace(ut,t,s)*laplace(ht,t,s),s,t);yt=simple(yt)
yt =
-(exp(-t)-exp(-t/T))/(T-1)
【 * 例 6.3.4 -3 】求函数
和
的卷积。
(1)在 5.2 版(配 Symbolic Math Toolbox 2.0.1 )中,采用以下指令。
syms t tao;ut=sym('Heaviside(t)-Heaviside(t-1)');ht=t*exp(-t);
yt=int(subs(ut,t,tao)*subs(ht,t,abs(t)-tao),tao,0,abs(t));
yt=collect(yt,'signum(abs(t)-1)'), yt=subs(yt,abs(t),t)%<3>
yt =
(-1/2+1/2*exp(-abs(t)+1)*abs(t))*signum(abs(t)-1)-exp(-abs(t))*abs(t)-exp(-abs(t))+1/2*exp(-abs(t)+1)*abs(t)+1/2
yt =
(-1/2+1/2*exp(-t+1)*t)*signum(t-1)-exp(-t)*t-exp(-t)+1/2*exp(-t+1)*t+1/2
(2) 5.3 版(配 Symbolic Math Toolbox 2.1 )中的运行指令和结果如下。
syms tao;t=sym('t','positive'); % 把 t 定义为限定性符号变量 <4>
ut=sym('Heaviside(t)-Heaviside(t-1)');ht=t*exp(-t);
yt53=int(subs(ut,t,tao)*subs(ht,t,t-tao),tao,0,t);
yt53=collect(yt53,'Heaviside(t-1)')
yt53 =
(-1+exp(1-t)*t)*Heaviside(t-1)+1+(-t-1)*exp(-t)
matlab 微分符号,Matlab 符号微积分相关推荐
- matlab第七章符号对象,MATLAB语言:第七章 MATLAB符号计算
<MATLAB语言:第七章 MATLAB符号计算>由会员分享,可在线阅读,更多相关<MATLAB语言:第七章 MATLAB符号计算(33页珍藏版)>请在人人文库网上搜索. 1. ...
- matlab中两个符号矩阵相加,MATLAB矩阵分析及符号运算.ppt
MATLAB矩阵分析及符号运算 第三讲 MATLAB的符号运算 -- matlab 不仅具有数值运算功能,还开发了在matlab环境下实现符号计算的工具包Symbolic Math Toolbox 符 ...
- matlab画图时特殊符号的输入
matlab画图时特殊符号的输入 (2018-03-08 16:00:43) 转自http://blog.sina.com.cn/s/blog_a7c0455c0102yaa6.html 1.在MAT ...
- matlab中的箭头符号怎么打开,MATLAB中上下标、斜体、箭头等符号的使用方法
在MA TLAB中,如何标注上标.下标.斜体.黑体.箭头.上圆圈.正负号等特殊符号?如何标注各种特殊数学符号,如约等于号.积分.乘号.无穷等?本文中详细介绍了这些方法 在使用MATLAB的过程中,经常 ...
- matlab中如何输入积分,在matlab中怎么输入特殊符号,Matlab的符号积分
一.文档中的Text Properties: 下标用 _(下划线) 上标用^ (尖号) 斜体 \it 黑体 \bf << \ll >> \gg 正负 \pm 左箭头 \left ...
- 第十一章:MATLAB:符号运算(符号与数值,符号矩阵)
第十一章:MATLAB符号运算 11.1. 符号与数值 11.1.1. 符号与数值间的转换 实例-数值与符号转换 11.1.2. 符号表达式与数值表达式的精度设置 实例-魔方矩阵的数值解 实例-稀疏矩 ...
- matlab求线性方程组的符号解,符号代数线性方程组的MATLAB求解.PPT
符号代数线性方程组的MATLAB求解 建模常见问题MATLAB求解 Content 符号代数线性方程(组)的MATLAB求解 最优化方法的MATLAB实现 一维搜索问题 线性规划 无约束非线性最优化问 ...
- matlab程序特殊符号,matlab中怎么输入特殊符号
matlab中怎么输入特殊符号: word文档中怎么插入心形特殊符号?,缺失:matlab中怎么输入特殊符号12628/11 以下是的一些我们精选的word文档中怎么插入心形特殊符号? word文档中 ...
- matlab 点符号,matlab 中在运算符号前加一点是什么意思?
matlab 中在运算符号前加一点是什么意思? mip版 关注:210 答案:2 悬赏:20 解决时间 2021-01-23 18:50 已解决 2021-01-23 08:24 例如,矩阵A. ...
- matlab 创建 符号矩阵,如何在MATLAB中动态创建符号矩阵?
我需要在MATLAB中创建一个符号矩阵.它可以静态地完成为如何在MATLAB中动态创建符号矩阵? syms a11 a12 a21 a22; A = [a11 a12; a21 a22]; 或使用紧凑 ...
最新文章
- 在小公司混水摸鱼3年,被面试官直接开怼:能力太差了,还不如应届生哇~
- IOC和DI 控制反转和依赖注入
- php redis删除所有key,php redis批量删除key的方法
- map reduce相关程序
- 小心使用IOCP完成端口
- 从零开始 - iOSRTMP推流篇(1)
- CDN缓存服务器现状
- IOS多选单选相册图片
- B - C语言实验——整数位
- for update行级锁的作用
- VMware网络连接模式——桥接模式、NAT模式以及仅主机模式的介绍和区别
- ^全^ 获取SQL SERVER2000/2005、MySql、Oracle元数据的SQL语句 [SQL语句来自CodeSmith]
- RPC failed; curl 56 GnuTLS recv error (-54): Error in the pull function.
- AntColony 磁力搜索BT引擎的核心
- 易基因 - 外泌体let-7d-3p和miR-30d-5p作为宫颈癌及其癌前病变无创筛查的诊断标志物|早期筛查
- google 安装去广告插件
- oracle to_char 进制转换_ORACLE函数TO_CHAR以及数字转换格式
- html easyui怎么实现折叠面板,Easyui 创建折叠面板_EasyUI 教程
- python输入一个分数、约分后输出_Python中分数的相关使用教程
- Filter过滤器介绍及使用
热门文章
- Theory for the information-based decomposition of stock price
- 工作笔记:如何用Django连接Kerberized甲骨文(Oracle)数据库
- 微信小程序分享生成海报
- 一寸等于多少厘米,一寸照片尺寸是多少?
- Android Compose 新闻App(二)ViewModel、Hlit、数据流
- Non-local Networks Meet Squeeze-Excitation Networks and Beyond 论文解读。
- windows录屏_录屏软件:如何录制游戏视频?简单方法了解下
- 函数参数缺省值/默认值
- pr模板.mogrt格式安装方法
- android 仿QQ音乐歌单效果