含答案《matlab实用教程》.doc

MATLAB实用教程3版版权所有翻版必究0第二章MATLAB语言及应用实验项目实验一MATLAB数值计算三、实验内容与步骤1创建矩阵987654321A(1)直接输入。(2)用FROMSTEPTO方式。(3)用LINSPACE函数。(4)使用特殊矩阵函数，并修改元素。ZEROS、ONES、RANDN、EYE。2矩阵的运算(1)利用矩阵除法解线性方程组。12247358431421XX将方程表示为AXB，计算XA\B。(2)利用矩阵的基本运算求解矩阵方程。已知矩阵A和B满足关系式A1BA6ABA，计算矩阵B。其中，7/104APSATRATA1/30001/40001/7BINVAINVINVAEYE36A(3)计算矩阵的特征值和特征向量。已知矩阵，计算其特征值和特征向量。10452X(4)PAGE322利用数学函数进行矩阵运算。已知传递函数GS1/2S1,计算幅频特性LW20LGX1202514101V,DEIGX验证特征值和牲向量与该矩阵的关系XVVDA123456789A134679ALINSPACE1,3,3LINSPACE4,6,3LINSPACE7,9,3ZEROS3,3A19147258369ONES3,3A19147258369RANDN3,3A19147258369EYE3,3A19147258369A2302152131114122B82712XA\BMATLAB实用教程3版版权所有翻版必究1和相频特性FWARCTAN2W,W的范围为001，10，按对数均匀分布。12W3多项式的运算(1)多项式的运算。已知表达式GXX4X5X26X9,展开多项式形式，并计算当X在0，20内变化时GX的值，计算出GX0的根。PAGE324(2)多项式的拟合与插值。将多项式GXX45X317X2129X180，当X在0，20多项式的值上下加上随机数的偏差构成Y1，对Y1进行拟合。对GX和Y1分别进行插值，计算在55处的值。PAGE325四、思考练习题1使用LOGSPACE函数创建0～4Π的行向量，有20个元素，查看其元素分布情况。PSLOGSPACELOG100,LOG104PI,202SORTC,2顺序排列31多项式1)FX2X23X5X8用向量表示该多项式，并计算F10值2)根据多项式的根0534I34I创建多项式。R05,34I,34IPPOLYRPOLY2STRP, X 3)把表达式(X22X2)X25X4展开为多项式形式，并求对应方程的根。PSA1122A2154ACONVA1,A24)求多项式3X612X54X47X38X1除以(X3)X35X的商和余式。32曲线拟合有一组测量数据如下表所示X1152253354455Y142735984122166188262使用曲线拟合，求出此数据对应的2阶和10阶多项式。33曲线插值假设某日气温的实际测量值为时间1357911温度141415182225时间131517192123温度313127242017使用三次样条插值方法，求出10点、12点的温度。X1357911131517192123Y141415182225313127242017PPOLYFITX,Y,5XILINSPACE0,25,100ZPOLYVALP,XI多项式求值PLOTX,Y, O ,XI,Z, K ,X,YHOLDONTINTERP1X,Y,12, SPLINE PLOT12,T, R TINTERP1X,Y,10, SPLINE PLOT10,T, R P131247081A113A21050ACONVA1,A2Q,RDECONVP1,AX1152253354455Y142735984122166188262POLYFITX,Y,2POLYFITX,Y,10F288POLYVALF,10MATLAB实用教程3版版权所有翻版必究2实验二MATLAB符号计算三、实验内容与步骤1符号表达式的代数运算与化简对符号表达式FX23X2和GX31进行运算。(1)符号表达式的代数运算。FSYM X23X2 GSYM X31 (2)符号表达式的化简。EXPAND、COLLECT、SIMPLIFY。EXPANDFCOLLECTFSIMLIFYFG函数同理(3)符号表达式与多项式的转换。SYM2PLOY、PLOY2SYM。PAGE952符号矩阵的操作2COSTX(1)创建符号矩阵。(2)符号矩阵的代数运算。(3)符号矩阵的微积分。(4)求特征值。(5)求积分。3符号方程的求解(1)用代数方程求解。对方程组进行求解。12247358431421XX(2)用符号微分方程求解。解方程组的通解。当Y01,Z05时，求特解。1COSYDXZ四、思考练习题1分解因式。4YX2化简表达式。2121SINCOSIN3用符号方法求下列极限、导数或积分。(1)1ARCSLIM1XXA2302152131114122B82712XA\BY,ZDSOLVE DYZCOSX,DZY1 , X 求通解Y,ZDSOLVE DYZCOSX,DZY1 , Y01,Z05 , X 求特解FSYM X4Y4 FACTORFFSYM SQRTPISQRTARCCOSX/SQRTX1 LIMITF, X ,1, RIGHT ASYM AXT3TCOSXLOGX DADTDIFFA, T D2ADT2DIFFA, T ,2DADXDIFFA, X D2ADXDTDIFFDADX, T MATLAB实用教程3版版权所有翻版必究3(2)已知，分别求、、。XTTAAXLNCOS3DTA2DXT(3)已知，求、。XYEYXF2,21,0YXF(4)D0421(5)已知,1求2求。2YXF2DF1FDX4已知方程63D(1)求通解(2)求满足的特解。30,2,10YYSYMSXYFX22XEXP1X2Y2XYPYPXDIFFY, X PFPXDIFFF, X P2FPXPYDIFFPFPX, Y X0,Y1F1SUBSP2FPXPYSYMSXFX21/X41INTF, X ,0,INFSYMSXYFX2Y2D2FDY2DIFFF, Y ,2INTF, X ,1,1YDSOLVE D3Y3DY2Y6 , X 求通解Y1DIFFY, X 对Y求一阶导数Y2DIFFY1, X 对Y求二阶导数YDSOLVE D3Y3DY2Y6 , Y03,Y102,Y201 , X 求特解MATLAB实用教程3版版权所有翻版必究4实验三MATLAB的计算可视化三、综合实例图31双调谐滤波器结构图图32单调谐滤波器结构图供电系统中，常存在谐波电流，所以设计电力系统滤波器对谐波电流进行滤除，电力系统无源滤波器中双调谐滤波器的结构如图31所示，其阻抗表达式为11321211/NZRJLRJRJNLNCN式中，为基波角频率。当设，，，，，，110200H014285H，，在08，20取值时，双调谐滤波器的阻抗与谐波次数的曲线如图33所示，02CF04285FNZ编写程序完成曲