MATLAB 代数运算

若要得到一天24小时的温度曲线，则：
xi=0:1/3600:24;
yi=interp1(x,y,xi, ‘spline’);
plot(x,y,’o’ ,xi,yi)
MATLAB的初等代数运算
factor(s) 因式分解
expand(s) 展开符号表达式s
collect(s,x) 对表达式s中的每一个函数按x的次数合并系数。
collect(s,v) 按指定变量v的次数合并系数
[n,d]=numden(s) 符号表达式s的通分
simple(s) 显示多种方法化简，最后给出表达式s的最简型。
simplify(s) 用一般化简法化简符号表达式s。 radsimp(s)
对根式表达式s化简 horner(s) 符号表达式s的嵌套形式

syms x y %就是定了符号变量x y以后x y就可以直接使用了，有他们运算出来的结果也是符号变量当然上面的也可以x=sym(‘x’),y=sym(‘y’)sys(‘a+b’)就是将a+b转化为符号表达式
阅读下列程序：
1）分解因式1x3

syms x
factor(x^3-1) ans =
(x-1)*(x^2+x+1)
2）将cos(3x)展开 >> syms x
expand(cos(3*x)) ans =
4*cos(x)^3-3*cos(x)
3）将5
2x展开 >> expand((x+2)^5) ans =
x^5+10*x^4+40*x^3+80*x^2+80*x+32
4）将xxyxyx22
2按x的次数合并 >> syms x y
collect(x^2*y+y*x-x^2-2*x)

代数运算sym-
symadd(f,g)符号表达式f加上符号表达式g
symsub(f,g) 符号表达式f减去符号表达式g
symmul(f,g) 符号表达式f乘上符号表达式g
symdiv(f,g) 符号表达式f除符号表达式g
sympow(f,p) 符号表达式f的p次幂,p可以是表达式

复合函数的建立
compose(f,g) 生成函数f(g(x))，注意：此时与手写方式完全不同！！！
compose(f,g,z) 生成复合函数f(g(z)), 复合函数以z为自变量

syms x z
f=2^x;g=sin(x); >> compose(f,g) ans =
2^sin(x)
compose(g,f) ans =
sin(2^x)
compose(g,f,z) ans =
sin(2^z)

多项式运算poly-
多项式按降幂排列,写出系数向量,一定要把缺项的系数补0
用poly2sym命令建立多项式。
格式：poly2sym(p) % p为多项式的系数向量
如：输入多项式：32,652234xxsxxxp

p=[1 ,2, 0 ,-5, 6];s=[1, 2, 3]; >>p1= poly2sym(p) p1 =
x^4+2*x^3-5*x+6 >>s1= poly2sym(s)
s1 =
x^2+2*x+3
syms是生成符号变量的命令，其格式为：syms 变量名1 变量名2 变量名3…… 注：变量之间用空格分隔，不能用逗号。
如：多项式 7432xxp >> syms x;
q=3*x^2+4*x+7 q =
3*x^2+4*x+7(因为syms定义变量，结果还为变量)

命令格式：polyval (p,a) % 求 x=a 时多项式p的值
x=[a,b,c]; polyval(p,x) % 求 x=a ,x=b,x=c 时多项式p的值
如：polyval(p,-2) % 把x=-2代入多项式p求值
ans =
16
多项式乘法除法运算
conv(p,q) 多项式p乘以 q ,显示乘积的系数向量
deconv(p,q) 多项式p除以 q , 只显示商的系数向量
[t,r]=deconv(p,q) 多项式p除以q，商的系数向量
为t，余式的系数向量为r

多项式的微分运算 poly-
polyder(p)求p的一阶导数

多项式拟合 polyfit
多项式拟合是将给定的数据点采用最小二乘法按多项式形式进行拟合，确定多项式的系数。多项式拟合有命令语句和图形窗口两种方法