z=peaks(x,y); %定义内差点
[xi,yi]=meshgrid(-5:0.25:5); %利用四种方法进行内插 %最邻近内插
zi1=interp2(x,y,z,xi,yi,'nearest'); %线性内插
zi2=interp2(x,y,z,xi,yi,'linear'); %三次样条内插
zi3=interp2(x,y,z,xi,yi,'spline'); %三次曲线内插
zi4=interp2(x,y,z,xi,yi,'cubic'); %绘出曲面图进行比较 figure(1) subplot(2,3,1)
surf(x,y,z),title('原始数据曲线') subplot(2,3,2)
surf(xi,yi,zi1),title('最邻近内插') subplot(2,3,3)
surf(xi,yi,zi2),title('线性内插') subplot(2,3,4)
surf(xi,yi,zi3),title('三次样条内插') subplot(2,3,5)
surf(xi,yi,zi4),title('三次曲线内插') %绘制轮廓图进行比较 figure(2) subplot(2,3,1)
contour(x,y,z),title('原始数据曲线') subplot(2,3,2)
contour(xi,yi,zi1),title('最邻近内插')
subplot(2,3,3)
contour(xi,yi,zi2),title('线性内插') subplot(2,3,4)
contour(xi,yi,zi3),title('三次样条内插') subplot(2,3,5)
contour(xi,yi,zi4),title('三次曲线内插') 16.%产生三个信号：
%x1＝sin(kt)+randn(size(t)) %x2＝cos(kt)+randn(size(t)) %x3＝sin(kt)+randn(size(t))
%试计算x1与x2，x1与x3之间的相关系数，从中可得出什么结论? %如果信号不含正余弦信号分量，结论有如何？ %程序如下： x0=2*pi*20; t=0:1/119:1;
x1=sin(x0*t)+randn(size(t)); x2=cos(x0*t)+randn(size(t)); x3=sin(x0*t)+randn(size(t)); y1=corrcoef(x1,x2) y2=corrcoef(x1,x3) %该程序尚未完成 20.微分方程组 %x1(t)=0.5-x1(t) %x2(t)=x1(t)-4*x2(t)
%当t＝0时,x1(0)=1,x2(0)=-0.5,求微分方程在t∈[0,25]上的解， %并画出x1－x2的系统轨迹。 %程序如下： %含有参变量函数 X0=[1;-0.5];
tspan=[0,25];
[T,X]=ode45('char5_200',tspan,X0); figure(1)
subplot(2,1,1),plot(T,X(:,1),'r'),title('x_{1}'),gridon subplot(2,1,2),plot(T,X(:,2),'k'),title('x_{2}'),gridon figure(2)
plot(X(:,1),X(:,2)),title('系统轨迹'),gridon xlabel('X_{1}'),ylabel('x_{2}') %该程序尚未完成 %微分方程组 %x1(t)=0.5-x1(t) %x2(t)=x1(t)-4*x2(t)
%当t＝0时,x1(0)=1,x2(0)=-0.5,求微分方程在t∈[0,25]上的解， %并画出x1－x2的系统轨迹。 %程序如下：
functiondy=char5_200(t,y) dy=[0.5-y(1);y(1)-4*y(2)];
实验五数字通信信道编译码
一、 实验目的
在前面熟悉matlab语言编程的基础上，接触学习matlab在电子通信领域的应用，特别是matlab语言在信道编译码算法仿真方面的强大优势
二、实验内容
1．数字通信系统信道编码AMI编译码 （一）编码规则
消息代码中的0传输码中的0
消息代码中的1传输码中的+1、-1交替 例如:
消息代码:1010100010111 AMI码:+10-10+1000-10+1-1+1 （二）AMI码的特点:
①由AMI码确定的基带信号中正负脉冲交替，而0电位保持不变;所以由AMI码确定的基带信号无直流分量，且只有很小的低频分量; ②不易提取定时信号，由于它可能出现长的连0串。 （三）解码规则
从收到的符号序列中将所有的-1变换成+1后，就可以得到原消息代码 2．数字通信系统信道编码HDB3编译码
HDB3码即三阶高密度双极性码。 （一）编码规则:
1，先将消息代码变换成AMI码，若AMI码中连0的个数小于4,此时的AMI码就是
HDB3码;
2，若AMI码中连0的个数大于4,则将每4个连0小段的第4个0变换成与前一个非0符号(+1或-1)同极性的符号,用
表示(+1
+
,-1
-);
3，为了不破坏极性交替反转，当相邻V符号之间有偶数个非0符号时，再将该小段
的第1个0变换成＋B或-B,B符号的极性与前一非符号的相反，并让后面的非零符号从V符号开始再交替变化。 相邻的v和v之间要不同极性 而v总和前面的非0同极性 包括B 所以很好找到v （破坏点） 例如:
消息代码:100001000011000011 AMI码: +10000-10000+1-10000+1-1 HDB3码：+1000+V-1000-V+1-1+B00+V-1+1
二、HDB3码的特点:
1，由HDB3码确定的基带信号无直流分量，且只有很小的低频分量; 2，HDB3中连0串的数目至多为3个，易于提取定时信号。 3，编码规则复杂，但译码较简单。
三、解码规则
1，从收到的符号序列中找到破坏极性交替的点
，可以断定
符号及其前面的3个
符号必是连0符号，从而恢复4个连码;（先把特殊规律解决，再解决非特殊规律） 2，再将所有的-1变换成+1后，就可以得到原消息代码。
三、实验要求
要求在实验前必须预习，将实验内容事先准备好，否则不允许上机。
上机过程中由指导老师检查结果后方可做其他内容。每次实验结束后完成实验报告并在下次实验之前由学委统一交给指导教师。