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

作者:冷哲
链接:https://www.zhihu.com/question/24499729/answer/28087804
来源:知乎
著作权归作者所有,转载请联系作者获得授权。

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,将一个矩阵重复扩展为更大的矩阵。
作者:芷菡
链接:https://www.zhihu.com/question/24499729/answer/121898170
来源:知乎
著作权归作者所有,转载请联系作者获得授权。

匿名用户

作者:tengmoon
链接:https://www.zhihu.com/question/24499729/answer/28105239
来源:知乎
著作权归作者所有,转载请联系作者获得授权。

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)
极其简洁!!!我看到后非常佩服

  • 4
    点赞
  • 17
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值