文档介绍：第5章 MATLAB程序设计
脚本文件和函数文件
M文件有两种形式:M脚本文件和M函数文件。
1、M文本编辑器
2、 M文件的基本形式
M脚本文件(Script File)
M函数文件(Function File)
3、 M脚本文件
脚本文件的特点:
(1) 脚本文件中的命令格式和前后位置,与在命令窗口中输入的没有任何区别。
(2) MATLAB在运行脚本文件时,只是简单地按顺序从文件中读取一条条命令,送到MATLAB命令窗口中去执行。
(3) 与在命令窗口中直接运行命令一样,脚本文件运行产生的变量都是驻留在MATLAB的工作空间(workspace)中,可以很方便地查看变量,除非用clear命令清除。
4、 M函数文件
函数文件的特点:
(1) 第一行总是以“function”引导的函数声明行;
 函数声明行的格式:
function [输出变量列表] = 函数名(输入变量列表) 
(2) 函数文件在运行过程中产生的变量都存放在函数本身的工作空间;
(3) 当文件执行完最后一条命令或遇到“return”命令时,就结束函数文件的运行,同时函数工作空间的变量就被清除;
(4) 函数的工作空间随具体的M函数文件调用而产生,随调用结束而删除,是独立的、临时的,在MATLAB运行过程中可以产生任意多个临时的函数空间。
【】根据不同阻尼系数,编写M脚本文件和M函数文件分别绘制二阶系统的时域曲线。欠阻尼系统的时域输出y与x的关系为:
程序流程控制
1、 for ... end循环结构
语法: for 循环变量= array (表达式)
循环体
end
说明:array可以是向量也可以是矩阵,循环变量依次取array的各列,每取一次循环体执行一次。循环体被循环执行,执行的次数就是array的列数。
【】使用for ... end循环的array向量编程求出 1+3+5...+99 的值。
% EX0501 使用向量for循环
sum=0;
for n=1:2:99
sum=sum+n;
end
sum
【】与【】相同,计算1+3+5...+99 的值。
% 使用while循环
sum=0;
n=1;
while n<=99
sum=sum+n;
n=n+2 ;
end
sum
n
sum =
2500
n =
101
3、 If…else…end条件转移结构
语法:
if 条件式1
   语句段1
elseif 条件式2
  语句段2
...
else
   语句段n+1 
end
说明:当有多个条件时,条件式1为假再判断elseif的条件式2,如果所有的条件式都不满足,则执行else的语句段n+1,当条件式为真则执行相应的语句段;If…else…end结构也可以是没有elseif和else的简单结构。
条件转移结构框图
if
elseif
elseif
else
【】用If结构执行二阶系统时域响应。根据阻尼系数0<ζ<1和ζ=1两种情况,得出不同的时域响应表达式:
function y=Ex0503(zeta)
% EX0503 使用if结构的二阶系统时域响应
x=0::20;
if (zeta>0)&(zeta<1)
y=1-1/sqrt(1-zeta^2)*exp(-zeta*x).*sin(sqrt... (1-zeta^2)*x+acos(zeta));
elseif zeta==1
y=1-exp(-x).*(1+x);
end
plot(x,y)