• 设为首页
  • 点击收藏
  • 手机版
    手机扫一扫访问
    迪恩网络手机版
  • 关注官方公众号
    微信扫一扫关注
    迪恩网络公众号

matlab符号运算的一些笔记

原作者: [db:作者] 来自: [db:来源] 收藏 邀请

--求一元二次方程ax^2+bx+c=0的根
solve(\'a*x^2+b*x+c=0\')
--求f(x)=(cosx)^2的一次导数
x=sym(\'x\'); x定义为符号变量
diff(cos(s)^2)
--计算f(x)=x^2在区间[a,b上的定积分
syms a b x;
int(x^2,a,b)
符号对象的建立 sym syms
S = sym(A)将非符号对象(如,数字,表达式,变量等)A转换为符号对象,并存储在符号变量S中。
sym函数用来建立单个符号变量 一般调用格式为
符号变量=sym(A)
A可以为一个数或数值矩阵 也可以是字符串
例如
a=sym(\'a\')
--符号运算操作
X=sym(\'[x11,x12,x13;x21,x22,x23]\');
Y=sym(\'[y11,y12,y13;y21,y22,y23]\');
Z=X*Y 这里我测试的报错了
--查找符号变量
findsym(expr) 按字符顺序列出符号表达式espr中的所有符号变量
--符号表达式的替换
用给定的数据替换符号表达式中指定的符号变量
subs(f,x,a)
用a 代替字符函数f中的字符变量x
a可以是数/数值变量/表达式或字符变量/表达式
--六类常见的符号运算
1.因式分解
factor(f)
f=sym(\'2*w-3*y+z^2+5*a\')
factor(f)
factor也可以用于正整数的分解
s=factor(100)
2.函数展开
expand(f)
多项式展开
syms x; f=(x+1)^6;
expand(f)
三角函数展开
syms x y; f=sin(x+y);
expand(f)
3.合并同类项
collect(f,v) 按指定变量v进行合并
collect(f) 按默认变量进行合并
syms x,y;
f=x^2*y+y*x-x^2+2*x;
collect(f);
collect(f,y)
4.函数简化 y=simplify(f); 对f尝试多种不同的算法进行简化 返回其中最简短的形式
[How,y]=simplify(f) y为f的最简形式 How为简化方法 可以为因式分解 展开 合并同类项等
y=simpify(f)对f进行简化
syms x;f=sin(x)^2+cos(x)^2;
simplify(f)
5.分式通分
[N,D]=numden(f) N为通分后的分子 D为通分后的分母
syms x,y;
f=x/y+y/x;
[n,d]=numden(f)
6.horner多项式 
syms x;
f=x^4+2*x^3+4*x^2+x+1;
g=horner(f) g =x*(x*(x*(x + 2) + 4) + 1) + 1
7.求导 diff
g=diff(f,v) 求f关于v的导数
g=diff(f);求f关于默认变量的导数
g=diff(f,v,n)f关于v的n阶导数
9.计算积分
int(f,v,a,b) 计算f(v)dv在[a,b]的积分
int(f,a,b)计算默认变量的定积分
int(f,v)f(v)的不定积分
10.符号求和
symsum(f,v,a,b)求和f(v) v从a到b
symsum(f,a,b)关于默认变量求和
11.代数方程求解solve(f,v)
12.微分方程求解 dsolve
y=dsolve(\'eq1\',\'eq2\',...,\'cond1\',\'cond2\',...,\'v\')
y为输出解 eq1 eq2 ...为微分方程
cond1 cond2 ...为初值条件 v为自变量

求微分方程dy/dx+2xy=xe^(-x^2)
y=dsolve(\'Dy+2*x*y=x*exp(-x^2)\',\'x\')
例 求微分方程xy\'+y-e^x=0满足初值条件 y(1)=2e的特解 并画出解函数的图像
y=dsolve(\'x*Dy+y-exp(x)=0\',\'y(1)=2*exp(1)\',\'x\')
ezplot(y) % plot是绘制二维图形,并且是x,y的表达式是已知的或者是形如 y=f(x)这样确切的表达式,而ezplot是画出隐函数图形,
13.反函数finverser(f,v)f对于v的反函数


鲜花

握手

雷人

路过

鸡蛋
该文章已有0人参与评论

请发表评论

全部评论

专题导读
上一篇:
MATLAB中的符号运算发布时间:2022-07-18
下一篇:
Delphi操作XML(三)发布时间:2022-07-18
热门推荐
阅读排行榜

扫描微信二维码

查看手机版网站

随时了解更新最新资讯

139-2527-9053

在线客服(服务时间 9:00~18:00)

在线QQ客服
地址:深圳市南山区西丽大学城创智工业园
电邮:jeky_zhao#qq.com
移动电话:139-2527-9053

Powered by 互联科技 X3.4© 2001-2213 极客世界.|Sitemap