第四讲MATLAB-程序设计.pptxVIP

第四讲. MATLAB 程序设计MATLAB 语言简洁高效，被称为第四代编程语言。第四讲MATLAB-程序设计全文共31页，当前为第1页。一、M 文件包含 MATLAB 代码的文件称为 M 文件，扩展名为 .m。第四讲MATLAB-程序设计全文共31页，当前为第2页。M 文件的形式M 文件有两种形式：独立的 M 文件 — 称脚本文件 (Script File)可调用 M 文件 — 称函数文件 (Function File) 第四讲MATLAB-程序设计全文共31页，当前为第3页。1. 命令文件命令文件是一串指令的集合，没有输入输出参数。与在命令窗口逐行执行文件中的所有指令等价。第四讲MATLAB-程序设计全文共31页，当前为第4页。【例】绘制“花瓣”图案% 一个用来绘制花瓣 (flower petal) 图案的命令文件theta = -pi:0.01:pi;rho(1,:) = 2*sin(5*theta).^2;rho(2,:) = cos(10*theta).^3;rho(3,:) = sin(theta).^2;rho(4,:) = 5*cos(3.5*theta).^3;for i = 1:4 polar(theta,rho(i,:),r) pauseend将该文本文件命名为 petal.m，并在指令窗口运行 petal，将看到下面的图案第四讲MATLAB-程序设计全文共31页，当前为第5页。第四讲MATLAB-程序设计全文共31页，当前为第6页。2. 函数文件如果 M 文件的第一行包含 function，此文件就是函数文件。函数文件区别于脚本文件之处是： 脚本文件中的变量在执行完后保存在内存中 函数文件中的变量在执行完后将被清除 函数文件有输入输出参数，脚本文件没有第四讲MATLAB-程序设计全文共31页，当前为第7页。【例1】计算向量元素的平均值function y = average(x)% 本函数计算向量元素的平均值y = sum(x) / length(x); % 计算过程第四讲MATLAB-程序设计全文共31页，当前为第8页。几点说明 第一行执行指令的作用：指明该文件是函数文件；定义函数名、输入参数和输出参数。 M 文件的文件名必须是函数名 .m在 M 文件前的几行注释可供 help 命令查询。第四讲MATLAB-程序设计全文共31页，当前为第9页。二、数据结构【例2】字符变量的输入和检查a = this is a stringisstr(a)ans = 1其他可能的取值还有： char（字符）、double（双精度）、sparse （稀疏矩阵）…… 第四讲MATLAB-程序设计全文共31页，当前为第10页。三、程序结构MATLAB 有顺序、循环和分支 3 种基本程序结构。■循环由 while 和 for 语句实现，■分支由 if 和 switch 语句实现。第四讲MATLAB-程序设计全文共31页，当前为第11页。1. for-end 循环语法为： for i = 表达式 （表达式为一个向量） {可执行语句} end【例3】 ：利用 for 循环求1+ 3 + ??? +9 的值 sum=0; for i = 1:2:9 sum=sum+i; endsum = 25第四讲MATLAB-程序设计全文共31页，当前为第12页。for 循环的嵌套【例4】 ：利用 for 循环求1！+2！+3！+ ??? +5！的值 sum = 0; for i = 1:5 pdr = 1; for k = 1:i pdr = pdr * k; end sum = sum + pdr; endsum = 153第四讲MATLAB-程序设计全文共31页，当前为第13页。for 循环的嵌套【例5】 ：利用 for 循环求给矩阵3行四列矩阵A赋值，A的元素为行标与列标的和 A= []; for i = 1:3 for j = 1:4 A(i,j)= i+j; end endA = 2 3 4 5 3 4 5 6 4 5 6 7第四讲MATLAB-程序设计全文共31页，当前为第14页。2. while - end 循环while 循环将循环体中的语句循环执行不定次数。 语法为： while 表达式 {循环体语句} end 表达式一般是由逻辑运算和关系运算以及一般 运算组成的，以判断循环的进行和停止； 只要表达式的值非 0，继续循环； 直到表达式值为 0，循环停止。第四讲MATLAB-程序设计全文共31页，当前为第15页。【例6 】：用 while 循环求 1~100 间整数的和 sum = 0; i = 1; while i = 100 s