线性分组码

#include<stdio.h>
#define N 100
int G[N][N];
int M[N][N];
int GM[N][N];
void XIE();
void write();
void main()
{
char a,b;
int i=0,j=0,n=0,m=0,k=0,x=0,y=0,z=0;
FILE *read;
read=fopen("G.txt","r");
while(!feof(read))
{
a=fgetc(read);
m++; //m为G的列数
if(a=='/n')
{
n++; //n为G的行数
i++;
j=-1;
k=m;
m=0;
}
G[i][j]=a-'0';
j++;
}
i=0;j=0;
read=fopen("M.txt","r");
while(!feof(read))
{
b=fgetc(read);
x++; //x为M的列数
if(b=='/n')
{
y++; //y为M的行数
i++;
j=-1;
z=x;
x=0;
}
M[i][j]=b-'0';
j++;
}
XIE(G,m,n,M,x,y);
write(m,n,x,y);
}
void XIE(int G[N][N],int m,int n,int M[N][N],int x,int y)
{
int i,j,k,gm=0;
for(i=0;i<=y;i++)
{
for(k=0;k<m-1;k++)
{
for(j=0;j<x-1;j++)
{
gm+=G[j][k]*M[i][j];
}
GM[i][k]=gm%2;
gm=0;
}
}
}
void write(int m,int n,int x,int y)
{
FILE *write;
int i,j;
write=fopen("结果.txt","w+");
fprintf(write,"G矩阵为：/n");
for(i=0;i<=n;i++)
{
for(j=0;j<m-1;j++)
{
fprintf(write,"%d",G[i][j]);
}
fprintf(write,"/n");
}
fprintf(write,"M矩阵为：/n");
for(i=0;i<=y;i++)
{
for(j=0;j<x-1;j++)
{
fprintf(write,"%d",M[i][j]);
}
fprintf(write,"/n");
}
fprintf(write,"生成码字矩阵为：/n");
for(i=0;i<=y;i++)
{
for(j=0;j<m-1;j++)
{
fprintf(write,"%d",GM[i][j]);
}
fprintf(write,"/n");
}
printf("结果已写入"结果"文件中/n");

}