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