信息论与编码学上机实习报告

for(i=0;i<X_num;i++)
{
for(j=0;j<Y_num;j++)
{
printf("P[%d][%d]=",i,j);
scanf("%lf",&P[i][j]);
}
}
printf("\n请输入迭代法精度误差（e）:");
scanf("%lf",&e);
for(i=0;i<X_num;i++)
double Calculate_cap_max(double a[])
{
int i;
double max_a = a[0];
for(i=0;i<X_num;i++)
{
if(a[i] > max_a)
{
max_a = a[i];
}
}
printf("较大值=%f",log(max_a));
return log(max_a);
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
int X_num,Y_num; // X_num为信源个数，Y_num为信宿个数
int n=1;
double e; //迭代法精度误差
double PXi[50]; //输入符号的概率P(xi)数组
【思路】：利用迭代法，使cap_result初步逼近cap_max，当误差小于e时，cap_result即为信道容量
**********************************************************************************************/
{
char c3[5];
DataType C2;
int i;
int flag=1,flag1=0,flag2=0,k=1;
SLNode *p;
printf("The message is:\n %s\n",c1);
printf("The codeword is:\n");
for(i=0;i<x1;i=i+flag2)
函数名：double Calculate_cap_max(double a[])
功能：计算并输出迭代法所需的参数cap_max
参数：
返回值：cap_max，也即a[k]的最大值*/
/**************************************************************************/
{
c3[0]=c1[i];
c3[1]='\0';
p=head;
while(p->next!=NULL)
{
if(strcmp(c3,p->ble)) //c3与链表中的值不相等
{
p=p->next;
flag1=0;
}
else //c3与链表中的值相等
{
c3[flag++]=c1[i+k];
double Calculate_a(int k,double PXi[])
{
int i,j;
double sum1 = 0;
for(j=0;j<Y_num;j++)
{
double sum = 0; //必须在此处令sum置0
for(i=0;i<X_num;i++) //求和PXi[i] * P[i][j]得到P(yj)
char c2[]="The cat cannot sit on the canopy of the car.";
char s1[]="\0",s2[]="\0";
ListInitiate(&head);
x1=strlen(c1);
ListInitiate(&head1);
x2=strlen(c2);
4.main()函数主要是实现初始字典的构建，即’a’-’z’和’A’-’Z’的存储。
五、程序代码
子函数：
#ifndef DICTIONARY_H_INCLUDED
#define DICTIONARY_H_INCLUDED
void dictionary(SLNode *head,int x1,char c1[],int num)
p=p->next;
}
}
#endif // DICTIONARY_H_INCLUDED
主函数：
#include <stdio.h>
#include <stdlib.h>
#include <malloc.h>
#include <string.h>
typedef struct
{
char lable[5];
2、输出：经过字典编码后生成的码字。
二、实验目的
1、理解和掌握字典编码的原理；
2、理解和掌握字典译码的原理。
三、算法设计
四、程序分析
1.用结构体存储字符串及其对应的字典编码；
2.用链表存储整个字典；
3.dictionary()子函数的作用是将当前字符同链表中的字典相比较，若存在与字典中，则向结构体的数组里存放发送信息的下一位，直到当前数组中的字符串和链表字典不相等则将该字符串存储到链表中，并将标号存储到加最后一位字符之前的字符串对应的结构体的code，更新字典。
信息论与编码上机实习
学生姓名：XX
班级：XXXXX
学号：XXXXXXXX
指导老师：XX
实验一 字典编码
一、实验题目
1、输入：
①The Lempel Ziv algorithm can compress theEnglish text by about fifty five percent.
②The cat cannot sit on the canopy of the car.
PXi[i] = 1 / ((double)X_num); //初始化输入符号的概率P(xi)，等概率预处理
函数名：double Calculate_a(int k,double PXi[])
功能：计算并输出迭代法所需的参数a[k]
参数：
返回值：a[k]，也就是P(yi|xj)*log(P(yi|xj)/P(yi))/log2对i从0到Y_num-1求和
**************************************************************************/
int code;
}DataType;
#include "list.h"
#include "dictionary.h"
int main()
{
SLNode *head,*head1;
DataType C2;
int x1,x2,num=0;
char c1[]="The Lempel Ziv algorithm can compress the English text by about fifty five percent.";
参数：
返回值：cap_result，也即是用来逼近最值的参数
**************************************************************************/
double Calculate_cap_result(double PXi[],double a[])
}
/**************************************************************************
函数名：double Calculate_cap_result(double PXi[],double a[])
功能：计算并输出迭代法所需的参数cap_result
}
for(s1[0]='a';s1[0]>='a'&&s1[0]<='z';s1[0]++,num++)
{
strcpy(ble,s1);
C2.code=num;
ListInsert(head,ListLength(head),C2);
ListInsert(head1,ListLength(head1),C2);
}
s1[0]='.';
strcpy(ble,s1);C2.code=num++;
ListInsert(head,ListLength(head1),C2);
dictionary(head,x1,c1,num);
dictionary(head1,x2,c2,num);
return 0;
}
六、
k++;
c3[flag]='\0';
flag2=flag-1;
flag1=1;
printf("%d ",p->data.code);
if(i+k>=x1) break;
continue;
}
}
if(flag1==0) //插入新元素
{
strcpy(ble,c3);