[求助]频谱分析

[ct8613]

本人菜鸟一个，刚开始学习MATLAB。编了一个有关频谱分析的程序，但是实现不了，还望大家帮忙。小弟在此谢过了。在线等。
Ts=0.0005;
t=0:Ts:0.1;
%采样信号频率是fs=2000；
fs=1/Ts;
%原信号频率是50；
x=60*sin(2*pi*50*t);
%N是总共的采样点数；
N=length(x);
Ix=fft(x);
w=2*pi*50,c=5.72e-5;
Ux=Ix/(i*w*c);
u=ifft(Ux);
Fx=fft(u^2);
%N个点有N-1个间隔；
df=fs/(N);
n=0:N-1;
f=n*df;
subplot(211);
plot(t,x);grid;
subplot(212);
plot(f,abs(Fx));grid;

[fanxing39]

报什么错误 ？ 或者你想得到的和实现出来的哪里有差距

[fanxing39]

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

[ct8613]

谢谢，正是这个错误。
谢谢！！

[fanxing39]

谢谢，正是这个错误。
谢谢！！

不客气，没事常来这个地方逛逛，只有多交流，自己提高的才更快:redface: