基本应用2012.pptVIP

Matlab基本应用;什么是数值运算？什么是符号运算？;;;基本应用一：数值计算;1、多项式（polynominals）;2、多项式运算的函数 1）p与s 加，减，乘，除 2）求p (x)=0 的根 3）求p(a)的值 4）微分 d p(x) / d x ;;举例： x=-5:0.5:5; y1=polyval(p,x); y2=polyval(s,0) r=roots(s); ss=poly(r);;2、曲线拟合与插值运算;polyfit 函数：求得拟合多项式的系数 调用格式为：P=polyfit(X,Y,m) 例 采样点 x=[-3 -1 0 2 5.5 7]； y=[3.3 4.5 2.0 1.5 2.5 -1.2]; p3=polyfit(x,y,3) p5=polyfit(x,y,5) xx=-3:0.1:7; yy3=polyval (p3,xx); yy5=polyval (p5,xx); plot(xx,yy3,b,xx,yy5, ‘g, x,y,or);2、多项式插值 插值的定义——找到一个解析函数（多项式）连 接自变量相邻的两个点(xi,xi+1)，由此得到中间点的 近似值。;例3.16：分别用分段线性插值，三次样条插值方法 对余弦曲线进行插值。 x=0:10; y=cos(x); xi=0:0.25:10; y1=interp1(x,y,xi,linear); y2=interp1(x,y,xi,spline); plot(x,y,ro); hold on;pause plot(xi,y1,b); hold on;pause plot(xi,y2,r); hold on;pause plot(xi,cos(xi),:k); hold off;二维插值 Z1=interp2(X,Y,Z,X1,Y1,method) 由(X,Y,Z) 的值，据相应的插值方法得到在X1,Y1的 插值结果;例: 考虑一个平面范围内的温度变化，在多个离散点上测量出当地的温度，数据如下 1) x0= 1:0.1:6 2) x0=1:0.1:6 y0= 2 y0=1:0.06:4;x= [1 2 3 4 5 6]; y= [1 2 3 4]; t= [12 10 11 11 13 15;16 22 28 35 27 20;18 21 26 32 28 25;20 25 30 33 32 30]; x0= 1:0.1:6; y0=2.5; z0= interp2 (x,y,t,x0,y0) mesh(x,y,t); hold on;pause plot3(x0,y0*ones(1,length(x0)),z0);hold on;pause x1=1:0.1:6; y1=1:0.06:4; y2=y1; z1= interp2 (x,y,t,x1,y2); mesh(x1,y1,z1);hold off;; 3、数值微积分;x=linspace(-4,4,50); y=sin(x); s=(1-x.^3)./(1+x.^4); dfy=diff(y)./diff(x); dfs=diff(s)./diff(x);;2、数值积分 ;;4、求解线性代数方程组;1. 直接法 象解一元一次方程那样简单地求解： x=A\b 当系数矩阵A为N*N的方阵时，MATLAB会自行 用高斯消去法求解线性代数方程组。 2. 矩阵求逆法 等号两侧各左乘A-1，有 A-1Ax=A-1b 由于A-1A=I，故得： x=A-1b = inv (A)*b;3. 迭代法 迭代法是指用某种极限过程去逼近线性方程组的 精确解的过程。 基本思想是：先任取一组近似解初值X(0),然后按 照某种迭代函数计算出新的X(1),类似得到X(2),X(3), X(4),...X(k),...,当{X(k)}收敛时，有 X*即为方程组的解。;例 用数值法求解线性代数方程组Ax=b。;基本应用二：符号计算;解方程;§1 创建符号对象;2) syms 可以定义多个符号变量 调用格式为： syms 符号变量名1 ( 符号变量名2 … 符号变量名n);3) 符号矩阵的创建 A=[1,2;3,4] — 数值矩阵 A=sym(‘[