juanjima.cpp

经常遇到的问题 计算方法里面的卷积码

#include<iostream.h>
#include<math.h>
void main()
{
	int i,j,s,t,a,b,c,m,n,k;
	int h[13],g[4][13],H[4][13],G[10][30],X[100],Z[100],q[10][5];
	for(i=0;i<13;i++)
	{
		for(j=0;j<=3;j++)
		{
			g[j][i]=0;
			H[j][i]=0;
		}
		for(j=0;j<10;j++)
			G[j][i]=0;
		h[i]=0;
	}
	cout<<"	****05计科2		20052840	胡浩	卷积马码****\n\n";
	cout<<"请确定卷积码的类型(n,k,m):";
	cin>>n>>k>>m;
	cout<<"请输入基本监督矩阵:\n";       //基本监督矩阵
	for(i=1;i<=12;i++)
		cin>>h[i];
	s=12;
	for(i=1;i<=3;i++)                  //基本生成矩阵
	{	
		s=s-4;
		g[i][i]=1;
		for(j=1;j<=3;j++)
			g[j][4*i]=h[s+j];
	}
	cout<<"基本生成矩阵为:"<<endl;
	for(i=1;i<=3;i++)
	{
		for(j=1;j<=12;j++)
			cout<<g[i][j]<<" ";
		cout<<endl;
	}
	for(i=1;i<=3;i++)                  //一致监督矩阵
		for(j=1,t=13-4*i;j<=4*i;j++,t++)
			H[i][j]=h[t];
	cout<<"一致监督矩阵为:"<<endl;
	for(i=1;i<=3;i++)
	{
		for(j=1;j<=12;j++)
			cout<<H[i][j]<<" ";
		cout<<endl;
	}
	for(i=1;i<=3;i++)                  //一致生成矩阵
		for(j=1;j<=12;j++)
			for(t=0;t<3;t++)
				G[i+3*t][j+4*t]=g[i][j];
	cout<<"一致生成矩阵为:"<<endl;
	for(i=1;i<=9;i++)	{
		for(j=1;j<=12;j++)
			cout<<G[i][j]<<" ";
		cout<<endl;
	}
	for(i=0;i<100;i++)	{
		X[i]=-1;
		Z[i]=0;
	}
	cout<<"请输入待编码序列起始码字:"<<endl;       //编码过程
	for(i=1;i<=9;i++)
		cin>>X[i];
	for(i=1;i<=12;i++)	{
		for(j=1;j<=9;j++)
			Z[i]=Z[i]+G[j][i]*X[j];
		Z[i]=Z[i]%2;
	}
	c=2;
	if(c==1)	{
		while(c==1)		{
			cout<<"请确定待编码序列由第二组开始的长度:";
			cin>>b;
			cout<<"请由第二组开始，依次输入码元:";
			for(i=4;i<100;i++)
				X[i]=0;
			for(i=4;i<4+b;i++)
				cin>>X[i];
			for(i=13;i<100;i++)
				Z[i]=0;
			for(i=0;i<=9;i++)
				for(j=0;j<=4;j++)
					q[i][j]=0;
				for(i=1;i<=9;i++)
					for(j=1;j<=4;j++)
						q[i][j]=G[i][j+8];
					a=1;
					do{
						for(i=1;i<=4;i++)			{
							for(j=1;j<=9;j++)
								Z[8+i+4*a]=Z[8+i+4*a]+q[j][i]*X[j+3*a];
							Z[8+i+4*a]=Z[8+i+4*a]%2;
						}
						a++;
					}while(a<=(b-6)/3);
				c=2;
		}
		cout<<"编码完成，结果为:"<<endl;
		for(i=1;i<=4*(b+3)/3;i++)
			cout<<Z[i]<<" ";
		cout<<endl;
	}
	else{
		cout<<"编码完成，结果为:"<<endl;
		for(i=1;i<=12;i++)
			cout<<Z[i]<<" ";
		cout<<endl;
	}
}