u=1+3* rand(1,2001);
s=cos(2*pi*2*t);
y=u.*s;
y=y(1:100);
figure(1);
plot(y);
t=0:0.01:20 實際上是有 2001 個點的, 而 u 是壹個長度為 2000 的向量, u.*s 當然是有問題的.
寫成 u=1+3* rand(1,2001); 即可.