2016第3章MATLAB数值计算1.pptVIP

例:用?一个8次多?项式在[0?,3]逼近?sin(x?) ?x=0:0?.1:3;? y=si?n(x);? p=po?lyfit?(x,y,?8); y?1=pol?yval(?p,x);? plot?(x,y1?,’k*’?,x,y,?’k-’)?; x?1=0:0?.1:6;? y1=s?in(x1?); y2?=poly?val(p?,x1);? plot?(x,y1?,’k*’?,x,y,?’k-’)?; 3.3?.3数值插?值 一维数值?插值: 格?式:Y1=?inter?p1(X,?Y,X1,?’meth?od’);? meth?od插值方?式:lin?ear(线?性插值),?neare?st(最近?插值),s?pline?(三次样条?插值),c?ubic(?三次多项式?插值)默认?插值是li?near(?线性插值)?. 例?: X=0:?0.2:p?i; y1?=inte?rp1(X?,Y,pi?/2); ?y2=in?terp1?(X,Y,?pi/2,?near?est)?; y3=?inter?p1(X,?Y,pi/?2,‘sp?line?); y4?=inte?rp1(X?,Y,pi?/2,‘c?ubic?); y1? = ? ? y2 =? ? y3=? ? y4=? 0.9?975 ? 0.9?996 ? ? 1 ? ? ?1 3.3?.4数值微?分 差?分函数 d?x=dif?f(X);? dxi=?X(i+1?)-X(i?) dx=?diff(?X,n);? diff?(X,2)?=diff?(diff?(X));? 3.3?.5数值积?分 积?分函数: ?quad(?f,a,b?,tol,?trace?); f在?[a,b]?的积分,t?ol计算精?度,默认为?0.001?,trac?e非0时画?出积分图形? 例:? f=in?line(?exp(?-x.^2?)); ? ? qua?d(f,0?,1); ?quad(?f,a,b?,c,d,?tol,t?race)?;函数f(?x,y)在?[a,b]?*[c,d?]上的积分?. 3?.3.6多?项式及其运?算 多?项式的建立? ?按降幂用向?量表示.如?: ?p(x)=?2x4-3?x-10 ? p=?[2 0 ?0 -3 ?-10];? 多项式的?求根 ? 命令:?x=roo?ts(p)?; x =? 1?.6538? ? ? -0?.1677? + 1.?5050i? -0?.1677? - 1.?5050i? -1?.3185? 多项?式求值 ? 函数?格式: p?olyva?l(p,x?);x为标?量或矩阵 ?例: x=?[1,2;?3,4];? poly?val(p?,x); ?ans =? -?11 ? 16 ? 143? 49?0 多?项式四则运?算 多?项式的加减?法 ?例: p1?(x)=5?x2-2x?+20, ?p2(x)?=x3-3?x2+x+?13 ? ? 则p1?=[0 5? -2 2?0]; p?2=[1 ?-3 1 ?13]; ? ? p?=p1+p?2 ? ? p= ? ? ? 1 ?2 -1 ?33 多项?式的乘除法? 函数?: con?v (乘)?,deco?nv(除)?. ?多项式的导?数 例:? a=[4?,1,2,?3]; ? ? p=?polyd?er(a)?;计算多项?式a的导数? ? ?p = ? ? ? 12 ? 2? ?2 ? p=? poly?der(a?,b); ?求a*b的?导数 ? [?p,q]=?polyd?er(a,?b);求a?/b的导数? 第三章? MATL?AB数值计?算 3.1矩?阵运算 3?.2稀疏矩?阵 3.3? 数据处理?与多项式运?算 3?.1矩阵运?算 3.1?.1矩阵变?换 3.1?.2矩阵的?逆与秩 3?.1.3矩?阵的特征值?与特征向量? 3.1?.1矩阵变?换 MATL?AB提供了?几个矩阵变?换函数及调?用方式: ?B=fli?pud(A?);A对于?水平轴沿上?下方向进行?列维翻转 ?B=fli?plr(A?); A对?于垂直轴沿?左右方向进?行列维翻转? B=fl?ipdim?(A,di?m); A?沿着特定维?进行翻转.? B=ro?t90(A?);矩阵A?逆时针方向?旋转90度? B=ro?t90(A?,k);矩?阵A逆时针