如何在 MATLAB 中使用带有系列的 for 循环

Question

我有一个在图像中的方程

我正在尝试绘制 y 轴是这个方程，x 轴是时间的图 - 只是一个向量。 我所有的缩写：

%Initials
Beta=[24 123 117 262 108 45]*10^-5; %pcm
Lambda=[0.0127 0.0317 0.1160 0.3106 1.4006 3.8760]; %1/s
LAMBDA=10^-4 ;   %s
W=[ 0.376 -0.0133 -0.0426 -0.153  -0.972 -3.38 -29.5]
Rho=400*10^-5
t=linspace(1,30,7)

这是我正在使用的代码：

for n=1:7
    for j=1:6
    S1=Rho*sum(exp(W(n)'.*t)/(W(n)'.*(LAMBDA+(sum(Beta(j).*Lambda(j)./(W(n)+Lambda(j))^2)))))
    end
end
semilogy(t,S1,'b','linewidth',2); 

而 S1 返回的答案太多了，据我所知，它应该只给出 7 个... 而且我是matlab和一般编码的新手，所以如果答案很明显，我仍然不知道如何使它工作：D

Answer 1

让我们先澄清一些事情。

为了绘制（任何类型的）二维绘图，您需要在 Matlab 中使用两个向量。它们的长度应该相等。一个用于所有 x 坐标。另一个用于所有 y 坐标。

您在t=linspace(1,30,7) 中获得 x 坐标。但是，您没有相应的 y 坐标。

在您的情况下，最好将您的公式表述为t 的函数。为了清楚起见，让我们分解总和。例如，

function num = oscillation_modes_of_sort(t)
    outer_sum = 0;
    for j=1:numel(W)        
        inner_sum = 0;
        for i=1:numel(Beta)
            inner_sum = inner_sum + Beta(i)*Lambda(i)/(W(j)+Lambda(i))^2;
        end
        outer_sum = out_sum + exp(W(j)*t)/(W(j)*(LAMBDA+inner_sum));
    end
    num = Rho * outer_sum;
end

现在，您的 y 坐标将为 oscillation_modes_of_sort(t)。

如果W 和Beta 更长，有一些方法可以使代码更简短或更友好。但是，让我们在以后做这些。