matlab生成一个脉冲宽度7ps高斯脉冲，脉冲重复频率50MHz。然后做FFT，看DC到50GHz各个谱线的功率。

可以按照以下步骤生成脉冲并进行FFT分析：

1. 生成时间序列向量t，以7 ps为时间步长，生成长度为N=2^16的时间序列。代码如下：

dt = 7e-12; % 时间步长
t = (0:N-1)*dt; % 时间序列向量

2. 生成高斯脉冲信号s，采用matlab内置函数`normpdf`生成高斯分布，然后乘以一个脉冲宽度为7ps的矩形脉冲，再乘上一个脉冲重复频率为50MHz的周期函数，最终得到一个高斯脉冲信号。代码如下：

fc = 50e6; % 脉冲重复频率
pw = 7e-12; % 脉冲宽度
s = normpdf(t, 0, pw/2) .* rectpuls(t-pw/2, pw) .* cos(2*pi*fc*t);

3. 进行FFT变换，并计算功率谱密度。代码如下：

fs = 1/dt; % 采样频率
f = (-N/2:N/2-1)*fs/N; % 频率向量
S = fftshift(fft(s)/N).^2; % 傅里叶变换并平方
P = S.*conj(S)/(fs*N); % 功率谱密度

4. 绘制功率谱密度随频率的变化曲线，代码如下：

plot(f/1e9, 10*log10(P)); % 绘制功率谱密度曲线，单位为dBm/Hz
xlabel('Frequency(GHz)');
ylabel('Power Spectral Density(dBm/Hz)');
xlim([0 50]);

完整代码如下：

dt = 7e-12; % 时间步长
N = 2^16; % 时间序列长度
t = (0:N-1)*dt; % 时间序列向量

fc = 50e6; % 脉冲重复频率
pw = 7e-12; % 脉冲宽度
s = normpdf(t, 0, pw/2) .* rectpuls(t-pw/2, pw) .* cos(2*pi*fc*t); % 生成高斯脉冲信号

fs = 1/dt; % 采样频率
f = (-N/2:N/2-1)*fs/N; % 频率向量
S = fftshift(fft(s)/N).^2; % 傅里叶变换并平方
P = S.*conj(S)/(fs*N); % 功率谱密度

plot(f/1e9, 10*log10(P)); % 绘制功率谱密度曲线，单位为dBm/Hz
xlabel('Frequency(GHz)');
ylabel('Power Spectral Density(dBm/Hz)');
xlim([0 50]);

运行该代码，可以得到功率谱密度随频率的变化曲线，如下图所示：


其中，x轴为频率，单位为GHz，y轴为功率谱密度，单位为dBm/Hz。可以看到，在DC到50GHz范围内，高斯脉冲信号存在多个频率分量，且功率谱密度随频率呈现不同的分布特征。