M函数的基本结构分析

studyboy

M函数由function 语句引导，基本结构为：
function [ out1, out2, out3,... ] = 函数名（in1, in2, in3, ... ）
%注释说明语句
输入、返回变量格式的检测
函数体语句

输入和返回的变量个数由 nargin 和 nargout 两个Matlab保留变量决定。返回变量多于一个，用方括号，若为1个，可以省略。
正规的变量个数检测非常有必要，特别是对于要做一个比较复杂的工程来说。
以薛定宇老师书中的生成Hilbert矩阵程序为例：

function A = myhilb(n, m)
% A=MYHILB(N,M)产生一个N行M列的Hilbert矩阵A；
% A=MYHILB(N)产生一个N行N列的Hilbert方阵A；
% MYHILB(N,M)调用格式只显示N行M列的矩阵，但不返回任何矩阵。
%
% Designed  by Pro. Dingyu XUE at 30 July,2001
% Modified by studyboy at 29 April, 2008

if nargout>1, error('Too many output arguments.');end
if nargin == 1,m=n;
    elseif  nargin ==0 | nargin>2
    error('Wrong number of input arguments.');
end

J=1:double(m);
J=J(ones(n,1),:);
I=(1:double(n))';
I=I(:,ones(1,m));
E=ones(n,m,'double');
A1=E./(I+J-1);

if nargout==1,
    A=A1;
    elseif nargout==0,
        disp(A1);
end

Tips：
1.应当养成一个良好的编程习惯，比如这里对输入变量和返回变量的检测和说明；
2.在Matlab编程中，应尽量使用向量编程，这样可以提高效率。

对比原程序：

A1=zeros(n,m);
for i=1:n
    for j=1:m
        A(i,j)=1/(i+j-1);
    end
end

若n=1000;n=500；比较计算时间：
原程序结果：Elapsed time is 3.469414 seconds.
修改后的结果：Elapsed time is 0.036210 seconds.

可以看出效率有很大的提高。
参考：
1.薛定宇 陈阳泉 基于MATLAB/Simulink的系统仿真技术与应用 清华大学出版社
2.Matlab Help

[ 本帖最后由 ChaChing 于 2009-4-29 20:28 编辑 ]

ljan765

好东西  谢谢楼主的分享