实验一 实验准备

1、信息论与编码实验刘美春应用数学系课程目的：1、加深对信息理论的核心内容的理解，并提高计算和解决实际问题的能力2、掌握编码算法的实现过程课程要求： 理解原理，独立编程 3、实验报告模板实验内容实验一 实验准备：matlab软件复习(4学时)实验二 构建信源数学模型 （2学时）实验三 构建传输信道 （2学时）实验四 无失真信源编码技术在数据压缩中的应用 （8学时）实验考查 （2学时）（1）窗口命令方式（2）运行以 .M 为扩展名磁盘文件%用plot函数画一个方波t=-1 0 0 1 1 3;x=0 0 1 1 0 0;plot(t,x);xlabel(t);ylabel(x(t);axis(-1

2、3 0 2);直接在命令窗口输入以上命令建一个名字为my_file.M的文件，然后在命令窗口输入文件名及回车。命令窗口输入： help+函数名例如 help plot 不需要变量的类型说明 变量名的第一个字符必须是字母 变量名长度：不超过31个字符 大写和小写的字母视为不同的字符 例如：num_students = 25 特殊变量： pi 表示圆周率，inf 表示无穷大，NaN(Not a Number)表示不定量，如0/0。 MATLAB用常规的十进制表示数值 用i或j作为后缀来表示复数的虚部 例 1.235e5表示1.235105，x=2+3jabs(x) 求复数x的模angle(x) 求

3、复数x的相角(弧度)real(x) 求复数x的实部imag(x) 求复数x的虚部conj(x) 求复数x的共轭u + 加u - 减u * 乘u / 除u 乘方u 矩阵的复共轭转置u A & B逻辑与(and)u A | B 逻辑或(or)u A 逻辑非(not) 值为0时表示逻辑假(F)，其它任何非零值表示逻辑真。u A B 大于uA = B 大于等于uA = B 等于uA = B不等于u表达式 1:10 表示产生一个行向量，它的值为 1 2 3 4 5 6 7 8 9 10u表达式 10:-2:1 表示产生一个递减的行向量，它的值为 10 8 6 4 2 用冒号:产生数组例 x=2:

4、5 产生一个数组，它的值为x(1)=2, x(2)=3, x(3)=4, x(4)=5x=linspace(0,2,11)将区间0，2均匀抽样11点作为数组x 给2维数组赋值时，用分号表示一行的结束，如：z=1 2; 3 4。用linspace产生数组MATLAB 提供了一些产生基本矩阵的函数 zeros 产生矩阵元素全为0的矩阵 ones 产生矩阵元素全为1的矩阵 rand 产生(0,1)均匀分布随机数矩阵 randn 产生正态分布随机数矩阵size(A) 返回值数组A的行数和列数（二维） length(B) 确定数组B的元素个数（一维）例：编一个绘制图示波形的函数。function y=t

5、ri(t)y= abs(t) eps s=s+1/(n*n); n=n+1; end+222312111s例例: : 计算 的值，且误差小于10-6fprintf(s=%.5fn,s)t=linspace(0,4*pi,512);plot(t,sin(t),t,cos(t),-.);title(my figure);xlabel(t);legend(sin(t),cos(t);24八、文件文件可分为两类： 文本文件 由字符和与字符的显示格式有关的控制符构成常见扩展名：“TXT” 、“BAT”、“HTM”二进制文件 二进制文件为非文本文件常见扩展名：“COM”、“EXE” 、“BMP” 、“WA

6、V”25MATLAB中基本的低级文件I/O指令：1打开和关闭文件fopen、fclose2格式读写fprintf、fscanf、fgetl、fgets3非格式读写fread、fwrite4文件定位和状态feof、fseek、ftell、ferror、frewind26文件的打开与关闭1fopen函数fidfopen(filename，permission)fid,message=fopen(filename,permssion)filename,permission,machineformat=fopen(fid)fid =fopen(all)2fclose函数status=fclose(fi

7、d) status=fclose(all)例：打开一个名为std.dat的数据文件并进行读操作，命令为： fidfopen(my.txt,r)27格式化数据的读写fprintf函数 count=fprintf(fid,format , A .)fscanf函数A, count=fscanf(fid, format, size)28例8-1计算当x=0 1时f(x)=ex的值，并将结果写入到文件my.txt中。程序 x=0:0.1:1;y=x;exp(x); %y有两行数据fid=fopen(my.txt,w);fprintf(fid,%6.2f %12.8fn,y);fclose(fid);29例8-2从上例中生成的文件my.txt中读取数据，并将结果输出到屏幕。程序 fid = fopen(my.txt,r);a,count = fscanf(fid,%f