上一篇画出图:
明明函数只有三个频率为啥六条线,其实是横坐标设置问题及频谱搬移问题,
因为plot(k)没有x,只有一个变量,那x应该怎么取呢,freq=(0:N/2)*fs/N;这样,取采样点的一半作图,画的时候这样plot(freq,abs(k(1:N/2+1)))
再看看图,代码改为下面这样
close all;
clear all;
clc;
%采样频率
fs=1024;
%采样间隔
Ts=1/fs;
%信号
N=2048;
t=(0:N-1)/fs;
S=sin(2*pi*50*t)+sin(2*pi*100*t)+sin(2*pi*200*t);
k=fft(S,N);
freq=(0:N/2)*fs/N;
figure(1)
plot(freq,abs(k(1:N/2+1)))
是不是明显多了。。。。。。。。。。。。。。。。。。。。。。。。。。