如何使用matlab画出带宽100Khz，时长100ms信号的LFM形式的时频图啊

Question

如何使用matlab画出带宽100Khz，时长100ms信号的LFM形式的时频图啊，还有16位Costas频率编码波形的图

Answer 1

该回答引用ChatGPT GPT-4

运行结果

代码

% 参数设置
M = 16; % 阶数
N = 100; % 码元数
Tc = 1e-3; % 码元周期1ms
fc = 100e3; % 载波频率100kHz
t_costas = 0:1/Fs:(Tc-1/Fs); % 码元时间向量

% 生成Costas编码
costas_code = randi([0, M-1], 1, N);

% 生成Costas编码波形
costas_signal = [];
for ii = 1:N
    costas_signal = [costas_signal, exp(1i * 2 * pi * (fc + costas_code(ii) * B / M) * t_costas)];
end


% 计算短时傅里叶变换（STFT）
[s_costas, f_costas, t_stft_costas] = spectrogram(costas_signal, window, noverlap, nfft, Fs, 'yaxis');

% 绘制时频图
figure;
surf(t_stft_costas, f_costas / 1e3, 10 * log10(abs(s_costas)), 'EdgeColor', 'none');
axis xy;
axis tight;
view(0, 90);
xlabel('Time (s)');
ylabel('Frequency (kHz)');
title('16-bit Costas-encoded Signal Time-Frequency Representation');
colorbar;
colormap(jet);

% 参数设置
Fs = 1e6; % 采样频率1MHz
T = 100e-3; % 信号时长100ms
t = 0:1/Fs:T-1/Fs; % 时间向量
B = 100e3; % 带宽100kHz
f0 = 100e3; % 起始频率100kHz
k = B/T; % 调频斜率

% 生成LFM信号
LFM_signal = exp(1i * 2 * pi * (f0 * t + 0.5 * k * t.^2));
% 计算短时傅里叶变换（STFT）
window = hamming(256); % 窗函数
noverlap = 128; % 重叠长度
nfft = 1024; % FFT长度
[s, f, t_stft] = spectrogram(LFM_signal, window, noverlap, nfft, Fs, 'yaxis');

% 绘制时频图
figure;
surf(t_stft, f / 1e3, 10 * log10(abs(s)), 'EdgeColor', 'none');
axis xy;
axis tight;
view(0, 90);
xlabel('Time (s)');
ylabel('Frequency (kHz)');
title('LFM Signal Time-Frequency Representation');
colorbar;
colormap(jet);