有哪些MATLAB代碼驚艷到了你？

02-01

希望大家分享一下自己看到過或者寫過的精彩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. 電腦都給你打開了就只看我報告!?

※MATLAB學習-2-啟動
※MATLAB教學視頻：詳解Excel文件的讀取和寫入操作
※LTE 調製與編碼（一）
※hdlcoder 靠譜嗎?
※基於Matlab HDL Coder的8x8實矩陣QR分解

TAG:MATLAB |