sunshine

% 方案1：“x = a*cos(2*pi*w*t)”的形式：
% -------------------------------------------------------------------------
% 注意：1.时域的持续时间范围应较大；
%    2.频率w与序列k的对应关系：w =   k * df;
%    3.采样频率 fs 应大于 w 的2倍
%    4.结果曲线的峰值的横坐标对应的就是 w 和 -w 值
fs = 10;       %采样频率
N = 1024;           %采样点数
t = (0:N-1)/fs;     %采样时间序列s

a = 0.75;
w = 4;
x = a*cos(2*pi*w*t);

subplot(2,1,1);
plot(t, x);
xlabel('t/s');

xf = fft(x,N)/N;
xf = fftshift(xf); %双边复数谱

df = fs/N;       %频率分辨率Hz

f = (-N/2+1:N/2)*df; %频域序列

subplot(2,1,2);
plot(f, abs(xf));
xlabel('f/Hz');



% 方案2：“x = a*cos(w*t)”的形式：
% -------------------------------------------------------------------------
% 注意：1.时域的持续时间范围应较大；
%    2.频率w与序列k的对应关系：w = 2 * pi* k * df;
%    3.采样频率 fs 应大于 w/(2*pi) 的2倍
%    4.结果曲线的峰值的横坐标对应的就是 w 和 -w 值
fs = 10;       %采样频率
N = 1024;           %采样点数
t = (0:N-1)/fs;     %采样时间序列s

a = 0.75;
w = 20;
x = a*cos(w*t);

subplot(2,1,1);
plot(t, x);
xlabel('t/s');

xf = fft(x,N)/N;
xf = fftshift(xf); %双边复数谱

df = fs/N;       %频率分辨率Hz

f = 2*pi*(-N/2:N/2-1)*df; %频域序列

subplot(2,1,2);
plot(f, abs(xf));
xlabel('f/Hz');