假设两个矩阵A、B,A是一个3x16的矩阵,B矩阵是一个3x3的矩阵,用B矩阵对 A矩阵做中值运算 ,得到一个1X14的矩阵C
中值运算就是用矩阵B作为掩膜,在矩阵A上滑动。每一次掩膜运算,就是矩阵A/B对应的、相同位置的数相乘,得到9个积。这九个积相加的和,再除以16,就得到一个C矩阵的对应掩膜起始位置的值。矩阵B每向右移动一列,就进行一次掩膜运算,在矩阵A上,共进行了14次掩膜运算。得到矩阵C的值。
原理:一次掩膜运算即9个数对应相乘,再相加。
本人思路:以矩阵A的列为基本单位,将一次掩膜运算所需要的9个数,对应相乘后,拆成3列,
每一列进行列加处理(一列的数全部相加),再将列加处理后的结果储存在数组c[n];最后,将数组c[n]中每三个的元素相邻作为一组,相加后除以16,再赋给数组C[t]
网上思路:按矩阵A的行处理
#include<stdio.h>
int main()
{
int A[3][16]={{16,16,16,8,18,8,15,15,15,20,22,23,17,17,17,17},
{15,15,15,7,17,8,15,15,15,25,23,21,16,16,16,16},
{15,15,15,8,18,8,15,15,15,27,25,24,16,17,17,17}};
int B[3][3]={{1,2,1},{2,4,2},{1,2,1}};
int i=0,n=-1,l=0;
float C[100];
int k,m,j,c[100],t;
for (m=1;m<15;m++)//m代表掩膜运算
{
for(k=0,j=m-1;k<3;k++,j++) //3*14=42,总共得到42次列加数据
{
n++;//第一次循环n=0
c[n]=A[i][j]*B[i][k]+A[i+1][j]*B[i+1][k]+A[i+2][j]*B[i+2][k];
printf("i=%d m=%d j=%d n=%d c[n]=%d\n",i,m,j,n,c[n]);
}
}
for(t=0;t<14;t++)
{
C[t]=(float)(c[l]+c[l+1]+c[l+2])/16;//加上强制类型转换float,使最终能够输出小数位部分的数据
l=l+3;
printf("t=%-d c[%d]=%7.4f\n",t,t,C[t]);
}
return 0;
}
下图是本人手动计算的最终正确答案,仅供参考对照。
看完觉得不错的话点个赞吧!
欢迎评论区交流!
注:该帖仅是提供思路和灵感!!
不足之处,请多指教。