fs=100; T=1/fs; N=2000; dur=N*T; t=0:T:dur-T; X=fft(yn1); Y=abs(X); A=Y/(N/2); width=fs/N; f=0:width:fs-width; subplot(412); plot(f(1:N/2),A(1:N/2)); title('原始信号频谱图');每段代码意思

这段代码是用来进行信号处理的，具体解释如下：

1. `fs=100;` 设置采样率为100Hz。
2. `T=1/fs;` 计算采样时间间隔。
3. `N=2000;` 设置采样点数目为2000个。
4. `dur=N*T;` 计算采样时间长度。
5. `t=0:T:dur-T;` 生成采样时间序列。
6. `X=fft(yn1);` 对信号进行傅里叶变换，得到频域信号。
7. `Y=abs(X);` 取频域信号的幅值。
8. `A=Y/(N/2);` 对幅值进行归一化处理。
9. `width=fs/N;` 计算每个频率间隔的宽度。
10. `f=0:width:fs-width;` 生成频率序列。
11. `subplot(412);` 在一个图中生成四个子图，这里选择第二个子图。
12. `plot(f(1:N/2),A(1:N/2));` 绘制频谱图。
13. `title('原始信号频谱图');` 设置子图标题为“原始信号频谱图”。