当前位置：文档之家 › matlab符号计算实例

matlab符号计算实例

  • 格式：pdf
  • 大小：94.33 KB
  • 文档页数：8

下载文档原格式

  / 8

合集下载

matlab符号运算 多项式

matlab符号运算 多项式

matlab符号运算多项式【提纲】1.MATLAB符号运算简介MATLAB是一款功能强大的数学软件，其中符号运算功能允许用户进行高级数学计算、分析和可视化。

符号运算可以帮助工程师、科学家和数学家在各种领域解决问题，如线性代数、微积分、概率论等。

2.多项式基本概念与MATLAB表示多项式是数学中一个重要的概念，它表示为一个无穷级数，其中包含常数、变量及其幂次。

在MATLAB中，多项式可以用符号表达式表示，如：f(x) = 2x^3 + 4x^2 - 3x + 1。

3.多项式运算实例以下是几个MATLAB中进行多项式运算的实例：- 多项式加法：将两个多项式相加，如f(x) + g(x)。

- 多项式减法：将两个多项式相减，如f(x) - g(x)。

- 多项式乘法：将两个多项式相乘，如f(x) * g(x)。

- 多项式除法：将一个多项式除以另一个多项式，如f(x) / g(x)。

- 多项式求导：对一个多项式求导，如diff(f(x))。

- 多项式积分：对一个多项式进行积分，如int(f(x))。

4.多项式函数与应用MATLAB提供了许多与多项式相关的函数，如：- polyfit：根据一组数据拟合多项式。

- polyval：根据多项式系数计算多项式的值。

- roots：求多项式的根。

- legendre：勒让德多项式。

- laguerre：拉格朗日多项式。

这些函数在信号处理、控制系统、优化等领域具有广泛的应用。

5.总结与建议MATLAB的符号运算功能为多项式计算提供了便捷的工具和函数。

掌握这些功能和函数可以帮助用户在各种应用场景中解决问题。

matlab符号计算练习题目

matlab符号计算练习题目

1 将(x -6)(x -3)(x -8)展开为系数多项式的形式。

2 求解多项式x 3-7x 2+2x +40的根。

3 求解在x =8时多项式(x -1)(x -2) (x -3)(x -4)的值。

4 计算多项式乘法(x 2+2x +2)(x 2+5x +4)。

5 计算多项式除法(3x 3+13x 2+6x +8)/(x +4)。

6 对下式进行部分分式展开：27243645232345234+++++++++x x x x x x x x x7 计算多项式9514124234++--x x x x 的微分和积分。

8 用符号函数法求解方程a t 2+b*t +c=0。

9用符号计算验证三角等式：sin(ϕ1)cos(ϕ2)-cos(ϕ1)sin(ϕ2) =sin(ϕ1-ϕ2) 10 因式分解：6555234-++-x x x x 11 ⎥⎥⎦⎤⎢⎢⎣⎡=)sin()log(12x x e x x a f ax ，用符号微分求df/dx 。

12 求代数方程组⎪⎩⎪⎨⎧=+=++002y x c by ax 关于x,y 的解。

13， 用符号函数法求解方程a t 2+b*t +c=0。

(应用solve)14， 因式分解：6555234-++-x x x x (应用syms, factor) 15， ⎥⎥⎦⎤⎢⎢⎣⎡=)sin()log(12x x e x x a f ax ，用符号微分求df/dx 。

(应用syms,diff)16计算极限求极限：n n h n x M hx h x L )1(lim )2(,)ln()ln(lim )1(0-=-+=∞→→ 17 计算导数22d d ,d d ,d d ,sin x y C a y B x y A ax y ====求．18计算不定积分、定积分、反常积分 x x x x I d )22(1222⎰+-+=，x x x x J d cos sin cos 2/0⎰π+=，⎰+∞-=0d e 2x K x ．19 符号求和求级数 ∑∞=121n n 的和S, 以及前十项的部分和S1． 20解代数方程和常微分方程例如：求一元二次方程a*x^2+b*x+c=0的根． 21求微分方程x y ='的通解．22求微分方程⎩⎨⎧==+=0)0(' ,1)0('"y y y x y 的特解． 23求微分方程组⎩⎨⎧=+=x y xy x 2''的通解．。

Matlab上机操作2——符号计算

Matlab上机操作2——符号计算

Matlab上机操作2——符号计算1、%创建数值常量和符号常量a1=2*sqrt(5)+pi %创建数值常量a2=sym('2*sqrt(5)+pi') %按符号表达式创建符号常量a3=sym(2*sqrt(5)+pi) %按最接近的两个正整数表示符号常量a4=sym(2*sqrt(5)+pi,'d') %按最接近的十进制浮点数表示符号常量a31=a3-a1 %数值常量和符号常量的计算a5='2*sqrt(5)+pi' %字符串常量2、%创建符号变量和符号表达式f1=sym('a*x^2+b*x+c') %创建符号表达式syms a b c x %创建符号变量f2=a*x^2+b*x+c %创建符号表达式syms('a','b','c','x') %注意这里是syms，不是sym>> sym('a','b','c','x')错误使用sym输入参数太多。

f3=a*x^2+b*x+c; %创建符号表达式syms x y real %创建实数符号变量z=x+i*y; %创建z为复数符号变量real(z) %复数z的实部是实数xsym('x','unreal'); %清除符号变量的实数特性real(z) %复数z的实部3、%比较符号矩阵与字符串矩阵A=sym('[a,b;c,d]') %创建符号矩阵B='[a,b;c,d]' %创建字符串矩阵C=[a,b;c,d] %创建数值矩阵C=sym(B) %转换为符号矩阵whos4、%计算符号矩阵的行列式值、非共轭转置和特征值syms a11 a12 a21 a22A=[a11 a12;a21 a22] %创建符号矩阵det(A) %计算行列式A.' %计算非共轭转置eig(A)%计算特征值5、%符号表达式的代数运算f=sym('2*x^2+3*x+4')g=sym('5*x+6')f+g %符号表达式相加f*g %符号表达式相乘6、%对符号表达式进行任意精度控制并用三种运算方式表示同一符号常量a=sym('2*sqrt(5)+pi')digits %显示默认的有效位数vpa(a) %用默认的位数计算并显示vpa(a,20) %按指定的精度计算并显示digits(15) %改变默认的有效位数vpa(a) %按digits指定的精度计算并显示a1 =2/3 %数值型a2 = sym(2/3) %有理数型digitsa3 =vpa('2/3',32) %VPA型format longa17、%符号变量与数值变量进行转换a1=sym('2*sqrt(5)+pi')b1=double(a1) %转换为数值变量a2=vpa(sym('2*sqrt(5)+pi'),32)b2=numeric(a2) %转换为数值变量b2=numeric(a2)未定义函数或变量'numeric'。

matlab符号运算(共57张)

matlab符号运算(共57张)

（3）limit（F）：计算符号函数F的极限值，变量为
findsym确定的默认变量；在未指定目标值时，默认变量趋 近于0； （4）limit（F,x,a,’left’）或limit（F,x,a,’right’）：分别计算 函数F的左极限和右极限。
17
第17页，共57页。
>> syms x t;
>> limit(sin(x)/x)
A= [ a, 2*b] [3*a, 0]
3、用函数syms创建矩阵
格式： syms arg1 arg2… 参数
A=[ ]
>> syms a b
A=
>>A=[a,2*b;3*a,0]
[ a, 2*b] [3*a, 0]
15
第15页，共57页。
4、将数值矩阵(jǔ zhèn)转化成符号矩阵
(jǔ zhèn)
ans =
11/15
3
第3页，共57页。
>> 3*sin(3)+pi/2
ans =
1.9942
>> sym('3*sin(3)+pi/2')
ans = 3*sin(3)+pi/2
>> sym(3*sin(3)+pi/2)
ans = 8980881799167258*2^(-52)
>> sym(3*sin(3)+pi/2,'d')
>> real(z) ans = x
>> syms x unreal >> real(z)
ans =
1/2*x+1/2*conj(x)
matlab符号运算(二)

matlab符号运算(二)

六大常见符号运算
因式分解、展开、合并、简化及通分等
计算极限 limit(f,x,a): 计算 lim f ( x )
xa
limit(f,a): 计算默认自变量趋向于a时f的极限 limit(f): 计算 a=0 时的极限 limit(f,x,a,’right’)：右极限 limit(f,x,a,’left’)：左极限
1 2 n 1 n

，以及其前10项的部分和。
>> syms n >> S=symsum(1/n^2,n,1,inf) >> S10=symsum(1/n^2,n,1,10)
x 2 n 1 n

S=1/6*pi^2 S10=1968329/1270080
例：求函数级数
S
>> syms n x >> S=symsum(x/n^2,n,1,inf)
符号矩阵中元素的引用和修改
>> A=sym(’[1+x, sin(x); 5, exp(x)]’) >> A(1,2) >> A(2,2)=sym(’cos(x)’)
Matlab 符号运算（二）
符号矩阵的基本运算
符号矩阵的基本运算与数值矩阵的基本运算相类似。
1) 基本运算符：+、-、*、\、/、
ans=10
ans=2*x+y
ans=10 ans=[2+y,4+y,6+y] ans=[7 10 13]
ans=3*a+b
?
Matlab 符号运算（二）
符号矩阵
使用sym函数直接生成
>> A=sym(’[1+x, sin(ห้องสมุดไป่ตู้); 5, exp(x)]’)
matlab符号求偏导和全微分及其实验

matlab符号求偏导和全微分及其实验

符号求导命令
功能
zx=diff(f(x,y),x)
求x的一阶偏导
Zy=diff(f(x,y),y)
求y的一阶偏导
Dz=zx*dx+zy*dy
Байду номын сангаас求z的全微分
Zxx=diff(zx,x)
求x的二阶偏导
Zxy=diff(z,x,y)
求z的xy混合二阶偏导
Zxn=diff(f(x,y),x,n)
求x的n阶偏导
Zyn=diff(f(x,y),y,n)
求y的n阶偏导
Ux=diff(f(x,y,z),x)
求三元u的x的偏导
Uy=diff(f(x,y,z),y)
求三元u的y的偏导
Uz=diff(f(x,y,z),z)
求三元u的z的偏导
Du=ux*dx+uy*dy+uz*dz
求u的全微分
Uyx=diff(u,y,x)
举例
求u的xy二阶混合偏微分
Uyxy=diff(u,y,x,y)
求u的xyx的三阶混合偏微分
Zx=-diff(F,x)/diff(F,z)
Zy=-diff(F,y)/diff(F,z)
求隐函数的xy的偏微分
Pretty(diff(f(x,y,z),x))
输出一个符合日常书写习惯的表达式
举例
拉格朗日插值多项式数值计算高阶导数MATLAB程序
end
C=Y*L1;L=Y*1;
syms x dyk
for k=1:n1
k
dyk=diff(L,x,k)
end
计算梯度和方向导数的MATLAB程序及其实验
符号求导命令
功能
[FX,FY]=gradient(F)
matlab符号计算实验报告

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.83382. 分解因式。

（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 +++64212575151 >> syms x >> factor(125*x^6+75*x^4+15*x^2+1) ans =(5*x^2+1)^33. 化简表达式（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+34、求下列极限，将完成实验的程序写到文件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 =NaN5、求下列函数的导数，将完成实验的程序写到文件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⎪⎭⎫ ⎝⎛-+=→31115lim()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 →∞=xx y 3cos cos 13-=xx x y ln sin 2=x xe y xsin 13-=cos x y e x =>> x = sym('x');y=cos(x).*exp(x); >> diff(y) ans =-sin(x)*exp(x)+cos(x)*exp(x)5、 >> x = sym('x');y=x^2.*sin(x); >> diff(y) ans = 2*x*sin(x)+x^2*cos(x)6、求下列函数的积分1、syms x a b c;int(sin(a*x).*sin(b*x).*sin(c*x)) ans =-1/4/(c+a-b)*cos((c+a-b)*x)+1/4/(-c+a-b)*cos((-c+a-b)*x)+1/4/(c+a+b)*cos((c+a+b)*x)-1/4/(-c +a+b)*cos((-c+a+b)*x) 2、>> syms x ;int(x^5+x^3-sqrt(x)/4) ans =1/6*x^6+1/4*x^4-1/6*x^(3/2)3、>> syms x ;int(x.*exp(x)/(x+1)^2,x,0,1) ans =1/2*exp(1)-1 4、 >> syms x y;F=int(int('x/(1+x*y)',x,0,1),y,0,1) F =2*log(2)-15、 由曲面22y x z +=，1=z ，2=z 所围成 >> syms x y z;F=int(int(int('x^2+y^2',x,-inf,inf),y,-inf,inf),z,1,2) F =Inf7、求下列级数的和（1） （2）（1）>> syms n;symsum(2*n-1/2^n,1,Inf) ans =Inf（2）syms n;symsum(1/n*(2*n-1),1,Inf) ans =Inf8、将函数 展开成2-x 的幂级数>> syms x;mtaylor(1/(x^2+5*x-3),x-2) 2sin 2y x x=⎰cxdx bx ax sin sin sin dx x x x )4(35⎰-+⎰+102)1(dx x xe x ⎰⎰+D dxdy xy x 1]1,0[]1,0[⨯=D ⎰⎰⎰Vzdxdydz 11212n n n I ∞=-=∑211(21)n I n n ∞==+∑21()53f x x x =+-。

matlab 求解符号方程

matlab 求解符号方程

一、背景介绍Matlab是一种强大的数学软件，常用于数学建模、仿真、数据分析等领域。

在工程和科学研究中，求解符号方程是一个常见的问题，Matlab提供了丰富的符号计算工具，可以帮助用户高效地求解符号方程。

二、Matlab符号计算工具1. 符号变量定义在Matlab中，我们可以通过syms命令定义符号变量，使用符号变量进行符号运算。

例如：```matlabsyms x y```2. 求解符号方程Matlab提供了solve函数，可以用来求解符号方程。

solve函数的基本语法如下：```matlabsol = solve(equations, variables)```其中，equations表示要求解的方程组，variables表示待求解的变量。

solve函数会返回符号方程的解。

三、示例接下来，我们通过一个示例来演示如何使用Matlab求解符号方程。

假设我们要求解如下的符号方程：```matlabsyms xeqn = x^2 - 4*x + 3 == 0;sol = solve(eqn, x);disp(sol);```运行以上代码，可以得到方程x^2 - 4*x + 3 = 0的解为x = 1或x = 3。

四、注意事项在使用Matlab求解符号方程时，有一些需要注意的事项：1. 可能存在多解或无解的情况，在求解后需要对解进行检查；2. 符号计算是一种复杂的运算，可能存在数值精度问题，需要注意数值的精确性；3. 在求解复杂的方程组时，可能需要对方程组进行化简或变形，以提高求解效率。

五、总结通过Matlab的符号计算工具，我们可以较为方便地求解符号方程，实现高效的符号计算。

在工程和科学研究中，这些工具能够帮助我们快速解决复杂的数学问题，提高工作效率。

希望本文的介绍和示例能够帮助读者更好地理解和应用Matlab的符号计算工具。

Matlab在求解符号方程方面具有广泛的应用。

通过利用Matlab的符号计算工具，用户可以轻松地进行符号方程的求解和符号计算，并获得高精度的结果。

matlab符号运算

matlab符号运算


2、创建符号变量和表达式
（1）使用sym命令创建符号变量和表达式； 语法：创建符号变量：sym(‘arg’,’参数’) 创建符号表达式：sym(‘表达式’) 说明：参数用来设置限定符号变量的数学 特性，可以选择 为：’positive’,’real’,’unreal’表示：正， 实和非实。不限定可以省略。
2、创建符号变量和表达式
例：>>x=sym(‘x’,’real’) >>y=sym(‘y’,’real’) >>z=x+i*y >>real(z) >>x=sym(‘x’,’unreal’) >>real(z) 比较一下两个real(z)命令的结果有什么不同？为 什么呢？
2、创建符号变量和表达式
x t sin( x ) e
>>syms t x
>>g=sym(‘[2*t t^2;tsin(x) exp(x)]’)
>>diff(f)
>>diff(f,’t’)
3、符号积分
函数int可以求符号表达式的积分 语法： int(f,’t’) 求符号变量t的不定积分 int(f,’t’,a,b) 求符号变量t的在[a,b]区间 的定积分，a,b是具体的数值 int(f,’t’,’m’,’n’) 求符号变量t的在[m,n]区间 的定积分，m,n是符号
第三讲 MATLAB的符号运算
—— matlab 不仅具有数值运算功 能，还开发了在matlab环境下实现 符号计算的工具包Symbolic Math Toolbox
一、符号运算入门
例 1 ：求解一元二次方程 x 2 5x 2 0
在matlab中利用命令 solve(‘a*x^2+b*x+c=0’) 可以计算出一元二次方程的解析解。 ans = 1/2/a*(-b+(b^2-4*a*c)^(1/2)) 1/2/a*(-b-(b^2-4*a*c)^(1/2)) 对于例1，可运行如下语句： solve('x^2+5*x+2=0') 即可得到方程的根： ans = -5/2+1/2*17^(1/2) -5/2-1/2*17^(1/2)
实验三 MATLAB符号计算

实验三 MATLAB符号计算

都可以在矩阵意义下进行。但应注意这些函数作用于符号矩阵 时，是分别作用于矩阵的每一个元素。 由于符号矩阵是一个矩阵，所以符号矩阵还能进行有关矩阵的运 算。MATLAB还有一些专用于符号矩阵的函数，这些函数作用于 单个的数据无意义。例如 transpose(s)：返回s矩阵的转置矩阵。 determ(s)：返回s矩阵的行列式值。 其实，曾介绍过的许多应用于数值矩阵的函数，如diag、triu、 tril、
expr1 =
x^3+2*exp(-t)*x^2+(1+exp(-t)^2)*x+exp(-t) expr2 = x*exp(-t)^2+(2*x^2+1)*exp(-t)+(x^2+1)*x
expand使用指令 y=0.14-(1.2e+002)*(-2.4005*(0.445-x)^7+4.2505*(0.445x)^6-2.2336*(0.445-x)^5+0.4993*(0.445-x)^40.0514*(0.445-x)^3+0.0025*(0.445-x)^2);
符号矩阵的生成
符号矩阵可通过函数sym来生成。符号矩阵中的元素是任何不带等号的符 号表达式，各符号表达式的长度可以不相同；符号矩阵中，以空格或逗号 分隔的元素指定的是不同列的元素而分号分隔的元素指定的是不同行的元 素。 例：
syms x; A=sym(‘[cos(x),sin(x),x;-x+1 x^2+x+1 tan(x)]’) A= [ cos(x), sin(x), x] [ -x+1, x^2+x+1, tan(x)] >> size(A) %求符号矩阵的大小 ans = 2 3 > a=[1 2 3 4;4 5 6 7]; >> b=sym(a) b= [ 1, 2, 3, 4] [ 4, 5, 6, 7]
关于matlab符号计算的综合实例

关于matlab符号计算的综合实例

关于matlab符号计算的综合实例在科学计算领域，matlab是一种非常常用的软件工具，它具有强大的符号计算功能，可以对符号表达式进行求导、积分、解方程等操作，为科学研究和工程设计提供了便利。

本文将通过一些实例，介绍matlab符号计算功能的应用和操作技巧。

一、符号变量的定义和基本运算我们需要明白matlab中符号计算需要先定义符号变量。

在matlab中，可以使用syms命令定义符号变量，比如：syms x y这样就定义了两个符号变量x和y。

接下来，我们就可以进行基本的符号运算，比如加法、减法、乘法和除法，示例代码如下：1. 加法：z = x + y2. 减法：w = x - y3. 乘法：4. 除法：v = x / y通过这些简单的示例，我们可以看到，matlab对于符号变量的基本运算操作和数学运算规则是一致的，只是使用符号变量进行运算，可以得到符号表达式作为结果。

二、符号函数的求导和积分在科学计算中，求导和积分是非常常见的操作，matlab可以对符号函数进行求导和积分操作，示例代码如下：1. 求导：f = x^2 + 3*x + 2df = diff(f, x)2. 积分：F = int(f, x)通过这些示例，我们可以看到，matlab可以对符号函数进行求导和积分操作，并得到相应的结果。

这对于解决一些数学问题和工程问题非三、符号方程的求解在科学研究和工程设计中，经常会遇到需要求解符号方程的情况，matlab提供了符号求解方程的功能，示例代码如下：1. 求解一元方程：syms xeqn = x^2 - 4*x + 3 == 0;sol = solve(eqn, x)2. 求解多元方程：syms x yeqn1 = x + y == 3;eqn2 = x - y == 1;sol = solve([eqn1,eqn2],[x,y])通过这些示例，我们可以看到，matlab可以对符号方程进行求解，并得到相应的结果。

matlab符号方程

在MATLAB 中，符号方程（symbolic equation ）是通过符号计算工具箱（Symbolic Math Toolbox ）来处理的。

符号计算允许您使用变量和符号来表示数学表达式，并执行各种代数和微积分操作。

以下是一些MATLAB 中处理符号方程的基本步骤和示例：创建符号变量：使用syms 函数创建符号变量。

例如，syms x y z 将创建三个符号变量x 、y 和z 。

定义符号方程：使用等号（=）将符号表达式定义为方程。

例如，eq = x^2 + y^2 - z^2 == 1将创建一个符号方程eq 。

求解符号方程：使用solve 函数求解符号方程。

例如，solve(eq, x)将求解方程eq 中的x 。

下面是一个具体的示例，展示如何在MATLAB 中创建和求解符号方程：matlab 复制代码%创建符号变量syms xy;%定义符号方程 equ ation =x^2+y^2==1;%求解符号方程solu% 创建符号变量 syms x y; % 定义符号方程 equation = x^2 + y^2 == 1; % 求解符号方程 solutions = solve(equation, x);在上述示例中，我们首先使用syms 函数创建了两个符号变量x 和y 。

然后，我们使用等号（=）创建了一个符号方程equation ，该方程表示一个单位圆。

最后，我们使用solve 函数求解该方程，将结果存储在solutions 变量中。

s =solve(equation, x);syms xy;%定义符号方程equation =x^2+ y^2 == 1; %求解符号方程solutions =solve(equation, x);syms x y; % 定义符号方程 equation = x^2 + y^2 == 1; % 求解符号方程 solutions = solve(equation, x);在上述示例中，我们首先使用syms 函数创建了两个符号变量x 和y 。

matlab阶乘符号

matlab阶乘符号
在MATLAB中，阶乘符号使用"!"表示。

阶乘是指将一个正整数n与小于等于n的所有正整数相乘的结果。

例如，n的阶乘可以表示为n!，计算公式为n! = n (n-1) (n-2) ... 3 2 1。

在MATLAB中，你可以使用阶乘符号来计算一个数的阶乘。

下面是一些示例代码：
1. 计算5的阶乘：
result = factorial(5);
disp(result);
运行以上代码，结果将会输出120，即5的阶乘。

2. 计算一个向量中每个元素的阶乘：
vector = [1, 2, 3, 4, 5];
result = factorial(vector);
disp(result);
运行以上代码，结果将会输出一个包含每个元素阶乘的向量。

需要注意的是，阶乘函数在计算大数的阶乘时可能会面临溢出
的问题。

MATLAB提供了`factorial`函数来计算整数的阶乘，但对
于非整数或较大的数值，可以考虑使用其他方法来计算阶乘，如使
用循环或递归。

matlab符号计算

logo>> randn('state',27)>> startprice=50;>> fracreturnsl=0.00015*randn(200,1)+0.0003;>> x1=[startprice;1+fracreturnsl];>> prices1=cumprod(x1);>> t=(1:length(prices1))';>> randn('state',27)>> fracreturns2=0.0015*randn(200,1)+0.0003;>> x2=[startprices;1+fracreturns2];??? Undefined function or variable 'startprices'.>> x2=[startprice;1+fracreturns2];>> prices=cumprod(x2);whosName Size Bytes ClassL 51x51 20808 double arrayfracreturns2 200x1 1600 double arrayfracreturnsl 200x1 1600 double arrayl1 1x1 8 double arrayl2 1x1 8 double arrayprices 201x1 1608 double arrayprices1 201x1 1608 double arrays 1x1 8 double arraystartprice 1x1 8 double arrayt 201x1 1608 double arrayx1 201x1 1608 double arrayx2 201x1 1608 double arrayGrand total is 4010 elements using 32080 bytes>> plot (fracreturnsl, 'DisplayName', 'fracreturnsl', 'YDataSource', 'fracreturnsl'); figure(gcf) >> hist (fracreturnsl); figure(gcf)>> bar (fracreturnsl, 'DisplayName', 'fracreturnsl', 'YDataSource', 'fracreturnsl'); figure(gcf) >> pie (fracreturnsl); figure(gcf)Warning: Ignoring non-positive data in pie chart.> In pie at 40>> stairs (fracreturnsl, 'DisplayName', 'fracreturnsl', 'YDataSource', 'fracreturnsl'); figure(gcf) >> plot (fracreturnsl, 'DisplayName', 'fracreturnsl', 'YDataSource', 'fracreturnsl'); figure(gcf) >> whosName Size Bytes ClassL 51x51 20808 double arrayfracreturns2 200x1 1600 double arrayfracreturnsl 200x1 1600 double arrayl1 1x1 8 double arrayl2 1x1 8 double arrayprices 201x1 1608 double arrayprices1 201x1 1608 double arrays 1x1 8 double arraystartprice 1x1 8 double arrayt 201x1 1608 double arrayx1 201x1 1608 double arrayx2 201x1 1608 double arrayGrand total is 4010 elements using 32080 bytes>> plot (prices1, 'DisplayName', 'prices1', 'YDataSource', 'prices1'); figure(gcf)>> plot (prices, 'DisplayName', 'prices', 'YDataSource', 'prices'); figure(gcf)>> plot (x1, 'DisplayName', 'x1', 'YDataSource', 'x1'); figure(gcf)>> hist (x1); figure(gcf)>> stem (x1, 'DisplayName', 'x1', 'YDataSource', 'x1'); figure(gcf)>> plot (prices2,'DisplayName','prices2','YDateSource','prices');figure(gef)??? Undefined function or variable 'prices2'.>>>> plot ('prices2','DisplayName','prices2','YDateSource','prices');figure(gef)??? Error using ==> plotMATLAB:Graphics:Plot.>> plot (prices2,'DisplayName','prices2','YDateSource','prices');figure(gcf)??? Undefined function or variable 'prices2'.>> plot (prices2,'DisplayName','prices2','YDateSource','prices2');figure(gcf)??? Undefined function or variable 'prices2'.>> figure; pcolor (x1); figure(gcf)??? Error using ==> surfaceZ must be a matrix, not a scalar or vector.Error in ==> pcolor at 53hh = surface(zeros(size(x)),x,'parent',cax);>> plottools>> plot (prices1, 'DisplayName', 'prices1', 'YDataSource', 'prices1'); figure(gcf)>> plot (fracreturns2, 'DisplayName', 'fracreturns2', 'YDataSource', 'fracreturns2'); figure(gcf) >> plot (prices, 'DisplayName', 'prices', 'YDataSource', 'prices'); figure(gcf)>> plot (fracreturns2, 'DisplayName', 'fracreturns2', 'YDataSource', 'fracreturns2'); figure(gcf)>> plot (fracreturns2, 'DisplayName', 'fracreturns2', 'YDataSource', 'fracreturns2'); figure(gcf) >> plot (t, 'DisplayName', 't', 'YDataSource', 't'); figure(gcf)>> plot (prices1, 'DisplayName', 'prices1', 'YDataSource', 'prices1'); figure(gcf)>> plot (fracreturnsl, 'DisplayName', 'fracreturnsl', 'YDataSource', 'fracreturnsl'); figure(gcf) >> plot (fracreturns2, 'DisplayName', 'fracreturns2', 'YDataSource', 'fracreturns2'); figure(gcf) >> plot (prices, 'DisplayName', 'prices', 'YDataSource', 'prices'); figure(gcf)>> plot (fracreturns2, 'DisplayName', 'fracreturns2', 'YDataSource', 'fracreturns2'); figure(gcf) >> plot (fracreturns2, 'DisplayName', 'fracreturns2', 'YDataSource', 'fracreturns2'); figure(gcf) >> plottoolsWarning: Plot empty.> In legend at 193In insertmenufcn at 53>> logo>> plot (prices, 'DisplayName', 'prices', 'YDataSource', 'prices'); figure(gcf)>> plot (prices1, 'DisplayName', 'prices1', 'YDataSource', 'prices1'); figure(gcf)>> plot (prices, 'DisplayName', 'prices', 'YDataSource', 'prices'); figure(gcf)>> plot (prices1, 'DisplayName', 'prices1', 'YDataSource', 'prices1'); figure(gcf)>> plot (prices1, 'DisplayName', 'prices1', 'YDataSource', 'prices1'); figure(gcf)>> plot (prices, 'DisplayName', 'prices', 'YDataSource', 'prices'); figure(gcf)>>。

matlab符号计算

Matlab符号计算
符号对象和数值对象的差异 x1=2; x2=sym('2'); y1=sin(y2 = sin(2)
符号对象的创建 两个基本函数：sym, syms 例：使用sym函数创建符号变量. y=sym( 'x^3+5*x^2+12*x+20'); syms a b c d ; y a y= x^3+5*x^2+12*x+20 a= a
符号表达式的化简 符号表达式的两个化简函数：simplify， simple , simplify：化简函数,可用于化简各种表达式 例:对表达式f=sin2(x)+cos2(x)进行化简. syms x f=sin(x)^2+cos(x)^2; simplify(f) ans = 1
例：最简表达式的获得。 syms x t f=cos(x)^2-sin(x)^2; simple(f)
符号代数方程组的求解(symbolic algebra equations set solution) 例：分别求解代数方程ax2+bx+c=0和cos(2x)+sin(x)=1 syms a b c x s=a*x^2+b*x+c; 例2： solve(s) [x,y] = solve('x^2 + x*y + solve('cos(2*x)+sin(x)=1') y = 3','x^2 - 4*x + 3 = 0') ans = 1/2/a*(-b+(b^2-4*a*c)^(1/2)) x= 1/2/a*(-b-(b^2-4*a*c)^(1/2)) 1 ans = pi 3 0 y= 1/6*pi 1 5/6*pi

第6章MATLAB符号计算

ans = -1/48*(x^3+x^4)^(1/2)*(16*(x^2+x)^(3/2)+12*(x^2+x)^(1/2)*x+6*(x^2+x)^(1/2)3*log(1/2+x+(x^2+x)^(1/2)))/x/((x+1)*x)^(1/2)
g= 1/3*((x+1)*x)^(1/2)*x^2+1/12*((x+1)*x)^(1/2)*x1/8*((x+1)*x)^(1/2)+1/16*log(1/2+x+((x+1)*x)^(1/2))
第6章MATLAB符号计算
❖ 6.1 符号计算基础 ❖ 6.2 符号导数及其应用 ❖ 6.3 符号积分 ❖ 6.4 级数 ❖ 6.5 代数方程的符号求解 ❖ 6.6 常微分方程的符号求解
例6.1考察符号变量和数值变量的差别。 在 MATLAB命令窗口，输入命令：
>>a=sym('a');b=sym('b');c=sym('c');d=sym('d'); %定义4个符号变量
%求极限(4)
-1/2*2^(1/2)/a^(1/2)
目录
6.2.2 符号函数求导及其应用 MATLAB中的求导的函数为： diff(f,x,n) diff函数求函数f对变量x的n阶导数。参数x的用
法同求极限函数limit，可以缺省，缺省值与 limit相同，n的缺省值是1。
目录
例6.11求函数的导数。
6.3.2 符号函数的定积分
定积分在实际工作中有广泛的应用。在MATLAB中，定 积分的计算使用函数：

第7章MATLAB符号计算

例7.20求数列 fn=e-n的Z变换及其逆变换。
命令如下： syms n z fn=exp(-n); Fz=ztrans(fn,n,z) %求fn的Z变换 f=iztrans(Fz,z,n) %求Fz的逆Z变换
第7章MATLAB符号计算
7.4 级数 7.4.1 级数的符号求和 级数符号求和函数symsum，调用格式为： symsum(a,n,n0,nn)
7.3.2 符号函数的定积分 定积分在实际工作中有广泛的应用。在MATLAB
中，定积分的计算使用函数：
int(f,x,a,b)
例7.6 求定积分。
命令如下：
x=sym('x');t=sym('t');
int(abs(1-x),1,2)
%求定积分(1)
f=1/(1+x^2); int(f,-inf,inf)
f=x^2+y^2+z^2-a^2;
zx=-diff(f,x)/diff(f,z) %求(5)。按隐函数求导公式求z对x的偏导数
zy=-diff(f,y)/diff(f,z) %求(第57)。章M按AT隐LAB函符数号计求算导公式求z对y的偏导数
例7.4 在曲线y=x3+3x-2上哪一点的切线与直 线y=4x-1平行。
第7章MATLAB符号计算
5．符号表达式与数值表达式之间的转换 利用函数sym可以将数值表达式变换成它的符
号表达式。 函数numeric或eval可以将符号表达式变换成
数值表达式。
第7章MATLAB符号计算
7.1.3 符号表达式中变量的确定 MATLAB中的符号可以表示符号变量和符号常数。
findsym可以帮助用户查找一个符号表达式中 的的符号变量。该函数的调用格式为： findsym(S,n) 函数返回符号表达式S中的n个符号变量，若没有 指定n，则返回S中的全部符号变量。 在求函数的极限、导数和积分时，如果用户没有 明确指定自变量，MATLAB将按缺省原则确定 主变量并对其进行相应微积分运算。可用 findsym(S,1)查找系统的缺省变量，事实上， MATLAB按离字符'x'最近原则确定缺省变量。
  1. 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

例题 1：
f
=3
1 x3
+
6 x2
+ 12 + 8 x
syms x f=(1/x^3+6/x^2+12/x+8)^(1/3) g1=simple(f) pretty(f) 例题 2： 公因子法简化表达 clear A=sym('[a b;c d]') [V,D]=eig(A) subexpr([V;D])%提取公因式默认 sigma Dw=subexpr(D,'w')%提取公因式 记为 w [RVD,w]=subexpr([V;D],'w') %与上相同 例题 3： 通用置换指令 （1） clear syms a b x; f=a*sin(x)+b （2）RES=subs(RS,old,new)或者 RES=subs(RS,new) f1=subs(f,sin(x),'log(y)') （3）数值代替 若其中还存在符号变量，表达式仍为符号 f2=subs(f,a,3.11) class(f2) （4）胞元代替 f3=subs(f,{a,b,x},{2,5,sym('pi/3')}) class(f3) （7）数组代替 t=0:pi/10:2*pi; f6=subs(f,{a,b,x},{2,3,t}) plot(t,f6) （8）两次替代
y2=subs(y(2),'C3',k); ezplot(y2,[-6,6,-4,8],1) end hold off box on legend('奇解','通解','Location','Best') ylabel('y') title(['\fontsize{14}微分方程',' (y '')^2 – xy '' + y = 0 ','的解']) 例题 13： 代数方程 S=solve('eq1','eq2',…,'eqn','v1','v2',…,'vn') 或者 S=solve(eq1,eq2,…，eqn,v1,v2,…,vn) 慎用要注意顺序 （1）符号解 eq1=sym('d+n/2+p/2-q'); eq2=sym('d+n-p+q-10'); eq3=sym('d-n/4-p+q'); eq4=sym('-8*d-n+p+q-1'); S=solve(eq1,eq2,eq3,eq4,'d','n','p','q'); disp([' d',' n',' p',' q']) disp([S.d,S.n,S.p,S.q]) （2）数值解 A=sym([1 1/2 1/2 -1;1 1 -1 1;1 -1/4 -1 1;-8 -1 1 1]); b=sym([0;10;0;1]); X1=A\b %X1=inv(A)*b （3）注意，无解析解则给出近视解 clear all,syms x; s=solve('(x+2)^x=2','x') xs=(s+2)^s
例题 12： 微分方程奇解和通解，注意既无显示解又无隐式解则会警告 S 为 空
clear all y=dsolve('(Dy)^2-x*Dy+y=0','x') clf,hold on hy1=ezplot(y(1),[-6,6,-4,8],1); set(hy1,'Color','r','LineWidth',5) for k=-2:0.5:2
（3）
dydx=solve(dgdx1,'dydx'2,….xn)=0， ∂xi ∂x j
=
∂ − ∂x j
∂ ∂ xi
f ( x1, x2 ,..., xn ) f ( x1, x2 ,..., xn )
syms x y;f=(x^2-2*x)*exp(-x^2-y^2-x*y) DyDX=-diff(f,x)/diff(f,y)
（1） clear all reset(symengine) Da=1.2;Dw=1/3; syms sa sw sx sy sz syms A B positive syms C real （2） evalin(symengine,'anames(Properties)') （3） clear A syms （4） evalin(symengine,'anames(Properties)') （5） syms B clear syms （6） evalin(symengine,'anames(Properties)')
4. 符号对象的清除： clear x clear all syms x clear 撤销假设 evalin(symengine,’getprop(x)’) 查看假设
evalin(symengine,’anames(Properties)’) 列出 MuPAD 假设 reset(symengine)
5. 符号数字的精度 digits 默认 32 位 digits(n)或 digits n xs=vpa(x) xs=vpa(x,n) 6. 符号表达式的操作 collect 合并同类项 expand 对指定项展开 factor 因式分解 numden 提取公因式 simplify 恒等式简化
pretty 习惯方式 coeffs 多项式系数 7. 例题
例题 6： 偏导 jacobian（f,v）
syms x1 x2; f=[x1*exp(x2) x2 cos(x1)*sin(x2)]; v=[x1;x2]; Jf=jacobian(f,v) 例题 7： 隐函数求导 （1） syms x g=sym('cos(x+sin(y(x)))=sin(y(x))') dgdx=diff(g,x) % （2） dgdx1=subs(dgdx,'diff(y(x),x)','dydx') % 替换
例题 5： 求导 diff(f,v,n) 默认一阶 r s t u v w x y z 先后关系，x,y,w,z... syms a t x; f=[a,t^3;t*cos(x),log(x)]; df=diff(f,t) dfdt2=diff(f,t,2)%2 阶导 dfdxdt=diff(diff(f,t),t)%2 阶导
符号计算实例
1. 符号计算的优点： 所谓符号计算是指解算数学表达式、方程时，不是在离散化的数值点上进行 的，而是凭借一系列的恒等式和数学定理，通过推理和演算获得的解析结果。这 种计算建立在数值完全准确表达和严格推演的基础之上，因而所得结果完全准 确。当然，也存在者不足，后文将会提到。 符号变量的优点是，使用符号变量运算得到的只是一个解析解，例如， 在符号变量运算过程中 pi 就用 pi 表示，而不是具体的近似数值 3.14 或 3.14159。使用符号变量进行运算能最大限度减少运算过程中因舍入造成的 误差。符号变量也便于进行运算过程的演示。 2. 符号对象的创建： 2.1 单个符号变量 S = sym(A)将非符号对象（如，数字，表达式，变量等）A 转换为符号 对象，并存储在符号变量 S 中。 x = sym('x') 创建符号变量 x，其名字是'x'。示例：alpha = sym('alpha') x = sym('x', 'real') 这 里 假 设 x 是 实 数 ， 因 此 有 x 的 共 轭 conj(x)等 于 x。 示 例 ： r = sym('Rho','real') k = sym('k', 'positive') %09 版不能用这个方法实现 这里创建一个正的（实数）符号变量。 x = sym('x', 'clear') 创建一个没有额外属性的纯形式上的符号变量 x（例如，创建符号变量 x，但是并没指定它是正的或它是一个实数）。为了兼容旧的 MATLAB 版本， x = sym('x','unreal')的功能和 x = sym('x', 'clear')一样。 S = sym(A, flag) 把一个数值标量或矩阵转换为符号型的对象。这里 flag 参数的值可以 是：'r', 'd', 'e', or 'f'，它指定了对浮点数进行转换时的规则： 'f':表示“floating-point”。这样，所有的数值都用 N*2^e 或-N*2^e 这种形 式表示。N 和 e 都为整数，且 N 不小于 0。例如，sym(1/10, 'f')被转换为 3602879701896397/36028797018963968。 'r ': 表 示 “r ati o nal ”（ 有 理 数 形 式 ）。 'e':表示“estimate error”。这种形式会在'r'的基础上添加一个由包含 eps
例题 8： taylor（f,n,v,a） 得到最高是 n-1 阶 求 8 阶展开 syms x r=taylor(x*exp(x),9,x,0) pretty(r) 例题 9： 级数求和 s=symsum(f,v,a,b)