matlab经典教案(土木工程专业).pptVIP

(iV). 多项式微分： polyder(p) polyfit(x,y,n) 其中x、y为拟合数据，n为拟和多项式的阶数。 例如： p=[1 –5 6 –1 9]; p4=poly2sym(p) Dp=polyder(p) P3=poly2sym(Dp) (V). 多项式拟合 例2.3:用5阶多项式对 上的正弦函数进行最小二乘拟合 x=0:pi/20:pi/2; y=sin(x); p=polyfit(x,y,5) x1=0:pi/30:2*pi; y1=sin(x1); y2=polyval(p,x1) plot(x1,y1,’-b’,x1,y2,’*r’) legend(‘原曲线‘，’拟合曲线‘) axis([-0.2,6.5,-1.2,1.2]) 演示5: poly01.m (5) 其它算例 执行字符串——用 eval 函数生成四阶的Hilbert函数 n=4; t=‘1/(i+j-1)’; a=zeros(n); for i=1:n for j=1:n a(i,j)=eval(t); end end (6). 上机练习题 1. 对于给定数据x=0.5,1,1.5,2,2.5,3,3.5,4,4.5,5,计算如下函数值 。 2. 给定矩阵 （1）.计算矩阵行列式的值 （2）.计算矩阵的逆 （3）.计算矩阵的秩、特方程 （4）.计算矩阵的特征值、特征向量 （5）.计算矩阵的条件数 （6）.对矩阵进行LU分解 3. 对于 解方程组Ax=b。 4. 对于离散数据 P=[1 2 3 2 4 –1 2],试拟合三次多项式， (1).求出该多项式的具体表达形式。 (2).求出该多项式的根。 (3). 求出该多项式在x=2.5处的值。 (4). 求出该多项式的微分。 5. 对于向量 a=[1 2 3 ], b=[2 4 –1],c=[2,1,3],试求 (1). a与b的向量积。 (2). a、b、c的混合积。 三、MATLAB的符号运算功能 1. 符号表达式的生成 (1). 创建符号函数 f=‘log(x)’ (2). 创建符号方程 equation=‘a*x^2+b*x+c=0’ (3). 创建符号微分方程 diffeq=‘Dy-y=x’ (4). 用sym命令来创建 f=sym(‘sin(x)’) f=sym(‘a*x^2+b*x+c=0’) (5). 用sym命令创建的另一种方式 syms x f=sin(x)+cos(x) 2. 符号与数值之间的转换 符号运算得到的是精确解,有时候要对得到的解析解 进行数值转换。主要有三个函数: digits, vpa, Subs： 变量替换函数. Digits(D): 函数设置有效数字个数为D的近似解精度. Vpa(s) : 符号表达式s在digits函数设置下的精度的 近似解. Vpa(s,D) : 符号表达式s在digits(D)精度下的 近似解. Subs(s,old,new): 将符号表达式中 old 的换成new变量. Numeric(s): 将不含自由变量的符号表达式转换成数 值形式. 例3. 1 符号运算 1. 解方程: 3*x-1=0 s=solve(‘3*x-1=0’) :解析解 digits(10) :设置10位有效数字 x=vpa(s) :10位有效数字解析解 or x=vpa(s,20) :20位有效数字解析解 numeric(s) :将解析解转换成数值解 2. 变量替换 x=sym(‘x’) f=x-cos(x)