fangbo.m

x=2*sin(4*pi*t)+5*cos(8*pi*t) 对其进行快速傅立叶变换 加入高斯噪声 再进行频谱分析

N=200;n=0:N-1;t=0.0001*n;q=n*0.02/N;
%N=64;n=0:N-1;t=0.01*n;q=n*2*pi/N;
y=square(2*pi*50*t);
figure(1);
plot(t,y);
%figure(2);
%plot(y);
%figure(2);
%bar(y);
%x=fft(y);


x=fft(y,200);
c=abs(2/200*x);
figure(2);
bar(c);
figure(3);
plot(q,c)