如何制作不带的二维矩阵？

Question

我现在正在为二维图编码矩阵。 虽然是这么简单的等式，但执行起来却需要很多时间。我认为它可以变得更快。 特别是，我认为“for - command term”可以简化。 我该如何简化？

q=1:1:30
x(q)=330+q*0.3
F=1:30:8970
T=x(1)-0.3:0.001:x(30)+0.3
n=size(T,2)
k=1:1:n
for a=1:1:30
I(a,k)=F(a)*exp(-2.*(T(:,k)))
end
happy=sum(I)
plot(k,I)

---

Answer 1

我会说时间是用来打印结果的。尝试在每行末尾使用;，它会加快计算速度。

您还可以通过以下逐个元素的计算来替换 for 循环：

a = (1:1:30).';
aux = repmat(exp(-2.*(T(:,k))), length(a), 1);
a = repmat(a, 1, length(k));
I = a.'.*aux.';