卷积码编码程序的c语言实现
213卷积码编码和译码
No.15 (2,1,3)卷积码的编码及译码摘要：本报告对于 (2,1,3)卷积码原理部分的论述主要参照啜刚教材和课件，编程仿真部分绝对原创，所有的程序都是在Codeblocks 8.02环境下用C语言编写的，编译运行都正常。
完成了卷积码的编码程序，译码程序，因为对于短于 3 组的卷积码，即2 bit 或4 bit 纠错是没有意义的，所以对正确的短序列直接译码，对长序列纠错后译码，都能得到正确的译码结果。
含仿真结果和程序源代码。
如果您不使用Codeblocks 运行程序，则可能不支持中文输出显示，但是所有的数码输出都是正确的。
一、卷积码编码原理卷积码编码器对输入的数据流每次1bit或k bit进行编码，输出n bit编码符号。
但是输出的分支码字的每个码元不仅于此时可输入的k个嘻嘻有关，业余前m个连续式可输入的信息有关，因此编码器应包含m级寄存器以记录这些信息。
k通常卷积码表示为（n,k,m）.编码率r —n当k=1时，卷积码编码器的结构包括一个由m个串接的寄存器构成的移位寄存器（成为m级移位寄存器、n个连接到指定寄存器的模二加法器以及把模二加法器的输出转化为穿行的转换开关。
本报告所讲的（2,1,3）卷积码是最简单的卷积码。
就是n 2，k 1，m 3的卷积码。
每次输入1 bit输入信息，经过3级移位寄存器，2个连接到指定寄存器的模二加法器，并把加法器输出转化为串行输出。
EncoderOuTi>iit编码器如题所示。
二、卷积码编码器程序仿真C语言编写的仿真程序。
为了简单起见，这里仅仅提供数组长度30 bit的仿真程序，当然如果需要可以修改数组大小。
为了更精练的实现算法，程序输入模块没有提供非法字符处理过程，如果需要也可以增加相应的功能。
进入程序后，先提示输入数据的长度，请用户输入int （整型数）程序默认用户输入的数据小于30，然后提示输入01数码，读入数码存储与input数组中，然后运算输出卷积码。
卷积码
研09电子 庆栋良 09110023
(1)在实际应用中怎样去建立网格图？ (2)怎样计算最小汉明距离？ (3)在找到最优路径后怎样去译出原始的 信息码？
内容简介
卷积码简单介绍 卷积码编码实现 卷积码译码实现 总结
卷积码简单介绍
卷积码(convolutional code)是由伊利亚斯(p.Elias)发 明的一种非分组码。在前向纠错系统中，卷积码在实际 应用中的性能优于分组码，并且运算较简单。 卷积码在编码时将k比特的信息段编成n个比特的码 组，监督码元不仅和当前的k比特信息段有关，而且还同 前面m=(N-1)个信息段有关。 通常将N称为编码约束长度，将nN称为编码约束长 度。一般来说，卷积码中k和n的值是比较小的整数。将 卷积码记作(n,k,N)。
图5 (2,1,4)码网格图
②寻找最优路径
网格图建立之后，根据接收码组和网格图中生成 的码组比较，判断最优路径。假设某一时刻的状态i， 首先判断前一时刻所有状态中，是哪两个状态指向当 前状态i；其次，根据这两个指向当前状态i的状态生成 的码组和前一时刻接收的码组比较，保留汉明距离最 小的那条路径以及到达状态i时的最小汉明距离。下一 时刻，同样操作，但是保留的最小汉明距离是前面最 小汉明距离累加。 在所有的接收码组处理完之后，会得到一组汉明 距离以及所对应的最优路径。比较选择出最小的汉明 距离，那么该最小汉明距离所对应的路径即为最优路 径。
③译码
如图5所示，其中红线表示计算得到的最优路 径。观察其变化规律，可以发现，假设i时刻的状 态和第(i+1)时刻比较得知，若第(i+1)时刻小于第 i时刻，那么对应代表实际信息中0；若大于，则 为1。然后逐次比较译出信息值。
总结
实验九-(2-1-5)卷积码编码译码技术
实验九 (2,1,5）卷积码编码译码技术一、实验目的1、掌握(2,1,5）卷积码编码译码技术2、了解纠错编码原理。
二、实验内容1、(2,1,5）卷积码编码。
2、(2,1,5）卷积码译码。
三、预备知识1、纠错编码原理。
2、(2,1,5）卷积码的工作原理。
四、实验原理/卷积码是将发送的信息序列通过一个线性的，有限状态的移位寄存器而产生的编码。
通常卷积码的编码器由K级（每级K比特）的移位寄存器和n个线性代数函数发生器（这里是模2加法器）组成。
若以（n,k,m）来描述卷积码，其中k为每次输入到卷积编码器的bit数，n 为每个k元组码字对应的卷积码输出n元组码字，m为编码存储度，也就是卷积编码器的k元组的级数，称m+1= K为编码约束度m称为约束长度。
卷积码将k 元组输入码元编成n元组输出码元，但k和n通常很小，特别适合以串行形式进行传输，时延小。
与分组码不同，卷积码编码生成的n元组元不仅与当前输入的k元组有关，还与前面m-1个输入的k元组有关，编码过程中互相关联的码元个数为n*m。
卷积码的纠错性能随m的增加而增大，而差错率随N的增加而指数下降。
在编码器复杂性相同的情况下，卷积码的性能优于分组码。
编码器随着信息序列不断输入，编码器就不断从一个状态转移到另一个状态并同时输出相应的码序列，所以图3所示状态图可以简单直观的描述编码器的编码过程。
因此通过状态图很容易给出输入信息序列的编码结果，假定输入序列为110100，首先从零状态开始即图示a状态，由于输入信息为“1”，所以下一状态为b并输出“11”，继续输入信息“1”，由图知下一状态为d、输出“01”……其它输入信息依次类推，按照状态转移路径a->b->d->c->b->c->a输出其对应的编码结果“”。
译码方法⒈代数代数译码是将卷积码的一个编码约束长度的码段看作是[n0(m+1），k0(m+1)]线性分组码，每次根据（m+1）分支长接收数字，对相应的最早的那个分支上的信息数字进行估计，然后向前推进一个分支。
matlab的conv的c源代码
matlab的conv的c源代码MATLAB的conv函数是一种用于进行卷积运算的函数。
卷积运算在信号处理、图像处理、语音识别等领域中具有广泛的应用。
本文将介绍MATLAB中conv函数的C源代码实现。
在MATLAB中，conv函数可以用于计算一维或二维信号的线性卷积。
其基本语法如下：```matlabC = conv(A, B)```其中A和B是需要进行卷积运算的输入向量或矩阵，C是卷积运算的结果。
接下来，我们将展示一个C源代码实现的例子。
请注意，为了简化代码，我们将仅考虑一维信号的卷积运算。
```c#include <stdio.h>#include <stdlib.h>void conv(double A[], int size_A, double B[], int size_B, double C[]) {int size_C = size_A + size_B - 1; // 计算卷积结果的长度int i, j;// 初始化结果数组Cfor(i = 0; i < size_C; i++){C[i] = 0;}// 进行卷积运算for(i = 0; i < size_C; i++){for(j = 0; j < size_A; j++){if(i - j >= 0 && i - j < size_B) // 检查是否越界 {C[i] += A[j] * B[i - j];}}}}int main(){double A[] = {1, 2, 3, 4, 5};double B[] = {0.5, 0.5};int size_A = sizeof(A) / sizeof(A[0]);int size_B = sizeof(B) / sizeof(B[0]);int size_C = size_A + size_B - 1;double C[size_C];int i;conv(A, size_A, B, size_B, C);// 输出卷积结果for(i = 0; i < size_C; i++){printf("%lf ", C[i]);}return 0;}```以上是使用C语言实现的MATLAB conv函数的简单示例。
实验九 (2,1,5)卷积码编码译码技术
实验九 (2,1,5）卷积码编码译码技术一、实验目的1、掌握(2,1,5）卷积码编码译码技术2、了解纠错编码原理。
二、实验内容1、(2,1,5）卷积码编码。
2、(2,1,5）卷积码译码。
三、预备知识1、纠错编码原理。
2、(2,1,5）卷积码的工作原理。
四、实验原理卷积码是将发送的信息序列通过一个线性的，有限状态的移位寄存器而产生的编码。
通常卷积码的编码器由K级（每级K比特）的移位寄存器和n个线性代数函数发生器（这里是模2加法器）组成。
若以（n,k,m）来描述卷积码，其中k为每次输入到卷积编码器的bit数，n 为每个k元组码字对应的卷积码输出n元组码字，m为编码存储度，也就是卷积编码器的k元组的级数，称m+1= K为编码约束度m称为约束长度。
卷积码将k 元组输入码元编成n元组输出码元，但k和n通常很小，特别适合以串行形式进行传输，时延小。
与分组码不同，卷积码编码生成的n元组元不仅与当前输入的k元组有关，还与前面m-1个输入的k元组有关，编码过程中互相关联的码元个数为n*m。
卷积码的纠错性能随m的增加而增大，而差错率随N的增加而指数下降。
在编码器复杂性相同的情况下，卷积码的性能优于分组码。
编码器随着信息序列不断输入，编码器就不断从一个状态转移到另一个状态并同时输出相应的码序列，所以图3所示状态图可以简单直观的描述编码器的编码过程。
因此通过状态图很容易给出输入信息序列的编码结果，假定输入序列为110100，首先从零状态开始即图示a状态，由于输入信息为“1”，所以下一状态为b并输出“11”，继续输入信息“1”，由图知下一状态为d、输出“01”……其它输入信息依次类推，按照状态转移路径a->b->d->c->b->c->a输出其对应的编码结果“110101001011”。
译码方法⒈代数代数译码是将卷积码的一个编码约束长度的码段看作是[n0(m+1），k0(m+1)]线性分组码，每次根据（m+1）分支长接收数字，对相应的最早的那个分支上的信息数字进行估计，然后向前推进一个分支。
C语言卷积码(附注释)
ret_dis(j) = dis_temp1;
else
ret_msg(j,1:chip+i) = msg_temp2;
ret_dis(j) = dis_temp2;%选择较小汉明距的状态储存并输出在ret_msg中，最小汉明距存于ret_dis
word1=word
msg_1 = decode_conv213(word1)
msg-msg_1
function[output,nextState] = state_machine(input,current_state)
output(1) = mod(current_state(1)+current_state(3),2);
output(2) = mod(input+current_state(2)+current_state(1),2);
nextState(1) = current_state(2);
nextState(2) = current_state(3);
nextState(3) = input;
functiondistance = hamming_distance(a,b)
fori = 1:2^chip
M(i,:) = de2bi(i-1,chip);%把所有可能性按二进制输出
W(i,:) = encode_conv213(M(i,:));%得到相应的二进制编译后的码字
D(i) = hamming_distance(W(i,:),word(1:chip*2));%与出错码字对比得到汉明距
ret_msg(i,1:chip) = M(index(i),:);
对卷积码的编译码实现进行仿真编程和仿真实验
实验二：卷积码编译器1120510217 吴宏晶一、实验目的对卷积码的编译码实现进行仿真编程和仿真实验二、仿真软件LabVIEW 2013版（前面版用于数据显示，后面板为程序源代码）三、实验流程1.卷积码的编码以（3,1,2）卷积码为例，下图为LabVIEW仿真前面板每一时刻都有一个信息比特进入第一个寄存器中，如图中D1。
编码器利用D1和D2这两个寄存器来产生3比特的输出码字如图中所示C1、C2、C3。
编码的规则是C1=M1，C2=M1+D1+D2，C3=M1+D2。
定义D1和D2的状态为S1和S2，S1和S2共有四种组合状态。
描述的在不同的输入比特状态下，编码器对应的输出和下一个时刻的状态变化。
它们的转化关系如下图所示2.卷积码的译码卷积码的译码采用维特比译码法。
维特比算法的思路：在每个节点处，留下到达该节点处局部度量最大的路径，丢弃其余的局部路径，这条留下来的路径叫做幸存路径。
如果所有幸存路径在某个支路处是重合的，就可以输出这个支路上数据的译码结果。
为了避免延时，在当前局部路径中，找出局部累积量最大的，沿着这个路径回退N步会有一个分支，输出这个分支的数据。
（3,1,2）卷积码t0~t5六个时刻，第一次选取幸存路径发生在从t2到t3转移时刻。
四、实验源程序说明1.卷积码的编码程序源代码如下图所示1.坐标1处是输入码元数组，它经过一个FOR循环依次进入循环中。
2.坐标2处是D1和D2两个移位寄存器，D1的结果在一次循环过后移入D2中。
3.坐标3处表示C1的值为每次直接进来的码元值。
4.坐标4处为C2的计算值，对局部进行放大，有如下图说明：C2=D1+D2+M1（模2加法），M1为码元值。
因为没有模2加法的函数，故调用除2函数，再对余数与0进行比较是否相等，若结果为TRUE，则表明C2为2的倍数，即为0，故将0送入C2中，否则将1送入C2中。
5.坐标5处为C3的计算值，对局部进行放大，有如下图说明：C3=M1+C2，故将两数相加，再对结果与1比较，因为只有01或10相加才为1，其他值为0或者2，故若判断结果为TRUE,则将1的值赋值给C3，如上图所示。
卷积码C语言
include <stdio.h>include "Conio.h"define N 7include "math.h"include <stdlib.h>include<time.h>define randomize srandunsignedtimeNULLencodeunsigned int symbols; /编码输出/unsigned int data; /编码输入/unsigned int nbytes; /nbytes=n/16;n为实际输入码字的数目/ unsigned int startstate /定义初始化状态/{int j;unsigned int input;a1=0;a2=0;forj=0;j<nbytes;j++{ input=data;data++;symbols = input^a1^a2;symbols++;symbols = input^a2;symbols++;a2=a1;a1=input;}return 0;}int trandistanceint m; /符号m与从state1到state2时输出符号的汉明距离;如果state1无法到state2则输出度量值为100/int state1;int state2{int c;int sym;sym1;sym2;sym1=state2>>1&1^state2&1^state1&1;sym2=state2>>1&1^state1&1;sym=sym1<<1 | sym2;if state1&2>>1==state2&1c=m&1^sym&1+m>> 1&1^sym >> 1&1;elsec=10000;returnc;}int traninputint a;int b /状态从a到b时输入卷积码的符号/{int c;c=b&2>>1;returnc;}int tranoutputint a;int b /状态从a到b时卷积码输出的符号/{int c;s1;s2;s1=a&1^a&2>>1^b&2>>1;s2=a&1^b&2>>1;c=s1<<1|s2;returnc;}void viterbiint initialstate; /定义解码器初始状态/int viterbiinput; /解码器输入码字序列/int viterbioutput /解码器输出码字序列/{struct sta /定义网格图中每一点为一个结构体;其元素包括/{int met; /转移到此状态累计的度量值/int value; /输入符号 /struct sta last; /及指向前一个状态的指针/};struct sta state4N;struct sta g;head;int i;j;p;q;t;r;u;l;fori=0;i<4;i++ / 初始化每个状态的度量值/forj=0;j<N;j++stateij.met=0;forl=0;l<4;l++{statel0.met=trandistanceviterbiinput;initialstate;l;statel0.value=traninputinitialstate;l;st=NULL;}viterbiinput++; /扩展第一步幸存路径/fort=1;t<N;t++{forp=0;p<4;p++{statept.met=state0t-1.met+trandistanceviterbiinput;0;p; statept.value=traninput0;p;st=&state0t-1;forq=0;q<4;q++{ifstateqt-1.met+trandistanceviterbiinput;q;p<statept.met{statept.met=stateqt-1.met+trandistanceviterbiinput;q;p;statept.value=traninputq;p;st=&stateqt-1;}}}viterbiinput++;} /计算出剩余的幸存路径/r=state0N-1.met; /找出n步后度量值最小的状态;准备回溯路由/g=&state0N-1;foru=N;u>0;u-- /向前递归的找出最大似然路径 /{viterbioutput+u-1=g->value;g=g->last;}/ foru=0;u<8;u++viterbioutput+u=stateu2.met; / /此行程序可用于检测第n 列的度量值/}void decodeunsigned int input; int output;int n{int viterbiinput100;int j;forj=0;j<n+2;j++{viterbiinputj=inputj2<<1|inputj2+1;//printf"%3d";viterbiinputj;}viterbi0;viterbiinput;output;}void main{unsigned int encodeinput100;wrong10={0;0;0;0;0;0;0;0;0;0};encodeoutput100;int n=5;i;m;j=0;decodeinput100;decodeoutput100;randomize;fori=0; i<n; i++encodeinputi=rand%2;encodeinputn= encodeinputn+1=0;encodeencodeoutput;encodeinput;n+2;0;printf"the input of encoder is :\n";fori=0;i<n; i++printf"%2d";encodeinputi;printf"\n";printf"the output of encoder is :\n";fori=0;i<n+22;i++{printf"%2d";encodeoutputi;ifi%20==19printf"\n";}printf"\n";printf"please input the number of the wrong bit\n";scanf"%d";&m;printf"please input the positions of the wrong bit0-9\n";fori=0;i<m;i++{ scanf"%d";&wrongm;ifencodeoutputwrongm==0encodeoutputwrongm=1;elseencodeoutputwrongm=0;}printf"the input of decoder is :\n";fori=0;i<n+22;i++{printf"%2d";encodeoutputi;ifi%20==19printf"\n";}printf"\n";decodeencodeoutput;decodeoutput;n+2;printf"the output of decoder is :\n";fori=0;i<n;i++printf"%2d";decodeoutputi;printf"\n";fori=0;i<n;i++{ ifencodeinputi=decodeoutputij++;}printf"the number of incorrect bit is:%d\n";j; }。
卷积码编码和维特比译码
卷积码编码维特比译码实验设计报告SUN一、实验目的掌握卷积码编码和维特比译码的基本原理，利用了卷积码的特性, 运用网格图和回溯以得到译码输出。
二、实验原理1.卷积码是由连续输入的信息序列得到连续输出的已编码序列。
其编码器将k个信息码元编为n个码元时，这n个码元不仅与当前段的k个信息有关，而且与前面的（m-1）段信息有关（m为编码的约束长度）。
2.一般地，最小距离d表明了卷积码在连续m段以内的距离特性，该码可以在m个连续码流内纠正(d-1)/2个错误。
卷积码的纠错能力不仅与约束长度有关，还与采用的译码方式有关。
3. 维特比译码算法基本原理是将接收到的信号序列和所有可能的发送信号序列比较，选择其中汉明距离最小的序列认为是当前发送序列。
卷积码的Viterbi 译码是根据接收码字序列寻找编码时通过网格图最佳路径的过程,找到最佳路径即完成了译码过程,并可以纠正接收码字中的错误比特。
4.所谓“最佳”, 是指最大后验条件概率:P( C/ R) = max [ P ( Cj/ R) ] , 一般来说, 信道模型并不使用后验条件概率,因此利用Beyes 公式、根据信道特性出结论:max[ P ( Cj/ R) ]与max[ P ( R/ Cj) ]等价。
考虑到在系统实现中往往采用对数形式的运算,以求降低运算量,并且为求运算值为整数加入了修正因子a1 、a2 。
令M ( R/ Cj) = log[ P ( R/ Cj) ] =Σa1 (log[ P( Rm/ Cmj ) ] + a2) 。
其中, M 是组成序列的码字的个数。
因此寻找最佳路径, 就变成寻找最大M( R/ Cj) , M( R/ Cj) 称为Cj 的分支路径量度,含义为发送Cj 而接收码元为R的似然度。
5.卷积码的viterbi译码是根据接收码字序列寻找编码时通过网格图最佳路径的过程，找到最佳路径即完成了译码过程并可以纠正接收码字中的错误比特。
三、实验代码#include<stdio.h>#include "Conio.h"#define N 7#include "math.h"#include <stdlib.h>#include<time.h>#define randomize() srand((unsigned)time(NULL))encode(unsigned int *symbols, /*编码输出*/unsigned int *data, /*编码输入*/unsigned int nbytes, /*nbytes=n/16,n为实际输入码字的数目*/unsigned int startstate /*定义初始化状态*/)////////////////////////////////////////////////////////////////////////////卷积码编码///////////////////////////////////////////////////////////////////////////////{unsigned int j;unsigned int input,a1=0,a2=0,a3=0,a4=0,a5=0,a6=0;for(j=0;j<nbytes;j++){input=*data;data++;*symbols = input^a1^a2^a3^a6; //c1(171)symbols++;*symbols = input^a2^a3^a5^a6; //c2(133)symbols++;a2=a1;a1=input;}return 0;}int trandistance(int m, int state1, int state2)/*符号m与从state1到state2时输出符号的汉明距离,如果state1无法到state2则输出度量值为100*/{int c;int sym,sym1,sym2;sym1=((state2>>1)&1)^(state2&1)^(state1&1);sym2=((state2>>1)&1)^(state1&1);sym=(sym1<<1) | sym2;if ( ((state1&2)>>1)==(state2&1))c=((m&1)^(sym&1))+(((m>> 1)&1)^((sym >> 1)&1));elsec=10000;return(c);}int traninput(int a,int b) /*状态从a到b时输入卷积码的符号*/{int c;c=((b&2)>>1);return(c);}int tranoutput(int a,int b) /*状态从a到b时卷积码输出的符号*/{int c,s1,s2;s1=(a&1)^((a&2)>>1)^((b&2)>>1);s2=(a&1)^((b&2)>>1);c=(s1<<1)|s2;return(c);}////////////////////////////////////////////////////////////////////////////维特比译码///////////////////////////////////////////////////////////////////////////////void viterbi(int initialstate, /*定义解码器初始状态*/int *viterbiinput, /*解码器输入码字序列*/int *viterbioutput /*解码器输出码字序列*/){struct sta /*定义网格图中每一点为一个结构体,其元素包括*/ {int met; /*转移到此状态累计的度量值*/int value; /*输入符号*/struct sta *last; /*及指向前一个状态的指针*/};struct sta state[4][N];struct sta *g,*head;int i,j,p,q,t,r,u,l;for(i=0;i<4;i++) /* 初始化每个状态的度量值*/for(j=0;j<N;j++)state[i][j].met=0;for(l=0;l<4;l++){state[l][0].met=trandistance(*viterbiinput,initialstate,l);state[l][0].value=traninput(initialstate,l);state[l][0].last=NULL;}viterbiinput++; /*扩展第一步幸存路径*/for(t=1;t<N;t++){for(p=0;p<4;p++){state[p][t].met=state[0][t-1].met+trandistance(*viterbiinput,0,p);state[p][t].value=traninput(0,p);state[p][t].last=&state[0][t-1];for(q=0;q<4;q++){if(state[q][t-1].met+trandistance(*viterbiinput,q,p)<state[p][t].met){state[p][t].met=state[q][t-1].met+trandistance(*viterbiinput,q,p);state[p][t].value=traninput(q,p);state[p][t].last=&state[q][t-1];}}}viterbiinput++;} /*计算出剩余的幸存路径*/r=state[0][N-1].met; /*找出n步后度量值最小的状态准备回溯路由*/g=&state[0][N-1];for(u=N;u>0;u--) /*向前递归的找出最大似然路径*/{*(viterbioutput+(u-1))=g->value;g=g->last;}/* for(u=0;u<8;u++)*(viterbioutput+u)=state[u][2].met; */ /*此行程序可用于检测第n列的度量值*/}void decode(unsigned int *input, int *output,int n){int viterbiinput[100];int j;for(j=0;j<n+2;j++){viterbiinput[j]=(input[j*2]<<1)|input[j*2+1];}viterbi(0,viterbiinput,output);}void main(){unsigned intencodeinput[100],wrong[10]={0,0,0,0,0,0,0,0,0,0},encodeoutput[100];int n=5,i,m,j=0,decodeinput[100],decodeoutput[100];randomize();for(i=0; i<n; i++)encodeinput[i]=rand()%2;encodeinput[n]= encodeinput[n+1]=0;encode(encodeoutput,encodeinput,n+2,0);printf("the input of encoder is :\n"); //信息源输入的信息码（随机产生）for(i=0;i<n; i++)printf("%2d",encodeinput[i]);printf("\n");printf("the output of encoder is :\n"); //编码之后产生的卷积码for(i=0;i<(n+2)*2;i++){printf("%2d",encodeoutput[i]);if(i%20==19)printf("\n");}printf("\n");printf("please input the number of the wrong bit\n"); //信道传输收到干扰而产生的错误码scanf("%d",&m);printf("please input the positions of the wrong bit(0-9)\n");for(i=0;i<m;i++){scanf("%d",&wrong[m]);if(encodeoutput[wrong[m]]==0)encodeoutput[wrong[m]]=1;elseencodeoutput[wrong[m]]=0;}printf("the input of decoder is :\n");for(i=0;i<(n+2)*2;i++){printf("%2d",encodeoutput[i]);if(i%20==19)printf("\n");}printf("\n");decode(encodeoutput,decodeoutput,n+2);printf("the output of decoder is :\n");for(i=0;i<n;i++)printf("%2d",decodeoutput[i]);printf("\n");for(i=0;i<n;i++){if(encodeinput[i]!=decodeoutput[i])j++;}printf("the number of incorrect bit is:%d\n",j);}四、实验总结（1）了解实验原理，分析实验所占数组变量很重要，也是相对考虑较多的；（2）对于读写文件，通过本实验更加熟悉；（3）记录实验程序最佳路径是本实验的难点；。
卷积码编码程序的c语言实现
#include<iostream>using namespace std;int table1[8]={1,2,4,8,16,32,64,128};int myn=0;int stalen=0;int stan0[256][2]={0};//输入0时个状态的输出int stan1[256][2]={0};//输入1时各状态的输出int stachn[256][2]={0};//状态装换表int myg1[10]={0};int myg2[10]={0};int myout[100]; //int myoutsym=0;void chartobits(char ch,int *bits);char bitstochar(int *bits);void convolution(void);void creatsta(void);void myinput(void);int main(){char exit_char;myinput();creatsta();convolution();cin>>exit_char;}void myinput(void){int i,j;cout<<"输入编码的约束长度N:(3<N<9)"<<endl;cin>>myn;stalen=int(pow(2、0,myn-1));cout<<"选择默认的编码矢量则输入1,输入2则可输入其她的编码矢量"<<endl;cin>>i;if(i==1){switch(myn){case 3: myg1[0]=1,myg1[1]=1,myg1[2]=1;myg2[0]=1,myg2[1]=0,myg2[2]=1;break;case 4: myg1[0]=1,myg1[1]=1,myg1[2]=1,myg1[3]=1;myg2[0]=1,myg2[1]=0,myg2[2]=1,myg2[3]=1;break;case 5: myg1[0]=1,myg1[1]=0,myg1[2]=1,myg1[3]=1,myg1[4]=1;myg2[0]=1,myg2[1]=1,myg2[2]=0,myg2[3]=1,myg2[4]=1;break;case 6: myg1[0]=1,myg1[1]=0,myg1[2]=1,myg1[3]=1,myg1[4]=1,myg1[5]=1;myg2[0]=1,myg2[1]=1,myg2[2]=0,myg2[3]=1,myg2[4]=0,myg2[5]=1;break;case 7: myg1[0]=1,myg1[1]=0,myg1[2]=0,myg1[3]=1,myg1[4]=1,myg1[5]=1,myg1[6]=1;myg2[0]=1,myg2[1]=1,myg2[2]=0,myg2[3]=1,myg2[4]=1,myg2[5]=0,myg2[6]=1;break;case 8: myg1[0]=1,myg1[1]=0,myg1[2]=0,myg1[3]=1,myg1[4]=1,myg1[5]=1,myg1[6]=1,myg1[7]=1;myg2[0]=1,myg2[1]=1,myg2[2]=1,myg2[3]=0,myg2[4]=0,myg2[5]=1,myg2[6]=0,myg2[7]=1;break;case 9: myg1[0]=1,myg1[1]=1,myg1[2]=0,myg1[3]=1,myg1[4]=0,myg1[5]=1,myg1[6]=1,myg1[7]=1,m yg1[8]=1;myg2[0]=1,myg2[1]=0,myg2[2]=0,myg2[3]=0,myg2[4]=1,myg2[5]=1,myg2[6]=1,myg2[7]=0,m yg2[8]=1;break;}}else{cout<<"输入g1"<<endl;for(j=0;j<myn;j++)cin>>myg1[j];cout<<"输入g2"<<endl;for(j=0;j<myn;j++)cin>>myg2[j];}cout<<"连接矢量1为"<<endl;for(j=0;j<myn;j++)cout<<myg1[j]<<" ";cout<<endl;cout<<"连接矢量2为"<<endl;for(j=0;j<myn;j++)cout<<myg2[j]<<" ";cout<<endl;cout<<"stalen: "<<stalen;cout<<endl;}void creatsta(void){int i,j,k,myi,myj;int tembits[10];for(i=0;i<stalen;i++){stan1[i][0]=0;stan1[i][1]=0;stan0[i][0]=0;stan0[i][1]=0;stachn[i][0]=i/2;myi=i;for(j=0;j<myn;j++){if(myi>=pow(2、0,myn-1-j)){tembits[j]=1;myi=myi-pow(2、0,myn-1-j);}else{tembits[j]=0;}}/*for(k=0;k<myn;k++)cout<<tembits[k]<<" ";cout<<endl;*/for(k=0;k<myn;k++){stan0[i][0]+=myg1[k]*tembits[k];stan0[i][1]+=myg2[k]*tembits[k];}stan0[i][0]=stan0[i][0]%2;stan0[i][1]=stan0[i][1]%2;myi=i+int(pow(2、0,myn-1));stachn[i][1]=myi/2;for(j=0;j<myn;j++){if(myi>=pow(2、0,myn-1-j)){tembits[j]=1;myi=myi-pow(2、0,myn-1-j);}else{tembits[j]=0;}}/*for(k=0;k<myn;k++)cout<<tembits[k]<<" ";cout<<endl;*/for(k=0;k<myn;k++){stan1[i][0]+=myg1[k]*tembits[k];stan1[i][1]+=myg2[k]*tembits[k];}stan1[i][0]=stan1[i][0]%2;stan1[i][1]=stan1[i][1]%2;}cout<<"状态转移出"<<endl;for(i=0;i<stalen;i++)cout<<stachn[i][0]<<","<<stachn[i][1]<<" ";cout<<endl;cout<<"输入0状态转移后的输出"<<endl;for(i=0;i<stalen;i++)cout<<stan0[i][0]<<","<<stan0[i][1]<<" ";cout<<endl;cout<<"输入1状态转移后的输出"<<endl;for(i=0;i<stalen;i++)cout<<stan1[i][0]<<","<<stan1[i][1]<<" ";cout<<endl;}void chartobits(char ch,int *bits){int i;for(i=0;i<8;i++){if(ch<0)bits[i]=1;elsebits[i]=0;ch=ch<<1;}}char bitstochar(int *bits){char temp=0;int i;for(i=0;i<8;i++){if(bits[i]==1)temp+=table1[7-i];}return temp;}void convolution(){FILE *fp_input,*fp_output;if(!(fp_input=fopen("D://input、txt","r"))==1){cout<<"failed to open input_file"<<endl;exit(0);}elsecout<<"we opened the input_file "<<endl;if(!(fp_output=fopen("D://output、txt","w+"))==1){ cout<<"failed to open output_file"<<endl;exit(0);}cout<<"we opened the output_file "<<endl;char ch;int i,j;int mybits[8],mytembits[8];int mysta=0;int wcout;char wch;for(ch=fgetc(fp_input);feof(fp_input)==0;ch=fgetc(fp_input)){ chartobits(ch,mybits);/*cout<<"输入数据为"<<endl;for(i=0;i<8;i++)cout<<mybits[i]<<" ";cout<<endl;*/for(i=0;i<8;i++){if(mybits[i]==0){myout[myoutsym++]=stan0[mysta][0];myout[myoutsym++]=stan0[mysta][1];mysta=stachn[mysta][0];}else{myout[myoutsym++]=stan1[mysta][0];myout[myoutsym++]=stan1[mysta][1];mysta=stachn[mysta][1];}}/*cout<<"输出数据1为"<<endl;for(temi=0;temi<myoutsym;temi++)cout<<myout[temi]<<" ";cout<<endl;*/if(myoutsym>=8){wcout=myoutsym/8;for(i=0;i<wcout;i++){for(j=0;j<8;j++)mytembits[j]=myout[8*i+j];wch=bitstochar(mytembits);fputc(wch,fp_output);/*cout<<"输出数据2为"<<endl;cout<<wch<<" ";*/}for(i=0;i<myoutsym-wcout*8;i++)myout[i]=myout[wcout*8+i];myoutsym=myoutsym-wcout*8;}for(i=0;i<myn-1;i++){myout[myoutsym++]=stan0[mysta][0];myout[myoutsym++]=stan0[mysta][1];mysta=stachn[mysta][0];}wcout=myoutsym/8;for(i=0;i<wcout;i++){for(j=0;j<8;j++)mytembits[j]=myout[8*i+j];wch=bitstochar(mytembits);fputc(wch,fp_output);}for(i=0;i<myoutsym-wcout*8;i++)myout[i]=myout[wcout*8+i];myoutsym=myoutsym-wcout*8;if(myoutsym!=0){for(i=0;i<8-myoutsym;i++)myout[myoutsym++]=0;}for(j=0;j<8;j++)mytembits[j]=myout[j];wch=bitstochar(mytembits);fputc(wch,fp_output);fclose(fp_input);fclose(fp_output);}。
通信原理实验19 卷积码的编解码实验
实验十九卷积码的编解码实验实验内容1. 熟悉卷积码编码实验2．熟悉卷积码译码实验一、实验目的1．了解卷积码的基本概念和原理2．加深对卷积码的编解码过程的理解3. 学习通过CPLD编程实现卷积码编译码实验二、实验电路工作原理卷积码又称连环码，是1955年提出来的一种纠错码，它和分组码有明显的区别，但在编码器复杂度相同的情况下, 卷积码的性能优于分组码,因此卷积码几乎被应用在所有无线通信的标准之中, 如GSM, IS95和CDMA 2000 的标准中。
1．卷积码编码方法：卷积码通常记作( n0 , k0 , m) ,它将k0 个信息比特编为n0 个比特, 其编码效率为k0/ n0 , m为约束长度。
( n0 , k0 , m )卷积码可用k0 个输入、n0 个输出、输入存储为m的线性有限状态移位寄存器及模2 加法计数器来实现。
实验中所选(2 ,1 ,6) 卷积编码器上图所示,其子生成元为: g(1 ,1) ( D) = 1 , g(1 ,2) ( D) = 1 + D2 + D5 +D6 ,生成矩阵G( D) = (1 ,1 + D2 + D5 + D6) 。
设输入信息序列M = (1111) ,即M( D) = 1 + D + D2 + D3 ,则编码器的输出C( D) = M( D) ·G( D) ,即:C( D) = (1+D+D2+D3)·(1,1+D2+D5+D6)= (1+D+D2+D3 ,1+D+D2+D3+D2+D3+D4+D5+D5+D6+D7+D8+D6+D7+D8+D9)= (1+D+D2+D3 ,1+D+D4+D9)= (11)+(11)D+(10)D2+(10)D3+(01)D4+(00)D5+(00)D6+(00)D7+(00)D8+(01)D9+⋯因此,编码器输出序列为11111010010000000001。
2．卷积码编码算法process(clk,clr)beginif(clr='1')thenif(clk'event and clk='1')thentemp(0)<=datain;temp(1)<=temp(0);temp(2)<=temp(1);temp(3)<=temp(2);temp(4)<=temp(3);end if;else temp<="00000";end if;end process;y2j<= (datain xor temp(2) xor temp(3) xor temp(4));y1j<=datain;3．大数逻辑解码器大数逻辑解码器是卷积码代数解码最主要的解码方法, 既可用于纠随机错误, 又可用于纠突发错误,但要求卷积码是自正交码或可正交码。
卷积 维特比译码 c语言
卷积维特比译码c语言卷积码和维特比（Viterbi）译码是两种在通信系统中常用的编码和解码技术。
卷积码是一种线性分组码，它通过将输入信息序列与一组预定的约束条件进行卷积运算来生成编码序列。
而维特比译码是一种高效的解码算法，用于从接收到的信号中恢复出原始的编码序列。
下面是一个简单的C语言示例，演示了如何实现卷积编码和维特比译码。
请注意，这只是一个基本的示例，实际的实现可能会更复杂，并且需要更多的错误处理和优化。
c复制代码:#include <stdio.h>#include <stdlib.h>#define N 4#define K 2#define G 2// 卷积编码函数void convolutional_encode(int input[K], int code[N]) {int i, j;for (i = 0; i < N; i++) {code[i] = 0;}for (i = 0; i < K; i++) {code[i] = input[i];}for (i = K; i < N; i++) {code[i] = (code[i - 1] ^ code[i - 2]) & G;}}// 维特比译码函数void viterbi_decode(int received[N], int output[K]) {int branch_metric[N][2];int path_metric[N];int max_metric, new_max_metric;int max_path, new_max_path;int i, j, k;for (i = 0; i < N; i++) {path_metric[i] = abs(received[i] - 0); // 初始化路径度量branch_metric[i][0] = abs(received[i] - 0); // 初始化分支度量branch_metric[i][1] = abs(received[i] - G); // 初始化分支}for (i = 0; i < N; i++) {output[i] = 0; // 初始化输出序列}for (i = 0; i < N; i++) {if (branch_metric[i][0] > branch_metric[i][1]) {output[i] = 0; // 选择分支0作为当前的最优路径path_metric[i] = branch_metric[i][0]; // 更新路径度量} else {output[i] = 1; // 选择分支1作为当前的最优路径path_metric[i] = branch_metric[i][1]; // 更新路径度量}if (path_metric[i] > max_metric) {max_metric = path_metric[i]; // 记录最大路径度量max_path = i; // 记录最大路径度量对应的路径值} else if (path_metric[i] == max_metric) { // 如果当前路径度量与最大路径度量相等，则选择路径值较小的路径作为最优路径new_max_metric = path_metric[i]; // 记录新的最大路径度量new_max_path = i; // 记录新的最大路径度量对应的路} else if (path_metric[i] < max_metric && path_metric[i] > new_max_metric) { // 如果当前路径度量比新记录的最大路径度量要小，但是比之前的最大路径度量要大，则更新新的最大路径度量和对应的路径值new_max_metric = path_metric[i]; // 更新新的最大路径度量new_max_path = i; // 更新新的最大路径度量对应的路径值} else if (path_metric[i] < max_metric && path_metric[i] < new_max_metric) { // 如果当前路径度量比新记录的最大路径度量要小，但是比之前的最大路径度量要小，则更新新的最大路径度量和对应的路径值，同时更新最优路径为新记录的最大路径对应的路径值和对应的分支值new_max_metric = new_max_metric; // 更新新的最大路径度量不变new_max_path = i; // 更新新的最大路径度量对应的路径值为当前路径。
c语言编写卷积计算的串行版本
文章标题：深入探讨C语言编写的卷积计算串行版本导言在计算机科学和工程领域，卷积计算是一种重要的数学运算，被广泛应用于信号处理、图像处理、神经网络等领域。
而在实际应用中，为了提高计算效率和节约资源，我们常常需要编写高效的卷积计算代码。
本文将深入探讨C语言编写的卷积计算串行版本，旨在帮助读者更深入地理解卷积计算及其在C语言中的实现方式。
一、卷积计算概述卷积计算是一种数学运算，用于描述两个函数之间的关系。
在信号处理中，卷积计算常用于滤波、特征提取等应用；在神经网络中，卷积计算则被广泛应用于图像识别、语音识别等领域。
其数学定义如下：\[ (f*g)[n] = \sum_{m=-\infty}^{\infty} f[m]g[n-m] \]其中，\( f[n] \) 和 \( g[n] \) 分别表示两个函数，\( (f*g)[n] \) 表示它们的卷积。
二、C语言中的卷积计算在C语言中，我们可以使用数组和循环来实现卷积计算。
我们需要定义两个数组来表示函数 \( f[n] \) 和 \( g[n] \)；我们可以通过嵌套循环来计算它们的卷积。
具体实现代码如下：```c// 定义函数f和g的数组int f[] = {1, 2, 3, 4, 5};int g[] = {1, 1, 1};// 定义卷积计算的结果数组int result[5] = {0};// 计算卷积for (int n = 0; n < 5; n++) {for (int m = 0; m < 3; m++) {if (n - m >= 0 && n - m < 5) {result[n] += f[n - m] * g[m];}}}```通过以上代码，我们就可以实现了一个简单的卷积计算。
但是，在实际应用中，为了提高计算效率，我们还可以进行一些优化。
三、优化卷积计算的串行版本我们可以注意到，卷积计算中存在大量的重复计算。
matlab的conv的c源代码
matlab的conv的c源代码MATLAB的conv函数是一个非常常用的函数，用于计算两个向量的卷积。
在MATLAB中，conv函数是通过C语言编写的，下面是conv函数的C源代码。
```c#include <stdio.h>#include <stdlib.h>void conv(double *x, int len_x, double *h, int len_h, double *y) {int len_y = len_x + len_h - 1; // 计算卷积结果的长度// 分配内存空间double *temp = (double *)malloc(len_y * sizeof(double));// 初始化卷积结果为0for (int i = 0; i < len_y; i++) {temp[i] = 0;}// 计算卷积for (int i = 0; i < len_x; i++) {for (int j = 0; j < len_h; j++) {temp[i + j] += x[i] * h[j];}}// 将卷积结果复制到输出数组for (int i = 0; i < len_y; i++) {y[i] = temp[i];}// 释放内存空间free(temp);}int main() {double x[] = {1, 2, 3, 4, 5};double h[] = {1, 1, 1};int len_x = sizeof(x) / sizeof(x[0]);int len_h = sizeof(h) / sizeof(h[0]);int len_y = len_x + len_h - 1;double *y = (double *)malloc(len_y * sizeof(double)); conv(x, len_x, h, len_h, y);// 打印卷积结果for (int i = 0; i < len_y; i++) {printf("%f ", y[i]);}printf("\n");// 释放内存空间free(y);return 0;}```以上是conv函数的C源代码。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
#include<iostream>using namespace std;int table1[8]={1,2,4,8,16,32,64,128};int myn=0;int stalen=0;int stan0[256][2]={0};//输入0时个状态的输出int stan1[256][2]={0};//输入1时各状态的输出int stachn[256][2]={0};//状态装换表int myg1[10]={0};int myg2[10]={0};int myout[100]; //int myoutsym=0;void chartobits(char ch,int *bits);char bitstochar(int *bits);void convolution(void);void creatsta(void);void myinput(void);int main(){char exit_char;myinput();creatsta();convolution();cin>>exit_char;}void myinput(void){int i,j;cout<<"输入编码的约束长度N:(3<N<9)"<<endl;cin>>myn;stalen=int(pow(2.0,myn-1));cout<<"选择默认的编码矢量则输入1，输入2则可输入其他的编码矢量"<<endl;cin>>i;if(i==1){switch(myn){case 3:myg1[0]=1,myg1[1]=1,myg1[2]=1;myg2[0]=1,myg2[1]=0,myg2[2]=1;break;case 4:myg1[0]=1,myg1[1]=1,myg1[2]=1,myg1[3]=1;myg2[0]=1,myg2[1]=0,myg2[2]=1,myg2[3]=1;break;case 5:myg1[0]=1,myg1[1]=0,myg1[2]=1,myg1[3]=1,myg1[4]=1;myg2[0]=1,myg2[1]=1,myg2[2]=0,myg2[3]=1,myg2[4]=1;break;case 6:myg1[0]=1,myg1[1]=0,myg1[2]=1,myg1[3]=1,myg1[4]=1,myg1[5]=1;myg2[0]=1,myg2[1]=1,myg2[2]=0,myg2[3]=1,myg2[4]=0,myg2[5]=1;break;case 7: myg1[0]=1,myg1[1]=0,myg1[2]=0,myg1[3]=1,myg1[4]=1,myg1[5]=1,myg1[6]=1;myg2[0]=1,myg2[1]=1,myg2[2]=0,myg2[3]=1,myg2[4]=1,myg2[5]=0,myg2[6]=1;break;case 8: myg1[0]=1,myg1[1]=0,myg1[2]=0,myg1[3]=1,myg1[4]=1,myg1[5]=1,myg1[6]=1,myg1[7]=1;myg2[0]=1,myg2[1]=1,myg2[2]=1,myg2[3]=0,myg2[4]=0,myg2[5]=1,myg2[6]=0,myg2[7]=1;break;case 9: myg1[0]=1,myg1[1]=1,myg1[2]=0,myg1[3]=1,myg1[4]=0,myg1[5]=1,myg1[6]=1,myg1[7]=1,m yg1[8]=1;myg2[0]=1,myg2[1]=0,myg2[2]=0,myg2[3]=0,myg2[4]=1,myg2[5]=1,myg2[6]=1,myg2[7]=0,m yg2[8]=1;break;}}else{cout<<"输入g1"<<endl;for(j=0;j<myn;j++)cin>>myg1[j];cout<<"输入g2"<<endl;for(j=0;j<myn;j++)cin>>myg2[j];}cout<<"连接矢量1为"<<endl;for(j=0;j<myn;j++)cout<<myg1[j]<<" ";cout<<endl;cout<<"连接矢量2为"<<endl;for(j=0;j<myn;j++)cout<<myg2[j]<<" ";cout<<endl;cout<<"stalen: "<<stalen;cout<<endl;}void creatsta(void){int i,j,k,myi,myj;int tembits[10];for(i=0;i<stalen;i++){stan1[i][0]=0;stan1[i][1]=0;stan0[i][0]=0;stan0[i][1]=0;stachn[i][0]=i/2;myi=i;for(j=0;j<myn;j++){if(myi>=pow(2.0,myn-1-j)){tembits[j]=1;myi=myi-pow(2.0,myn-1-j);}else{tembits[j]=0;}}/*for(k=0;k<myn;k++)cout<<tembits[k]<<" ";cout<<endl;*/for(k=0;k<myn;k++){stan0[i][0]+=myg1[k]*tembits[k];stan0[i][1]+=myg2[k]*tembits[k];}stan0[i][0]=stan0[i][0]%2;stan0[i][1]=stan0[i][1]%2;myi=i+int(pow(2.0,myn-1));stachn[i][1]=myi/2;for(j=0;j<myn;j++){if(myi>=pow(2.0,myn-1-j)){tembits[j]=1;myi=myi-pow(2.0,myn-1-j);}else{tembits[j]=0;}}/*for(k=0;k<myn;k++)cout<<tembits[k]<<" ";cout<<endl;*/for(k=0;k<myn;k++){stan1[i][0]+=myg1[k]*tembits[k];stan1[i][1]+=myg2[k]*tembits[k];}stan1[i][0]=stan1[i][0]%2;stan1[i][1]=stan1[i][1]%2;}cout<<"状态转移出"<<endl;for(i=0;i<stalen;i++)cout<<stachn[i][0]<<","<<stachn[i][1]<<" ";cout<<endl;cout<<"输入0状态转移后的输出"<<endl;for(i=0;i<stalen;i++)cout<<stan0[i][0]<<","<<stan0[i][1]<<" ";cout<<endl;cout<<"输入1状态转移后的输出"<<endl;for(i=0;i<stalen;i++)cout<<stan1[i][0]<<","<<stan1[i][1]<<" ";cout<<endl;}void chartobits(char ch,int *bits){int i;for(i=0;i<8;i++){if(ch<0)bits[i]=1;elsebits[i]=0;ch=ch<<1;}}char bitstochar(int *bits){char temp=0;int i;for(i=0;i<8;i++){if(bits[i]==1)temp+=table1[7-i];}return temp;}void convolution(){FILE *fp_input,*fp_output;if(!(fp_input=fopen("D://input.txt","r"))==1){cout<<"failed to open input_file"<<endl;exit(0);}elsecout<<"we opened the input_file "<<endl;if(!(fp_output=fopen("D://output.txt","w+"))==1){cout<<"failed to open output_file"<<endl;exit(0);}elsecout<<"we opened the output_file "<<endl;char ch;int i,j;int mybits[8],mytembits[8];int mysta=0;int wcout;char wch;for(ch=fgetc(fp_input);feof(fp_input)==0;ch=fgetc(fp_input)){ chartobits(ch,mybits);/*cout<<"输入数据为"<<endl;for(i=0;i<8;i++)cout<<mybits[i]<<" ";cout<<endl;*/for(i=0;i<8;i++){if(mybits[i]==0){myout[myoutsym++]=stan0[mysta][0];myout[myoutsym++]=stan0[mysta][1];mysta=stachn[mysta][0];}else{myout[myoutsym++]=stan1[mysta][0];myout[myoutsym++]=stan1[mysta][1];mysta=stachn[mysta][1];}}/*cout<<"输出数据1为"<<endl;for(temi=0;temi<myoutsym;temi++)cout<<myout[temi]<<" ";cout<<endl;*/if(myoutsym>=8){wcout=myoutsym/8;for(i=0;i<wcout;i++){for(j=0;j<8;j++)mytembits[j]=myout[8*i+j];wch=bitstochar(mytembits);fputc(wch,fp_output);/*cout<<"输出数据2为"<<endl;cout<<wch<<" ";*/}for(i=0;i<myoutsym-wcout*8;i++)myout[i]=myout[wcout*8+i];myoutsym=myoutsym-wcout*8;}}for(i=0;i<myn-1;i++){myout[myoutsym++]=stan0[mysta][0];myout[myoutsym++]=stan0[mysta][1];mysta=stachn[mysta][0];}wcout=myoutsym/8;for(i=0;i<wcout;i++){for(j=0;j<8;j++)mytembits[j]=myout[8*i+j];wch=bitstochar(mytembits);fputc(wch,fp_output);}for(i=0;i<myoutsym-wcout*8;i++)myout[i]=myout[wcout*8+i];myoutsym=myoutsym-wcout*8;if(myoutsym!=0){for(i=0;i<8-myoutsym;i++)myout[myoutsym++]=0;}for(j=0;j<8;j++)mytembits[j]=myout[j];wch=bitstochar(mytembits);fputc(wch,fp_output);fclose(fp_input);fclose(fp_output);}。