MatLab常见函数和运算符号解读
MatLab常见函数和运算符号
基本运算
convhull :凸壳函数
cumprod :累计积
cumsum :累计和
cumtrapz :累计梯形数值积分
delaunay :Delaunay三角化
dsearch :求最近点(这是两个有趣的函数 factor :质数分解inpolygon :搜索多边形内的点
max :最大元素
mean :平均值
median :数组的中间值
min :最小值
perms :向量所有排列组成矩阵
polyarea :多边形的面积
primes :生成质数列表
prod :数组元素积
sort :元素按升序排列
sortrows :将行按升序排列
std :标准差
sum :元素和
trapz :梯形数值积分
tsearch :搜索Delaunay三角形var :方差
voronoi :Voronoi图
del2 :Laplacian离散
diff :差分和近似微分gradient:数值梯度
corrcoef :相关系数
cov :协方差矩阵
xcorr :互相关系数
xcov :互协方差矩阵
xcorr2 :二维互相关
conv :卷积和多项式相乘conv2 :二维卷积
deconv :反卷积
filter :滤波
filter2 :二维数字滤波
傅立叶变换
abs :绝对值和模
angle :相角
cplxpair :按复共扼把复数分类
fft :一维快速傅立叶变换
fft2 :二维快速傅立叶变换
fftshit :将快速傅立叶变换的DC分量移到谱中央ifft :以为逆快速傅立叶变换
ifft2 :二维逆快速傅立叶变换
ifftn :多维逆快速傅立叶变换
ifftshift :逆fft平移
nextpow2 :最相邻的2的幂
unwrap :修正相角
cross :向量叉积
intersect:集合交集
ismember :是否集合中元素
setdiff :集合差集
setxor :集合异或(不在交集中的元素
union :两个集合的并
unique :返回向量作为一个集合所有元素(去掉相同元素基本数学函数
abs :绝对值
acos :反余弦
acosh :反双曲余弦函数
acot :反余切
acoth :反双曲线余切
acsc :反余割
acsch :反双曲线余割
angle :相位角
asec :反正割
asech :反双曲线正割
asin :反正弦
asinh :反双曲线正弦
atan :反正切
atanh :反双曲线正切
atan2 :四象限反正切
ceil :向正无穷方向舍入:ceil(4.3=5
complex : 复数
conj :求共扼
cos :余弦
cosh :双曲余弦
cot :余切
coth :双曲线余切
csc :余割
csch :双曲线余割
exp :指数
fix :向零舍入
floor :向负无穷大舍入
gcd :最大公约数
imag :复数的虚部
lcm :最小公倍数
log :自然对数
log2 :以2为底的对数
log10 :以10为底的对数
mod :模除
nchoosek :二项式系数nchoosek(n,k=n!/{k!(n-k!} real :复数实部rem :余数
round :四舍五入
sec :正割
sech :双曲线正割
sign :符号函数
sin :正弦
sinh :双曲线正弦
sqrt :平方根
tan :正切
tanh :双曲线正切
特殊函数
airy :airy函数
besselh:第三类贝塞尔函数; besseli
besselk:改良型besselh函数besselj
bessely:贝塞尔函数
beta
betainc
betaln :贝塔函数
ellipj :雅克比椭圆函数
ellipke:完全椭圆积分
erf
erfc
erfcx
erfinv :误差函数
expint :指数积分
factorial:阶乘函数
gamma
gammaln
gammainc:伽马函数
legendre:勒让德函数
pow2 :2的幂次
rat
rats :有理逼近
坐标变换
cart2pol :笛卡儿坐标变换为极坐标或圆柱坐标cart2sph :笛卡儿坐标变换为球坐标
pol2cart :极坐标变换为笛卡儿坐标
sph2cart :球坐标变换为笛卡儿坐标
矩阵和数组基础
blkdiag :构造一个分块对角矩阵
eye :创建单位矩阵
linespace :产生线性间隔向量logspace :产生对数间隔向量
ones :产生元素全为1的数组
rand :产生均匀分布随机数组
randn :产生正态分布随机数组
zeros :产生全零数组
特殊变量和常量
ans :对最近输入的反应
computer :当前计算机类型
eps :浮点精度
flops :计算浮点操作次数,现已不再常用i :虚部单位
inf :无穷大
inputname :输入参数名
j :虚部单位
nan :非数值
nargin :输入参数的数目
nargout :输出参数的数目(用户定义函数
pi :圆周率
realmax :最大正浮点数
realmin :最小正浮点数
varargin,varargout:返回参数数目(matlab函数时间和日期
calendar :返回日历
clock :当前时间
cputime :从启动经过的cpu时间
date :日期
datenum :转换成连续日期数
datestr :日期字符串
datevec :日期分隔部分
eomday :一个月的日期
etime :经过的时间
now :当前的日期和时间
tic, toc :秒表定时器
weekday :星期几
矩阵操作
cat :把矩阵按行或列连接起来
diag :给定向量,构造对角矩阵
fliplr :矩阵左右翻转
flipud :矩阵上下翻转
repmat :复制数组,repmat(A,m,n表示把A复制m行n列组成新数组reshape :按逐列来的方式重新整形数组
rot :按逆时针方向旋转90度
tril :返回一个矩阵的下三角矩阵
triu :返回一个矩阵的上三角矩阵
特殊函数矩阵
compan :返回向量的伴随矩阵
eig :矩阵特征值
gallery :测试矩阵,或者说大约50个矩阵模版
hadamard :哈达马得矩阵
hankel :汉克尔矩阵
hilb :希尔波特矩阵
invhilb :逆希尔波特矩阵
magic :魔术方阵
pascal :帕斯卡矩阵
toeplitz :托普利茨矩阵
wilkinson :维尔金森特征值测试矩阵算数运算符
+:矩阵加
- :减
* :乘
.* :数组乘
^ :矩阵乘方
.^ :数组乘方
\ :矩阵左除
/ :矩阵右除
.\ :数组左除
./ :数组右除
kron :克罗内克张量积
关系运算符
& :逻辑与
| :逻辑或 ~ :逻辑非 xor :逻辑异或∶ :冒号,用于创建数组和下表 [] :方括号，构成向量或矩阵 ( :圆括号，表示算术表达式优先级、放置函数参数、放置矩阵下标{} :单元数组(cell array专用 . :句点，小数点表示、数组运算符组成部分、字段(属性访问 .. :父目录，与 cd 一起使用 ... :连续三点，一行末尾表示续行 , :逗号，格开参数 ; :分号，表示矩阵一行结束或者语句末尾使得执行结果不显示 % :注释符 ' :单引号，矩阵后表示转置或者成对出现括起字符串 .' :数组转置 = :赋值逻辑函数 all :向量中是否所有分量均非零 any :是否有元素非零 exist :指定变量或文件是否存在find :返回矩阵中非零元素的索引和置 is* :is 系列，大部分根据名称可以知道功能iscell( iscellstr( isempty( :是否空数组 isequal( :数组是否相等 isfield( :是否结构数组中的字段 isfinite( :数组中元素是否有限 isglobal( :是否全局变量 ishandle( :是否有效图形句柄 ishold( :hold 命令处于 on 状态与否 isieee :是否使用 IEEE 算法 isinf( :是否无穷大 isletter( :数组元素是否是字母表中字母 islogical( :是否逻辑数组 isnan( :是否非数值数 isnumeric( :是否数值数组 isobject( :是否对象 isprime( :是否质数 isreal( :是否实数 isspace( :是否空字符 issparse( :是否是按稀疏类别存储 isstruct( :是否结构类isstudent :matlab 版本是否学生版
isunix :是否 unix 版本 isa( :指定对象是否属于指定类 logical( :把数值数组转变为逻辑数组 A(B :A 为数值数组,B 为逻辑数组时，按 B 索引的非零元素返回 A 中相应位置元素 mislocked( :当前文件是否锁定仿真命令： sim ---仿真运行一个simulink 模块 sldebug ---调试一个 simulink 模块 simset ---设置仿真参数 simget ---获取仿真参数线性化和整理命令： linmod ---从连续时间系统中获取线性模型linmod2 ---也是获取线性模型，采用高级方法 dinmod ---从离散时间系统中获取线性模型 trim ---为一个仿真系统寻找稳定的状态参数构建模型命令： open_system --打开已有的模型 close_system --关闭打开的模型或模块 new_system --创建一个新的空模型窗口 load_system --加载已有的模型并使模型不可见 save_system --保存一个打开的模型 add_block --添加一个新的模块 add_line --添加一条线（两个模块之间的连线） delete_block --删除一个模块 delete_line --删除一根线 find_system --查找一个模块 hilite_system --使一个模块醒目显示 replace_block --用一个新模块代替已有的模块 set_param --为模型或模块设置参数 get_param --获取模块或模型的参数
add_param --为一个模型添加用户自定义的字符串参数 delete_param --从一个模型中删除一个用户自定义的参数 bdclose --关闭一个 simulink 窗口 bdroot --根层次下的模块名字 gcb --获取当前模块的名字 gcbh --获取当前模块的句柄 gcs --获取当前系统
的名字 getfullname --获取一个模块的完全路径名 slupdate --将 1.x 的模块升级为 3.x 的模块 addterms --为未连接的端口添加 terminators 模块
boolean slhelp --将数值数组转化为布尔值 --simulink 挠没虻蓟蛘吣？榘镏?nbsp; 封装命令: hasmask --检查已有模块是否封装 hasmaskdlg --检查已有模块是否有封装的对话框 hasmaskicon --检查已有模块是否有封装的图标 iconedit --使用 ginput 函数来设计模块图标 maskpopups --返回并改变封装模块的弹出菜单项 movemask --重建内置封装模块为封装的子模块诊断命令： sllastdiagnostic --上一次诊断信息sllasterror --上一次错误信息 sllastwarning --上一次警告信息 sldiagnostics --为一个模型获取模块的数目和编译状态硬拷贝和打印命令： frameedit --编辑打印画面 print --将 simulink 系统打印成图片，或将图片保存为 m 文件 printopt --打印机默认设置orient --设置纸张的方向 help rtw 看看，rtw 相关的命令有三个命令： rtwgen -从一个模型中创建一个 rtw 文件(model.rtw，用此函数可以指定一些 rtw 的属性设置 tlc -调用目标语言编译器 rtwbuild -对一个模型调用 rtw 的 build 程序相关的命令好像
还有 make_rtw,rtw_c,tlc_c 以上摘自：水木清华站
实验四 MATLAB符号运算
实验四 MATLAB 符号运算 一、实验目的 掌握符号变量和符号表达式的创建，掌握MATLAB 的symbol 工具箱的一些基本应用。 二、实验内容 (1) 符号变量、表达式、方程及函数的表示。 (2) 符号微积分运算。 (3) 符号表达式的操作和转换。 (4) 符号微分方程求解。 三、实验步骤 1. 符号运算的引入 在数值运算中如果求x x x πsin lim 0→，则可以不断地让x 接近于0，以求得表达式接近什么数，但是终究不能令0=x ，因为在数值运算中0是不能作除数的。MATLAB 的符号运算能解决这类问题。输入如下命令： >>f=sym('sin(pi*x)/x') >>limit(f,'x',0) >> f=sym('sin(pi*x)/x') f = sin(pi*x)/x >> limit(f,'x',0) ans = Pi 2. 符号常量、符号变量、符号表达式的创建 1) 使用sym( )创建 输入以下命令，观察Workspace 中A 、B 、f 是什么类型的数据，占用多少字节的内存空间。 >> A=sym('1') >> B=sym('x') >> f=sym('2*x^2+3*y-1') >> clear >> f1=sym('1+2') >> f2=sym(1+2) >> f3=sym('2*x+3') >> f4=sym(2*x+3) >> x=1 >> f4=sym(2*x+3) > A=sym('1') A = 1
>> B=sym('x') B = x >> f=sym('2*x^2+3*y-1') f = 2*x^2+3*y-1 >> clear >> f1=sym('1+2') f1 = 1+2 >> f2=sym(1+2) f2 = 3 >> f3=sym('2*x+3') f3 = 2*x+3 >> f4=sym(2*x+3) ??? Undefined function or variable 'x'. >> x=1 x = >> f4=sym(2*x+3) f4 =
matlab符号运算
MATLAB程序设计教程(9)——MATLAB符号计算 by：ysuncn（欢迎转载，请注明原创信息） 第9章MATLAB符号计算 9.1 符号对象 9.2 符号微积分 9.3 级数 9.4 符号方程求解 9.1 符号对象 9.1.1 建立符号对象 1．建立符号变量和符号常量 MATLAB提供了两个建立符号对象的函数：sym和syms，两个函数的用法不同。 (1) sym函数 sym函数用来建立单个符号量，一般调用格式为： 符号量名=sym('符号字符串') 该函数可以建立一个符号量，符号字符串可以是常量、变量、函数或表达式。 应用sym函数还可以定义符号常量，使用符号常量进行代数运算时和数值常量进行的运算不同。
下面的命令用于比较符号常量与数值常量在代数运算时的差别。 (2) syms函数 函数sym一次只能定义一个符号变量，使用不方便。MATLAB提供了另一个函数syms，一次可以定义多个符号变量。syms函数的一般调用格式为： syms 符号变量名1 符号变量名2 … 符号变量名n 用这种格式定义符号变量时不要在变量名上加字符串分界符(‘)，变量间用空格而不要用逗号分隔。 2．建立符号表达式 含有符号对象的表达式称为符号表达式。建立符号表达式有以下3种方法： (1)利用单引号来生成符号表达式。 (2)用sym函数建立符号表达式。 (3) 使用已经定义的符号变量组成符号表达式。 9.1.2 符号表达式运算 1．符号表达式的四则运算 符号表达式的加、减、乘、除运算可分别由函数symadd、symsub、symmul和symdiv来实现，幂运算可以由sympow来实现。
实验MATLAB符号计算
实验四符号计算 符号计算的特点：一，运算以推理解析的方式进行，因此不受计算误差积累问题困扰；二，符号计算，或给出完全正确的封闭解，或给出任意精度的数值解（当封闭解不存在时）；三，符号计算指令的调用比较简单，经典教科书公式相近；四，计算所需时间较长，有时难以忍受。 在MATLAB中，符号计算虽以数值计算的补充身份出现，但涉及符号计算的指令使用、运算符操作、计算结果可视化、程序编制以及在线帮助系统都是十分完整、便捷的。 MATLAB的升级和符号计算内核Maple的升级，决定着符号计算工具包的升级。但从用户使用角度看，这些升级所引起的变化相当细微。即使这样，本章还是及时作了相应的更新和说明。如MATLAB 6.5+ 版开始启用Maple VIII的计算引擎，从而克服了Maple V计算“广义Fourier变换”时的错误（详见第5.4.1节）。 5.1符号对象和符号表达式 5.1.1符号对象的生成和使用 【例5.1.1-1】符号常数形成中的差异 a1=[1/3,pi/7,sqrt(5),pi+sqrt(5)] % <1> a2=sym([1/3,pi/7,sqrt(5),pi+sqrt(5)]) % <2> a3=sym([1/3,pi/7,sqrt(5),pi+sqrt(5)],'e') % <3> a4=sym('[1/3,pi/7,sqrt(5),pi+sqrt(5)]') % <4> a24=a2-a4 a1 = 0.3333 0.4488 2.2361 5.3777 a2 = [ 1/3, pi/7, sqrt(5), 6054707603575008*2^(-50)] a3 = [ 1/3-eps/12, pi/7-13*eps/165, sqrt(5)+137*eps/280, 6054707603575008*2^(-50)] a4 = [ 1/3, pi/7, sqrt(5), pi+sqrt(5)] a24 = [ 0, 0, 0, 189209612611719/35184372088832-pi-5^(1/2)] 【例5.1.1-2】演示：几种输入下产生矩阵的异同。 a1=sym([1/3,0.2+sqrt(2),pi]) % <1> a2=sym('[1/3,0.2+sqrt(2),pi]') % <2> a3=sym('[1/3 0.2+sqrt(2) pi]') % <3> a1_a2=a1-a2 % a1 = [ 1/3, 7269771597999872*2^(-52), pi] a2 = [ 1/3, 0.2+sqrt(2), pi] a3 = [ 1/3, 0.2+sqrt(2), pi] a1_a2 = [ 0, 1.4142135623730951010657008737326-2^(1/2), 0]
matlab符号运算函数大全
m a t l a b符号运算函数大 全 The Standardization Office was revised on the afternoon of December 13, 2020
算术符号操作 命令 +、-、*、.*、\、.\、/、./、^、.^、’、.’ 功能符号矩阵的算术操作 用法如下： A+B、A-B 符号阵列的加法与减法。 若A与B为同型阵列时，A+B、A-B分别对对应分量进行加减；若A与B中至少有一个为标量，则把标量扩大为与另外一个同型的阵列，再按对应的分量进行加减。 A*B 符号矩阵乘法。 A*B为线性代数中定义的矩阵乘法。按乘法定义要求必须有矩阵 A的列数等于矩阵B的行数。即：若 A n*k* B k*m=(a ij)n*k.*(b ij)k*m= C n*m=(c ij)n*m，则，i=1,2,…,n； j=1,2,…,m。或者至少有一个为标量时，方可进行乘法操作，否则 将返回一出错信息。 A.*B 符号数组的乘法。 A.*B为按参量A与B对应的分量进行相乘。A与B必须为同型 阵列，或至少有一个为标量。即： A n*m.* B n*m=(a ij)n*m.*(b ij)n*m= C n*m=(c ij)n*m，则c ij= a ij* b ij， i=1,2,…,n；j=1,2,…,m。 A\B 矩阵的左除法。 X=A\B为符号线性方程组A*X=B的解。我们指出的是，A\B近 似地等于inv(A)*B。若X不存在或者不唯一，则产生一警告信 息。矩阵A可以是矩形矩阵（即非正方形矩阵），但此时要求方 程组必须是相容的。 A.\B 数组的左除法。 A.\B为按对应的分量进行相除。若A与B为同型阵列时， A n*m.\ B n*m=(a ij)n*m.\(b ij)n*m= C n*m=(c ij)n*m，则c ij= a ij\ b ij，i=1,2,…,n； j=1,2,…,m。若若A与B中至少有一个为标量，则把标量扩大为 与另外一个同型的阵列，再按对应的分量进行操作。 A/B 矩阵的右除法。 X=B/A为符号线性方程组X*A=B的解。我们指出的是，B/A粗 略地等于B*inv(A)。若X不存在或者不唯一，则产生一警告信 息。矩阵A可以是矩形矩阵（即非正方形矩阵），但此时要求方 程组必须是相容的。 A./B 数组的右除法。 A./B为按对应的分量进行相除。若A与B为同型阵列时， A n*m./ B n*m=(a ij)n*m./(b ij)n*m= C n*m=(c ij)n*m，则c ij= a ij/b ij，i=1,2,…,n； j=1,2,…,m。若A与B中至少有一个为标量，则把标量扩大为与 另外一个同型的阵列，再按对应的分量进行操作。 A^B 矩阵的方幂。
MATLAB符号计算实验报告
实验六符号计算 学院：数计学院班级：1003班姓名：黄晓丹学号：1051020144 一、实验目的 1、了解富符号对象和数值对象之间的差别，以及它们之间的互相转换 2、了解符号运算和数值运算的特点、区别和优缺点 3、掌握符号对象的基本操作和运算，以及符号运算的基本运用 二、实验内容 1、符号常数形成和使用 （1）符号常数形成中的差异 >> a1=[1/3,pi/7,sqrt(5),pi+sqrt(5)] a1 = 0.3333 0.4488 2.2361 5.3777 >> a2=sym([1/3,pi/7,sqrt(5),pi+sqrt(5)]) a2 = [ 1/3, pi/7, sqrt(5),
6054707603575008*2^(-50)] >> a3=sym('[1/3,pi/7,sqrt(5),pi+sqrt(5)]') a3 = [ 1/3, pi/7, sqrt(5), pi+sqrt(5)] >> a24=a2-a3 a24 = [ 0, 0, 0, 189209612611719/35184372088832-pi-5^(1/2)] （2）把字符表达式转化为符号变量 >> y=sym('2*sin(x)*cos(x)') y = 2*sin(x)*cos(x) >> y=simple(y)
y = sin(2*x) （3）用符号计算验证三角等式 >> syms fai1 fai2;y=simple(sin(fai1)*cos(fai2)-cos(fai1)*sin(fai2)) y = sin(fai1-fai2) （4）求矩阵的行列式值、逆和特征值 >> syms a11 a12 a21 a22;A=[a11,a12;a21,a22] A = [ a11, a12] [ a21, a22] >> DA=det(A),IA=inv(A),EA=eig(A) DA =
完整word版,MATLAB符号运算
符号运算 科学计算包括数值计算和符号计算两种计算，数值计算是近似计算；而符号计算则是绝对精确的计算。 符号变量的生成和使用 1、符号变量、符号表达式和符号方程的生成 （1）、使用sym函数定义符号变量和符号表达式 单个符号变量 sqrt(2) sym(sqrt(2)) %显示精确结果 a=sqrt(sym(2)) %显示精确结果 double(a) sym(2)/sym(3) %显示精确结果 2/5+1/3 sym(2/5+1/3) %显示精确结果 sym(2)/sym(5)+sym(1)/sym(3) %显示精确结果 sym函数定义符号表达式：单个变量定义法，整体定义法 单个变量定义法 a=sym('a') b=sym('b') c=sym('c') x=sym('x') f=a*x^2+b*x+c 整体定义法 f=sym('a*x^2+b*x+c') g=f^2+4*f-2 （2）、使用syms函数定义符号变量和符号表达式 一次可以创建任意多个符号变量syms var1 var2 var3… syms a b c x f=a*x^2+b*x+c g=f^2+4*f-2 （3）、符号方程的生成 函数：数字和变量组陈的代数式 方程：函数和等号组成的等式 用sym函数生成符号方程： equation1=sym('sin(x)+cos(x)=1') 2、符号变量的基本操作 （1）、findsym函数用于寻找符号变量 findsym(f)：找出f表达式中的符号变量 findsym(s,n)：找出表达式s中n个与x接近的变量 syms a alpha b x1 y findsym(alpha+a+b)
matlab符号运算符
Matlab符号运算符的使用 一、&&/||/&/| |：数组逻辑或 ||：先决逻辑或 &：数组逻辑与 &&：先决逻辑与 &&和||被称为&和|的short circuit形式。 先决逻辑符号含义： 先判断左边是否为真；若为真，则不再判断右边；若为假，才继续进行或运算 先判断左边是否为假；若为假，则不再判断右边；若为真，才继续进行与运算两种运算符号的区别： 先决逻辑运算的运算对象只能是标量 数组逻辑运算可为任何维数组，运算符两边维数要相同 举例分析： A&B ：首先判断A的逻辑值，然后判断B的值，然后进行逻辑与的计算。 A&&B：首先判断A的逻辑值，如果A的值为假，就可以判断整个表达式的值为假， 就可以判断整个表达式的值为假，就不需要再判断B的值。这种用法非常有用， 如果A是一个计算量较小的函数，B是一个计算量较大的函数，那么首先判断A 对减少计算量是有好处的。 另外这也可以防止类似被0除的错误。 Matlab中的if和while语句中的逻辑与和逻辑或都是默认使用short-circuit形式。// 这可能就是有时候用&和| 会报错的原因。
二、系统结构体内的变量 一般都是小写。 matlab区分大小写。 三、== 表示逻辑相等，返回结果，相等为1，不等为0。 四、.*（times）点乘 times Array multiply 数组乘 Syntax c = a.*b c = times(a,b) Description c = a.*b multiplies arrays a an d b element-by-element and returns th e result in c. Inputs a and b must have the same size unless one is a scalar. 注释：a、b要同尺寸，或其中一个为标量。 c = times(a,b) is calle d for th e syntax a.*b when a or b is an object. Example a = [1 2 3]'; b = [5 6 7]'; c = a.*b; 五、矩阵或向量共轭转置“’”和转置“.’” 若矩阵由实数构成，二者作用一样；
matlab实验五多项式和符号运算
实验五：Matlab多项式和符号运算 一、实验目的 1．掌握Matlab多项式的运算。 2．了解符号运算。 二、实验内容 1．将多项式()(2)(3)(7)(1) =-+-+化为x的降幂排列。 P x x x x x syms x; y=(x-2)*(x+3)*(x-7)*(x+1); expand(y) ans = x^4-5*x^3-19*x^2+29*x+42 2．求一元高次方程的根。 98765432 --++--++= 53015027313658204100576-28800 x x x x x x x x x syms x y; y=x^9-5*x^8-30*x^7+150*x^6-1365*x^4-820*x^3+410 0*x^2+576*x-2880; solve(y,x) ans = 6.81947687944124431946 1.42761488953013276419+.8192491831*i 2.865487219+2.49263348244446271927*i
-1.887673354+1.812452594*i -.9583509633 -5.922730991 -1.887673354-1.812452594*i 2.865487219-2.49263348244446271927*i 1.42761488953013276419-.8192491831*i 3．求一元高次方程的根，并画出左边多项式函数在[2,2] x∈-区间内的曲线。 42 -+= x x 210 a=[1 0 -2 0 1]; r=roots(a) syms x; x=-2:2; y=[1 0 -2 0 1]; plot(x,y) r = 1.0000 + 0.0000i 1.0000 - 0.0000i -1.0000 -1.0000
MatLab常见函数和运算符号解读
MatLab常见函数和运算符号 基本运算 convhull :凸壳函数 cumprod :累计积 cumsum :累计和 cumtrapz :累计梯形数值积分 delaunay :Delaunay三角化 dsearch :求最近点(这是两个有趣的函数 factor :质数分解inpolygon :搜索多边形内的点 max :最大元素 mean :平均值 median :数组的中间值 min :最小值 perms :向量所有排列组成矩阵 polyarea :多边形的面积 primes :生成质数列表 prod :数组元素积 sort :元素按升序排列 sortrows :将行按升序排列
std :标准差 sum :元素和 trapz :梯形数值积分 tsearch :搜索Delaunay三角形var :方差 voronoi :Voronoi图 del2 :Laplacian离散 diff :差分和近似微分gradient:数值梯度 corrcoef :相关系数 cov :协方差矩阵 xcorr :互相关系数 xcov :互协方差矩阵 xcorr2 :二维互相关 conv :卷积和多项式相乘conv2 :二维卷积 deconv :反卷积 filter :滤波 filter2 :二维数字滤波
傅立叶变换 abs :绝对值和模 angle :相角 cplxpair :按复共扼把复数分类 fft :一维快速傅立叶变换 fft2 :二维快速傅立叶变换 fftshit :将快速傅立叶变换的DC分量移到谱中央ifft :以为逆快速傅立叶变换 ifft2 :二维逆快速傅立叶变换 ifftn :多维逆快速傅立叶变换 ifftshift :逆fft平移 nextpow2 :最相邻的2的幂 unwrap :修正相角 cross :向量叉积 intersect:集合交集 ismember :是否集合中元素 setdiff :集合差集 setxor :集合异或(不在交集中的元素 union :两个集合的并
matlab符号计算实验报告
1. 已知x=6,y=5, 利用符号表达式求z =>> syms x >> z=(x+1)/(sqrt(x+3)-sqrt(y)); >> subs(z,x,5) ans =6/(8^(1/2)-y^(1/2)) >> subs(ans,6) ans = 15.8338 2. 分解因式。 （1）x y -44； >> syms x y >> factor(x^4-y^4) ans =(x-y)*(x+y)*(x^2+y^2) （2）x x x +++642 12575151 >> syms x >> factor(125*x^6+75*x^4+15*x^2+1) ans =(5*x^2+1)^3 3. 化简表达式 （1）sin cos cos sin ββββ-1212； >> syms x y >> f=sin(x).*cos(y)-cos(x).*sin(y)； >> sfy1=simple(f) 结果：sfy1 =sin(x-y) （2）x x x +++248321 >> syms x >> f=(4*x^2+8*x+3)/(2*x+1);sfy1=simplify(f) sfy1 =2*x+3 4、求下列极限，将完成实验的程序写到文件sy1.m 中： （1） （2） （3） （4） （5） （1）>> syms x >> F1=atan(x)/(x); >> w=limit(F1) w =1 （2）>> syms x F2=((1+x)/(1-x))^(1/x); >> w=limit(F2) w =exp(2) （3）>> syms x F3=(x.*log(1+x))/(sin(x^2)); >> w=limit(F3) w =1 （4）>> syms x F4=atan(x)/(x); >> w=limit(F4,x,inf) w =0 （5）>> syms x F5=(1/(1-x)-1/(1-x^3)); >> w=limit(F5,x,1) w =NaN 5、求下列函数的导数，将完成实验的程序写到文件sy2.m 中： 1、 >> x = sym('x'); >> y1=(cos(x))^3-cos(3*x); >> diff(y1)ans =-3*cos(x)^2*sin(x)+3*sin(3*x) 2、 >> x = sym('x'); >> y2=x.*sin(x).*(log(x)); >> diff(y2)ans =sin(x)*log(x)+x*cos(x)*log(x)+sin(x) 3、 >> x = sym('x'); >> y3=(x.*exp(x)-1)/sin(x); >> diff(y3) ans =(exp(x)+x*exp(x))/sin(x)-(x*exp(x)-1)/sin(x)^2*cos(x) 4、 x x x x F 1011lim 2??? ??-+=→3 1115lim()11x F x x →=---20sin )1ln(lim 3x x x F x +=→x x F x arctan lim 10→=arctan 4lim x x F x →∞=x x y 3cos cos 13-=x x x y ln sin 2=x xe y x sin 13-=cos x y e x =
MATLAB实验——符号运算讲解
实验一符号运算 班级:电气4班姓名:叶元亮学号:B2012052409 一、实验目的 1、了解符号、数值、字符等数据类型的差别 2、了解符号运算的特点、优缺点 3、掌握符号变量的创建和运算，以及其运算的基本应用 4、掌握基本的符号绘图指令 二、实验内容 1、指出下面的 M1，M2，M3 分别是什么，并上机验证。 取a=1、b=2、c=3、d=4，M1=[a,b;c,d]，M2='[a,b;c,d]'，M3=sym('[a,b;c,d]'); >> a=1,b=2,c=3,d=4 a = 1 b = 2 c = 3 d = 4 >> M1=[a,b;c,d] M1 =
1 2 3 4 >> M2='[a,b;c,d]' M2 = [a,b;c,d] >> M3=sym('[a,b;c,d]') M3 = [ a, b] [ c, d] 结论:M1是矩阵，2是字符串，M3是字符变量。 2、下面2种取值情况下，计算b a b a- + 并赋给相应情况下的c1、c2，问c1、c2相等吗，为什么？上机验证。 （1） a1=1010; b1=10-10; （2）将a1、a2作为符号变量赋给a2、b2; >> a1=1e10; b1=1e-10; >> c1=(a1+b1-a1)/b1 c1 = >> a2=sym(a1); b2=sym(b1); >> c2=(a2+b2-a2)/b2 c2 = 1
结果:c1~=c2，因为c1=0,c2=1,a1、b1是具体的数值，a2、b2是符号变量。 3、符号表达式中自由变量的确定生成符号变量a 、b 、x 、X 、Y 、 k=3、z=a y w c sin +，表达式为 Y k bx azX f )(2++=。 （1）找出f 中的全部自由符号变量 （2）在f 中确定最优先的自由符号变量 （3）在f 中确定2个和3个自由变量时的执行情况 （4）试通过对各符号变量与x 的ASCII 值做绝对差值，分析自 由变量优秀顺序，能得出什么结论？ >> syms a b x X Y k=sym('3'); z=sym('c*sqrt(w)+y*sin(a)'); f=a*z*X+(b*x^2+k)*Y; >> findsym(f) ans = X, Y, a, b, c, w, x, y >> findsym(f,1) ans = x >> findsym(f,2) ans = x,y
实验3 Matlab 符号运算及求函数极值
实验3 Matlab 符号运算及求函数极值一、实验目的和要求 掌握用Matlab软件进行符号运算以及求函数的极值。 二、实验环境 Windows系列操作系统，Matlab软件。 三、实验内容 1.用MATLAB进行符号运算； 2.编程求函数的极值。 四、实验步骤 3.开启软件平台——Matlab，开启Matlab编辑窗口； 4.根据求解步骤编写M文件； 5.保存文件并运行； 6.观察运行结果(数值或图形)； 7.根据观察到的结果和体会写出实验报告。 五、示例 1．计算一元函数的极值 例1求 2 2 344 1 x x y x x ++ = ++ 的极值 解首先建立函数关系： s yms x y=(3*x^2+4*x+4)/( x^2+x+1); 然后求函数的驻点： dy=diff(y); xz=solve(dy) xz= [0] [-2] 知道函数有两个驻点x 1=0和x 2 =-2， 接下来我们通过考察函数的图形，则它的极值情况和许多其它特性是一目了然的。而借助MATLAB的作图功能，我们很容易做到这一点。 例2 画出上例中函数的图形
解 syms x y=(3*x^2+4*x+4)/( x^2+x+1); 得到如下图形 ezplot(y) 2．计算二元函数的极值 MATLAB 中主要用diff 求函数的偏导数,用jacobian 求Jacobian 矩阵。 例1 求函数42823z x xy y =-+-的极值点和极值. 首先用diff 命令求z 关于x,y 的偏导数 >>clear; syms x y; >>z=x^4-8*x*y+2*y^2-3; >>diff(z,x) >>diff(z,y) 结果为 ans =4*x^3-8*y ans =-8*x+4*y 即348,84z z x y x y x y ??=-=-+??再求解方程，求得各驻点的坐标。一般方程组的符号解用solve 命令，当方程组不存在符号解时，solve 将给出数值解。求解方程的MA TLAB 代码为：
实验二2MATLAB地符号计算与可视化
实验二MATLAB的符号计算与可视化 1:完成教材实验三第1节“1.创建符号表达式和符号表达式的操作”中（1）-（5）部分的内容，分别用sym和syms创建符号表达式f和g，并对它们进行相关操作，思考每一条命令的作用是什么，并提交命令行和结果； （1）创建符号变量。 ①使用sym命令创建符号表达式： >> f=sym('sin(x)') f = sin(x) >> g=sym('y/exp(-2*t)') g = y*exp(2*t) ②使用syms命令创建符号表达式： >> syms x y t >> f=sym(sin(x)) f = sin(x) >> g=sym(y/exp(-2*t)) g = y*exp(2*t) （2）：自由变量的确定：
>> symvar(g) ans = [ t, y] >> symvar(g,1) ans = y >> findsym(g,2) ans = y,t （3）：用常数替换符号变量： >> x=0:10; >> y=subs(f,x) y = Columns 1 through 8 0 0.8415 0.9093 0.1411 -0.7568 -0.9589 -0.2794 0.6570 Columns 9 through 11 0.9894 0.4121 -0.5440 练习：用y替换x，查看结果及其数据类型。 z=subs(f,y) z = Columns 1 through 8
0 0.7456 0.7891 0.1407 -0.6866 -0.8186 -0.2758 0.6107 Columns 9 through 11 0.8357 0.4006 -0.5176 >> class(z) ans = double （4）：符号对象与数值的转换和任意精度控制： >> f1=subs(f,'5') f1 = sin(5) >> y1=double(f1) y1 = -0.9589 >> y2=eval(f1) y2 = -0.9589 练习：将y1用sym函数转换为符号对象，并用’d’,’f’,’e’,’r’4种格式表示。>> y2=sym(y1,'d') y2 = -0.95892427466313845396683746002964
第9章MATLAB符号计算_习题答案
第9章MATLAB符号计算 习题9 一、选择题 1 .设有a=sym(4)。则1/a+1/a 的值是( A . 0.5 B . 1/2 2 .函数factor(sym(15))的值是( A . '15' B. 15 3 .在命令行窗口输入下列命令： >> f=sym(1); >> eval(i nt(f,1,4)) 则命令执行后的输 出结果是 A . 3 4 . MATLAB A . tailor 5. MATLAB A . solve 二、填空题 1. 在进行符号运算之前首先要 建立符号对象，sym, syms 2. 对于“没有定义”的极限， 大的极限，MATLAB给出的结果为 3. 在命令行窗口输入下列命 令： >> syms n; >> s=symsu m(n ,1,10) 命令执行后s的 值是 ________________________ ， 4. 在MATLAB 中，函数 )。B C . 1/4+1/4 D . 2/a )。D C . [ 1, 3, 5] D . [ 3, 5] ,所使用的函数或命令有__________ 和 ________________________________ 代表________ 。符号代数方程，求解变量 5. 在MATLAB符号计算中 三、应用题 1 .分解因式。 (1) x9-1 (3) 125X6+75X4+15X2+1)。A B . 4 C . 5 D . 1将函数展开为幕级数，所使用的函数是( )。D B . tayler C . diff 用于符号常微分方程求解的函数是( )。C B . solver C . dsolve D . taylor D . dsolver MATLAB给出的结果为 _________ ;对于极限值为无穷 _______ 。 NaN, Inf 55 solve(s,v)用于代数方程符号求解，其中s代表________ , v y的二阶导数表示为__________ 。D2y (2) X4+X3+2X2+X+1 / 、 2 2 2 (4) X +y +z +2(xy+yz+zx) (1):
Matlab符号变量
Matlab的符号运算功能强大，看了些资料，都比较啰嗦，然后再次总结为一个m 文件测试大部分符号运算功能%% 符号变量与符号表达式%%%%%%%%%%%%%%%%%%%%%%%%%%% %1.符号变量与符号表达式 %%%%%%%%%%%%%%%%%%%%%%%%%%% clear all ; clc; close all; % f =sym( 'sin(x)+5x') % f ——符号变量名 % sin(x)+5x——符号表达式 % ' '——符号标识 % 符号表达式一定要用' ' 单引号括起来matlab才能识别 % ' ' 的内容可以是符号表达式，也可以是符号方程。 % 例： % f1=sym('a*x^2+b*x+c') ——二次三项式 % f2=sym('a*x^2+b*x+c=0' )——方程 % f3=sym('Dy+y^2=1') ——微分方程 % 符号表达式或符号方程可以赋给符号变量，以后调用方便；也可以不赋给符号变量直接参与运算 % syms 命令用来建立多个符号量，一般调用格式为： % syms 变量1 变量2 ... 变量n %% 符号矩阵的创建 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %2.符号矩阵的创建 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % 数值矩阵A=[1,2;3,4] % A=[a,b;c,d] ——不识别 % @1.用matlab函数sym创建矩阵（symbolic的缩写) % 命令格式：A=sym('[ ]') % ※ 符号矩阵内容同数值矩阵 % ※ 需用sym指令定义 % ※ 需用' '标识 % 例如： A = sym('[a , 2*b ; 3*a , 0]') % A = % [ a, 2*b] % [3*a, 0] % 这就完成了一个符号矩阵的创建。 % 注意：符号矩阵的每一行的两端都有方括号，这是与 matlab数值矩阵的一个重要区别。%@2.用字符串直接创建矩阵(这种方法创建的没有什么用处)
MATLAB符号计算函数用法总结
MATLAB符号计算函数用法总结 符号计算是对未赋值的符号对象（可以是常数、变量、表达式）进行运算和处理。MTALAB具有符号数学工具箱（Symbolic Math toolbox），将符号运算结合到MATLAB的属具运算环境。符号数学工具箱是建立在Maple软件基础上的。 算术符号操作： 命令有：+、-、*、.*、\、.\、/、./、^、.^、’、.’ 用法如下： A+B、A-B符号阵列的加法和减法。 若A与B为同型阵列时，A+B、A-B分别对对应分量进行加减；若A与B中至少有一个为标量，则把标量扩大为与另外一个同型的阵列，再按对应的分量进行加减。 A*B符号矩阵乘法。 A*B为线性代数中定义的矩阵乘法。按乘法定义要求必须有矩阵A的列数等于矩阵B的行数。即：若 An*k*Bk*m=(aij)n*k.*(bij)k*m=Cn*m=(cij)n*m，则，i=1,2,…,n；j=1,2,…,m。 或者至少有一个为标量时，方可进行乘法操作，否则将返回一出错 信息。 A.*B符号数组的乘法。 A.*B为按参量A与B对应的分量进行相乘。A与B必须为同型阵列，或至少有一个为标量。即： An*m.*Bn*m=(aij)n*m.*(bij)n*m=Cn*m=(cij)n*m，则cij= aij* bij，i=1,2,…,n； j=1,2,…,m。 A\B矩阵的左除法。 X=A\B为符号线性方程组A*X=B的解。我们指出的是，A\B近似地等于inv(A)*B。若X不存在或者不唯一，则产生一警告信息。矩阵A可以是矩形矩阵（即非正方形矩阵），但此时要
求方程组必须是相容的。 A.\B数组的左除法。 A.\B为按对应的分量进行相除。若A与B为同型阵列时， An*m.\Bn*m=(aij)n*m.\(bij)n*m=Cn*m=(cij)n*m，则cij= aij\ bij，i=1,2,…,n； j=1,2,…,m。若若A与B中至少有一个为标量，则把标量扩大为与另外一个同型的阵列，再按对应的分量进行操作。 A/B矩阵的右除法。 X=B/A为符号线性方程组X*A=B的解。我们指出的是，B/A粗略地等于B*inv(A)。若X不存在或者不唯一，则产生一警告信息。矩阵A可以是矩形矩阵（即非正方形矩阵），但此时要求方程组必须是相容的。 A./B数组的右除法。 A./B为按对应的分量进行相除。若A与B为同型阵列时， An*m./Bn*m=(aij)n*m./(bij)n*m=Cn*m=(cij)n*m，则cij= aij/bij，i=1,2,…,n； j=1,2,…,m。若A与B中至少有一个为标量，则把标量扩大为与另外一个同型的阵列，再按对应的分量进行操作。 A^B矩阵的方幂。 计算矩阵A的整数B次方幂。若A为标量而B为方阵，A^B用方阵B的特征值与特征向量计算数值。若A与B同时为矩阵，则返回一错误信息。 A.^B数组的方幂。 A.^B为按A与B对应的分量进行方幂计算。若A与B为同型阵列时， An*m..^Bn*m=(aij)n*m..^(bij)n*m=Cn*m=(cij)n*m，则cij= aij^bij，i=1,2,…,n； j=1,2,…,m。若A与B中至少有一个为标量，则把标量扩大为与另外一个同型的阵列，再按对应的分量进行操作。 A'矩阵的Hermition转置。 若A为复数矩阵，则A'为复数矩阵的共轭转置。即，若A=(aij)=(xij+i*yij)，则 。
实验四 MATLAB符号运算
实验四MATLAB符号运算 一、实验目的： 1、掌握定义符号对象的方法； 2、掌握符号表达式的运算法则以及符号矩阵运算。 3、掌握求符号函数极限及导数的方法。 4、掌握求符号函数定积分和不定积分的方法。 二、实验原理 1、符号常量、符号变量、符号表达式的创建 (1) 使用sym( )创建 输入以下命令，观察Workspace 中A、B、f是什么类型的数据，占用多少字节的内存空间。 >>A=sym('1') %符号常量 >>B=sym('x') %符号变量 >>f=sym('2*x^2+3y-1') %符号表达式 >>clear >>f1=sym('1+2') %有单引号，表示字符串 >>f2=sym(1+2) %无单引号 >>f3=sym('2*x+3') >>f4=sym(2*x+3) %为什么会出错 >>x=1 >>f4=sym(2*x+3) 通过看MATLAB 的帮助可知，sym( )的参数可以是字符串或数值类型，无论是哪种类型都会生成符号类型数据。 (2) 使用syms 创建 >>clear >>syms x y z %注意观察x,y,z都是什么类型的，它们的内容是什么 >>x,y,z >>f1=x^2+2*x+1 >>f2=exp(y)+exp(z)^2 >>f3=f1+f2 通过以上实验，知道生成符号表达式的第二种方法：由符号类型的变量经过运算(加减乘除等)得到。又如： >>f1=sym('x^2+y +sin(2)') >>syms x y >>f2=x^2+y+sin(2) >>x=sym('2') , y=sym('1') >>f3=x^2+y+sin(2)
matlab符号运算函数大全
2.1 算术符号操作 命令+、-、*、.*、\、.\、/、./、^、.^、’、.’ 功能符号矩阵的算术操作 用法如下： A+B、A-B 符号阵列的加法与减法。 若A与B为同型阵列时，A+B、A-B分别对对应分量进行加减；若A与B中至少有一个为标量，则把标量扩大为与另外一个同型的阵列，再按对应的分量进行加减。 A*B 符号矩阵乘法。 A*B为线性代数中定义的矩阵乘法。按乘法定义要求必须有矩阵A的列数等于矩阵B的 行数。即：若A n*k*B k*m=(a ij)n*k.*(b ij)k*m=C n*m=(c ij)n*m，则，i=1,2,…,n；j=1,2,…,m。或者 至少有一个为标量时，方可进行乘法操作，否则将返回一出错信息。 A.*B 符号数组的乘法。 A.*B为按参量A与B对应的分量进行相乘。A与B必须为同型阵列，或至少有一个为 标量。即：A n*m.*B n*m=(a ij)n*m.*(b ij)n*m=C n*m=(c ij)n*m，则c ij= a ij* b ij，i=1,2,…,n；j=1,2,…,m。A\B 矩阵的左除法。 X=A\B为符号线性方程组A*X=B的解。我们指出的是，A\B近似地等于inv(A)*B。若X 不存在或者不唯一，则产生一警告信息。矩阵A可以是矩形矩阵（即非正方形矩阵）， 但此时要求方程组必须是相容的。 A.\B 数组的左除法。 A.\B为按对应的分量进行相除。若A与B为同型阵列时， A n*m.\ B n*m=(a ij)n*m.\(b ij)n*m= C n*m=(c ij)n*m，则c ij= a ij\ b ij，i=1,2,…,n；j=1,2,…,m。若若A与B 中至少有一个为标量，则把标量扩大为与另外一个同型的阵列，再按对应的分量进行操 作。 A/B 矩阵的右除法。 X=B/A为符号线性方程组X*A=B的解。我们指出的是，B/A粗略地等于B*inv(A)。若X 不存在或者不唯一，则产生一警告信息。矩阵A可以是矩形矩阵（即非正方形矩阵）， 但此时要求方程组必须是相容的。 A./B 数组的右除法。 A./B为按对应的分量进行相除。若A与B为同型阵列时， A n*m./ B n*m=(a ij)n*m./(b ij)n*m= C n*m=(c ij)n*m，则c ij= a ij/b ij，i=1,2,…,n；j=1,2,…,m。若A与B 中至少有一个为标量，则把标量扩大为与另外一个同型的阵列，再按对应的分量进行操 作。 A^B 矩阵的方幂。 计算矩阵A的整数B次方幂。若A为标量而B为方阵，A^B用方阵B的特征值与特征 向量计算数值。若A与B同时为矩阵，则返回一错误信息。 A.^B 数组的方幂。 A.^B为按A与B对应的分量进行方幂计算。若A与B为同型阵列时， A n*m..^ B n*m=(a ij)n*m..^(b ij)n*m= C n*m=(c ij)n*m，则c ij= a ij^b ij，i=1,2,…,n；j=1,2,…,m。若A与B 中至少有一个为标量，则把标量扩大为与另外一个同型的阵列，再按对应的分量进行操 作。 A' 矩阵的Hermition转置。 若A为复数矩阵，则A'为复数矩阵的共轭转置。即，若A=(a ij)=(x ij+i*y ij)，则。 A.' 数组转置。 A.'为真正的矩阵转置，其没有进行共轭转置。 例2-1
2014秋实验四_MATLAB的符号计算二
2014年10月26日星期日 湖北科技学院 MATLAB 专业：电气工程及其自动化班级：电气一班 姓名：尹锥（133521015）指导教师：刘芳华 电子与信息工程学院
实验四MATLAB的符号计算（二） 一、实验目的 1.掌握MATLAB7.0 subs、simple、simplify、finverse等函数的应用； 2.掌握利用MATLAB7.0计算极限和级数的方法，计算复合、反函数的方法； 3.了解利用MATLAB7.0计算taylor级数。 4.掌握利用MATLAB7.0来求解常微分方程的方法。 5.掌握利用MATLAB7.0符号计算中的ezplot及ezplot3绘图方法。 二、实验内容 1.设x为符号变量， 42 ()21 f x x x =++，32 ()635 g x x x x =+++，试进行如下运算。 （1） ()() f x g x +， （2） ()() f x g x ?， （3）对 () f x进行因式分解， （4）求 () g x的反函数。 2
2.(1)指出下面程序中的f1、f2、f3、f4、f5的值。f5=subs(f,{a,x},{0:6,0:pi/6:pi}) 3
(2)指出下面替换的结果。 3. 用符号计算验证三角等式： sin(?1)cos(?2)-cos(?1)sin(?2) =sin(?1-?2) 4
4.设,求 7.求微分方程的解：y'''-y''=x，y''(2)=4，y'(1)=7，y(1)=8 5
8.计算函数级数 2 1 n x S n ∞ = =∑ 9. f(x)=e x分别求5阶、6阶泰勒展开式 10.符号函数绘图法绘制函数x=sin(3t)cos(t)，y=sin(3t)sin(t)的图形，t的变化范围为[0,2π]。 6