OFDM功率谱密度绘制matlab，采样周期Ts=0.01，频率范围是0~400HZ

OFDM信号的功率谱密度可以通过对其进行离散傅里叶变换（DFT）得到。在MATLAB中，可以使用fft函数来进行DFT计算，并使用pwelch函数来计算功率谱密度。

假设OFDM信号的时域采样点数为N，可以先生成一个长度为N的随机OFDM信号：

N = 1024; % 时域采样点数
Ts = 0.01; % 采样周期
t = 0:Ts:(N-1)*Ts; % 时域采样点
f = 0:1/Ts/N:1/Ts-1/Ts/N; % 频域采样点
x = randn(1,N); % 随机OFDM信号

然后对x进行DFT计算，得到频域信号X：

X = fft(x);

接着使用pwelch函数计算功率谱密度：

[Pxx, F] = pwelch(x, [], [], [], 1/Ts);

其中，第一个参数是要计算功率谱密度的信号，第二个参数是窗口长度（默认为256），第三个参数是重叠长度（默认为50%），第四个参数是FFT点数（默认为窗口长度），最后一个参数是采样频率（这里为1/Ts）。

最后，可以使用plot函数绘制功率谱密度曲线：

plot(F, Pxx);
xlabel('Frequency (Hz)');
ylabel('Power Spectral Density');

完整代码如下：

N = 1024; % 时域采样点数
Ts = 0.01; % 采样周期
t = 0:Ts:(N-1)*Ts; % 时域采样点
f = 0:1/Ts/N:1/Ts-1/Ts/N; % 频域采样点
x = randn(1,N); % 随机OFDM信号
X = fft(x); % DFT计算
[Pxx, F] = pwelch(x, [], [], [], 1/Ts); % 计算功率谱密度
plot(F, Pxx); % 绘制功率谱密度曲线
xlabel('Frequency (Hz)');
ylabel('Power Spectral Density');