MATLAB 中有哪些命令,让人相见恨晚?

https://www.zhihu.com/question/24499729

谢不邀

提问都说了是命令,大家回答那么多函数干什么... 我来给一个超级大杀器
在命令行敲入  dbstop if error

如果运行出现错误,matlab会自动停在出错的那行,并且保存所有相关变量。再也不用设断点了有没有!!!

bsxfun
强大的、万能的、不同维数的矩阵扩展混合运算,从此告别矩阵运算中的for循环
另,matlab里所有以fun为后缀的命令都很好用,arrayfun,cellfun,structfun,等等

@() 匿名函数
使用函数式编程,在编写以数学公式为主的程序中比传统编程方法好得多

set(gca,...)等
所有与画图相关的命令都包含了千奇百怪的属性和参数,非常实用,如著名的 set(gca,...)

$ $
画图中的标题、注释等文字内容是可以直接使用latex格式的,只须在包含latex命令的部分前后用$框起来即可

slice
高维切片,要在三维空间中表现三个以上变量的函数关系,必备此命令


Matlab里面应该多用向量运算,把循环语句转变为向量运算会省很多时间,程序也更简洁易读。
比方说,一个名叫array数组里面,你要将里面大于1的都变成0,就不必用到循环:
array(array>1)=0;
把大于1小于3的变成0。
array(and(array>1,array<3))=0;
当然,还可以使用find,这个也很好用。
—————————————————————————————
另一个重要的技巧是delete/clf-plot-pause
用plot可以画图(注意记录句柄),然后用delete删掉特定图象,或用clf清图,再绘制,这可以在figure窗口产生动画。但是如果只plot,往往只会在全部程序执行结束时显示,这时候需要用pause让figure完成图像的更新。drawnow貌似也可以,但是我比较喜欢用pause,能够简单地控制动画的速度。
这会方便调试和展示。这个技巧尤其适合使用matlab的图形用户界面设计功能时构造一个显示运行状态等信息的figure。
—————————————————————————————
mathworks 有一个fig函数(不是系统自带的,是别人编写的),可以很容易地调整字体、尺寸以及绘图是否有边框等等,不必画出来再自己手动调整。这对于写论文的人来说会很方便。
——————————————————————————————
善用eval,可以让你的程序的灵活度大大加强。尤其是在变量名的问题上。当然,这可能会对代码维护和调试带来麻烦。很多情况可以用其他方法代替。 @王备 指出了一个链接,值得一读 Alternatives to the eval Function
—————————————————————————————
save、load可以将部分或全部变量、结构体等存入mat文件或从mat导入workspace
global可以将变量变为全局变量,在各函数之间共享。不过这不太好用,尽量慎用吧。
——————————————————————————————
exist可以检测某目标(如变量)是否存在,减少一些麻烦。
——————————————————————————————
surf、mesh都很漂亮,不过surf之后记得用shading interp,看起来更漂亮。
——————————————————————————————
对于一些重复性的矩阵赋值,比如:1、2、3、4、1、2、3、4
可以使用repmat,将一个矩阵重复扩展为更大的矩阵。
——————————————————————————————
很多函数都有高级的用法,当使用到了,但又觉得有点麻烦的时候,不妨help一下,看看其他的用法。
总而言之,还是多上网搜索,一般问题总有很好的答案。


1. 转置但不共轭
b = a.'
2. 得到列向量
b = a(:)
相应的,得到行向量就是
b = a(:).'
3. 清除并更新有 class 的定义
clear classes
当你修改了某些 class 的定义时,需要调用此命令——仅仅 clear all 是不够的。
4. 忽略某些返回参数
[~, idx] = max(a)
忽略最大值,只返回最大值的索引。
【注】
这种语法是在 Matlab 7.9 (2009b) 中引入的,在此之前的版本不支持这种用法。如果你使用早期版本,可以使用如下写法:
[tmp, idx] = max(a); clear tmp
或者一句话搞定
[idx, idx] = max(a)
不过第二种写法不见得每个人都会喜欢。

5. 去除长度为 1 的维度(常用于 3D 矩阵转为 2D)
b = squeeze(a)
Simulink 的 save to workspace 模块经常得到一个 [m x 1 x n] 的三维矩阵,直接使用非常不便。使用 squeeze 后,会将其转换为 [m x n] 的二维矩阵。

6. 代码运行耗时显示
tic; some_code_to_run(); toc
tic/toc 会在 Matlab 命令行中打印出两条命令之间的代码的运行耗时。更详细的代码效率分析,可以使用 profile 命令。
7. 电脑跑分
bench
只能用于相同版本的 Matlab 下电脑之间的比较。
8. 记录命令行命令
diary on
% Lots of my Matlab commands here.
diary off
使用 diary 好处是运行结果也会被记录下来,并且以纯文本方式存储,方便编辑。更正式的报告,则推荐使用 publish。

9. 反转向量
b = fliplr(a)    % For row-vector a.
b = flipud(a)    % For column-vector a.
b = wrev(a)      % For any vector a.
b = a(end:-1:1); % This is the implementation of function wrev.

10. 去除数据中的 NaN
a(isnan(a)) = []
对于 inf 和 -inf,使用 isinf 命令。
答个冷门的,未必人人都相见恨晚,但对于合适的项目还是有点用。

其实也很简单,就是sound函数,能发出声音,可以在程序跑完的时候来两下。因为我的一个项目需要反复跑程序调参数,程序跑一次从10分钟到几小时不等都有可能。很多时候电脑开着跑程序,我在旁边干别的事,如果为了不浪费时间就要隔一会看看电脑跑完程序没,觉得很烦。后来想到能不能让电脑跑完主动告诉我?就搜到了这个函数。当然网上还有更高端的,程序运行到指定位置可以给自己的邮箱发个邮件,这样只要在有网的地方就能通过手机之类随时远程获知了。。。

然后觉得自己之前真是sb,为什么早没想到。。。

另外再补充一条:shutdown命令,可以自动关机,并且可以设置延时。很多时候为了赶时间,我都在睡前最后一次跑程序之前在程序末尾加上shutdown(当然之前还要加save),然后就可以呼呼睡大觉了,程序跑完会自动关机的。另外windows的命令框里也可以shutdown(可能就是从那里来的),不过估计很少有人用。

1. 关于向量化
大家都知道,如果用一个函数对一个矩阵中的每一个元素进行求值,那么就要用到 .*, ./之类
function r = fmat(x)
    r = x.^2 + 1./x;
end 

A = [1 2 3 4];
fmat(A)
上面这样就可以做到用fmat对A每一个元素分别求值。但不是所有的函数都像fmat这样,下面这种情况就不行:
function r = fmat2(x)
if x > 0
    r = x.^2;
else 
    r = 1./x;
end

fmat2(A)    % 会出错
可以看到,因为上面fmat2 里面有判断语句,就不能把整个矩阵作为输入参数。我知道的解决方法有3种:
% 法1
for i = 1:length(A)
    r(i) = fmat2(A(i));
end

% 法2 用arrayfun
r = arrayfun(@fmat2, A) 

% 法3 用逻辑矩阵

r = zeros(1,length(A));
l1 = (A>0);
l2 = ~l1;
r(l1) = A(l1).^2;
r(l2) = 1./A(l2);
值得说一说的是第二种和第三种。 arrayfun只是比用for更加简洁,速度上没有太大差别。但是,用 逻辑矩阵的话,就比前两种快10倍。逻辑矩阵是向量化的利器。比如 A>0 就会返回一个逻辑矩阵,里面全是1, 0 。并且, 逻辑矩阵可以当作索引, A(A>0) 就把所有A>0 的数返回。

2. 类型转换
matlab中一共有4种类型:func(函数句柄), num(数值), sym(符号), str(字符)
- 在数值运算中,使用func, num
- 需要输入输出出,要用到str
- 进行符号运算(符号微分,积分)时,用sym

我们使用matlab,一般有这么个过程:
- 先进行公式推导,这时要用sym, symfun类型
- 再进行数值计算,这里要用普通的func, 和num 类型
- 最后与GUI界面交互的时候,要用str类型
这里介绍几个用于转换的函数
str2func('@(x,y)sin(x*y)') %  str -> @func 返回一个函数句柄

syms x y
fs(x,y) = x^2+sin(x*y);
fh=matlabFunction(fs); % symfun -> @func 返回一个函数句柄 !!!强烈推荐

fh = @(x)x.^2+sin(x);
fun2str(fh)    % @func ->str 函数句柄变为字符

3. 三维画图
大家都晓得,作图第一步,对 x, y 划分网格,用的是 meshgrid 。这样画出来的图在一个立方体的范围内,但是有时候我们要显示一个 圆柱体内的三维图形这时候画网格可以用 cylinder
[x, y] = cylinder(linspace(0,10,100),200)   % linspace(0,10,100)指在半径[0, 10]上划分为100份,参数200指的是在圆周方向上 200等分
z = sin(x) + cos(y);
mesh(x, y, z);

4. 学无止境
学matlab,还是要去论坛,看看大神们的说法,大神们的做法。比如说1中的向量化
function r = fmat2(x)
if x > 0
    r = x.^2;
else 
    r = 1./x;
end

fmat2(A)    % 会出错
论坛中的大神的做法是这样的
fmat3 = @(x)x.^2.*(x>0) + 1./x.*(x<=0);
fmat3(A)
极其简洁!!!我看到后非常佩服


先这些吧





  • 15
    点赞
  • 59
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值