2

我正在尝试使用正在运行的函数中的值填充数组。在此函数中，我以from到z2的增量更改变量。同时我想用这些值填充一个数组。我尝试了一个while循环和一个for循环。如何只用一个 for 循环来写这个？0.0100.99

我的代码：

Ys = y(2000);
Mp = ((max(y)-Ys)/Ys)*100;
[max, i] = max(y);
tp = t(i);
z = sqrt(((log(Mp))^2)/((pi*pi)+(log(Mp))^2));


%Given Equations
wd = pi/tp;
wn = wd/(sqrt(1-(z*z)));
ts = 4.6/(z*wn);
tr = 1.8/wn;

%for loop for z from 0 to 0.99
for i = 1:100;
   for z2 = 0:0.01:0.99
      %tr*wn = Fa
      Fa(i) = 2.917*z2^2 - 0.4167*z2 + 1;
      i = i + 1;
      if i >= 100;
          break;
      end
    end
end

find_zeta = interp1(Fa, 1.8);

disp(Fa);
disp(find_zeta);

我遇到的错误：我1的数组中只有 s。

4

2 回答 2

4

您可以在“矢量化”matlab 语句中执行此操作而无需循环

z2 = 0:0.01:0.99;
Fa = 2.917*z2.^2 - 0.4167.*z2 + 1;

首先，计算z2您需要的值。接下来，逐个元素（注释 .^语法）对向量进行操作，计算 Fa。^和运算符之间存在根本区别.^，您真的很想知道区别。通常，您可以将 a 添加.到运算符，这意味着它将对数组的各个标量元素起作用。在此处阅读有关矩阵和数组运算符的信息。例如：

A = rand(10);
B = rand(10);
% matrix multiplication
C=A*B; 
% element by element multiplication
C=A.*B

您收到的投诉仅意味着您的工作效率低下。由于 matlab 事先不知道它的大小Fa，因此每次附加到它时都需要重新分配数组。你仍然可以使用它，它只是效率低下。

于 2012-10-08T13:31:32.703 回答
1

matlab 中的索引从 1 开始，而不是像在其他编程语言中那样从 0 开始

您得到的错误是关于预分配数组。添加Fa = zeros(1,100);到脚本的开头，它将消失。

Fa = zeros(1,100);
i = 1;
for z2 = 0:0.01:0.99
   Fa(i) = 2.917*z2^2 - 0.4167*z2 + 1;
   i = i +1;
end

正如另一个答案中所建议的那样，使用“矢量化”是更可行的方法。

于 2012-10-08T13:37:38.827 回答