因为一只蝙蝠，我又熬夜写下自己的第二篇博客
(っ•̀ω•́)っ✎⁾⁾ 开始：

问题描述：

希望在MATLAB里面构造这样一个实际信号并画出频谱图
$x=0.5(sin2\pi f_{1}t)+2(sin2\pi f_{2}t)$
$f_1=10Hz,f_2=40Hz$

构造信号

fs=100;%采样频率
N=100;%数据长度，比如N=100,就是取1秒的信号
n=0:N-1;
t=n/fs;%时间序列
f1=10;
f2=40;
x=0.5*sin(2*pi*f1*t)+2*sin(2*pi*f2*t);

此处我最开始的疑惑：
时间t应该如何取呢？为什么是t=n/fs;
fs是采样频率，也就是一秒钟采样多少个点，换言之：每1/fs秒，采样一次；fs=100Hz,采样时间间隔就是0.01s，100个点的时间长度就是1秒；因此我的时间序列就是t=n*1/fs

打印频谱图

y=abs(fft(x));
y=y(1:end);
figure;
plot(y);

此处我最开始的疑惑：
为什么频谱的位置都+1了？原本应该是在10Hz，40Hz处取最大值，现在是在11Hz，41Hz处取的最大值；
原来是fft后的结果第一位是直流分量，所以画图应该plot（2：end）

y=abs(fft(x));
y=y(2:end);
figure;
plot(y);

修改后的图就正确了：(๑◡๑)

如果数据点数不是100，比如取128（方便进行fft运算）：
此时频谱图为：
此时我想直接在横轴显示实际频率，那么如何更改呢？
代码如下：

set(gca, 'XTick', [0:N/fs*10:N]) ;%X坐标轴刻度数据点位置；
set(gca,'XTickLabel',{
   '0','10','20','30','40','50','60','70','80','90','100'});
xlabel('频率轴/Hz');
ylabel('幅度轴/V');
axis([0 128 0 120])