Matlab实验指导书河北大学(郑晓昆)

1、Matlab 实验指导书河北大学电子信息工程学院2004 年 1 月目录MATLAB实验教学计划2实验一 MATLAB基本操作.3实验二MATLA图形系统.5实验三MATLAB程序设计.6实验四MATLAB基本应用领域 .7实验五设计性综合实验1-数字信道编译码 .14实验六设计性综合实验2-fir滤波器设计.16MATLAB实验教学计划指导教师：郑晓昆 薛文玲 王竹毅学时数：12学时 周4学时2次实验，共3周6次实验，第79教学周, 每次实验2学时所用仪器设备：MATLAB7.0实验软件系统实验指导书：Matlab实验指导书自编实验参考书：VMATLAB7.X程序设计语言 >,楼顺天等

2、编著，西安电子科大出版社，06年5月第二版 实验项目：A, MATLAB基本操作内容：矩阵操作，基本数学函数，逻辑函数操作等；要求：循序渐进完成P83练习题1 10B, MATLAB图形系统内容：图形绘制，图形标注，对数和极坐标，坐标轴控制，颜色控制等要求：循序渐进完成P146练习题1 10C, MATLAB程序设计内容：脚本script和函数function认识，流程控制，参数交互输入，基本 程序设计技巧练习，程序调试 DEBUG等要求：循序渐进完成P184练习题1 10D, MATLAB基本应用领域内容：线性代数，多项式与内插，曲线拟合，数据分析与统计，泛函基 础等要求：循序渐进完成P14

3、6练习题1 4，619E, 设计性综合实验-数字信道编译码内容：1数字通信系统信道编码AMI编译码2数字通信系统信道编码HDB3编译码F, 设计性综合实验fir滤波器设计内容：设计一个由县冲击相应数字滤波器FIR是该滤波器能够滤出规定频率以上的信号，而该频率以下的信号不受影响。、 实验目的1. 熟悉MATLAB的使用方法及特点；学会建立MATLAB搜索路径；熟悉 MATLAB工作空间、MATLAB集成环境、命令窗口；掌握MATLAB勺通用命令、管理命令和函数、管理变量和工作空间的使用方法；2. 掌握MATLAB基本操作及矩阵基础知识，包括：输入矩阵、矩阵的转置、矩阵元素求和、矩阵下际、矩阵连接

4、、矩阵行列删除、矩阵产生和操作、逻辑和关 系运算、操作符和特殊字符、基本矩阵和矩阵操作、基本矩阵和阵列；掌握特 殊变量和常数；掌握基本数学函数。二、实验内容利用基本矩阵产生 3x3和15x8的单位阵，全1阵，全0阵, 机阵(均匀分布的随2.利用 diag()-1,1之间)，正态分布随机阵(方差4,均值1)函数和rot90()产生下列矩阵：a= 0 -7 5b= 0然后求解a阵的逆矩阵aa及b阵的特征值和对应特征向量，并利用reshape 将aa阵变换成行向量。3. 产生一均匀分布在(-5,5)随机阵(50x2)，精确到小数点后一位。4. 编程实现当a - n, n ，间隔为1°时，求

5、解正弦和余弦的值，并利用plot() 函数绘制正弦，余弦曲线。5. 利用rand函数产生(0,1)间均匀分布的10x10随机矩阵a,然后统计a中大于等于 0.6的元素个数。6. 利用randn函数产生均值为 0,方差为1的10x10正态分布随机阵，然后统计其 中大于-0.5，小于0.5的元素个数。7. 编程实现下表功能：b<=0.5b>0.5a<1显示case 1 '显示case 2 'a>=1显示case 3 '显示error '&有一矩阵a,找出矩阵中其值大于 1的元素，并将他们重新排列成列向量bo存在9. 在一保定市区9月份

6、平均气温变化测量矩阵temp_Baoding_sep中(48x30), 有奇异值(大于42。C,小于0° C),编程实现删除奇异值所在的行。10. 在给定的100x100矩阵中，删除整行内容全为0的行，删除整列内容全为 列。三、实验要求要求在实验前必须预习，将实验内容事先准备好，否则不允许上机。上机过程中由指导老师检查结果后方可做其他内容。每次实验结束后完成实验报告并在下次实验之前由学委统一交给指导教师。一、实验目的和要求Matlab 提供了强大的图形处理功能， 本次实验旨在使学生熟悉和掌握应用 Matlab 实现二维图形和三维图形的绘制和控制与表现方法。二、实验内容 1，画出对数和

7、指数函数曲线，并分别加上标题、轴标记和曲线说明。 2，将图形窗口分为两格，分别绘制正割和余割函数曲线，并加上适当的标注。 3，根据教材 3.3 节内容，循序渐进的绘制对数和极坐标系图形。 4，根据教材 3.4 节内容，绘制多峰函数和三角函数的多条曲线。 5，将图形窗口分为两个窗格分别绘制函数：y=2x+5 和 y=x 2-3x+1在 -3 ， 3 区间上的曲线，并利用 axis 调整轴刻度，使他们具有相同的缩放尺寸。6，按图 3.19 的方式显示出 autumn、bone、cool 、hot 、hsv 、 gray 、等颜色条形图。 7，有一位研究生，在一年中平均每月的费用为生活费190元，资

8、料费 33 元，电话费45元，购买衣服 42 元，其他费用 45元。请以一饼图表示出他每月的消费比例，并 分离出表示资料费用的切片。8，参照图 3.26 和图 3.27 ，画出下列函数的三维曲线和网格曲线：22z=(x-2)2+(y-1.2) 29, 参照图 3.28 ，画出下列函数的曲面及等高线：22Z=x +y +sin (xy)10, 参照图 3.29 、图 3.30 、图 3.31 画出各种大小和形状的球、柱体。 三、实验要求要求在实验前必须预习，将实验内容事先准备好，否则不允许上机。 上机过程中由指导老师检查结果后方可做其他内容。每次实验结束后完成实验 报告并在下次实验之前由学委统一

9、交给指导教师。一、实验目的熟练掌握MATLAB勺三种程序结构一顺序结构、循环结构和分支结构。掌握并善于利用MATLAB勺控制流语句命令。学会 MATLAB勺M文件的编写方法，包括命令文件和 函数文件。掌握 M文件的调试命令和调试方法。本章重点是掌握 MATLAB勺三种程序设计结构：顺序结构、循环结构和分支结构。 难点是学习MATLAB言的编程机巧。二、实验内容1， 编写 M 函数实现求一个数是否为素数，再编写一主程序（脚本文件） ，要求通过 键盘输入一个整数，然后调用判断素数函数，从而确定它是否素数。2, 编写M函数统计一数值中零的个数，然后编写脚本文件，实现统计从1 2007中零勺总个数。3

10、,编写程序计算x -3,3,字长0.01 ：(x24*x3)/23 x 1yx211x1(x24* x3)/2 1 x 34， p158 t65， p158 t86， p158 t97,有n个人围成一圈，按顺序编号。从第1个人开始报数，数到 m时该人退出，并且下一个从 1 重新开始报数。求出出圈人勺顺序。三、实验要求要求在实验前必须预习，将实验内容事先准备好，否则不允许上机。 上机过程中由指导老师检查结果后方可做其他内容。每次实验结束后完成实验报 告并在下次实验之前由学委统一交给指导教师。实验四 MATLAB 的应用一、实验目的从工程教学的角度，详细并系统地学习MATLAB在高等数学、线性代数

11、以及数据处理、数字信号处理等方面的应用，这一章是实验课最重要、最核心的部分。通过本实 验的练习，应该重点掌握如下内容：能对矩阵作多种变换和运算，包括求解矩阵的特 征值、特征向量和矩阵的对角化等，熟练掌握各类方程组的多种解法。在解方程的过 程中，注意数组运算和符号运算之间的区别和联系。了解和掌握多项式的创建和基本 运算，熟练掌握多项式的各种化简、提取和替换命令，掌握多项式因式分解和展开。 初步掌握曲线拟合的方法，学会多项式拟合和非线性最小二乘估计。在插值和样条方 面，要掌握一维插值、二维函数插值和样条函数插值的方法。熟练掌握一重和多重数 值积分的命令，以及用多项式求导法求数值微分和用 diff

12、计算差分法求数值微分。熟 练掌握本章中符号微积分应用的内容，包括符号自变量的确定、求函数的极限、对符 号表达式（符号数组和多元向量函数）求导数和微分、符号积分、符号求和等，同时 熟练掌握通过调用 taylor 命令求函数的泰勒级数展开式。熟练掌握各类常微分方程的各种求解方法和函数命令，包括数值解和符号解，了解MATLAB的ODE文件模板及其使用方法。熟练掌握数据分析函数的基础运算和有限差分，包括拉普拉斯微分算子、数值梯度、向量运算、协方差矩阵和相关阵；根据需要， 熟练掌握傅立叶变换及其逆变换。二、实验内容理论课教材 p222 课后习题 1-4、 6-11 、 13-194. 设 A=11.9

13、12.81 15.66;15.2 4.18 13.61;10.59 7.59 9.22 ,B=16.00 4.41 -10.37 -21.61;0.88 -20.04 12.86 8.56;-1.43 10.7118.815.99;-12.48 24.35 -23.9 10.34分别求出这两个矩阵的LU和QR的分解。8. 将下列多项式进行因式分解，也即计算出多项式的根。（1）: p1（x）=xA4-2*xA3-3*xA2+4*x+2（2）: p2（x）=xA4-7*xA3+5*xA2+31*x-30(3) : p3(x)=xA3-xA2-25*x+25(4) : p4(x)=-2*xA5+3*

14、xA4+xA3+5*xA2+8*x12.对函数y=10*exp(-|x|) 取x-5,-4,- 3,，3, 4, 5点的值作为粗值， 分别采 用最邻近内插，线性内插，三次样条内插和三次曲线内插方法，对-5,5 内的点进行内插,比较其结果。16. 产生三个信号：x1 = sin( kt)+ra ndn( size(t)x2 = cos(kt)+ra ndn( size(t)x3 = sin( kt)+ra ndn( size(t)试计算 x1 与 x2, x1 与 x3 之间的相关系数,从中可得出什么结论 ? 如果信号不含正余弦信号分量,结论有如何？20. 微分方程组x1(t)=0.5-x1(t

15、)x2(t)=x1(t)-4*x2(t)当t = 0时,x1(0)=1,x2(0)=-0.5,求微分方程在t 0,25上的解，并画出 x1 x2 的系统轨迹。三、实验要求 要求在实验前必须预习,将实验内容事先准备好,否则不允许上机。 上机过程中由指导老师检查结果后方可做其他内容。每次实验结束后完成实验报告并在下次实验之前由学委统一交给指导教师。四、部分参考答案4. % 设 A=11.9 12.81 15.66;15.2 4.18 13.61;10.59 7.59 9.22 ,% B=16.004.41-10.37-21.61;0.88-20.0412.868.56;-1.4310.7118.8

16、1-5.99;-12.48 24.35 -23.9 10.34%分别求出这两个矩阵的 LU和QR的分解。%程序如下：% A的LU分解A=11.9 12.81 15.66;15.2 4.18 13.61;10.59 7.59 9.22 L,U=lu(A)% B的LU分解B=16.00 4.41 -10.37 -21.61;0.88 -20.04 12.86 8.56;-1.43 10.71 18.81 -5.99;-12.48 24.35 -23.9 10.34L,U=lu(B)% A的QR分解A=11.9 12.81 15.66;15.2 4.18 13.61;10.59 7.59 9.22

17、Q,R=qr(A)% B的QR分解B=16.00 4.41 -10.37 -21.61;0.88 -20.04 12.86 8.56;-1.43 10.71 18.81 -5.99; -12.48 24.35 -23.9 10.34Q,R=qr(B)8. % 将下列多项式进行因式分解，也即计算出多项式的根。%（1）：p1(x)=xA4-2*xA3-3*xA2+4*x+2%（2）：p2(x)=xA4-7*xA3+5*xA2+31*x-30%（3）：p3(x)=xA3-xA2-25*x+25%（4）：p4(x)=-2*xA5+3*xA4+xA3+5*xA2+8*x%程序如下： p1=1 -2 -3

18、 4 2;r1=roots(p1) p2=1 -7 5 31 -30;r2=roots(p2) p3=1 -1 -25 25;r3=roots(p3) p4=-2 3 1 5 8 0;r4=roots(p4)12.% 对函数 y=10*exp（-|x|） 取 x-5,-4,-3,，3, 4, 5点的值作为粗值,% 分别采用最邻近内插，线性内插，三次样条内插和三次曲线内插方法，% 对-5,5 内的点进行内插，比较其结果。% 程序如下： %产生低分辨率峰值函数 x,y=meshgrid(-5:1:5); z=peaks(x,y);%定义内差点 xi,yi=meshgrid(-5:0.25:5);

19、%利用四种方法进行内插 %最邻近内插 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),

20、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)原始数据曲线 ')最邻近内插 ') 线性内插 ') 三次样条内插 ') 三次曲线内插 ')最邻近内插 ')线性内插 ')三次样条内插 ')三次曲线内插 ')subplot(2,

21、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)+ra ndn (size(t)% x2 = cos(kt)+ra ndn (size(t)% x3 = sin (kt)+ra ndn (size(t)% 试计算 x1 与 x2，x1 与 x3 之间的相关系数

22、，从中可得出什么结论% 如果信号不含正余弦信号分量，结论有如何？% 程序如下： 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 的系统轨迹

23、。%程序如下：%含有参变量函数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'),grid on subplot(2,1,2),plot(T,X(:,2),'k'),title('x_2'),grid onfigure(2)plot(X(:,1),X(:,2),title(' 系统轨迹 '),grid on xlabel('

24、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 的系统轨迹。%程序如下：function dy=char5_200(t,y) dy=0.5-y(1);y(1)-4*y(2);实验五 数字通信信道编译码一、 实验目的在前面熟悉 matlab 语言编程的基础上， 接触学习 matlab 在电子通信领域的应用， 特别是 matlab 语言在信道编译码算法仿真方面的强大优

25、势二、实验内容1数字通信系统信道编码AMI 编译码（一）编码规则消息代码中的 0 传输码中的 0消息代码中的 1 传输码中的 +1、 -1 交替例如:消息代码 :1 0 1 0 1 0 0 0 1 0 1 1 1AMI 码:+1 0 -1 0 +1 0 0 0 -1 0 +1 -1 +1（二）AMI 码的特点 : 由 AMI 码确定的基带信号中正负脉冲交替，而 0 电位保持不变 ; 所以由 AMI 码 确定的基带信号无直流分量，且只有很小的低频分量 ; 不易提取定时信号，由于它可能出现长的连0串。（三）解码规则从收到的符号序列中将所有的 -1 变换成 +1 后，就可以得到原消息代码2. 数字通

26、信系统信道编码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（破坏点）例如：消息代码:1 0 0 0 0 1 0 0 0 0 1 1 0 0 0 0 1 1AMI 码:+1 0 0 0 0 -1 0 0 0 0 +1 -1 0 0 0 0 +1 -1HDB3码：+1 0 0 0 +V -1 0 0 0 -V +1 -1 +B 0 0 +V -1 +1二、HDE3码的特点：1, 由HDB3码确定的基带信号无直流分量，且只有很小的低频分量；2, HDB3中连0串的数目至多为3个，易于提取定时信号。3, 编码规则复杂，但译码较简单。三、解码规则1, 从收到的符号序列中找到破坏极性交替的点，可以断定厂符号及其前面的 3个符号必是连0符号，