IDL矩阵运算.docVIP

阵就是数组。 矩阵初始化：用中括号[]而不是圆括号() 矩阵元素取值：可用中括号或者圆括号 例： IDL A=[[1,2],[3,4]] IDL print,a 1 2 3 4 IDL print,a[0] 1 IDL print,a(1) 2 IDL print,a[1,1] 4 获取矩阵的第一行： print,A[*,0] IDL print,a[*,0] 1 2 获取矩阵的第二行： print,A[*,1] 获取矩阵的第一列： print,A[0,*] IDL print,a[0,*] 1 3 获取i行j列元素,A[j,i](注意列在前行在后) IDL PRINT,A[0,1] 3 A#B相当于 B乘A，而A*B没有什么意义。 如果不可乘，也会有一个结果，自动转换比如A=(1,2),B=(2,3,4) A#B= 2 4 3 6 4 8 B#A= 2 3 4 4 6 8 矩阵转置： TRANSPOSE(A) 或者ROTATE(A,1) 如果仅仅是打印，可以使用PM,A PM, array 就是 PRINT, TRANSPOSE(array) 例子： ; Define an array arr arr = [[1.0, 3.0], [0.0, 4.0], [2.0, 1.0]] ; Print using PM and PRINT PM, arr PRINT PRINT, arr IDL prints: 1.00000 0.000000 2.00000 3.00000 4.00000 1.00000 1.00000 3.00000 0.000000 4.00000 2.00000 1.00000 求逆矩阵： MATRIX_POWER(A,-1) 该函数同样可以求矩阵的其他次方 MATRIX_POWER(A,-2) 矩阵元素乘方： A^2 IDL print,a 1 2 3 IDL print,a^2 1 4 9 要注意： a^-1语法错误，要写成a^(-1) a^-1不是a的逆阵，而是内部元素的倒数构成的矩阵，但是有一个问题：都是整数，比如2,倒数不是0.5，而是0 解决办法是，A的其中一个元素必须是浮点数，比如 [1,2],[2,4] 可以写成[1,2.0],[2,4]这样就对了。 如何对不同的列/行增加不同的数值?比如A的第一列加2，第二列加1 IDL A=[[1,2,3],[4,5,6],[7,8,9]] IDL PRINT,A 1 2 3 4 5 6 7 8 9 IDL B=[TRANSPOSE(replicate(2,3)),TRANSPOSE(replicate(1,3)),TRANSPOSE(replicate(0,3))] IDL PRINT,B 2 1 0 2 1 0 2 1 0 IDL PRINT,A+B 3 3 3 6 6 6 9 9 9 更简便的方法，就是直接对某行某列进行运算 第一列加2：A(0,*)= A(0,*)+2 如何交换一行、列？ 交换第一列和第三列过程如下： IDL A=[[1,2,3],[4,5,6],[7,8,9]] IDL PRINT,A 1 2 3 4 5 6 7 8 9 IDL TEMP = A[0,*] IDL A[0,*]=A[2,*] IDL A[2,*]=TEMP IDL PRINT,A 3 2 1 6 5 4 9 8 7 行列求和 IDL A=[[1,2,3