-2

我有时间(第 1 列)和电流幅度(第 2 列)的 csv 数据。我想绘制电流的 FFT。这实际上是超过 30ns 的仿真数据,数据步长为 1ps。我可以在 MATLAB 中绘制电流与时间的关系。但是在执行 FFT 功能时,它根本没有像它所说的那样绘图

Error using plot
Vectors must be the same length.

Error in FFT_Ideal_current_maxstep_1ps (line 25).
plot(f,Y)

谁能帮我?我还附上了 MATLAB 代码和CSV文件。

我还想绘制功率谱密度。如果有人可以提供帮助,那就太好了。我想获得超过 2GHz 或更多频谱范围的 FFT 和 psd

MATLAB代码1:

% open data file
fid = fopen('current_maxstep_1ps.csv');

% Read data in from csv file
readData = textscan(fid,'%f %f','Headerlines',1,'Delimiter',',');

% Extract data from readData
t = readData{1,1}(:,1);
x = readData{1,2}(:,1);

N = length(x);
ts = 0.000000000001;
Fs = 1/ts;
tmax = (N-1)*ts;
tm = 0:ts:tmax;
f = 0:Fs/(N-1):Fs/2;
y = fftshift(fft(x));
Y = abs(y);
plot(f,Y)

我还尝试了另一个 MATLAB 代码,其中的绘图(这里是图片:代码 2 的 FFT 图片)但在时域中显示,我想要频谱,如沿频谱的振幅尖峰。

MATLAB 代码 2:

% open data file
fid = fopen('Ideal_current_maxstep_1ps.csv');

% Read data in from csv file
readData = textscan(fid,'%f %f','Headerlines',1,'Delimiter',',');

% Extract data from readData
xData = readData{1,1}(:,1);
yData = readData{1,2}(:,1);

Ts = 1e-12;
Fs = 1/Ts;
%Fs = 1000;
%Ts = 1/Fs;

X = fft(yData);
plot(xData, abs(X))
4

1 回答 1

1

问题是和的长度fY一样。您可以使用length(f)和进行检查length(Y)。原因是它fft也计算了负频率。因此,您应该定义f如下:

f = -Fs/2:Fs/(N-1):Fs/2;

请注意,fft 是共轭对称的,因为输入数据是真实的。

您可以使用以下命令限制绘制的频率范围xlim

xlim([0 3*10^9]) % limit x range between 0Hz and 3GHz

在此处输入图像描述

于 2017-06-06T11:28:19.523 回答