这是一个办法。
clc,清晰
%的基本频率
f1 = 1/3;
f2 = 5/6;
%采样频率
fs = 50 *马克斯((f1、f2));
%采样周期
Ts = 1 / f;
%时间跨度(5期)
T = 5 *马克斯([1 / f1, 1 / f2]);
%时间轴
t = 0: Ts: T-Ts;
%频率轴
f = - f / 2: fs /长度(t): fs / 2 - f /长度(t);
%的信号
x = 2 + cos(2 *π* f1 * t) + 4 *罪(2 *π* f2 * t);
% fft
X = fftshift (fft (X));
%图
无花果=图(“单位”,“归一化”,“位置”(0.15 - 0.35 0.75 - 0.55));
%绘制信号
次要情节(1、3、1)、情节(t, x)
包含(的时间(秒));ylabel (“振幅”);
标题(“信号”);网格“小”;
%的阴谋fft级
次要情节(1、3、2),图(f, abs (X)),
包含(的频率(赫兹)),ylabel (“级”),
标题(谱,大小);网格“小”;
%的阴谋fft相
次要情节(1、3、3),情节(f,打开(角(X))),
包含(的频率(赫兹));ylabel (“阶段”);
网格“小”;
标题(的频谱,打开相);
