MATLAB实践教程 教学课件 作者 王玉顺 第1 5章 第5章.pptVIP

5.1 MATLAB程序结构及其运行机制 5.2 MATLAB程序的流程控制 5.3 编写MATLAB函数 5.4 编程案例 　　MATLAB程序由主程序和函数(亦称子程序)两部分组成，如图5-1所示。编程采用模块化程序设计方法，用户任务分解为多个子任务，子任务编写成用户自行命名的MATLAB函数。主程序设计成若干个模块，每个模块实现特定的数据处理目标，模块调用MATLAB函数完成子任务，多个模块顺序衔接构成主程序，在Command Window中键入主程序的存盘文件名则可执行和完成全部任务。 图5-1 MATLAB程序的结构与数据流 5.2.1 for循环　　(1) 利用循环向量化技术设计for循环，计算数列｛1 3 5 …｝的前15项累加和与累乘积。程序如下： clc; close all; clear all;x=0; y=1; %赋初值for k=1:2:2*15-1; %创建循环变量x=x+k; y=y*k; %计算累加和与累乘积end;format short e; Results=[x y] %计算结果　　程序执行结果如下：Results = 2.2500e+002 6.1903e+015 　　(2) 利用阵列预分配技术设计for循环，计算下面离散系统在正弦输入下的响应yk。 　　程序如下:clc; close all; clear all; tic %初始化和启动秒表计时器space=0.001; t=0:space:16; u=sin(2*pi*t); %计算正弦输入n= length(u); %计算u的数据个数y=zeros(n, 1); %按u的数据个数n、实施阵列y的预分配y(1)=2*u(1); y(2)=0.75*y(1)+2*u(2); %计算前2个不能递推计算的响应for k=3:n; %创建循环变量y(k)=0.75*y(k-1)-0.125*y(k-2)+2*u(k); %用递推公式计算响应yend;Time=toc %测定程序运行时间plot(t,y) %响应y的计算结果展示　　程序执行结果如下，图形如图5-2所示。Time = 5.1431e-003 图5-2 离散系统的时间响应y(t) 5.2.2 while循环　　while循环程序设计应注意两点，一是循环变量累加，二是循环变量关系表达式，两者构成循环控制逻辑。 用关系表达式设计while循环，计算函数 　 　　在区间　　　上的面积。程序如下：clc; close all; clear all; tic %初始化和启动秒表计时器space=pi/30; theta=0:space:pi; y=sin(theta); %生成被积函数数据x=0; k=1; %给面积变量x和循环变量k赋初值while k=length(y); %用下标k作循环变量，设计关系表达式%while klength(y)+1; %可选的循环变量关系表达式 x=x+y(k)*space; %累计计算面积和k=k+1; %循环变量增1累加end;time＝tocdisp(k=), disp(k)disp(aera=), disp(x)　　程序执行结果：time = 0.0065k = 32aera= 1.9982 5.2.3 if-else-end结构　　if-else-end结构关系表达式根据检验结果有条件地执行程序。　　(1) 编写一个MATLAB函数用于判断一个小于50的数是否为素数，并存盘为prime.m。程序如下：function y=prime(x)if rem(x,2)==0x2; y=不是素数,;elseif rem(x,3)==0x3; y=不是素数,;elseif rem(x,5)==0x5; y=不是素数,;elseif rem(x,7)==0x7; y=不是素数,;else; y=是素数，;end; 　　编写主程序指定要判断的多个数字(向量x)，存盘名为myprime01.m。程序如下：clc; close all; clear all;x=[1 2 3 4 5 6 7 8 9 10 11 12 13 14 15]; %指定要判别的数字n=length(x); p=[]; %判别结果变量p赋初值for i=1:n; %设计循环逐个判别p=[p [num2s