[求助]请问这个程序错在什么地方？

恩，已经改了，可以运行了，改成(sin(800*t))^2就可以了
随便在问一下下面这个程序是我自己写的，但是老师说不对，红色字是老师写的，我想问下那个相关迭代程序改怎样写，我把我自己画的流程图也附上，另外我还想问下，我要在y(t)中加入一个信噪比为65dB的高斯白噪声，该怎样编程序，请帮忙看下，谢谢
u=input('输入数值u：')
p(1)=0
r=0
a(2)=0
for t=2:512
p(t)=0.99*p(t-1)+(sin(800*t))^2;
e(t)=5*sin(800*(t+pi/3))-r;
y(t)=10*sin(800*t);
a(t+1)=a(t)+u*e(t)*y(t)/p(t);
end
for k=1:511
H(k)=a(k)*10*sin(800*(512-k));
end
r=sum(H(k));
e(t)=5*sin(800*(t+pi/3))-r;//这里没有进行相关的迭代
a(t+1)=a(t)+u*e(t)*y(t)/p(t);//这里没有进行相关的迭代
for k=1:511
H(k)=a(k)*10*sin(800*(512-k));;//这里进行相关的迭代与前面一样
end
r=sum(H(k));
e(t)=5*sin(800*(t+pi/3))-r;
plot(e)

[此贴子已经被作者于2007-6-2 15:02:13编辑过]