[fanxing39]

Fx=fft(u^2); % 这个地方错了
改为：Fx=fft(u.^2); % ^ 是矩阵的幂， .^是向量的幂，你这里是向量，所以应该用后者
这是个新手经常碰到的问题，比如 * / 也有这样的问题