<pre name="code" class="cpp">double x[3][3] = {{1, 50, -100}, {1, -110, 220}, {-1, 150, 150}};
double y[3][3] = {{1, 0, 0}, {0, 1, 0}, {0, 0, 1}};void showMatdate(Mat tmpMat){int i, j;CvScalar s1;Width = tmpMat.rows;Height = tmpMat.cols;IplImage tmp;tmp = tmpMat;for(i=0; i< Width; i++){for(j=0; j<Height; j++){s1 = cvGet2D(&tmp, i, j);printf("%0.1lf  ", s1.val[0]);}printf("\n");}printf("\n\n");
}int main(int argc, char *argv[]){/*************初始化矩阵*****************************/mat1 = Mat(3, 3, CV_64FC1, x);src1 = mat1;mat2 = Mat(3, 3, CV_64FC1, y);src2 = mat2;
/*************显示矩阵数据***************************/printf("mat1:\n");showMatdate(mat1);printf("mat2:\n");showMatdate(mat2);return 0;
}

<pre name="code" class="cpp">/*************两个矩阵绝对值之差*********************/mat3 = Mat(3, 3, CV_64FC1);src3 = mat3;cvAbsDiff(&src1, &src2, &src3);printf("cvAbsDiff(mat1, mat2):\n");showMatdate(mat3);

<pre name="code" class="cpp">/*************矩阵与数值差的绝对值*******************/mat3 = Mat(3, 3, CV_64FC1);src3 = mat3;cvAbsDiffS(&src1, &src3, Scalar(5));printf("cvAbsDiffS(mat1, 5):\n");showMatdate(mat3);

<pre name="code" class="cpp">/************两个矩阵相加***************************/mat3 = Mat(3, 3, CV_64FC1);src3 = mat3;cvAdd(&src1, &src2, &src3);printf("cvAdd(mat1, mat2):\n");showMatdate(mat3);

<pre name="code" class="cpp">/**********两个矩阵以权重形式相加******************/mat3 = Mat(3, 3, CV_64FC1);src3 = mat3;cvAddWeighted(&src1, 0.5, &src2, 0.5, 3, &src3);printf("cvAddWeighted(mat1, 0.5,mat2, 0.5, 3):\n");showMatdate(mat3);

<pre name="code" class="cpp">/*******************矩阵与计算**********************
*******可以使用掩码来决定计算被的数据元素*********
*******对应的还有cvNot cvOr cvXor *****************/mat3 = Mat(3, 3, CV_64FC1);src3 = mat3;cvAnd(&src1, &src2, &src3);printf("cvAnd(mat1, mat2):\n");showMatdate(mat3);

<pre name="code" class="cpp">/******************两个矩阵比较********************/int tmp = CMP_EQ;mat3 = Mat(3, 3, CV_8UC1);src3 = mat3;mat1.convertTo(mat1, CV_8UC1);mat2.convertTo(mat2, CV_8UC1);cvCmp(&src1, &src2, &src3, tmp);printf("cvCmp(mat1, mat2):\n");showMatdate(mat3);

<pre name="code" class="cpp">CMP_EQ src1 is equal to src2.
CMP_GT src1 is greater than src2.
CMP_GE src1 is greater than or equal to src2.
CMP_LT src1 is less than src2.
CMP_LE src1 is less than or equal to src2.
CMP_NE src1 is unequal to src2.

<pre name="code" class="cpp">/******************计算绝对值************************/mat3 = Mat(3, 3, CV_8UC1);src3 = mat3;cvConvertScaleAbs(&src1, &src3, 0.5, 4);printf("cvConvertScaleAbs(mat1, 0.5, 4):\n");showMatdate(mat3);注意:mat3 = | mat1 * 0.5 + 4 |;

<pre name="code" class="cpp">/****************计算矩阵非0元素的个数***************/int sum = 0;showMatdate(mat2);sum = cvCountNonZero(&src2);printf("sum=%d\n", sum);

<pre name="code" class="cpp">double x[3][3] = {{1, 50, -100}, {1, -110, 220}, {-1, 150, 150}};
double y[3][3] = {{1, 0, 0}, {0, 2, 0}, {0, 0, 3}};/*****************数组除法运算************************/mat3 = Mat(3, 3, CV_64FC1);src3 = mat3;cvDiv(&src1, &src2, &src3, 2);printf("cvDiv(mat1, mat2, 2)\n");showMatdate(mat3);注意：1、mat3 = mat1 * 2 / mat2。2、如果矩阵中mat2某一位为0，则mat3对应那一位的结果直接为0。

<pre name="code" class="cpp">/*****************矩阵行列式计算**********************/double sum = 0;sum = cvDet(&src1);showMatdate(mat1);printf("cvDet(mat1)=%.1lf\n", sum);

<pre name="code" class="cpp">/*****************矩阵元素指数计算**********************/mat3 = Mat(3, 3, CV_64FC1);src3 = mat3;cvExp(&src1, &src3);printf("cvExp(mat1)\n");showMatdate(mat3);计算公式如下：

<pre name="code" class="cpp">/*****************矩阵水平或者垂直翻转**********************/mat3 = Mat(3, 3, CV_64FC1);src3 = mat3;cvFlip(&src1, &src3, -1);printf("cvFlip(mat1, -1)\n");showMatdate(mat3);cvFlip(&src1, &src3, 0);printf("cvFlip(mat1, 0)\n");showMatdate(mat3);cvFlip(&src1, &src3, 1);printf("cvFlip(mat1, 1)\n");showMatdate(mat3);