实验3－MATLAB程序设计.docx

《实验3－MATLAB程序设计.docx》由会员分享，可在线阅读，更多相关《实验3－MATLAB程序设计.docx（6页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、授课班次与时间:实验 3 MATLAB程序设计班次时间课题名称:教学重点、难点和教学方法设计：重点：1 .M文件与函数的编程方法2 .程序控制结构难点：1 .无学时：2学时授课方法：实验课，以MATLAB程序演示为主，配合基本知识点讲授, 布置学生上机练习，并对学生上机练习进行指导。说明：一、教案还应包含教具、幻灯、电化教学（若使用）使用手段的说明；新课内容小结; 作业布置；后记二、课时授课计划（教案）以一次课（2学时）为单元编写，每一单元有一首页 三、教学内容、小结、作业布置、后记等书写在竖直线左边，其它内容书写右边 四、青年教师需提供板书设计（最后）课时授课计划(教案)教学主要内容：1 .

2、文件存储路径我们要将编写的M文件放置自己设定的路径下，如何存储及操作。看一个简单的求向量平均值的例子。命令窗口中输入edit,启动编辑器。x=input ( 1 input a vector:x=1);m,n=size (x);if (ml)& (nl) | ( (m=l)& (n=l)error ( 1 you must input a vector! 1);endaverage=sum(x)/length(x)学会下面这些用法/信息量有点大吧。input;size;error;if/end;=;&;errorp23/24表2-2矩阵算术运算符运算符名 称示例法则或使用说明JjlC=A+B矩

3、阵加法法则，即。(5=1(0广即力减C=.4-B矩阵减法法则，即*乘C=AB矩阵乘法法则/右除C=A B定义为线性方程组X*3=.4的解，即(X4左除CAB定义为线性方程组-4*升=5的解，即c=a =a-i*bAA. 5其中一个为标量时有定义1共枕转置B=.45是的共扰转置矩阵表2-3数组算术运算符运算符名称示例法则或使用说明数组乘C=vl.*BJ数组右除C=A. BX数组左除C=A B3 MBA数组乘解C=A B入即力1转置A：将数组的行投放成列，豆数元素不做共扼课时授课计划(教案)表2-4关系运算符运算符名称示例法则或使用说明小于AB1. A. 3都是标量，结果是或为1(真)或为0(假)

4、的标量2. 3若一个为标量，另一个为数组，标量将与数组各元素逐一比较，结果为与运算数组行列相同的数组，其中各兀素取值 或1或03. A. 8均为数组时,必须行、列数分别相同，4与3各对应元 素相比较，结果为与d或3行列相同的数组,其中各元素取值 或1或04. 一和=运算对参与比较的量同时比较实部和虚部,其他运 算只比较实部AB大于等于A恒等于A=B不等于A =B编辑完成后，保存到指定路径。单击运行，弹出对应话框如下:读一下这个对话框的内容。(当出现提示信息的时候，要认真阅读，就能明白解决问题的办法)我们是要更改文件夹还是添加路径？更改文件夹吧，简单些。2函数文件1)函数的基本概念函数的概念大家

5、应该不陌生，几乎所有的语言都会大量使用这种编程模式。除了 FPGA中的HDL语言，这个语言虽然有函数的设计功能，但一般不用，而是用模块来写。因为它不是程序，它是电路。直接看一个例子，将前面的求平均值的文件改写为函数。function y=e3_2_average(x) % get the average of x% this is my first matlab function.m, n =size (x);if (ml)&(nl)|( (m=l)&(n=l) error(1 you must input a vector! 1);endy=sum(x)/length(x);课时授课计划(教

6、案)编辑完成后，直接运行：出现错误提示信息。采用函数调用的方式调用。E3_2_average (Ii nspace (1,10)注意函数文件最后一条语句加了分号。反同值仍然会出现在命令窗口中。以上，就是函数的基本要素：函数的关键词function函数要求与文件名相同函数可以有输入参数输出参数，也可以都没有。输入help e3_2_average,看出我们在文件中写的注释语句(紧显示function下面的注释)2)输入输出参数个数Nargin, nargout分别表示函数的输出及输入参数个数。有些函数可以输入3个参数，也可以输入2个参数，就是使用了这个nargin变量看看例子：我们主要使用nar

7、gin这个变量。function y=e3_3_max(a,bz c);%get the maxim abs value of a,b,c;if nargin=2c=0;endy=abs (a),abs (b),abs (c);y=max(y);程序控制结构l)for循环求和1-10s=o；For i=l:10;S=s+i;End;步长可设置S=0;For i=l:0.5:10;课时授课计划（教案）S=s+i;End2) while 循环num=O;EPS=1;while (1+EPS)1EPS=EPS/2;num=num+l;endnum=num-lEPS=EPS*23)选择结构if exp

8、ressionstatementselseif expressionstatementselsestatementsend分段统计随机数据的个数；function y=e3_5_statics(x);y=zeros(1,4);for i=l:length(x)if x (i)3y(1)=y (1)+1;elseif x (i)2y(2)=y (2)+1;elseif x (i)1y(3)=y (3)+1;elsey(4)=y(4)+1;endend掌握这3种程序控制结构，足够数字信号处理这门课程应用了。4课堂练习课时授课计划（教案）要求：1）编写一个函数文件e3_sumpro.m,求输入数据的

9、乘积或和。Function y=e3_sumpro（x,n,c）;X为输入的数组；如果输入不是单个数据或数组，则输出提示信息：please input a vector!c为运算方法：为0时求积，否则求和。n为求取数据的乘积或和的个数，当n大于数组长度时，仅求取全部数组个数的 运算。当n小于数组长度时，仅求前n个元素的运算。当不输入n时，默认n取全部数组元素参与运算；当不输入c时，默认求数组元素的乘积。2）编写M文件e3_sort.m,利用rand产生10个随机数，完成从大到小排序， 并输出排序前后的数据。3）提交资料：A: e3_sumpro . m文件B： e3_sort.m文件C： Dl_E3.doc文件，包括两个程序清单及运行结果展示。D:文件夹名称：DI E3学号中文姓名