（1）利用matlab产生载频fc=20Mhz，带宽B=1MHz，脉宽T=200us的LFM时域信号，并分析其频谱； （2）将该基带信号按照公式： st=y(t)*cos(2*pi*f*t) 进行调制，其中调制频率f为2GHz，利用matlab分析调制以后信号的时域和频域波形； （3）对s(t)进行解调，调制频率f为2GHz，利用matlab分析解调后的时域和频域波形；

以下是matlab代码实现：

（1）产生LFM信号并分析频谱：

% 信号参数
fc = 20e6; % 载频
B = 1e6; % 带宽
T = 200e-6; % 脉宽
fs = 100e6; % 采样频率
t = -T/2:1/fs:T/2; % 时间序列

% 产生LFM信号
y = chirp(t, -B/2, T/2, B/2, 'linear', 90);

% 分析频谱
Y = fftshift(fft(y));
f = linspace(-fs/2, fs/2, length(Y));
plot(f, abs(Y));
xlabel('频率（Hz）');
ylabel('幅度');

（2）将LFM信号调制并分析时域和频域波形：

% 调制参数
f = 2e9; % 调制频率

% 调制信号
st = y.*cos(2*pi*f*t);

% 分析时域波形
subplot(2,1,1);
plot(t, st);
xlabel('时间（s）');
ylabel('幅度');

% 分析频域波形
subplot(2,1,2);
S = fftshift(fft(st));
f_S = linspace(-fs/2, fs/2, length(S));
plot(f_S, abs(S));
xlabel('频率（Hz）');
ylabel('幅度');

（3）进行解调并分析时域和频域波形：

% 解调参数
f = 2e9; % 调制频率

% 解调信号
dt = t(2)-t(1);
s_I = st.*cos(2*pi*f*t)*2*dt;
s_Q = st.*sin(2*pi*f*t)*2*dt;
s = sqrt(s_I.^2 + s_Q.^2);

% 分析时域波形
subplot(2,1,1);
plot(t, s);
xlabel('时间（s）');
ylabel('幅度');

% 分析频域波形
subplot(2,1,2);
S = fftshift(fft(s));
f_S = linspace(-fs/2, fs/2, length(S));
plot(f_S, abs(S));
xlabel('频率（Hz）');
ylabel('幅度');

注意：以上代码仅供参考，具体实现可能需要根据实际情况进行调整。