快忘光了MATLAB,把不会的函数存一下
★matlab中的下标都是从1开始的
★原文:https://blog.csdn.net/caomin1hao/article/details/78510298
Matlab函数定义格式,用关键字function定义函数
function [输出变量] = 函数名称(输入变量)
% 注释
% 函数体
简单的函数如下:
function [a,b,c]=test(d,e,f)
a=d+e;
b=e+f;
c=f+d;
end
函数的调用:
[q,w,e]=test(1,2,3)
★单引号是转置
vector=vector(:)';%输入向量转换为行向量
find()函数
复制的 https://www.cnblogs.com/anzhiwu815/p/5907033.html
[1]find()函数基本功能
find()函数的基本功能是返回向量或者矩阵中不为0的元素的位置索引。
>> X = [1 0 4 -3 0 0 0 8 6];
>> ind = find(X)
ind =
1 3 4 8 9
其有多种用法,比如返回前k的不为0的元素:
>> ind = find(X,2)
ind =
1 3
也可以写成:
>> ind = find(X,2,'first')
ind =
1 3
返回后k个不为0的元素:
>> ind = find(X,2,'last')
ind =
8 9
若X是一个矩阵,索引该如何返回呢?
>> X = [1 -3 0;0 0 8;4 0 6]
X =
1 -3 0
0 0 8
4 0 6
>> ind = find(X)
ind =
1
3
4
8
9
这是因为在Matlab在存储矩阵时,是一列一列存储的,我们可以做一下验证:
>> X(4)
ans =
-3
假如你需要按照行列的信息输出该怎么办呢?
用法:[a1,a2]=find(a),找出a矩阵中非零元素所在行和列,并存在a1,a2中
>> [r,c] = find(X)
r =
1
3
1
2
3
c =
1
1
2
3
3
如果你还需要输出具体的元素值:
>> [r,c,v] = find(X)
r =
1
3
1
2
3
c =
1
1
2
3
3
v =
1
4
-3
8
6
[2] 进阶技巧
(1)find()函数的功能是找到向量或者矩阵中不为0的元素,那如果需要找到其中满足一定条件的元素,比如,等于4的元素该怎么办呢?
>> X = [1 0 4 -3 0 0 0 8 6];
>> ind = find(X == 4)
ind =
3
在Matlab中,有一个logical数据类型,和C++中的bool值相同。find()函数将logical值0也视为数值0,logical值1视为数值1,因此可以用上面的方法找到矩阵中满足一定条件的元素。
(2)在程序中,我们经常会以向量中是否包含某一元素为判断条件,比如X中是否有等于9的元素,这时也可以应用find()函数:
>> if isempty(find(X == 9))
log = 0
else
log = 1
end
log =
0
这里需要补充说明下,当输入中没有非0元素时,findf()返回一个空数组。
[3] 高级技巧
在进阶技巧(2)中,假如向量X的阶数很高,比如1000万维,如果我们用上面的方法,运行速度会很慢。仔细分析这个语句,我们其实并不希望找到X中等于9的全部元素,我们只想找到其中是否有等于9的元素,所以,我们只需作如下更改:
isempty(find(X == 9,1))
我们做一个小实验:
X = randi(100,1,10000000);
tic
for i = 1:100
isempty(find(X == 9));
end
toc
插播:【简单地说,tic和toc是用来记录matlab命令执行的时间
tic用来保存当前时间,而后使用toc来记录程序完成时间。
两者往往结合使用,用法如下:
tic
operations
toc
显示时间单位: 秒
Tic和toc函数可以计算运行一段时间的代码】
输出是:
Elapsed time is 1.954537 seconds.
更改后:
X = randi(100,1,10000000);
tic
for i = 1:100
isempty(find(X == 9,1));
end
toc
输出为:
Elapsed time is 0.757994 seconds.
运行效率约是原来的2.5倍。