对实验数据进行频谱分析,使用matlab中的fft函数。
Y=fft(X,N),其中X为需要分析的信号,N为数据点数,一般不指定N,简化为Y=fft(b)。
*****************************************************************************************
*****************************************************************************************
针对实验项目具体情况,代码如下:
vds_test; %长度为1251,周期为400kHz,一个周期,滤波之后。
fs=5e8; %采样间隔2e-9,采样频率500 MHz
N=1251-1; %使N为偶数
Y=fft(vds_test);
magY=abs(Y(1:1:N/2))*2/N;
logmagY=20*log10(magY)+120; %将幅频变为对数幅频,加120是变为dBμV
f=(0:N/2-1)*fs/N; %对数幅频的横坐标
semilogx(f,logmagY)
*****************************************************************************************
*****************************************************************************************