频率和残留扭矩的IDEXing，可变W和S从1到1000开始，并重新注射2至11。
因此，我认为索引不正确。
尝试修改索引如下：
为了i = 1：1000
w（i）= 0.2*（i-1）;
a（1，i）= 1;％振幅（假设）
t（1，i）= j（1）*a（1，i）*（w（i）^2）;％扭矩
s（1，i）= t（1，i）;第一成员的残余扭矩％
为了n = 2：11％成员2至11
a（n，i）= a（n-1，i） - （j（n-1）*e（n-1）*a（n-1，i）*w（i）^2）;
t（n，i）= j（n）*a（n，i）*w（i）^2;
s（n，i）= s（n-1，i） + t（n，i）;
a（n，i）= a（n-1，i）-s（n-1，i）*e（n-1）;％a（n）=（先前的res torques的总和）*灵活性
结尾
结尾
情节（w，s（11，:)）;
Xlabel（“频率（rad/s）”）；
ylabel（“残留扭矩”）；
