1、理论依据:
复数点乘:a=x1+y1i,b=x2+y2;a*b=(x1*x2-y1*y2)+(x1*y2+x2*y1)i
2、代码:
bool matSplit(Mat &src,Mat &dst_real, Mat &dst_img)
{
Mat planes[2];
split(src, planes);
dst_real = planes[0];//实部
dst_img = planes[1];//虚部
return true;
}
bool matDotMultiply(Mat src1, Mat src2, Mat dst)
{
//src1_real(x1),src1_imag(y1),src2_real(x2),src2_imag(y2)
Mat src1_real, src1_imag, src2_real, src2_imag;
Mat channel[2];
//将输入复数Mat数据拆分成虚部与实部
matSplit(src1, src1_real, src1_imag);
matSplit(src2, src2_real, src2_imag);
/*double x1, x2, x3, x4;*/
//计算src1.*src2
for (int i = 0; i<src1.rows; i++) {//用时125ms
for (int j = 0; j<src1.cols; j++)
{
//x1*x2-y1*y2
channel[0].at<double>(i, j) = src1_real.at<dou