matlab 编程学习,matlab编程入门学习（4）

之前把matlab的一些基本知识点讲解了下，下面继续讲函数的部分

第五章、 自定义函数

5.1、简单介绍

好的编程习惯把大的程序分解成函数，有很多的好处，例如，程序部分的独立检测，代码的可复用性，避免意想不到的错误。

MATLAB 函数的基本形式如下：

function [outarg1, outarg2, ...] = fname(inarg1, inarg2, ...)

%H1 comment line

%Other comment lines

...

(Executable code)

...

(return)

简单介绍 在一个函数中的初始注释行有特定的目的。在function 语句的第一个行注释被称为H1注释行。它应当是对本函数功能的总结。这一行的重要性在于，通过lookfor 命令它能被搜索到并显示出来。从H1 注释行到第一个空行或第一个可执行性语句可以通过help 命令或帮助窗口搜索到。它们则应包含如何使用这个函数的简单总结。

5.2、 在MATLAB 中传递变量：按值传递机制

在C语言中只有按址传递，值传递的并不能改变实参，但是在matlab中却不是这样的，看下面的代码

function aa = ssort(a)

%this is a program about how to sort number

disp('Before sorting:');

fprintf('\n');

for i = 1:size(a,2)

fprintf('%d\t',a(i));

end

for i = 1:size(a,2)-1

for j=i:size(a,2)

if a(i)>a(j)

temp = a(i);

a(i) = a(j);

a(j) = temp;

end

end

end

aa = a;这个若在C中，最后的aa=a这个语句绝对不能够对数组a产生什么影响，但是在matlab中却起到了排序的作用，看test函数

sort_number = input('please input number to be sorted!!\n');

a = zeros(1,sort_number);

for i=1:sort_number

a(i) = input('please input a number:\n');

end

ss = ssort(a);

fprintf('\n');

disp('After sorting:\n');

fprintf('\n');

for i = 1:size(a,2)

fprintf('%d\t',ss(i));

end

5.3 选择性参数

模拟噪声一般由随机数发生器产生。随机数发生器是一个函数，当它每一次被调用时，将会返回一个不同的随机出现的数。事实上，这些数是由一个确定性算法产生的，它们只是表现为随机。

代码如下：

function ran = random0(n,m)

tic;

global ISEED

%this is a program about how to create random number

msg = nargchk(1,2,nargin);

error(msg);

if nargin<2

m = n;

end

ran = zeros(n,m);

for i = 1:n

for j = 1:m

ISEED = mod(8121*ISEED + 28411, 134456);

ran(i,j) = ISEED / 134456;

end

end

toc;种子代码是

function seed(new_seed)

global ISEED

msg = nargchk(1,1,nargin);

error(msg);

new_seed = round(new_seed);

ISEED = abs(new_seed);

理论上应该没有matlab自带的rand函数快，可是运行出来的时间却不是这样的，我的random0函数的时间是random0(1,200000);

Elapsed time is 0.033057 seconds.

tic

rand(1,200000);

toc

Elapsed time is 5.417915 seconds.

如果有大神知道的话，可以给你个解释，我不是很清楚，不应该比它更快才对啊！

5.4、持久内存

持久内存(Persistent Memory )是内存的一种类型，在函数上一次调用之后，这一步调用之前，本地变量被保存在持久内存，值不变。形式如下：

persistent var1 var2 var3 ...

看一下例子：

function [ave,std] = runstatx(x)

%RUNSTATS Generate running ave / std deviation

%declare persistent calues

persistent n

persistent sum_x

persistent sum_x2

%check for a legel number of input arguments

msg = nargchk(1,1,nargin);

error(msg);

%if the argument is 'reset',reset the running sums

if x == 'reset'

n = 0;

sum_x = 0;

sum_x2 = 0;

else

n = n + 1;

sum_x = sum_x + x;

sum_x2 = sum_x2 + x^2;

end

%calc ave and sd

if n==0

ave = 0;

std = 0;

elseif n==1

ave = sum_x;

std = 0;

else

ave = sum_x/n;

std = sqrt((n*sum_x2-sum_x^2)/(n*(n-1)));

end这个给出了求一系列数据的平均值和方差

脚本文件为：

[ave,std] = runstatx('reset');

navls = input('Enter number of values in data set:');

for i=1:navls

string = ['enter value' int2str(i) ':'];

x = input(string);

[ave,std] = runstatx(x);

%display running statistics

fprintf('Average = %6.2f\n; std dev = %6.2f\n',ave,std);

end

5.5、函数的函数以及其它 函数的函数是指函数的输入参数中含有其它的函数，传递给函数的函数的变量名一般情况应用于这个函数执行的过程中。

函数的函数操作的关键字有两个专门的maltab 函数，eval 和feval。。函数eval 对一个字符串进行求值，就如它在命令窗口中已经键入了一样。函数feval 用一个特定的输入值对命名的函数进行求值。

function quickplot(fun,xlim)

%QUICKPLOT Generate quick plot of a function

msg = nargchk(2,2,nargin);

error(msg);

if(size(xlim,1)==1 && size(xlim,2)==2) || (size(xlim,1)==2 && size(xlim,2)==1)

n_steps = 100;

step_size = (xlim(2)-xlim(1))/n_steps;

x = xlim(1):step_size:xlim(2);

y = feval(fun,x);

plot(x,y);

title(['\bfPlot of function ' fun '(x)']);

xlabel('\bfx');

ylabel(['\bf' fun '(x)']);

else

% Else wrong number of elements in xlim.

error('Incorrect number of elements in xlim.');

end

子函数这个用法和C语言的还是有很大的差别的，C语言中并没有子函数这一叫法，这里的子函数和C语言中的被调用函数(也叫子函数)用法差不多，但是C语言中子函数和主函数之间并不能嵌套的定义，只能嵌套的调用。Matlab并没有这个限制！！！！！！！！！

function [avg,med] = myfun(x)

%this is a program about how to show function and sub_function

n = length(x);

avg = mean(x,n);

med = median(x,n);

function a = mean(x,n)

%sub_function to calc average

a = sum(x)/n;

function b = median(x,n)

%sub_function to calc median

w = sort(x);

if rem(n,2)==1

b = w((n+1)/2);

else

b = (w(n/2)+w(n/2+1))/2;

end

OK！！基本上这一章差不多学习结束了。

总结

下面是对MATLAB 函数和命令的总结，并带有简短的描述。

nargin 这个函数返回调用这个函数时所需要的实际输入参数的个数 nargout 这个函数返回调用这个函数时所需要的实际输出参数的个数 nargchk 如要一个函数调用被调用时参数过多或过少，那么nargchk 函数将返回一个标准错误信息 error 显示错误信息，并中止函数以免它产生这个错误。如果参数错误是致命的，这个函数将会被调用。 warning 显示警告信息并继续执行函数，如果参数错误不是致命的，执行还能继续，则这个将会被调用。 inputname 这个函数将会返回对于特定参数个数的实际变量名。 rand 产生一个随机数 rand(n) 产生一个n×n 的随机数数组 rand(n,m) 产生一个n×m 的随机数数组 rand 用于产生等可能的随机数 randn 用于产生普通的随机数