※ 引述《k0939620064 (skan)》之铭言： : 标题: [问题]是否有更快的做法? : 时间: Wed Apr 30 17:59:24 2014 : : 各位大大 : 小弟目前遇上一个问题就是用for回圈做讯号处理时,往往都要跑个1.20分钟以上 : 想请问各位大大是否有更快的方法? : : 以下是小弟的程式码 : : for m=1:4 : for n=1:500000 在回圈中动态的扩大矩阵会有数量未知的记忆体配置动作， 所以回圈中的变数需要再回圈外先配置好确定大小的记忆体。 p.s. matlab editor 右边有一堆橘色的提示，看一下好吗.... : signal_square(n,m)=sum(sig((n-1)+1:n,m).^2); ~~~~~~~~~ 啊不就是n，你多做了一次加法与建立矩阵的动作，而且一个数字不需sum : % 将收到的讯号平方 : phase_noise(n,m) = angle(signal_square(n,m))./2; : % 将平方的讯号取角度後除2 = phase noise : phase_sig((n-1)+1:n,m) = (((angle(sig((n-1)+1:n,m)))- phase_noise(n,m))); : % 原始讯号= arg(收下来的讯号)-phase noise : phase_comp((n-1)+1:n,m) = exp(j*phase_sig((n-1)+1:n,m)); : % 原始讯号 : end : end : 整个程式可写成： signal_square = sig.^2; phase_noise = angle(signal_square)/2; phase_sig = angle(sig) - phase_noise; phase_comp = exp(j*phase_sig); 其实完全不用回圈的， 用matlab就要用matlab的精神来写阿。 p.s. 你若是想把原来的 sig 用大小和相位来表示的话， signal_square 不是应该为 abs(sig).^2 吗？ : -- :

※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 140.117.163.74 : ※ 文章网址: http://webptt.com/cn.aspx?n=bbs/MATLAB/M.1398851967.A.95F.html : → celestialgod:写MEX or 用GPU 04/30 18:58 : 推 sunev:(n-1)+1:n ? 04/30 19:29 --

※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 111.67.62.134 ※ 文章网址: http://webptt.com/cn.aspx?n=bbs/MATLAB/M.1398859031.A.B76.html