MATLAB 软件使用简介 轮廓线提取 实验2 图像轮廓线提取技术 实验3 RGB向量空间中的图像分割技术 实

解：x=-pi:pi/2:pi; ↙ %定义自变量x

y=sin(x) ↙ %求自变量x的每一个元素对应的正弦值

八、M文件

M文件有两种形式：命令文件和函数文件。它们都是由若干MATLAB语句或命令组成的文件。两种文件的扩展名都是.m。若程序为命令文件，则程序执行完以后，中间变量仍予以保留；若程序为函数文件，则程序执行完以后，中间变量被全部删除。

? 文件的操作

在MATLAB中新建M文件的操作是在命令窗口中选择File|New|M-File，然后用鼠标单击M-File，就可以打开MATLAB自带的M文件编辑器。也可以对已经存在的M文件进行编辑、存储、修改和读取。

? 命令文件

M文件有两种运行方式：一是在命令窗口直接写文件名，按Enter键; 二是在编辑窗口打开菜单Tools，再单击Run。M文件保存的路径一定要在搜索路径上,否则M文件不能运行。以下例题中如果不做特别说明，都是以第一种方式运行的。

例15：用M命令文件画出衰减振荡曲线y?et的取值范围是[0,4?]。

?t3sin3t及其它的包络线y0?e?t3。

解： 1.打开M文件编辑窗口； 2.在编辑窗口逐行写下列语句；

t=0:pi/50:4*pi; y0=exp(-t/3);

y=exp(-t/3).*sin(3*t);

plot(t,y,'-r',t,y0,':b',t,-y0,':b')

3.保存M文件，并且保存在搜索路径上，文件名为a1.m; 4.运行M文件。

? M函数文件

M函数文件的一般形式为：

function <因变量>=<函数名>(<自变量>)

M函数文件可以有多个因变量和多个自变量，当有多个因变量时用[]括起来。

例16：设可逆方阵为A，编写同时求 解：1.打开M文件编辑窗口； 2.在编辑窗口逐行写下列语句；

function [da,a2,inva,traa]=comp4(x) %M函数文件comp4.m同时求矩阵x的四个值 ú为矩阵x的行列式 ￠为矩阵x的平方 %inva为矩阵x的逆矩阵 %traa为矩阵x的转置 da=det(x) a2=x^2 inva=inv(x) traa=x'

3. 保存M函数文件，并且保存在搜索路径上，文件名为comp4.m; 4．命令窗口执行下列语句：

A=[1,2;5,8]；↙ %输入矩阵A。

comp4(A)↙ %调用comp4.m函数计算矩阵A的

A,A,A,A'2?1A,A,A,A'2?1的M函数文件。

。

九、程序设计语句

? for循环

for i=范围表达式

可执行语句

end

例17：利用for循环求1?100的整数之和。 解:(1). 建立命令文件exam1.m。

%利用for循环求1~100的整数之和 sum=0; for i=1:100

sum=sum+i; end sum

(2) 执行命令文件exam1.m。 exam1↙

例18：利用for循环找出100~200 之间的所有素数。 解: (1).建立命令文件exam2.m。

%利用for循环找出100~200 之间的所有素数 disp('100~200 之间的所有素数为：') for m=100:200

k=fix(sqrt(m)); %求m的算术平方根然后取整. for i=2:k+1

if rem(m,i)==0 %求整数m与i的余数 break; end end if i>=k+1

disp(int2str(m)) %以字符串的形式显示素数. end end

(2) 执行命令文件exam2.m。

? while循环

while 逻辑表达式 循环体语句 end

例19：利用while循环来计算1!+2!+?+50！的值。 解: (1).建立命令文件exam3.m

%利用while循环来计算1!+2!+...+50！的值 sum=0; i=1; while i<51 prd=1; j=1;

while j<=i prd=prd*j; j=j+1; end

sum=sum+prd; i=i+1; end

disp(‘1!+2!+...+50！的和为：’) sum

(2) 执行命令文件exam3.m

? if-else-end分支 (1)

if 逻辑表达式

执行语句

end

功能: 如果表达式的值为真,就执行语句,否则执行end后面的语句.

(2) if 逻辑表达式