数值计算功能>>258*369ans=95202例：计算水在温度为0、20、40、60和80度时的黏度。已知水的黏度随温度的变化公式为：其中为零度水的黏度，值为1.785*10-3(P22)几个基本的算术运算四则运算“+”、“-”、“*”“/”“^”“sqrt”常驻变量PI()I或j(代表复数单位)inf代表无穷大数1/0NaN代表不定值0/0例：SIN(PI/3)J^2二、向量1.Ｘ=[5，7，8，9]Ｘ=[5，7;8,9]X(2)Ans=72.利用冒号表达式生成向量x=x0:step:xn例a=1:2:12a=13579113.利用函数linspace来生成线性等向量y=linspace(x1,x2)生成（1*100）维的行向量使得y(1)=x1,y(100)=x2y=linspace(x1,x2,n)生成（1*n）维的行向量使得y(1)=x1,y(100)=x24.向量的基本运算加、减、数乘、点积、叉积、混合积a1=linspace(1,100,6)a1-1a1*2a=[1,2,3]b=[34,5]dot(a,b)&&&点积ans=26C=CROSS(a,b)&&&叉积Dot(a,cross(b,c))&&&混合积三、矩阵直接输入小矩阵(1)a=[123;1,11;4,5,6]a=123111456(2)b=[sin(pi/3),cos(pi/4);log(9),tanh(6)]b=0.86600.70712.19721.0000b=[sin(pi/3),cos(pi/4);log(9),tanh(6)];屏幕上没有显示，但已在内存创建M文件输入大矩阵%example.mexm=[456468873257955216785448881365456588335567457687678956589456510953377]example;>>size(exm)ans=6矩阵的基本数学运算a=[123;235；345]b=[111;222;333]c=a+bc=234457678(2)e=[b,[555]’]&&&b为上例中bf=a*ef=141414302323235026262660(3)左除(\)和右除(/)右除比左除计算的速度慢，且精度还要差。但只是在条件数很大时才表现得明显，对于一般矩阵，几乎无差别例若线性方程的标准型为Az=,则A=(x2,1)为矩阵的系数，z=(a,b)’是未知数，=y解：x=[1925313844]’y=[1932.34973.397.8]’a=[x.^2,ones(5,1)]a=3611625196111444119361b=y;ab=a\b&&&相当于Pinv(a)*b,,,,,,,,Ax=yx=A’*yPinv求逆矩阵inv(a)a是方阵ab=0.05009726x1=19:0.1:44;y1=ab(2)+x1.^2*ab(1);plot(x,y,’o’);hold;plot(x1,y1)&&作图阵的逆运算（inv）a=[21–3–1;3107;-124–2;10–15];inv(a)ans=-0.04710.5882-0.2706-0.94120.3882-0.35290.48240.7647-0.22350.2941-0.0353-0.4706-0.0353-0.05880.04710.2941（5）.矩阵的行列式运算a1=det(a);a2=det(inv(a));a1*a2and=1（6）矩阵的幂运算^（7）矩阵的指数运算(expmexpm1expm2expm3)分别由padeTaylor和特征值法计算矩阵指数(expm同expm1)例：计算3阶魔方矩阵的指数，并比较不同函数的结果b=magic(3)b=816357492>>expm(b)ans=1.0e+006*1.08981.08961.08971.08961.08971.08971.08961.08971.0897矩阵的对数运算(logm(b))矩阵的开方运算（sqrtm）矩阵的函数运算（最为实用的部分）特征值（eigeigs）eig(a)值为特征根例A=[73–2;34–1;-2–13];[x,y]=eig(A)>>[x,y]=eig(A)x=&&&x为特征向量0.5774-0.0988-0.8105-0.57740.6525-0.49080.57740.75130.3197y=&&以对角线为Y