初学Matlab第三天（结构化程式与自定义函数）

1、MATLAB语言

  matlab语言与C语言很像，但是matlab语言是解释型脚本，c语言是编译执行的语言，matlab的解释程序是看一句代码执行一次，而C语言的编译器必须编译好源文件生成目标代码，再执行整个目标代码。matlab底层是C/C++实现的。.m文件可以使用matlab内置的mcc命令编译成exe文件。

2、编写matlab语言

2.1Script Editor

新建一个脚本，在编辑器里面写下matlab语言，可以直接点运行代码（或者按下F5），在运行前会将这个脚本存储起来（文件名不能以数字开头，大小写含义不一样）。

2.2Script flow

结构化程序设计（structured programming）

Matlab提供的结构化程序

if 逻辑

while 逻辑

prod(1:n)=n!  1e100=1*10^100

练习：用while循环实现1+2+3+...+999

n=1;
s=0;
while(n<1000)
    s=s+n;
    n=n+1;
end

for 逻辑

在使用同一符号进行写计算时，别忘记清除符号里储存的数据。

去除矩阵中的零元素：

                                     a(find(a==0))=[];
预分配
每次经过 for 和 while 循环时，这些循环都会递增数据结构体的大小，这会对性能和内存的使用产生不利影响。反复重新调整数组大小往往需要MATLAB花费额外的时间来寻找更大的连续内存块，然后将数组移入这些块中。通常可以通过预分配数组所需的最大空间量来缩短代码的执行时间。

下面的代码显示了创建标量变量 x，然后在 for 循环中逐步增加 x 大小所需的时间量:

tic
x = 0;
for k = 2:1000000
   x(k) = x(k-1) + 5;
end
toc
Elapsed time is 0.301528 seconds.
如果为 x 预分配一个 1×1,000,000 的内存块并将其初始化为零，则代码的运行速度更快，这是因为无需反复为不断增长的数据结构体重新分配内存。

tic
x = zeros(1,1000000);
for k = 2:1000000
   x(k) = x(k-1) + 5;
end
toc
Elapsed time is 0.011938 seconds.
对要初始化的数组类型使用适当的预分配函数：

对于数值数组，使用 zeros

对于字符串数组，使用 strings

对于元胞数组，使用 cell

对于表数组，使用 table

练习

1、

A=[0 -1 4;9 -14 25;-34 49 64];
B=zeros(3,3);
for i=1:9
    B(i)=A(i);
end
disp(B)

2、

for i=1:9
    while B(i)<0
        B(i)=0;
    end
end
disp(B)

break 逻辑 与while逻辑相关

2.3Tips for Writing

1、clear all：清除所有变量  close all：关闭所有图像。

2、在命令后面加上分号，输出结果不会显示出来。

3、换行符号用...

A=[1 2 3 4 5 6;...

6 5 4 3 2 1];

4、

2.4 Functions

2.4.1 Functions&&Script

Function用来定义函数，一般一个函数放在一个.m文件里。

Script与Function最大的区别就是Function有keywords，可以通过keywords在command window中直接使用。

可以使用edit(which('mean.m'))来检索mean是不是keyword。

2.4.2 自定义Function

Function [y1,...,yN] = myfun(x1,...,xM) 声明名为 myfun的函数，该函数接受输入 x1,...,xM 并返回输出 y1,...,yN。此声明语句必须是函数的第一个可执行代码行。有效的函数名称以字母字符开头，并且可以包含字母、数字或下划线。

可以将函数保存在以下位置：

只包含函数定义的函数文件中。文件的名称须与文件中第一个函数的名称一致。

包含命令和函数定义的脚本文件中。函数必须位于该文件的末尾。脚本文件不能与文件中的函数具有相同的名称。

文件可以包含多个局部函数或嵌套函数。为提高可读性，可使用 end 关键字来表示文件中每个函数的末尾。

以下情况下需要 end 关键字：

文件中有任意函数包含嵌套函数。

该函数是函数文件中的局部函数，并且文件中有局部函数使用 end 关键字。

该函数是脚本文件内的局部函数。

例如：

在名为 stat.m 的文件中定义一个函数，返回输入向量的均值和标准差。

function [m,s] = stat(x)
    n = length(x);
    m = sum(x)/n;
    s = sqrt(sum((x-m).^2/n));
end

从命令行调用该函数。

values=[14.5,19.8,16.7,15,16.9];
>> [ave,stdev]=stat(values)

ave =

   16.5800

stdev =

    1.8605