MAtlab与科学计算(word版).docVIP

1.定义二维数组：TpTs=[ 13.0249 20.0511; 15.1966 23.2570; 17.0544 26.4523; 19.2718 30.1574; 21.0893 33.2776; 23.1923 36.4374; 25.0391 39.7976; 27.4119 43.2971; 29.0323 46.3749]; 2.将二维数组分配给X/Y: y=[TpTs(:,1)]; X=? 3.求取平均值x_ave=mean(x); 4.平方取和：C=sum((Z-Z_ave).^2); 5.向量的点积、差积、混合积：dot(X,Y),cross(X,Y),dot(X,cross(Y,Z)) 6.（1）分号；……区分行以及取消运行显示等。例： A=[1,2;3,4]与A=[1,2;3,4]；的区别。 (2) 逗号，……区分列及函数参数分隔符等。例：=[1,2;3,4] ，B=[1,4,3;3,2,1;4,5,6] 7. 续行号…：例： y=sin(x)/(2+cos(x)) 也可写为 y=sin(x)… /(2+cos(x)) 8. 常量：MATLAB中有些预定义的变量，这些特殊的变量称为常量。常用到的有： i,j: 虚数单位； pi: π； NaN: 表示不定值,比如0/0； inf: 无穷大（infinit），比如1/0。 9. 算术操作符： +、-：加，减；可以通用。 *，^,\ , / .*, .^ , .\ , ./ ：分别为数组乘，乘方，左除，右除；此时向量的运算不会满足矩阵的运算法则。 一．矩阵及其运算 注意矩阵的加点运算结果。 miu0=1.785e-3; a=0.03368; b=0.000221; t=0:20:80 miu=miu0./(1+a*t+b*t.^2) 运行后的结果为 ： miu = 0.0018 0.0010 0.0007 0.0005 0.0003 10.miu0=1.785e-3;即是定义miu0= 11. 字符串： 字符串的约定（1）字符串用单引号括起来；（2）字符串的每个每个字符（包括空格）都是字符数组的一个元素. 例 s=‘xingtai college’ 　　f=‘sin(x)’　 是字符串（char array） 12. 向量的生成： （1）直接输入：如a=[1,2,5,3] (2) 利用冒号表达式生成：如：b=[2:2:10],此时[ ]可省略，步长为1时，步长可省略。第一个数为首元素的值，第2个数为步长或差值，第三个数为尾元素的限值，不能超过这个值。如b=2:2:11等价于b=[2:2:10] （3）线性等份向量生成：y=linspace(x1,x2,n),生成n维向量，使得y(1)=x1,y(n)=x2。如：y=linspace(1,100,6)。 13.解一元矩阵方程： （1）如:A=[1,2,3,2;3,2, 4,1;3,1,5,6;2,5,3,2], b=[1;3;2;1],求方程组Ax=b的解.由于rank(A)=rank(B)=4(B为增广矩阵),所以有唯一解,x=A\b,或x=inv(A)*b. （2）A=[361,625,961,1444,1936;1,1,1,1,1], b=[1;1]‘,求方程组Ax=b的解.由于rank(A)=rank(B)=2(B为增广矩阵),所以有无穷多个解,MATLAB中用除法解方程组时所得到的解是所有解中范数最小的一个x=A\b。 14.矩阵定义与运算及特殊矩阵的生成、矩阵特殊操作： 矩阵求逆：inv(A)为A的逆(inverse). 求转置矩阵:A. 求矩阵的行列式:det(A) ,(determinant是行列式)。 矩阵幂运算：用^.如A^3,表示A*A*A。 矩阵指数运算：expm(A),A为方阵。 矩阵对数运算:logm(A),A为方阵。 矩阵开方：sqrtm(a). 求矩阵呢的秩：rank(a). 特殊矩阵的生成： （1）zeros(n):生成 n×n阶0矩阵。 （2）zeros(m,n):生成 m×n阶0矩阵。 （3）zeros(size(a)):生成与a阶数相同的0矩阵。 （4） eye(n):生成 n阶单位矩阵。 （5） eye(m,n):生成 m×n阶单位矩阵。 （6） eye(size(a)):生成与a阶数相同的单位矩阵。 （7） ones(n):生成 n阶全1矩阵。 （8） ones