How to speed the code? It does not want to end processing!!!!
1 回表示 (過去 30 日間)
古いコメントを表示
R = length(lambda) ;
T = size(X,2) ;
time_dmd = zeros(T-1,R);
for iter = 1:T-1
for p = 1:R
time_dmd(iter,p) = b(p)*(exp(omega(p)*t(iter)));
Xdm(:,iter,p) = real(Phi(:,p)*(b(p).*exp(omega(p)*t(iter))));
end
end
1 件のコメント
回答 (1 件)
DGM
2023 年 6 月 29 日
編集済み: DGM
2023 年 6 月 29 日
Here's my guess.
% assuming b,omega,t are all row vectors
% width of Phi is equal to length of b,omega
time_dmd = b.*exp(omega.*reshape(t(1:end-1),[],1))
Xdm = permute(Phi,[3 2 1]).*time_dmd;
Xdm = real(permute(Xdm,[3 1 2]));
2 件のコメント
Steven Lord
2023 年 6 月 29 日
What is the size and class of the variable x that you load from your MAT-file? You can display this information using the whos function.
x = magic(4);
whos x
参考
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!