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!