0

我有一个包含两个向量的变量。我想将这些数量从 i = 1 加到 n，其中 n = 100。当程序运行时，我可以单击工作空间中的变量来查看每个列位置的向量的元素值。执行其位置的变量时，我得到相应的值，（例如，z(1) = 73.2360, z(2) = 63.3701, .... 等等）。但我似乎无法弄清楚如何为这些向量编制索引，以使总和贯穿两个向量的元素。例如在 for 循环中。我也看了看symsum()，但那只需要泛型变量。我试图从 i = 1 到 100 求和 epsilon，其中 z_i 和 w_i 是感兴趣的向量。我希望 MWE 能展示更多我想要完成的事情。

MWE：

 a = 0.6;
 b = sqrt(16-a^2);
 c = 6.5;
 d = sqrt(225-c^2);
 xh = randn(1,100);
 yw = randn(1,100);
  z = a*xh + b*yw + 68;
  w = c*xh + d*yw + 160;

     n = 100;                              % # of data entries.
  zbar = ((1/n)*sum(z));                   % Height data.
  wbar = ((1/n)*sum(w));                   % Weight data.
 Zbarv = zbar*ones(1,100);                 % Height data vector.
 Wbarv = wbar*ones(1,100);                 % Weight data vector.

sz = sqrt((1/n)*((z-Zbarv)*(transpose(z)-transpose(Zbarv))));              
sw = sqrt((1/n)*((w-Wbarv)*(transpose(w)-transpose(Wbarv)))); 

czw = ((1/n)*(w-Wbarv)*(transpose(z)-transpose(Zbarv)));
rzw = czw/(sz*sw);

 b = wbar - ((czw/sz^2)*zbar);             % y-intercept
 m = czw/sz^2;                             % Slope
 epsilon = (1/(n-1))*sum((w-(m*z+b)).^2);  % Error

epsilon 应该对此求和，但应该是来自 的总和i=[1,100]。我不知道如何索引w_i和z_i向量的值以使总和通过它们。我尝试为他们制作一个输入变量z(i)，w(i)但结果很混乱。

有什么建议么？

4

1 回答 1

1

听起来你只是想要：

epsilon = (1/(n-1))*sum((w(1:100)-(m*z(1:100)+b)).^2);

虽然你w和z只有 100 个元素长，但我认为这与你所做的相同。所以我认为你的代码实际上是正确的。

例如 if a = [1 2 3], then sum(a)is 6，您实际上不必索引a即可获得总和。因此，您的代码正在查找所有 100 个元素的总和。

您应该阅读Matlab 冒号运算符

基本思想是在matlab中类似的东西1:10只是简写[1 2 3 4 5 6 7 8 9 10]，你可以用它来索引一个向量。因此，如果我们有一个向量v = rand(100,1)，那么要获取v它的前 10 个元素v(1:10)

于 2013-04-30T09:21:02.977 回答