分别产生一段时长10us、带宽10MHz的LEM实信号和复信号采样率为230MHz。绘制对应波形的时域波形与频谱。假定4=1，=100MHz matlab代码

在MATLAB中，你可以使用`awgn`函数生成噪声，并结合`sin`函数创建正弦信号来模拟LEMM序列（一种线性相位移键控调制）。假设你想要创建的是离散时间信号，因为实际采样率远大于带宽，所以可以忽略信号的频率响应，主要关注时域波形和周期性的表示。

首先，我们需要设置一些变量：

% 参数设置
duration = 1e-7; % 时长 (秒)
bandwidth = 10e6; % 带宽 (Hz)
sample_rate = 230e6; % 采样率 (Hz)
noise_power = 0; % 如果需要添加噪声，将此设置为非零值
freq_offset = 0; % 频偏 (一般不需要)

% 创建时长为duration，采样率为sample_rate的向量
t = 0:1/sample_rate:duration;

% 正弦信号 (LEMM信号通常包含载波和数据)
carrier_freq = bandwidth / 2 + freq_offset;
data_signal = sin(2 * pi * carrier_freq * t); % 线性相位调制信号

% 如果需要，加上白噪声
if noise_power > 0
    data_signal = data_signal + sqrt(noise_power) * randn(size(t));
end

% 由于是10us时长，我们可以取前10us的数据（近似无损）
t_sub = t(1:round(sample_rate * duration)); % 时间轴截取
data_signal_sub = data_signal(1:round(sample_rate * duration));

% 绘制时域波形
figure;
plot(t_sub, data_signal_sub);
xlabel('Time (s)');
ylabel('Amplitude');
title('10us LEMM信号时域波形');

% 对于10MHz带宽，我们无法直接计算频率响应，因为我们通常处理连续信号，而这里的样本点有限
% 但是为了展示频谱，我们可以使用fft并仅考虑部分频率范围
fft_len = length(data_signal_sub);
freq_domain = fftshift(fft(data_signal_sub)) / fft_len;
freq_range = (-sample_rate/2:sample_rate/fft_len:sample_rate/2 - 1)';

% 绘制频谱
figure;
plot(freq_range, abs(fft_domain(1:(length(freq_range)/2+1))));
xlim([-bandwidth/2 bandwidth/2]);
xlabel('Frequency (Hz)');
ylabel('Magnitude');
title(['10MHz LEMM信号频谱 (230MHz采样)']);

%