《卷积码编码程序的c语言实现》由会员分享,可在线阅读,更多相关《卷积码编码程序的c语言实现(6页珍藏版)》请在金锄头文库上搜索。
1、卷积码编码程序的c 语言实现#includeusing namespace std;int table18=1,2,4,8,16,32,64,128;int myn=0;int stalen=0;int stan02562=0;/输入 0 时个状态的输出int stan12562=0;/输入 1 时各状态的输出int stachn2562=0;/状态装换表int myg110=0;int myg210=0;int myout100; /int myoutsym=0;void chartobits(char ch,int *bits);char bitstochar(int *bits);voi
2、d convolution(void);void creatsta(void);void myinput(void);int main()char exit_char;myinput();creatsta();convolution();cinexit_char;void myinput(void)int i,j;cout 输入编码的约束长度N:(3N9)myn;stalen=int(pow(2 、0,myn-1);cout 选择默认的编码矢量则输入1,输入 2 则可输入其她的编码矢量i;if(i=1)switch(myn)case 3:myg10=1,myg11=1,myg12=1;myg2
3、0=1,myg21=0,myg22=1;break;case 4:myg10=1,myg11=1,myg12=1,myg13=1;myg20=1,myg21=0,myg22=1,myg23=1;break;case 5:myg10=1,myg11=0,myg12=1,myg13=1,myg14=1;myg20=1,myg21=1,myg22=0,myg23=1,myg24=1;break;case 6:myg10=1,myg11=0,myg12=1,myg13=1,myg14=1,myg15=1;myg20=1,myg21=1,myg22=0,myg23=1,myg24=0,myg25=1;卷
4、积码编码程序的c 语言实现break;case7:myg10=1,myg11=0,myg12=0,myg13=1,myg14=1,myg15=1,myg16=1;myg20=1,myg21=1,myg22=0,myg23=1,myg24=1,myg25=0,myg26=1;break;case8:myg10=1,myg11=0,myg12=0,myg13=1,myg14=1,myg15=1,myg16=1,myg17=1;myg20=1,myg21=1,myg22=1,myg23=0,myg24=0,myg25=1,myg26=0,myg27=1;break;case9:myg10=1,myg
5、11=1,myg12=0,myg13=1,myg14=0,myg15=1,myg16=1,myg17=1,myg18=1;myg20=1,myg21=0,myg22=0,myg23=0,myg24=1,myg25=1,myg26=1,myg27=0,m yg28=1;break;elsecout 输入 g1endl;for(j=0;jmyg1j;cout 输入 g2endl;for(j=0;jmyg2j;cout 连接矢量1 为 endl;for(j=0;jmyn;j+)coutmyg1j;coutendl;cout 连接矢量2 为 endl;for(j=0;jmyn;j+)coutmyg2j
6、;coutendl;coutstalen:stalen;coutendl;void creatsta(void)int i,j,k,myi,myj;int tembits10;for(i=0;istalen;i+)stan1i0=0;卷积码编码程序的c 语言实现stan1i1=0;stan0i0=0;stan0i1=0;stachni0=i/ 2;myi=i;for(j=0;j=pow(2 、0,myn-1-j)tembitsj=1;myi=myi-pow(2 、 0,myn-1-j);elsetembitsj=0;/*for(k=0;kmyn;k+)couttembitsk;coutendl
7、;*/for(k=0;kmyn;k+)stan0i0+=myg1k*tembitsk;stan0i1+=myg2k*tembitsk;stan0i0=stan0i0%2;stan0i1=stan0i1%2;myi=i+int(pow(2 、 0,myn-1);stachni1=myi/ 2;for(j=0;j=pow(2 、0,myn-1-j)tembitsj=1;myi=myi-pow(2 、 0,myn-1-j);elsetembitsj=0;/*for(k=0;kmyn;k+)couttembitsk;coutendl;*/for(k=0;kmyn;k+)stan1i0+=myg1k*t
8、embitsk;stan1i1+=myg2k*tembitsk;stan1i0=stan1i0%2;stan1i1=stan1i1%2;卷积码编码程序的c 语言实现cout 状态转移出 endl;for(i=0;istalen;i+)coutstachni0,stachni1;coutendl;cout 输入 0 状态转移后的输出endl;for(i=0;istalen;i+)coutstan0i0,stan0i1;coutendl;cout 输入 1 状态转移后的输出endl;for(i=0;istalen;i+)coutstan1i0,stan1i1;coutendl;void chart
9、obits(char ch,int *bits)int i;for(i=0;i8;i+)if(ch0)bitsi=1;elsebitsi=0;ch=ch1;char bitstochar(int *bits)char temp=0;int i;for(i=0;i8;i+)if(bitsi=1)temp+=table17-i;return temp;void convolution()FILE *fp_input,*fp_output;if(!(fp_input=fopen(D:/input、 txt,r)=1)coutfailed to open input_fileendl;exit(0);
10、elsecoutwe opened the input_file endl;if(!(fp_output=fopen(D:/output、 txt,w+)=1)coutfailed to open output_fileendl;exit(0);卷积码编码程序的c 语言实现elsecoutwe opened the output_file endl;char ch;int i,j;int mybits8,mytembits8;int mysta=0;int wcout;char wch;for(ch=fgetc(fp_input);feof(fp_input)=0;ch=fgetc(fp_in
11、put)chartobits(ch,mybits);/*cout 输入数据为 endl;for(i=0;i8;i+)coutmybitsi;coutendl;*/for(i=0;i8;i+)if(mybitsi=0)myoutmyoutsym+=stan0mysta0;myoutmyoutsym+=stan0mysta1;mysta=stachnmysta0;elsemyoutmyoutsym+=stan1mysta0;myoutmyoutsym+=stan1mysta1;mysta=stachnmysta1;/*cout 输出数据 1 为 endl;for(temi=0;temimyoutsym;temi+)coutmyouttemi;coutendl;*/if