用matlab写一下

当然，我们可以使用MATLAB来编写代码以实现这个需求。首先，我们需要创建两个部分的信号，然后进行采样，最后绘制时域波形图和频谱图。假设我们已知正弦波的幅度、起始时间、结束时间和频率，以及采样率。以下是一个简单的示例：

% 定义参数
A = 1; % 正弦波幅度
T0 = 0; % 开始时间 (0 ms)
T1 = 100; % 结束时间 (100 ms)
f0 = 10.23e3; % 第一频率 (10.23 kHz)
f1 = 12.23e3; % 第二频率 (12.23 kHz)
fs_high = 50e3; % 高采样率 (50 Ksps)
fs_low = 5e3;   % 低采样率 (5 Ksps)

% 计算总时间长度和对应的采样点数
t_total = T1 - T0 + 0.001; % 加入一个小数值防止整数除法
num_samples_high = round(fs_high * t_total);
num_samples_low = round(fs_low * t_total);

% 创建时间向量
time_vec = linspace(T0, T0 + t_total, num_samples_high);

% 根据时间创建正弦波
signal_high = zeros(num_samples_high, 1);
signal_high(1:(num_samples_high*T0/fs_high)) = A*sin(2*pi*f0*fs_high*(0:T0/fs_high));
signal_high((num_samples_high*T0+1):end) = A*sin(2*pi*f1*fs_high*((T0+0.001):t_total));

% 对低采样率信号做同样的操作
signal_low = zeros(num_samples_low, 1);
signal_low(1:(num_samples_low*T0/fs_low)) = A*sin(2*pi*f0*fs_low*(0:T0/fs_low));
signal_low((num_samples_low*T0+1):end) = A*sin(2*pi*f1*fs_low*((T0+0.001):t_total));

% 绘制时域波形
subplot(2,1,1);
plot(time_vec, signal_high, 'b', time_vec, signal_low, 'g');
xlabel('Time (ms)');
ylabel('Amplitude');
legend('High Sampling Rate', 'Low Sampling Rate');

% 计算并绘制频谱图
[H_high, freq_high] = periodogram(signal_high, [], [], fs_high);
[H_low, freq_low] = periodogram(signal_low, [], [], fs_low);
subplot(2,1,2);
plot(freq_high, abs(H_high), 'b', freq_low, abs(H_low), 'g');
xlabel('Frequency (kHz)');
ylabel('Magnitude Spectrum');
title(['High Sampling Rate: ' num2str(fs_high) ' Hz, Low Sampling Rate: ' num2str(fs_low) ' Hz']);

% 显示图形
grid on;

这段代码首先创建了一个时间向量，然后按照给定频率生成正弦波，接着分别绘制了高和低采样率的时域波形图和频谱图。你可以运行这段代码并在MATLAB环境中查看结果。