如果您无法下载资料,请参考说明:
1、部分资料下载需要金币,请确保您的账户上有足够的金币
2、已购买过的文档,再次下载不重复扣费
3、资料包下载后请先用软件解压,在使用对应软件打开
实验二符号运算与曲线拟合1.符号矩阵的运算1.符号矩阵的四则运算运算符为:加(+),减(-)、乘(×)、除(/、\)等或:符号矩阵的和(symadd),差(symsub),乘(symmul)。例1:采用字符串直接创建矩阵的方法创建矩阵a:例2:分别将下列数值矩阵与符号矩阵互换例3;C=B-AD=a\b则显示:C=x-1/x1-1/(x+1)x+2-1/(x+2)-1/(x+3)D=-6*x-2*x^3-7*x^21/2*x^3+x+3/2*x^26+2*x^3+10*x^2+14*x-2*x^2-3/2*x-1/2*x^32.其他基本运算符号矩阵的其他一些基本运算包括转置(')、行列式(det)、逆(inv)、秩(rank)、幂(^)和指数(exp和expm)等都与数值矩阵相同3.符号矩阵的简化符号工具箱中提供了符号矩阵因式分解、展开、合并、简化及通分等符号操作函数。(1)因式分解命令:factor符号表达式因式分解函数格式:factor(s)说明:s为符号矩阵或符号表达式。常用于多项式的因式分解例4将x9-1分解因式在Matlab命令窗口建入symsxfactor(x^9-1)则显示:ans=(x-1)*(x^2+x+1)*(x6+x^3+1)(2)符号矩阵的展开命令expand符号表达式展开函数格式:expand(s)说明:s为符号矩阵或表达式。常用在多项式的因式分解中,也常用于三角函数,指数函数和对数函数的展开中例5将(x+1)3、sin(x+y)展开symsxyp=expand((x+1)^3)q=expand(sin(x+y))则结果显示为p=x^3+3*x^2+3*x+1q=sin(x)*cos(y)+cos(x)*sin(y)例6解线性方程组5*x^2+2*y+5*z=1、x-y^3+z=10、2*x-5*y-3*z=52.曲线拟合例7:多项式拟合x0=0:0.1:1;y0=[-.4471.9783.115.255.024.664.014.583.455.359.22];p=polyfit(x0,y0,3)p=56.6915-87.117440.0070-0.9043xx=0:0.01:1;yy=polyval(p,xx);plot(xx,yy,'-b',x0,y0,'or')上机练习(二):练习符号计算与曲线拟合的相关命令,将前面的例子在命令窗口中执行通过;创建符号表达式f=sin(x)+cos(x)-tan(x);使用sym函数直接生成3行3列的符号矩阵;随机生成一数值矩阵,使用sym函数将数值矩阵转换为符号矩阵;已知:,,请对符号矩阵x、y进行四则运算设x为符号变量,f=x^4+2x^2+1,g=x^3+6x^2+3x+5,试进行如下运算。(1)f+g;(2)f*g;(3)对f进行因式分解。假设现在有下表所列的实验数据,试选择不同阶数对数据进行拟合,并画出拟合曲线;X12345678910Y251218232629343856