矩阵转换，高手请进！！！

[bestajian]

16 * 16矩阵，（供占用32字节）。初始运算如下：
BYTE convert_matrix_dot(unsigned char *pbuf, unsigned char buf_len)
{
char i, j, m, n, count;
unsigned char buf[32];

count = buf_len/8;

memcpy(buf, pbuf, buf_len);

if (count != 4)
return 0;

for(i=0; i<2; i++)
{
for(j=0; j<2; j++)
{
m=i*16+j*8;
n=i*16+j;
pbuf[0+m]= ((buf[0+n]>>7)&0x01)|((buf[2+n]>>6)&0x02)|((buf[4+n]>>5)&0x04)|((buf[6+n]>>4)&0x08)|((buf[8+n]>>3)&0x10)|((buf[10+n]>>2)&0x20)|((buf[12+n]>>1)&0x40)|((buf[14+n]<<0)&0x80);
pbuf[1+m]= ((buf[0+n]>>6)&0x01)|((buf[2+n]>>5)&0x02)|((buf[4+n]>>4)&0x04)|((buf[6+n]>>3)&0x08)|((buf[8+n]>>2)&0x10)|((buf[10+n]>>1)&0x20)|((buf[12+n]>>0)&0x40)|((buf[14+n]<<1)&0x80);
pbuf[2+m]= ((buf[0+n]>>5)&0x01)|((buf[2+n]>>4)&0x02)|((buf[4+n]>>3)&0x04)|((buf[6+n]>>2)&0x08)|((buf[8+n]>>1)&0x10)|((buf[10+n]<<0)&0x20)|((buf[12+n]<<1)&0x40)|((buf[14+n]<<2)&0x80);
pbuf[3+m]= ((buf[0+n]>>4)&0x01)|((buf[2+n]>>3)&0x02)|((buf[4+n]>>2)&0x04)|((buf[6+n]>>1)&0x08)|((buf[8+n]>>0)&0x10)|((buf[10+n]<<1)&0x20)|((buf[12+n]<<2)&0x40)|((buf[14+n]<<3)&0x80);
pbuf[4+m]= ((buf[0+n]>>3)&0x01)|((buf[2+n]>>2)&0x02)|((buf[4+n]>>1)&0x04)|((buf[6+n]<<0)&0x08)|((buf[8+n]<<1)&0x10)|((buf[10+n]<<2)&0x20)|((buf[12+n]<<3)&0x40)|((buf[14+n]<<4)&0x80);
pbuf[5+m]= ((buf[0+n]>>2)&0x01)|((buf[2+n]>>1)&0x02)|((buf[4+n]<<0)&0x04)|((buf[6+n]<<1)&0x08)|((buf[8+n]<<2)&0x10)|((buf[10+n]<<3)&0x20)|((buf[12+n]<<4)&0x40)|((buf[14+n]<<5)&0x80);
pbuf[6+m]= ((buf[0+n]>>1)&0x01)|((buf[2+n]>>0)&0x02)|((buf[4+n]<<1)&0x04)|((buf[6+n]<<2)&0x08)|((buf[8+n]<<3)&0x10)|((buf[10+n]<<4)&0x20)|((buf[12+n]<<5)&0x40)|((buf[14+n]<<6)&0x80);
pbuf[7+m]= ((buf[0+n]>>0)&0x01)|((buf[2+n]<<1)&0x02)|((buf[4+n]<<2)&0x04)|((buf[6+n]<<3)&0x08)|((buf[8+n]<<4)&0x10)|((buf[10+n]<<5)&0x20)|((buf[12+n]<<6)&0x40)|((buf[14+n]<<7)&0x80);
}
}

return buf_len;
}

请问我怎样根据转换后的矩阵获取到转换前的矩阵呢？（亦即上述运算的逆运算算法！）

请各位大侠帮忙！

[你妹的特盗不]

贴代码要象这样

//这里面是你的代码
//你那们贴，读得太累
//帮你顶

[sea_spray]

vc可以调用matlab的，matlab可以实现这种操作，lz可以试一下在matlab中进行转换啊。

[bestajian]

BYTE convert_matrix_dot(unsigned char *pbuf, unsigned char buf_len)
{
	char i, j, m, n, count;
	unsigned char buf[32];
	
	count = buf_len/8;
	
	memcpy(buf, pbuf, buf_len);
	if (count != 4)
		return 0;
	
	for(i=0; i<2; i++)
	{
		for(j=0; j<2; j++)
		{
			m=i*16+j*8;
			n=i*16+j;
			pbuf[0+m]= ((buf[0+n]>>7)&0x01)|((buf[2+n]>>6)&0x02)|((buf[4+n]>>5)&0x04)|((buf[6+n]>>4)&0x08)|((buf[8+n]>>3)&0x10)|((buf[10+n]>>2)&0x20)|((buf[12+n]>>1)&0x40)|((buf[14+n]<<0)&0x80);
			pbuf[1+m]= ((buf[0+n]>>6)&0x01)|((buf[2+n]>>5)&0x02)|((buf[4+n]>>4)&0x04)|((buf[6+n]>>3)&0x08)|((buf[8+n]>>2)&0x10)|((buf[10+n]>>1)&0x20)|((buf[12+n]>>0)&0x40)|((buf[14+n]<<1)&0x80);
			pbuf[2+m]= ((buf[0+n]>>5)&0x01)|((buf[2+n]>>4)&0x02)|((buf[4+n]>>3)&0x04)|((buf[6+n]>>2)&0x08)|((buf[8+n]>>1)&0x10)|((buf[10+n]<<0)&0x20)|((buf[12+n]<<1)&0x40)|((buf[14+n]<<2)&0x80);
			pbuf[3+m]= ((buf[0+n]>>4)&0x01)|((buf[2+n]>>3)&0x02)|((buf[4+n]>>2)&0x04)|((buf[6+n]>>1)&0x08)|((buf[8+n]>>0)&0x10)|((buf[10+n]<<1)&0x20)|((buf[12+n]<<2)&0x40)|((buf[14+n]<<3)&0x80);
			pbuf[4+m]= ((buf[0+n]>>3)&0x01)|((buf[2+n]>>2)&0x02)|((buf[4+n]>>1)&0x04)|((buf[6+n]<<0)&0x08)|((buf[8+n]<<1)&0x10)|((buf[10+n]<<2)&0x20)|((buf[12+n]<<3)&0x40)|((buf[14+n]<<4)&0x80);
			pbuf[5+m]= ((buf[0+n]>>2)&0x01)|((buf[2+n]>>1)&0x02)|((buf[4+n]<<0)&0x04)|((buf[6+n]<<1)&0x08)|((buf[8+n]<<2)&0x10)|((buf[10+n]<<3)&0x20)|((buf[12+n]<<4)&0x40)|((buf[14+n]<<5)&0x80);
			pbuf[6+m]= ((buf[0+n]>>1)&0x01)|((buf[2+n]>>0)&0x02)|((buf[4+n]<<1)&0x04)|((buf[6+n]<<2)&0x08)|((buf[8+n]<<3)&0x10)|((buf[10+n]<<4)&0x20)|((buf[12+n]<<5)&0x40)|((buf[14+n]<<6)&0x80);
			pbuf[7+m]= ((buf[0+n]>>0)&0x01)|((buf[2+n]<<1)&0x02)|((buf[4+n]<<2)&0x04)|((buf[6+n]<<3)&0x08)|((buf[8+n]<<4)&0x10)|((buf[10+n]<<5)&0x20)|((buf[12+n]<<6)&0x40)|((buf[14+n]<<7)&0x80);
		}
	}
	
	return buf_len;
}

[bestajian]

谁指教一下啊

[SullenSun]

太多了，我也没看明白。。

[bestajian]

自己顶。