有哪些MATLAB代码惊艳到了你？

希望大家分享一下自己看到过或者写过的精彩MATLAB代码。

向量化编程的利器：reshape、bsxfun、arrayfun、repmat，各种秒杀for循环。不过自从推出Compatible Array Sizes之后，bsxfun的作用有所下降（但我并不喜欢Compatible Array Sizes，导致有时候矩阵大小出错了还不会报错）。

下面是一个我之前遇到过的例子。

A = randn(37,1e5); B = randn(41,1e5)+1i*randn(41,1e5); C = randn( 1,1e5); %% 使用for循环 tic peak = 0; for i = 1:size(A,1) A_fft = conj(fft(A(i,:))); for j = 1:size(B,1) temp = B(j,:).*C; temp_fft = fft(temp); result = abs(ifft(A_fft.*temp_fft)); peak = max(peak,max(result)); end end toc

时间已过 6.287102 秒。

%% 使用bsxfun tic A1 = reshape(A,1e5,37, 1); B1 = reshape(B.,1e5, 1,41); C1 = reshape(C,1e5, 1, 1); result1 = abs(ifft(bsxfun(@times,fft(bsxfun(@times,B1,C1)),conj(fft(A1))))); peak1 = max(result1(:)); toc

时间已过 2.359991 秒。

%% 验证前后两种方法结果是否相同 if abs(peak-peak1)&<1e-10 disp(true); end

true。

上述结果，CPU为8700，matlab版本为R2018b，这个版本对for循环已经有了一定的优化，但通过向量化编程、列优先的手段，仍然可以有效地提高程序执行效率。

[x,y,z]=meshgrid(linspace(-1.3,1.3)); val=(x.^2 + (9/4)*y.^2 + z.^2 - 1).^3 - x.^2.*z.^3 - (9/80)*y.^2.*z.^3; isosurface(x,y,z,val,0) axis equal view(-10,24) colormap flag camlight lighting phong

看到这个问题马上打开程序粘上去了, 做课设的时候本来打算, 把这个代码隐藏起来等老师运行检查的时候调戏一下老师? 狗头jpg. 电脑都给你打开了就只看我报告!?