shiftdim

我有點疑問，妳這樣循環，distance每被循環壹次，它的值就被更新壹次，最後記錄的，不就只是最後壹個lig，col位置對應的值麽？

[col, lig] = find(fsm == nm); % 確定需要進行循環的坐標

for id = 1:length(col)

for area = 1:narea

dist(area) = ...

trace(cohinv(:,:,area)*shiftdim(T3mean(col(id),lig(id),:,:)));

end

end

% 這樣寫，原則上是沒有加快多少的。。。。

% 建議用matlab的c/c++接口，利用mex文件來加速。。。