基于MATLAB的循环码实验报告
课程名称：信息论与编码
课程设计题目：循环码的编码和译码程序设计
指导教师：
系别：专业：
学号：姓名：
合作者
完成时间：
成绩：评阅人：
一、实验目的：
1、通过实验了解循环码的工作原理。
2、深刻理解RS 码构造、RS 编译码等相关概念和算法。
二、实验原理
1、RS 循环码编译码原理与特点
设C 使某线性分组码的码字集合，如果对任C c c c C n n ∈=--),,,(021Λ，它的循环
移位),,,(1032)
1(---=n n n c c c c C
Λ也属于C ，则称该
码为循环码。
该码在结构上有另外的限制，即一个码字任意循环移位的结果仍是一个有效码字。其特点是：（1）可以用反馈移位寄存器很容易实现编码和伴随式的计算；（2）由于循环码有很多固有的代数结构，从而可以找到各种简单使用的译码办法。
如果一个
线性码具有以下的属性，则称为循环码：如果n 元组}
,,,{110-=n c c c c Λ是子空间S 的一个码字，则经过循环移位得到的},,,{201)
1(--=n n c c c c Λ也同样是S 中的一个
码字；或者，一般来说，经过j 次循环移位后得到的},,,,,,,{11011)
(---+--=j n n j n j n j c c c c c c c ΛΛ也是
S 中的一个码字。
RS 码的编码系统是建立在比特组基础上的，即字节，而不是单个的0和1，因此它是非二进制BCH 码，这使得它处理突发错误的能力特别强。 码长：12-=m
n
信息段：t n k 2-= （t 为纠错符号数） 监督段：k n t -=2 最小码段：12+=t d
最小距离为d 的本原RS 码的生成多项式为：g(x)=(x-α)(x -α2)(x -α3)…(x -αd -2) 信息元多项式为：：m(x)=m0+m1x+m2x2+…+mk -1xk-1
循环码特点有：
1）循环码是线性分组码的一种，所以它具有线性分组的码的一般特性，且具有循环性，纠错能力强。 2）循环码是一种无权码，循环码编排的特点为相邻的两个数码之间符合卡诺中的邻接条件，即相邻数码间只有一位码元不同，因此它具有一个很好的优点是它满足邻接条件，没有瞬时错误（在数码变
换过程中，在速度上会有快有慢，中间经过其他一些数码形式，即为瞬时错误）。
3）码字的循环特性，循环码中任一许用码经过牡环移位后，所得到的码组仍然是许用码组。
对所有的i=0,1,2,……k -1，用生成多项式g(x)除n k i
x -+，有：
)()()(x b x g x a x i i i k n +=+- （2—7） 式中)(x b i 是余式，表示为：
0,1,11,)(i i k n k n i i b x b x b x b ++=----Λ （2—8） 因此，)(x b x i k n ++-是g(x)的倍式，即)(1x b x i k n ++-是码多项式，由此得到系统形式的生成矩阵为：
(2—9)
它是一个k ?n 阶的矩阵。
同样，由G ?T H =0可以得到系统形式的一致校验矩阵为：
（2—10）
已知（7，4）循环码的生成多项式和校验多项式分别为：1)(3++=x x x g ，
1)(24+++=x x x x h 。写得其生成矩阵和校验矩阵分别为：
2、编码原理：
有信息码构成信息多项式
01
1)(m x m x m k k ++=--Λ，其中最高幂次为k-1； 用k
n x -乘以信息多项式m(x)，得到的)(x m x
k
n -，最高幂次为n-1，该过程相当于把信息码
?
?
??
?
????
???=1101000011010000110100001101G ??
??
?
?????=101110001011100010111H ???????????
?=------------0,01,01,00,1,21,20,11,11,10000010001b b b b b b b b b G k n x k k k n k k k k n k ΛM M O M M O M M ΛΛΛΛ?????
???????=
------------10001
00010
,00
,20
,11,01
,21,11
,01
,21,1Λ
Λ
ΛΛM M O M M O M M ΛΛb b b b b b b b b H k k k k k n k n k k n k
（
1
-k m ，
2
-k m ，……，1m ，
m ）移位到了码字德前k 个信息位，其后是r 个全为零的监
督位；
用g(x)除)(x m x k n -得到余式r(x),其次数必小于g(x)的次数，即小于（n-k ），将此r(x)加于信息位后做监督位，即将r(x)于
)(x m x k
n -相加，得到的多项式必为一码多项式。
1）有信息码构成信息多项式m(x)=m k-1x k-1+``````m 0 其中高幂次为k-1。
2）用x n-k 乘上信息多项式m(x)，得最高幂次为n-1，做移位。 3）用g(x)除x n-k m(x)和到余式r(x)。 编码过程流程图：
3、译码原理：
1） 有接收到的y(x)计算伴了随式s(x)。 2） 根据伴随式s(x)找出对应的估值错误图样。
否
3） 计算c^(x)=y(x)+e^(x)，得估计码字。若c^(x)= c(x)，则译码正确，否则错误。
由于g(x) 的次数为n - k 次，g(x) 除E(x) 后得余式（即伴随式）的最高次数为n-k-1次，故S(x) 共有2n-k 个可能的表达式，每一个表达式对应一个错误格式。可以知道(7,4）循环码的S(x) 共有2(7-4) = 8个可能的表达式，可根据错误图样表来纠正(7,4）循环码中的一位错误。
解码过程流程图：
4、纠错能力：
由于循环码是一种线性分组码，所以其纠检错能力与线性分组码相当。而线性分组码的最小距离可用来衡量码的抗干扰能力，那么一个码的最小距离就与它的纠检错能力有关。
初始化
由R(x)确定S(x)：
S(x)=0，无误码误
由S(x)确定错误图样E(x) 纠错)()()(X R X E x c +=
存储c(x)
定理： 对于任一个),(k n 线性分组码，若要在码字内 （1） 检测个错误，要求码的最小距离1+≥e d ；
（2） 纠正个错误，要求码的最小距离12+≥t d ；
（3） 纠正个错误同时检测
个错误，则要求1++≥e t d ；
循环码的译码分检错译码与纠错译码两类。在无记忆信道上，对码字c ，差错图案e 和接收向量
r 的多项式描述为
)()()
(x e x c x r +=
定义)(x r 的伴随多项式为)(x s
1
12210))
((mod )()(--++++==r r x s x s x s s x g x r x s Λ
由于)),((mod 0)()()
(x g x g x a x c ==所以
))()(mod ()(x g x e x s =
由此可见，0)
(≠x s 则一定有差错产生，或说满足0))()(mod (≠x g x e 的差错图样
)(x e 产生，它满足0))()(mod (=x g x e 。
循环码的检错译码即是计算)(x s 并判断是否为0
三、实验分析
1、实验测试结果，包括译码结果、误码率与信噪比之间的关系、生成多项式
理想状态下，对信号随机的提取，编码器输入为000，
通过encode函数后，因为加入了监督码，信号变得复杂密集，
编码输出为00
通过译码输出为000，与编码输入一致。说明循环码的检错和纠错能力性能好。输出多项式为：g(x)=(x+a) (x+a2) (x+a3)=a6+a5x+a4x+a3x2+ a3x +a2x+x3
以randint函数重新做一个输入信号并进行编码，结果与上例相似，输入与输出一致。
由上面所有的图可以发现，编码器输入信号并不完全相同，因为对信号的提取是随机的，所以码元也是随机的，信号经过编码器后，因为要加入监督码，所以波形变得更加密集了。信号经过译码后，波形和编码器输入信号大致相同，说明循环码的检错和纠错能力可以。
信噪比与误码率的关系比较，从图中可看出，当信噪比在20以内时，误码率相对比较大，最高达到以上，而当SNR大于２０后，信噪比保持很稳定。下面是误码率的数字显示：
2、实验过程遇到的问题及解决方法
刚开始并没有注意到运用循环码时可用上简便的Matlab自带函数，一直苦恼怎么进行纠错编码及解码，然后查找资料，收集了与循环码相关的函数（部分如下：）
1）encode函数
功能：编码函数
语法：code=encode(msg,N,K,method,opt)
说明：用method指定的方法完成纠错编码。其中msg代表信息码元，是一个K列矩阵，N是编码后的码字长度；K是信息位的长度；opt是有些编码方式需要的参数。
2）decode函数
功能：译码函数
语法：msg=decode（code，N，K，method，opt1,opt2,opt3,opt4）;
说明：这个函数对接收到的码字进行译码，恢复出原始的信息，译码参数和方式必须和编码时采用的
严格相同。它对接收到的码字，按method指定的方式进行译码；opt1，…，opt4是可选项的参数。
3）cyclpoly函数
功能：生成循环码的生成多项式。
语法：p=cyclpoly（N,K）;
p=cyclpoly（N,K,fd_flag）;
说明：从p=cyclpoly（N,K）中可找到一个给定码长N和信息位长度K生成多项式p，注意不是任意给定一个多项式都可以作为生成多项式。
4）randint函数
功能：引起一致地分布的任意整数矩阵
语法：out = randint(m)
out = randint(m,n)
out = randint(m,n,rg)
out = randint(m,n,rg,state)
在进行误码率与信噪比之间的关系编程后，Matlab一直显示
其后才发现原来没有对加噪后的信号进行整形输出，于是加入for循环以四舍五入对信号进行整形设定。
for i=1:100
for a=1:k+1
if noisycode(i,a)<
noisycode(i,a) = 0;
else
noisycode(i,a) = 1;
end
end
end
四、RS码在现代通信系统中的应用
RS编码起源于1960年MIT Lincoln实验室，经历了数十年的发展，RS码成为了研究最详尽，分析最透彻，应用最广泛，研究成果最多的码类之一。
１、在井下通信中的应用：由于井下空间小，供电系统布置受到很大局限，电磁干扰现象也很严重，通过选择合适的信道编码方式是抗干扰措施之一，可有效消除干扰，而由于RSce纠正t个m位的二进制错误会符号，而不管这t\个错误会符号是连续出现的还是离散出现，因此RS很适用于存在突发错误信道中，如井下工作，对系统码率的降低和可靠性的提高起重要作用。
井下采煤机与通信系统的硬件框图
２、RS码在PDS水声通信技术的应用
水声通信技术信道是随机时变空变的，其多途扩展产生的码间干扰直接影响到了水声通信的质量。而RS是一种扩展的非二进制BCH码，具有与PDS通信相结合的优势。
３、RS 码在无线高保真音频传输系统中的应用
在无线音频传输中, 如果传输的数据出现误码, 则会在播放时出现噪声或者啸叫声, 所以想要确保得到高保真的音频, 控制无线传输的误码率是必须的, 通常采用重传和纠错编码两种方法。重传机制对无线信道的带宽要求更高, 这里采用RS 纠错编码的方法来控制误码率。
五、实验程序
1、循环码编码与解码Matlab源程序（实验以（7，4）循环码进行分析）
m = 3;
n = 2^m-1; %定义码长
k = n-m; %信息位长
msg = randint(k*4,1,2); %随机提取信号，引起一致地分布的任意整数矩阵subplot(2,2,1)
stem(msg)
title('编码器输入信号')
p=cyclpoly(n,k) %循环码生成多项式，n=7，k=4
code = encode(msg,n,k,'cyclic',p); %编码函数，对信号进行差错编码
subplot(2,2,2)
stem(code)
title('编码器输出信号')
recode=decode(code,n,k,'cyclic',p) %对信号进行译码，对接收到的码字进行译码，恢复
出原始的信息，译码参数和方式必须和编码时采用
的严格相同
subplot(2,2,3)
stem(recode)
title('译码器输出信号')
t=-1::1;
x=recode; %将recode赋值给x，并进行长度与fft设定
N=length(x);
fx=fft(x);
df=100/N;
n=0:N/2;
f=n*df;
subplot(2,2,4);
plot(f,abs(fx(n+1))*2/N); grid;
title('频谱图')
2、误码率与信噪比之间的关系程序（以（3，2）循环码进行测试）m = 2;
n = 2^m-1; %定义码长
k = n-m; %信息位长
Fs=40; %系统采样频率
Fd=1; %码速率
N=Fs/Fd;
M=2; %进制数
for SNRpBit=1:100;%信噪比
SNR=SNRpBit/log2(M);
%制造100个信息组，每组k位
msg = randint(100,k,[0,1]);
code = encode(msg,n,k,'cyclic/binary');
%加入噪声
%在已调信号中加入高斯白噪声
noisycode=awgn(code,SNR-10*log10-10*log10(N),'measured',[],'dB');
%将浮点数转化为二进制，波形整形过程
for i=1:100
for a=1:k+1
if noisycode(i,a)<
noisycode(i,a) = 0;
else
noisycode(i,a) = 1;
end
end
end
%译码
newmsg = decode(noisycode,n,k,'cyclic');
%计算误码率
[number,ratio]=biterr(newmsg,msg);
result(SNRpBit)=ratio;
disp(['The bit error rate is',num2str(ratio)]);
end
%不同信噪比下循环码经过加性高斯白噪声信道的误码率
figure(1)
stem(result);
title('循环码在不同信噪比下的误码率')
legend('误码率','*')
xlabel('信噪比');
ylabel('在加性高斯白噪声下的误码率');
六、实验感想
好像在做这个实验时，一直拖一直拖，到最后都快没时间了才查资料看例子，在这次基于MATLAB 的通信原理课程设计中，通过查找了大量的循环码理论知识和MATLAB仿真应用的书籍之后，开始调试程序，并进一步了解MATLAB编程的关键字的运用，查阅资料之后进行不断的修正和更改，得出最终的结果。
这次MATLAB课程设计不仅仅加深了对课程理论知识的了解，并更加熟悉了计算机语言，软件的应用。并让我了解，学习知识不仅仅只在课本上，还有更多更广的渠道获得更宽广的知识。
七、参考文献
[1]徐明远，邵玉斌MATLAB仿真在通信与电子工程中的应用.西安：西安电子科技大学学出版社，
[2]唐向宏，岳恒立，郑雪峰MATLAB及在电子信息类课程中的应用(第二版）.北京：电子工业出版社，
MATLAB实验报告(1-4)
信号与系统MATLAB第一次实验报告 一、实验目的 1.熟悉MATLAB软件并会简单的使用运算和简单二维图的绘制。 2.学会运用MATLAB表示常用连续时间信号的方法 3.观察并熟悉一些信号的波形和特性。 4.学会运用MATLAB进行连续信号时移、反折和尺度变换。 5.学会运用MATLAB进行连续时间微分、积分运算。 6.学会运用MATLAB进行连续信号相加、相乘运算。 7.学会运用MATLAB进行连续信号的奇偶分解。 二、实验任务 将实验书中的例题和解析看懂，并在MATLAB软件中练习例题，最终将作业完成。 三、实验内容 1.MATLAB软件基本运算入门。 1). MATLAB软件的数值计算： 算数运算 向量运算：1.向量元素要用”[ ]”括起来，元素之间可用空格、逗号分隔生成行向量，用分号分隔生成列向量。2.x=x0:step:xn.其中x0位初始值，step表示步长或者增量，xn为结束值。 矩阵运算：1.矩阵”[ ]”括起来；矩阵每一行的各个元素必须用”,”或者空格分开； 矩阵的不同行之间必须用分号”;”或者ENTER分开。2.矩阵的加法或者减法运算是将矩阵的对应元素分别进行加法或者减法的运算。3.常用的点运算包括”.*”、”./”、”.\”、”.^”等等。
举例：计算一个函数并绘制出在对应区间上对应的值。 2).MATLAB软件的符号运算：定义符号变量的语句格式为”syms 变量名” 2.MATLAB软件简单二维图形绘制 1).函数y=f(x)关于变量x的曲线绘制用语：>>plot(x,y) 2).输出多个图像表顺序：例如m和n表示在一个窗口中显示m行n列个图像，p 表示第p个区域，表达为subplot(mnp)或者subplot(m,n,p) 3).表示输出表格横轴纵轴表达范围：axis([xmax,xmin,ymax,ymin]) 4).标上横轴纵轴的字母：xlabel(‘x’),ylabel(‘y’) 5).命名图像就在subplot写在同一行或者在下一个subplot前：title(‘……’) 6).输出：grid on 举例1：
基于MATLAB的循环码实验报告
课程名称：信息论与编码 课程设计题目：循环码的编码和译码程序设计指导教师： 系别：专业： 学号：姓名： 合作者 完成时间： 成绩：评阅人：
一、实验目的： 1、通过实验了解循环码的工作原理。 2、深刻理解RS 码构造、RS 编译码等相关概念和算法。 二、实验原理 1、RS 循环码编译码原理与特点 设C 使某 线性分组码的码字集合，如果对任C c c c C n n ∈=--),,,(021 ，它的循环 移位),,,(1032)1(---=n n n c c c c C 也属于C ，则称该 码为循环码。 该码在结构上有另外的限制，即一个码字任意循环移位的结果仍是一个有效码字。其特点是：（1）可以用反馈移位寄存器很容易实现编码和伴随式的计算；（2）由于循环码有很多固有的代数结构，从而可以找到各种简单使用的译码办法。 如果一个 线性码具有以下的属性，则称为循环码：如果n 元组 },,,{110-=n c c c c 是子空间S 的一个码字，则经过循环移位得到的},,,{201)1(--=n n c c c c 也 同样是S 中的一个码字；或者，一般来说，经过j 次循环移位后得到的 },,,,,,,{11011)(---+--=j n n j n j n j c c c c c c c 也是S 中的一个码字。 RS 码的编码系统是建立在比特组基础上的，即字节，而不是单个的0和1，因此它是非二进制BCH 码，这使得它处理突发错误的能力特别强。 码长：12-=m n 信息段：t n k 2-= （t 为纠错符号数） 监督段：k n t -=2 最小码段：12+=t d 最小距离为d 的本原RS 码的生成多项式为：g(x)=(x-α)(x -α2)(x -α3)…(x -αd -2) 信息元多项式为：：m(x)=m0+m1x+m2x2+…+mk -1xk-1 循环码特点有： 1）循环码是线性分组码的一种，所以它具有线性分组的码的一般特性，且具有循环性，纠错能力强。 2）循环码是一种无权码，循环码编排的特点为相邻的两个数码之间符合卡诺中的邻接条件，即相邻数码间只有一位码元不同，因此它具有一个很好的优点是它满足邻接条件，没有瞬时错误（在数码变换过程中，在速度上会有快有慢，中间经过其他一些数码形式，即为瞬时错误）。 3）码字的循环特性，循环码中任一许用码经过牡环移位后，所得到的码组仍然是许用码组。
《MATLAB与数值分析》第一次上机实验报告
电子科技大学电子工程学院标准实验报告（实验）课程名称MATLAB与数值分析 学生姓名：李培睿 学号：2013020904026 指导教师：程建
一、实验名称 《MATLAB与数值分析》第一次上机实验 二、实验目的 1. 熟练掌握矩阵的生成、加、减、乘、除、转置、行列式、逆、范数等运算 操作。（用.m文件和Matlab函数编写一个对给定矩阵进行运算操作的程序） 2. 熟练掌握算术符号操作和基本运算操作，包括矩阵合并、向量合并、符号 转换、展开符号表达式、符号因式分解、符号表达式的化简、代数方程的符号解析解、特征多项式、函数的反函数、函数计算器、微积分、常微分方程的符号解、符号函数的画图等。（用.m文件编写进行符号因式分解和函数求反的程序） 3. 掌握Matlab函数的编写规范。 4、掌握Matlab常用的绘图处理操作，包括：基本平面图、图形注释命令、 三维曲线和面的填充、三维等高线等。（用.m文件编写在一个图形窗口上绘制正弦和余弦函数的图形，并给出充分的图形注释） 5. 熟练操作MATLAB软件平台，能利用M文件完成MATLAB的程序设计。 三、实验内容 1. 编程实现以下数列的图像，用户能输入不同的初始值以及系数。并以x， y为坐标显示图像 x(n+1) = a*x(n)-b*(y(n)-x(n)^2); y(n+1) = b*x(n)+a*(y(n)-x(n)^2) 2. 编程实现奥运5环图，允许用户输入环的直径。 3. 实现对输入任意长度向量元素的冒泡排序的升序排列。不允许使用sort 函数。 四、实验数据及结果分析 题目一： ①在Editor窗口编写函数代码如下：
实验4循环伏安法测定电极反应参数实验报告
华南师范大学实验报告 学生姓名学号2014 专业新能源材料与器件年级、班级2014 课程名称电化学实验实验项目循环伏安法测定电极反应参数实验类型□√验证□设计□综合实验时间2016年4月25日 实验指导老师吕东生实验评分
一、实验目的 1.了解循环伏安法的基本原理及应用 2. 掌握循环伏安法的实验技术和有关参数的测定方法。 二、实验原理 循环伏安法(Cyclic Voltammetry)是一种常用的电化学研究方法。该法控制电极电势以不同的速率，随时间以三角波形一次或多次反复扫描，电势范围是使电极上能交替发生不同的还原和氧化反应，并记录电流-电势曲线。根据曲线形状可以判断电极反应的可逆程度，中间体、相界吸附或新相形成的可能性，以及偶联化学反应的性质等。常用来测量电极反应参数，判断其控制步骤和反应机理，并观察整个电势扫描范围内可发生哪些反应，及其性质如何。对于一个新的电化学体系，首选的研究方法往往就是循环伏安法。该方法使用的仪器简单，操作方便，图谱解析直观，在电化学、无机化学、有机化学、生物化学等许多研究领域被广泛使用。循环伏安法通常采用三电极系统，一支工作电极（被研究物质起反应的电极），，一支参比电极，一支对电极。外加电压在工作电极和辅助电极之间，反应电流通过工作电极与辅助电极。 图1 循环伏安法测得的氧化还原曲线 正向扫描的峰电流i p 与v^0.5和C都成线性关系，对研究电极过程具有重要意义。标准 电极电势为：EΘ=（E pa +E pc ）/2。所以对可逆过程，循环伏安法是一个方便的测量标准电极 电位的方法。 三、实验器材 CHI电化学工作站；玻碳电极；铂电极；Hg/Hg2SO4电极；0.1 mol/L VO2+ + 0.1 mol/L VO2+ +3 mol/L H2SO4溶液 四、实验步骤 1. 预处理电极
MATLAB实验报告实验二
实验二 MATLAB矩阵及其运算 学号：3121003104 姓名：刘艳琳专业：电子信息工程1班日期：2014.9.20 一实验目的 1、掌握Matlab数据对象的特点以及数据的运算规则。 2、掌握Matlab中建立矩阵的方法以及矩阵处理的方法。 3、掌握Matlab分析的方法。 二实验环境 PC_Windows 7旗舰版、MATLAB 7.10 三实验内容 4、1. （1）新建一个.m文件，验证书本第15页例2-1； （2）用命令方式查看和保存代码中的所有变量；
（3）用命令方式删除所有变量； （4）用命令方式载入变量z。 2. 将x=[4/3 1.2345e-6]在以下格式符下输出：短格式、短格式e方式、长格式、长格式e方式、银行格式、十六进制格式、+格式。 短格式 短格式e 长格式
长格式e方式 银行格式 十六进制格式 3.计算下列表达式的值 （1）w=sqrt(2)*(1+0.34245*10^(-6)) （2）x=(2*pi*a+(b+c)/(pi+a*b*c)-exp(2))/(tan(b+c)+a) a=3.5;b=5;c=-9.8; （3）y=2*pi*a^2*((1-pi/4)*b-(0.8333-pi/4)*a) a=3.32;b=-7.9; （4）z=0.5*exp(2*t)*log(t+sqrt(1+t*t)) t=[2,1-3i;5,-0.65];
4. 已知A=[1 2 3 4 5 ;6 7 8 9 10;11 12 13 14 15;16 17 18 19 20]，对其进行如下操作：（1）输出A在[ 7, 10]范围内的全部元素； （2）取出A的第2，4行和第1，3，5列； （3）对矩阵A变换成向量B，B=[1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20]； （4）删除A的第2，3，4行元素； (1) (2)
电化学实验报告
电化学分析实验报告 院系：化学化工学院专业班级：学号：姓名： 同组者： 实验日期： 指导老师： 实验一：铁氰化钾在玻碳电极上的氧化还原 一、实验目的 1.掌握循环伏安扫描法。 2.学习测量峰电流和峰电位的方法。 二、实验原理 循环伏安法也是在电极上快速施加线性扫描电压，起始电压从ei开始，沿某一方向变化， 当达到某设定的终止电压em后，再反向回扫至某设定的起始电压，形成一个三角波，电压扫 描速率可以从每秒数毫伏到1v。 当溶液中存在氧化态物质ox时，它在电极上可逆地还原生成还原态物质，即 ox + ne → red;反向回扫时，在电极表面生成的还原态red则可逆地氧化成ox，即 red → ox + ne.由 此可得循环伏安法极化曲线。 在一定的溶液组成和实验条件下，峰电流与被测物质的浓度成正比。从循环伏安法图中 可以确定氧化峰峰电流ipa、还原峰峰电流ipc、氧化峰峰电位φ pa 和还原峰峰电位φpc。 对于可逆体系，氧化峰峰电流与还原峰峰电流比为:ipa/ipc =1 25℃时，氧化峰峰电位 与还原峰峰电位差为:△φ条件电位为:φ=（φpa+ φpc）/2 由这些数值可判断一个电极过程的可逆性。 =φ pa - φpc≈56/z (mv) 三、仪器与试剂 仪器:：电化学分析仪va2020, 玻碳电极、甘汞电极、铂电极。试剂：铁氰化钾标准溶 液，0.5mol/l氯化钾溶液，蒸馏水。 四、实验步骤 1、溶液的配制 移取铁氰化钾标准溶液（10mol/l）5ml于50ml的塑料杯中，加入0.5mol/l氯化钾溶液， 使溶液达到30ml 。 2、调试 （1）打开仪器、电脑，准备好玻璃电极、甘汞电极和铂电极并清洗干净。（2）双击桌 面上的valab图标。 3、选择实验方法：循环伏安法 设置参数：低电位:-100mv；高电位600mv；初始电位-100mv； 扫描速度：50mv/s；取样间隔：2mv；静止时间：1s；扫描次数：1；量程: 200μa。 4. 开始扫描：点击绿色的“三角形”。 5. 将上述体系改变扫描速度分别为10mv/s、50mv/s、100mv/s、160mv/s、200mv/s，其 他条件不变，作不同速度下的铁氰化钾溶液的循环伏安曲线，其峰值电流与扫描速度的平方 根成正比关系。 -3 五、实验数据及处理 1. 找到循环伏安曲线上对应的氧化与还原峰，然后手动做切线。
实验二 MATLAB程序设计 含实验报告
实验二 MATLAB 程序设计 一、 实验目的 1.掌握利用if 语句实现选择结构的方法。 2.掌握利用switch 语句实现多分支选择结构的方法。 3.掌握利用for 语句实现循环结构的方法。 4.掌握利用while 语句实现循环结构的方法。 5.掌握MATLAB 函数的编写及调试方法。 二、 实验的设备及条件 计算机一台（带有MATLAB7.0以上的软件环境）。 M 文件的编写： 启动MATLAB 后，点击File|New|M-File ，启动MATLAB 的程序编辑及调试器（Editor/Debugger ），编辑以下程序，点击File|Save 保存程序，注意文件名最好用英文字符。点击Debug|Run 运行程序，在命令窗口查看运行结果，程序如有错误则改正 三、 实验内容 1.编写求解方程02=++c bx ax 的根的函数（这个方程不一定为一元二次方程，因 c b a 、、的不同取值而定） ，这里应根据c b a 、、的不同取值分别处理，有输入参数提示，当0~,0,0===c b a 时应提示“为恒不等式!”。并输入几组典型值加以检验。 （提示：提示输入使用input 函数） 2.输入一个百分制成绩，要求输出成绩等级A+、A 、B 、C 、D 、E 。其中100分为A+，90分～99分为A ，80分～89分为B ，70分～79分为C ，60分～69分为D ，60分以下为E 。 要求：（1）用switch 语句实现。 （2）输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。 （提示：注意单元矩阵的用法） 3.数论中一个有趣的题目：任意一个正整数，若为偶数，则用2除之，若为奇数，则与3相乘再加上1。重复此过程，最终得到的结果为1。如： 2?1 3?10?5?16?8?4?2?1 6?3?10?5?16?8?4?2?1 运行下面的程序，按程序提示输入n=1,2,3,5,7等数来验证这一结论。 请为关键的Matlab 语句填写上相关注释，说明其含义或功能。 4. 的值，调用该函数后，
matlab第一次实验报告
Matlab第一次实验报告 2012029010010 尹康 1. 编程实现以下数列的图像，用户能输入不同的初始值以及系数。并以x，y为坐标显示图像 x(n+1) = a*x(n)-b*(y(n)-x(n)^2); y(n+1) = b*x(n)+a*(y(n)-x(n)^2) 程序代码： n=input('input the number of pionts:'); a=input('input a:'); b=input('input b:'); x=[]; y=[]; x(1)=input('input x1:'); y(1)=input('input y1:'); %输入点数、初始值以及系数for i=2:n x(i)=a*x(i-1)-b*(y(i-1)-x(i-1)^2); y(i)=a*x(i-1)+b*(y(i-1)-x(i-1)^2); %根据已输入的数据进行迭代end figure;plot(x,y,'linewidth',2) axis equal %横纵坐标等比例 text(x(1),y(1),'1st point') %标记初始点 运行结果：
心得体会及改进：在输入某些数据时，所绘曲线可能是一条折线（如：n=5，a=b=x1=1,y1=2）甚至只有一个点（如：n=5,a=b=x1=y1=1),此时可能出现曲线与坐标轴重合或无法看到点的情况，为了更清晰地展现曲线，可以使线宽适当加宽并标记初始点。 2.编程实现奥运5环图，允许用户输入环的直径。 程序代码： 函数circle： %在指定的圆心坐标处，用指定颜色、宽度的线条绘出指定半径、圆心角的弧 function f=circle(r,x,y,color,linw,alp1,alp2) alp=linspace(alp1,alp2); X=r*cos(alp)+x; Y=r*sin(alp)+y; plot(X,Y,color,'linewidth',linw) end 主程序代码： r=input('input r:');
循环伏安法实验报告(有测定电极有效面积)
循环伏安法实验 【实验目的】 学习和掌握循环伏安法的原理和实验技术。 了解可逆波的循环伏安图的特性以及测算玻碳电极的有效面积的方法。 【实验原理】 循环伏安法是在固定面积的工作电极和参比电极之间加上对称的三角波扫 描电压（如图1），记录工作电极上得到的电流与施加电位的关系曲线（如图2），即循环伏安图。从伏安图的波形、氧化还原峰电流的数值及其比值、峰电位等可以判断电极反应机理。 与汞电极相比，物质在固体电极上伏安行为的重现性差，其原因与固体电极的表面状态直接有关，因而了解固体电极表面处理的方法和衡量电极表面被净化的程度，以及测算电极有效表面积的方法，是十分重要的。一般对这类问题要根据固体电极材料不同而采取适当的方法。 对于碳电极，一般以Fe(CN) 63-/4- 的氧化还原行为作电化学探针。首先，固体 电极表面的第一步处理是进行机械研磨、抛光至镜面程度。通常用于抛光电极的 材料有金钢砂、CeO 2、ZrO 2 、MgO和α-Al 2 O 3 粉及其抛光液。抛光时总是按抛 光剂粒度降低的顺序依次进行研磨，如对新的电极表面先经金钢砂纸粗研和细磨 后，再用一定粒度的α-Al 2O 3 粉在抛光布上进行抛光。抛光后先洗去表面污物， 再移入超声水浴中清洗，每次2～3分钟，重复三次，直至清洗干净。最后用乙 醇、稀酸和水彻底洗涤，得到一个平滑光洁的、新鲜的电极表面。将处理好的碳 图2：循环伏安曲线（i—E曲线）
电极放入含一定浓度的K 3Fe(CN)6和支持电解质的水溶液中，观察其伏安曲线。如得到如图2所示的曲线，其阴、阳极峰对称，两峰的电流值相等（i pc / i pa =1），峰峰电位差ΔE p 约为70 mV （理论值约59/n mV ），即说明电极表面已处理好，否则需重新抛光，直到达到要求。 有关电极有效表面积的计算，可根据Randles-Sevcik 公式： 在25°C 时，i p =(2.69×105 )n 3/2 AD o 1/2ν1/2 C o 其中A 为电极的有效面积（cm 2 ），D o 为反应物的扩散系数(cm 2 /s)，n 为电极反应的电子转移数，ν为扫速（V/s ），C o 为反应物的浓度（mol/cm 3 ），i p 为峰电流（A ）。 【仪器和试剂】 1. CHI 660D 电化学系统，玻碳电极（d = 4mm ） 为工作电极，银/氯化银电极为参比电极，铂片电极为辅助电极； 2. 固体铁氰化钾、H 2SO 4 溶液、高纯水； 3. 100 mL 容量瓶、50 mL 烧杯、玻棒。 【实验内容】 1. 配制5 mM K 3Fe(CN)6 溶液（含0.5 M H 2SO 4），倒适量溶液至电解杯中； 2. 将玻碳电极在麂皮上用抛光粉抛光后，再用蒸馏水清洗干净； 3. 依次接上工作电极(绿)、参比电极(白)和辅助电极(红)； 4. 开启电化学系统及计算机电源开关，启动电化学程序，在菜单中依次选择Setup 、Technique 、CV 、Parameter ，输入以下参数： 5. 点击Run 开始扫描，将实验图存盘后，记录氧化还原峰电位E pc 、E pa 及峰电流I pc 、I pa ； 6. 改变扫速为0.05、0.1 和0.2 V/s ，分别作循环伏安图； 7. 将4个循环伏安图叠加比较； Init E (V) 0.8 V Segment 2 High E (V) 0.8 V Smpl Interval (V) 0.001 Low E (V) ?0.2 V Quiet Time (s) 2 Scan Rate (V/s) 0.02 V Sensitivity (A/V) 5e?5
matlab实验二实验报告及程序
昆明理工大学信息工程与自动化学院学生实验报告 （ 2011——2012 学年 第一学期 ） 课程名称：控制系统计算机辅助设计 开课实验室：信自楼 234 2011年10月28日 年级、专业、班 学号 姓名 成绩 实验项目名称 实验二 控制系统分析 指导教师 胡蓉 教师评 语 该同学是否熟悉实验内容： A.熟悉□ B.比较熟悉□ C.不熟悉□ 该同学的实验能力： A.强 □ B.中等 □ C.差 □ 该同学的实验是否达到要求 ： A.达到□ B.基本达到□ C.未达到□ 实验报告是否规范： A.规范□ B.基本规范□ C.不规范□ 实验过程是否详细记录： A.详细□ B.一般 □ C.没有 □ 注：5个A 为优，5个B 为中，介于二者间为良，5个C 为不及格，3个Ｂ以上为及格。 教师签名： 年 月 日 实验二 控制系统分析 一、 实验目的 1. 掌握如何使用Matlab 进行系统的时域分析。 2. 掌握如何使用Matlab 进行系统的频域分析。 3. 掌握如何使用Matlab 进行系统的根轨迹分析。 4. 掌握如何使用Matlab 进行系统的稳定性分析。 二、 实验内容 1.时域分析 （1）典型二阶系统传递函数为：
当ζ=0.7, ω 取2、4、6、8、10、12的单位阶跃响应。 n 程序为：>> num1=4;den1=[1,2.8,4];sys1=tf(num1,den1); >> num2=16;den2=[1,5.6,16];sys2=tf(num2,den2); >> num3=36;den3=[1,8.4,36];sys3=tf(num3,den3); >> num4=64;den4=[1,11.2,64];sys4=tf(num4,den4); >> num5=100;den5=[1,14,100];sys5=tf(num5,den5); >> num6=144;den6=[1,16.8,144];sys6=tf(num6,den6); >> step(sys1,sys2,sys3,sys4,sys5,sys6); 运行单位阶跃响应结果图为： （2）典型二阶系统传递函数为： 当ω =6, ζ取0.2、0.4、0.6、0.8、1.0、1.5、2.0的单位阶跃响应。 n
MATLAB实验报告
MATLAB程序设计语言 实 验 报 告 专业及班级：电子信息工程 姓名：王伟 学号：1107050322 日期 2013年6月20日
实验一 MATLAB 的基本使用 【一】 实验目的 1.了解MATALB 程序设计语言的基本特点，熟悉MATLAB 软件的运行环境； 2.掌握变量、函数等有关概念，掌握M 文件的创建、保存、打开的方法，初步具备将一般数学问题转化为对应计算机模型处理的能力； 3.掌握二维图形绘制的方法，并能用这些方法实现计算结果的可视化。 【二】 MATLAB 的基础知识 通过本课程的学习，应基本掌握以下的基础知识： 一. MATLAB 简介 二. MATLAB 的启动和退出 三. MATLAB 使用界面简介 四. 帮助信息的获取 五. MATLAB 的数值计算功能 六. 程序流程控制 七. M 文件 八. 函数文件 九. MATLAB 的可视化 【三】上机练习 1. 仔细预习第二部分内容，关于MATLAB 的基础知识。 2. 熟悉MATLAB 环境，将第二部分所有的例子在计算机上练习一遍 3. 已知矩阵???? ??????=??????????=123456789,987654321B A 。求A*B ，A .* B ，比较二者结果是否相同。并利用MATLAB 的内部函数求矩阵A 的大小、元素和、长度以 及最大值。 程序代码： >> A=[1 2 3;4 5 6;7 8 9]; >> B=[9 8 7;6 5 4;3 2 1]; >> A*B ans =
30 24 18 84 69 54 138 114 90 >> A.*B ans = 9 16 21 24 25 24 21 16 9 两者结果不同 >> [m,n]=size(A) m = 3 n = 3 >> b=sum(A) b = 12 15 18 >> a=length(A) a = 3 >>max(A) ans =
超级电容器实验报告
实验报告 题目C,MnO2的电化学电容特性实验姓名许树茂 学号20104016005 所在学院化学与环境学院 年级专业新能源材料与器件创新班 指导教师舒东老师 完成时间2012 年 4 月
1.【实验目的】 1. 了解超级电容器的原理； 2. 了解超级电容器的比电容的测试原理及方法； 3. 了解超级电容器双电层储能机理的特点； 4. 掌握超级电容器电极材料的制备方法； 5. 掌握利用循环伏安法及恒流充放电的测定材料比电容的测试方法。 2. 【实验原理】 超级电容器的原理 超级电容器是由两个电极插入电解质中构成。超级电容与电解电容相比,具有非常高的功率密度和实质的能量密度。尽管超级电容器储存电荷的能力比普通电容器高,但是超级电容与电解电容或者电池的结构非常相似。 图1 超级电容器的结构图 从图中可看出,超级电容器与电解电容或者电池的结构非常相似,主要差别是用到的电极材料不一样。在超级电容器里,电极基于碳材料技术,可提供非常大的表面面积。表面面积大且电荷间隔很小,使超级电容器具有很高的能量密度。大多数超级电容器的容量用法拉(F)标定,通常在1F到5,000F之间。 (1) 双电层超级电容器的工作原理 双电层电容是在电极/溶液界面通过电子或离子的定向排列造成电荷的对峙所产生的。对一个电极/溶液体系,会在电子导电的电极和离子导电的电解质溶液界面上形成双电层。当在两个电极上施加电场后,溶液中的阴、阳离子分别向正、负电极迁移,在电极表面形成双电层;撤消电场后,电极上的正负电荷与溶液中的相反电荷离子相吸引而使双电层稳定,在正负极间产生相对稳定的电位差。这时对某一电极而言,会在一定距离内（分散层）产生与电极上的
MATLAB实验报告(二) 西安邮电大学
西安邮电学院 《Matlab》 实验报告 （二） 2011- 2012 学年第 1 学期 自动化 专业： 自动0903 班级： 学号： 姓名： 2011 年10月15日
实验二 MATLAB 的基本计算 一、实验目的 1．掌握建立矩阵的方法。 2．掌握MATLAB 各种表达式的书写规则以及常用函数的使用。 3．能用MATLAB 进行基本的数组、矩阵及符号运算。 4．掌握矩阵分析的方法以及能用矩阵求逆法解线性方程组。 二、实验设备及条件 计算机一台（带有MATLAB6.5或以上版本的软件环境）。 三、实验内容 1．利用diag 等函数产生下列矩阵。 ??????????-=032570800a ???? ? ?????=8040 507 2 b 2．利用reshape 函数将1题中的a 和b 变换成行向量。
3．产生一个均匀分布在（-5，5）之间的随即矩阵（10×2），要求精确到小数点后一位。 4．已知： ???????? ??-=765 3877 34434 12A ???? ? ?? ???--=7312 033 2 1 B 求下列表达式的值：
(1) B A K *611+=和I B A K +-=12（其中I 为单位矩阵） (2) B A K *21=和B A K *.22= (3) 331^A K =和3.32^A K = (4) B A K /41=和A B K \42= (5) ],[51B A K =和]2:);],3,1([[52^B A K =
5．下面是一个线性方程组： ???? ? ?????=????????????????????52.067.095.03216/15 /14 /15/14/13/14/13/12/1x x x (1) 求方程的解 (2) 将方程右边向量元素3b 改为0.53，再求解，并比较3b 的变化和解的相对变 6．利用randn 函数产生均值为0，方差为1的6×6正态分布随机矩阵C ，然后统计C 中大于-0.3，小于0.3的元素个数t
扫描电镜实验报告
HUNAN UNIVERSITY 姓名：扫描电镜实验报告 姓名：高子琪 学号： 2
一.实验目的 1.了解扫描电镜的基本结构与原理； 2.掌握扫描电镜样品的准备与制备方法； 3.掌握扫描电镜的基本操作并上机操作拍摄二次电子像； 4.了解扫描电镜图片的分析与描述方法。 二．实验设备及样品 1.实验仪器：D5000-X衍射仪 基本组成：1）电子光学系统：电子枪、聚光镜、物镜光阑、样品室等 2）偏转系统：扫描信号发生器、扫描放大控制器、扫描偏转线圈 3）信号探测放大系统 4）图象显示和记录系统 5）真空系统 2.样品：块状铝合金 三.实验原理 1.扫描电镜成像原理 从电子枪阴极发出的电子束，经聚光镜及物镜会聚成极细的电子束（0.00025微米-25微米），在扫描线圈的作用下，电子束在样品表面作扫描，激发出二次电子和背散射电子等信号，被二次电子检测器或背散射电子检测器接收处理后在显象管上形成衬度图象。二次电子像和背反射电子反映样品表面微观形貌特征。而利用特征X射线则可以分析样品微区化学成分。 扫描电镜成像原理与闭路电视非常相似，显像管上图像的形成是靠信息的传送完成的。电子束在样品表面逐点逐行扫描，依次记录每个点的二次电子、背散射电子或X射线等信号强度，经放大后调制显像管上对应位置的光点亮度，扫描发生器所产生的同一信号又被用于驱动显像管电子束实现同步扫描，样品表面与显像管上图像保持逐点逐行一一对应的几何关系。因此，扫描电子图像所包含的信息能很好地反映样品的表面形貌。 2.X射线能谱分析原理 X射线能谱定性分析的理论基础是Moseley定律，即各元素的特征X射线频率ν的平方根与原子序数Z成线性关系。同种元素，不论其所处的物理状态或化学状态如何，所发射的特征X射线均应具有相同的能量。
MATLAB全实验报告
《数学实验》报告 实验名称 Matlab 基础知识 学院 专业班级 姓名 学号 2014年 6月
一、【实验目的】 1.认识熟悉Matlab这一软件，并在此基础上学会基本操作。 2.掌握Matlab基本操作和常用命令。 3.了解Matlab常用函数，运算符和表达式。 4.掌握Matlab工作方式和M文件的相关知识。 5.学会Matlab中矩阵和数组的运算。 二、【实验任务】 P16 第4题 编写函数文件，计算 1! n k k = ∑，并求出当k=20时表达式的值。P27第2题 矩阵A= 123 456 789 ?? ?? ?? ?? ?? ，B= 468 556 322 ?? ?? ?? ?? ?? ，计算A*B,A.*B,并比较两者的区别。 P27第3题 已知矩阵A= 52 91 ?? ?? ?? ，B= 12 92 ?? ?? ?? ，做简单的关系运算A>B,A==B,A
P27第3题 >> A=[5 2;9 1];B=[1 2;9 2]; >>A>B >>A==B >>A> (A==B)&(A
湖南大学材料化学电化学实验报告汇总
实验一 线性极化法测定金属Fe 在稀H 2SO 4中的腐蚀速度 一、基本要求 1. 掌握动电位扫描法测定电极极化曲线的原理和实验技术。通过测定Fe 在 0.1M 硫酸溶液中的极化曲线，求算Fe 的自腐蚀电位，自腐蚀电流。 2. 讨论极化曲线在金属腐蚀与防护中的应用。 二、实验原理 当金属浸于腐蚀介质时，如果金属的平衡电极电位低于介质中去极化剂（如H +或氧分子）的平衡电极电位，则金属和介质构成一个腐蚀体系，称为共轭体系。此时，金属发生阳极溶解，去极化剂发生还原。以金属铁在盐酸体系中为例： 阳极反应： Fe-2e=Fe 2+ 阴极反应： H ++2e=H 2 阳极反应的电流密度以 i a 表示， 阴极反应的速度以 i k 表示， 当体系达到稳定时，即金属处于自腐蚀状态时，i a =i k =i corr （i corr 为腐蚀电流），体系不会有净的电流积累，体系处于一稳定电位c ?。根据法拉第定律，体系通过的电流和电极上发生反应的物质的量存在严格的一一对应关系，故可阴阳极反应的电流密度代表阴阳极反应的腐蚀速度。金属自腐蚀状态的腐蚀电流密度即代表了金属的腐蚀速度。因此求得金属腐蚀电流即代表了金属的腐蚀速度。 金属处于自腐蚀状态时，外测电流为零。 极化电位与极化电流或极化电流密度之间的关系曲线称为极化曲线。极化曲线在金属腐蚀研究中有重要的意义。测量腐蚀体系的阴阳极极化曲线可以揭示腐蚀的控制因素及缓蚀剂的作用机理。在腐蚀点位附近积弱极化区的举行集会测量可以可以快速求得腐蚀速度。还可以通过极化曲线的测量获得阴极保护和阳极保护的主要参数。 在活化极化控制下，金属腐蚀速度的一般方程式为： 其中 I 为外测电流密度，i a 为金属阳极溶解的速度，i k 为去极化剂还原的速度，βa 、βk 分别为金属阳极溶解的自然对数塔菲尔斜率和去极化剂还原的自然对数 )]ex p()[ex p( k c a c corr k a i i i I β? ?β??---=-=
实验二MATLAB程序设计含实验报告
实验二 M A T L A B 程序设计 一、 实验目的 1.掌握利用if 语句实现选择结构的方法。 2.掌握利用switch 语句实现多分支选择结构的方法。 3.掌握利用for 语句实现循环结构的方法。 4.掌握利用while 语句实现循环结构的方法。 5.掌握MATLAB 函数的编写及调试方法。 二、 实验的设备及条件 计算机一台（带有MATLAB7.0以上的软件环境）。 M 文件的编写： 启动MATLAB 后，点击File|New|M-File ，启动MATLAB 的程序编辑及调试器 （Editor/Debugger ），编辑以下程序，点击File|Save 保存程序，注意文件名最好用英文字符。点击Debug|Run 运行程序，在命令窗口查看运行结果，程序如有错误则改正 三、 实验内容 1.编写求解方程02=++c bx ax 的根的函数（这个方程不一定为一元二次方程，因c b a 、、的不同取值而定），这里应根据c b a 、、的不同取值分别处理，有输入参数提示，当0~,0,0===c b a 时应提示“为恒不等式!”。并输入几组典型值加以检验。 （提示：提示输入使用input 函数） 2.输入一个百分制成绩，要求输出成绩等级A+、A 、B 、C 、D 、E 。其中100分为A+，90分～99分为A ，80分～89分为B ，70分～79分为C ，60分～69分为D ，60分以下为E 。 要求：（1）用switch 语句实现。 （2）输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。 （提示：注意单元矩阵的用法） 3.数论中一个有趣的题目：任意一个正整数，若为偶数，则用2除之，若为奇数，则与3相乘再加上1。重复此过程，最终得到的结果为1。如： 21 3105168421 63105168421 运行下面的程序，按程序提示输入n=1,2,3,5,7等数来验证这一结论。 请为关键的Matlab 语句填写上相关注释，说明其含义或功能。 4. 编写一个函数，计算下面函数的值，给出标量x 的值，调用该函数后，返回y 的值。 function [y]=myfun1(x) 选择一些数据测试你编写的函数。 5. 编写一个函数求向量x 中元素的平均值、最大值、最小值、均方根值。 %classic "3n+1" problem from number theory. while 1 n=input('Enter n,negative quits:'); if n<=0 break end a=n; while n>1 if rem(n,2)==0
浅析Matlab数学实验报告
数学实验报告 姓名: 班级: 学号: 第一次实验任务 过程: a=1+3i; b=2-i; 结果: a+b =3.0000 + 2.0000i a-b =-1.0000 + 4.0000i a*b = 5.0000 + 5.0000i a/b = -0.2000 + 1.4000i 过程: x=-4.5*pi/180; y=7.6*pi/180; 结果: sin(abs(x)+y)/sqrt(cos(abs(x+y))) =0.2098 心得:对于matlab 中的角度计算应转为弧度。 (1)过程: x=0:0.01:2*pi; y1=sin(x); y2=cos(x); y3=exp(x); y4=log(x); plot(x,y1,x,y2,x,y3,x,y4) plot(x,y1,x,y2,x,y3,x,y4) 结果: (2)过程:>> subplot(2,2,1) >> plot(x,y1) >> subplot(2,2,2) >> plot(x,y2) ./,,,,2,311b a b a b a b a i b i a ?-+-=+=计算、设有两个复数 6,7,5.4)
cos()sin(2=-=++y x y x y x ,其中、计算的图形。 下分别绘制)同一页面四个坐标系)同一坐标系下(、在( x y e y x y x y x ln ,,cos ,sin 213==== >> subplot(2,2,3) >> plot(x,y3) >> subplot(2.2.4) >> subplot(2,2,4) >> plot(x,y4) 结果: 心得:在matlab中,用subplot能够实现在同一页面输出多个坐标系的图像,应注意将它与hold on进行区别,后者为在同一坐标系中划出多条曲线。 5、随机生成一个3x3矩阵A及3x2矩阵B,计算(1)AB,(2)对B中每个元素平方后得到的矩阵C,(3)sinB,(4)A的行列式,(5)判断A是否可逆,若可逆,计算A的逆矩阵,(6)解矩阵方程AX=B,(7)矩阵A中第二行元素加1,其余元素不变,得到矩阵D,计算D。 过程:A=fix(rand(3,3).*10) ; B=fix(rand(3,3).*10);
实验报告-循环伏安法测定亚铁氰化钾
循环伏安法测定亚铁氰化钾 实验目的 (1) 学习固体电极表面的处理方法； (2) 掌握循环伏安仪的使用技术； (3) 了解扫描速率和浓度对循环伏安图的影响 实验原理 铁氰化钾离子[Fe(CN)6]3--亚铁氰化钾离子[Fe(CN)6]4-氧化还原电对的标准电极电位为 [Fe(CN)6]3- + e -= [Fe(CN)6]4- φθ= 0.36V(vs.NHE) 电极电位与电极表面活度的Nernst 方程式为 φ=φθ+ RT/Fln(C Ox /C Red ) -0.2 0.00.20.4 0.60.8 -0.0005 -0.0004-0.0003-0.0002-0.00010.0000 0.00010.00020.0003i pa i pc I /m A E /V vs.Hg 2Cl 2/Hg,Cl - 起始电位:(-0.20V) 终止电位：(0.80 V) 溶液中的溶解氧具有电活性，用通入惰性气体除去。 仪器与试剂 MEC-16多功能电化学分析仪（配有电脑机打印机）；金电极；铂丝电极；饱和甘汞电极； 容量瓶：250 mL 、100mL 各2个，25 mL 7个。 移液管：2、5、10mL 、20mL 各一支。 NaCl 溶液、K 4[Fe(CN)6]、、Al 2O 3粉末（粒径0.05 μm ） 实验步骤
1、指示电极的预处理 金电极用金相砂纸细心打磨，超声波超声清洗，蒸馏水冲洗备用。 2、溶液的配制 配制0.20 mol/L NaCl溶液250mL，再用此溶液配制0.10 mol/L的K4[Fe(CN)6]溶液100mL备用。 3、支持电解质的循环伏安图 在电解池中，放入25mL 0.2 mol·L-1 NaCl溶液，插入电极，以新处理的铂电极为工作电极，铂丝电极为辅助电极，饱和甘汞电极为参比电极，进行循环伏安仪设定，扫描速率为0.1V/s；起始电位为-0.20V，终止电位为0.80V。开始循环伏安扫描. 4、K4 [Fe(CN)6]溶液的循环伏安图 在-0.20至0.80V电位范围内，以0.1V/s的扫描速度分别作0.01 mol·L-1、0.02 mol·L-1、0.04 mol·L-1、0.06 mol·L-1、0.08 mol·L-1的K4 [Fe(CN)6]溶液(均含支持电解质NaCl浓度为0.20mol·L-1)循环伏安图 5、不同扫描速率K4 [Fe(CN)6]溶液的循环伏安图 在0.08 mol·L-1 K4 [Fe(CN)6]溶液中，以0.1V/s、0.15 V/s、0.2V/s、0.25 V/s、0.3V/s、0.35V/s，在-0.20至0.80V电位范围内扫描，做循环伏安图 数据处理 1、从K4[Fe(CN)6]溶液的循环伏安图，测量i pa、i pc值。 -1；起始电位为-0.20V，终止电位为0.80V) 2、分别以i pa和i pc对K4[Fe(CN)6]溶液浓度c作图，说明峰电流与浓度的关系。
循环码
实验、循环码编译码系统 一、 实验目的： 1、熟悉循环码的编译码原理； 2、掌握Quartus Ⅱ开发软件的运用，在该软件下熟练的运用多种输入方式完成各种电路设计的要求； 3、初步掌握VHDL 语言，能够运用该语言编写简单的程序，完成设计要求； 4、熟悉对PLD 的下载和仿真，学会观察测试结果的正确性； 5、学会运用各方面知识，设计并实现一个系统。 二、 实验要求： 使用Quartus Ⅱ软件，用m 序列发生器作为信号源设计循环码编译码，速率可自定,并在实验箱上调试出编码和译码波形，比较信号源和译码后的信号波形。 三、实验设备： Quartus II 软件、Modelsim 软件、FPGA 实验箱、微机1台、示波器1台 四、实验原理： 1、 循环码的编码 循环码最大的特点就是码字的循环特性，所谓循环特性是指：循环码中任一许用码组经过循环移位后，所得到的码组仍然是许用码组。若（1n a - 2n a -…… 1a 0a ）为一循环码组，则(2n a - 3n a -……0a 1n a -)、（3n a - 4n a -……1n a - 2n a -）、……还是许用码组。也就是说，不论是左移还是右移，也不论移多少位，仍然是许用的循环码组。表1-2给出了一种（7，3）循环码的全部码字。 可以将循环码码组用代数多项是来表示，这个多项式被称为码多项式，对于表1-2中的任一码组可以表示为： 654326543210()A x a x a x a x a x a x a x a =++++++ （1-4） 表1-2一种（7，3）循环码的全部码字
在码多项式运算中采用按模运算法则。若一任意多项式F (x )被一个n 次多项式N (x )除，得到商式Q (x )和一个次数小于n 的余式R (x )，也就是： ()() ()()() F x R x Q x N x N x =+ （1-5） 则可以写为：F (x )≡R (x )（模N (x )）。 这时，码多项式系数仍按模2运算，即只取值0和1，假设：计算x 4+x 2+1除以x 3+1的值可得： 42233 11 11 x x x x x x x ++++=+++ （1-6） 循环码的生成多项式和生成矩阵：（全0码字除外）称为生成多项式，用g (x )表示。 可以证明生成多项式g (x )具有以下特性： （1）g (x )是一个常数项为1的r=n-k 次多项式； （2）g (x )是1n x +的一个因式； （3）该循环码中其它码多项式都是g (x )的倍式。 一旦生成多项式g (x )确定以后，该循环码的生成矩阵就可以确定，进而该循环码的所有码字就可以确定。 以表1-2的（7，3）循环码为例，来构造它的生成矩阵和生成多项式，这个循环码主要参数为，n ＝7，k ＝3，r ＝4。从表中可以看到，其生成多项式可以用第1码字构造： 421()()1g x A x x x x ==+++ （1-7） 2643253242()()()()1x g x x x x x G x xg x x x x x g x x x x ???? +++???? ==+++????????+++???? （1-8） 一个较简单的系统循环码编码方法：设要产生（n ，,k ）循环码，m (x )表示信息多项式，则其次数必小于k ，而()n k x m x -?的次数必小于n ，用()n k x m x -?除以g (x )， 可得余数r (x )，r (x )的次数必小于（n-k ），将r (x )加到信息位后作监督位，就得到了系统 循环码。下面就将以上各步处理加以解释。 （1）用n k x -这一运算实际上是把信息码后附加上（n-k ）个“0”。例如，信息码为110， 它相当于2 ()m x x x =+。当n-k ＝7-3＝4时，65()n k x m x x x -?=+，它相当于1100000。而希望的到得系统循环码多项式应当是()()()n k A x x m x r x -=?+。 （2）求r (x )。由于循环码多项式A (x )都可以被g (x )整除，也就是：