[高等教育]matlab2_24.pptVIP

MATLAB即Matrix Laboratory 矩阵实验室的意思 MATLAB的一个重要功能是能够进行向量和矩阵运算， MATLAB中的多数功能也是基于向量和矩阵。因此矩阵在MATLAB中具有非常重要的位置。在MATLAB中，向量和矩阵主要由数组来表示，数组是MATLAB的核心数据结构。 从数组结构的角度讲，向量是一个一维数组，而矩阵是一个二维数组。 MATLAB界面 熟悉，试用下面函数 sin，cos，tan，cot，sec，csc，asin，acos，atan，acot，asec，acsc，sinh，cosh，tanh sqrt （正的平方根），exp （e的指数），pow2（2的指数），log（自然对数，以e为底），log10，log2，abs（绝对值或复数的模） real（取实部），imag（取虚部），angle（取辐角），sign（符号函数） round（四舍五入取整），floor（向-inf取整），ceil（向+inf取整） ，fix（向0取整） 0.5^a ans = 1.7310 0.0910 -0.5489 -0.1606 0.7135 -0.4123 -1.0523 -0.6640 0.7242 例:a=[1 2 3;4 5 6;7 8 0]; p=poly(a) p =1.00 -6.00 -72.00 -27.00 p是多项式p(x)=x3-6x2-72x-27的matlab描述方法，我们可用： p1=poly2str(p,‘x’) — 函数文件，显示 数学多项式的形式 p1 =x^3 - 6 x^2 - 72 x - 27 2.roots —— 求多项式的根 a=[1 2 3;4 5 6;7 8 0];p=poly(a) p = 1.00 -6.00 -72.00 -27.00 r=roots(p) r = 12.12 -5.73 ——显然 r是矩阵a的特征值 -0.39 当然我们可用poly令其返回多项式形式 p2=poly(r) p2 = 1.00 -6.00 -72.00 -27.00 matlab规定多项式系数向量用行向量表示，一组根用列向量表示。 3.conv(卷积）多项式乘运算 例:a(x)=x2+2x+3; b(x)=4x2+5x+6; c = (x2+2x+3)(4x2+5x+6) a=[1 2 3];b=[4 5 6]; c=conv(a,b)=conv([1 2 3],[4 5 6]) c = 4.00 13.00 28.00 27.00 18.00 p=poly2str(c,x) p = 4 x^4 + 13 x^3 + 28 x^2 + 27 x + 18 4.deconv多项式除运算 a=[1 2 3]; c = [4.00 13.00 28.00 27.00 18.00] d=deconv(c,a) d =4.00 5.00 6.00 [d,r]=deconv(c,a) 余数 c除a后的整数 5. 多项式微分 matlab提供了polyder函数多项式的微分。 命令格式： polyder(p): 求p的微分 polyder(a,b): 求多项式a, b乘积的微分 [p,q]=polyder(a,b): 求多项式a, b商的微分 例：a=[1 2 3 4 5]; poly2str(a,x) ans = x^4 + 2 x^3 + 3 x^2 + 4 x + 5 b=polyder(a) b = 4 6 6 4 poly2str(b,x) ans =4 x^3 + 6 x^2 + 6 x + 4 四、代数方程组求解 matlab中有两种除运算左除和右除，避免矩阵求逆反操 作。 对于方程AX=B，A 为An×m矩阵，有三种情况： ? 当n=m时，此方程成为“恰定”方程，方阵，给出精确解。 ? 当nm时，此方程成为“超定”方程,MATLAB将给出最小二乘解。 ? 当nm时，此方程成为“欠定”方程，将给出一组基解。 matlab定义的除运算可以很方便地解上述三种方程 1.恰定方程组的解 方程AX=B(A为非奇异) X=A-1 B 矩阵逆 两种解: X=inv(A)?B — 采用求逆运算解方程 X=A\B — 采用左除运算解方程 方程ax=b a=[1 2;2 3];b=[8;13]; ?x=inv(a)*b