参考-atLab3(程序设计).doc

学****文档 仅供参考
MatLab & 数学建模


第三讲 MatLab程序设计
M文件
我们前面所介绍在MATLAB所做的运算，是适合于所要计算的)
首先建立一个M函数:
function Y= fun1(x)
Y=(x^3 - 2*x^2 + x - )/(x^2 + *x - );
在指令窗口运行以下指令:
fun1(1)*fun1(2)+fun1(3)*fun1(3)
学****文档 仅供参考
ans =
例:计算阶乘函数
function output = fact(n)
if n == 1
output = 1;
return;
end
output = n*fact(n-1);
return：通常情况下，当被调函数执行完后，MATLAB会自动地把控制转至主调函数或者指定窗口。如果在被调函数中插入return命令后，可以强制matlab结束执行该函数并把控制转出。Return命令可使正在运行的函数正常退出，返回调用它的函数继续运行，经常用于函数的末尾用来正常结束函数的运行。
【例】M函数文件例如。
[]
function sa = circle(r,s)
%CIRCLE plot a circle of radii r in the line specified by s.
% r 指定半径的数值
% s 指定线色的字符串
% sa 圆面积
%
% circle(r) 利用蓝实线画半径为 r 的圆周线.
% circle(r,s) 利用串 s 指定的线色画半径为 r 的圆周线.
% sa=circle(r) 计算圆面积，并画半径为 r 的蓝色圆面.
% sa=circle(r,s) 计算圆面积，并画半径为 r 的 s 色圆面.
if nargin>2
error('输入宗量太多。');
end;
if nargin==1
s='b';
end;
clf;
t=0:pi/100:2*pi;
学****文档 仅供参考
x=r*exp(i*t);
if nargout==0
plot(x,s);
else
sa=pi*r*r;
fill(real(x),imag(x),s)
end
axis('square')
三 MATLAB控制流
for循环结构
For循环允许一组命令以固定的和预定的次数重复。For循环的一般形式是：
for x = array
{commands}
end
在for和end语句之间的{commands}按数组中的每一列执行一次。
【例】 for n=1:10
x(n)=sin(n*pi/10);
end
» x
x =
Columns 1 through 7
Columns 8 through 10
For循环可按需要嵌套。
n=0:1:10;
for i=1:11
for j=1:11
y(i)=sin(n(i));
n(j)=n(j)*10;
学****文档 仅供参考
end
end
y =
Columns 1 through 7
Columns 8 through 11
n =
+012 *
Columns 1 through 7
Columns 8 through 11
可以利用break命令跳出for循环
【例】一个简单的for循环例如。
for i=1:10; %i依次取1,2,…10,.
x(i)=i; %对每个i值，重复执行由该指令构成的循环体，
end;
x %要求显示运行后数组x的值。
x =
1 2 3 4 5 6 7 8 9 10
while循环结构
与For循环以固定次数求一组命令的值相反，While 循环以不定的次数求一组语句的值。While循环的一般形式是：
while expression
{commands}
end
只要在表达式里的所有元素