请运用 MATLAB 语言，完成 2FSK 通信系统的仿真， 仿真它的误码性能，并与理论误码曲线相比较。 2 1 Tb f = f +1/ ，检测器 为平方律检测器，信道相移为 0。信道中的加性噪声可认为是统计独 立，均值为 0 的高斯分量。分别作出 2FSK 通信系统仿真得到的误码 率曲线和理论上的误码率曲线。

好的，以下是一个简单的 MATLAB 代码，用于模拟2FSK调制的通信系统，并绘制误码率曲线：

% 2FSK调制通信系统的误码率仿真

% 参数设置
Tb = 1; % 每个符号的时间长度
f1 = 1/Tb; % 第一个频率
f2 = 2/Tb; % 第二个频率
Eb = 1; % 符号能量
N = 100000; % 发送的符号数量
SNR_dB = 0:2:20; % 信噪比范围

% 产生随机二进制数据
data = randi([0 1], 1, N);

% 将二进制数据转换为符号序列
symbols = 2*data-1;

% 生成调制信号
t = 0:Tb/1000:(N*Tb)-Tb/1000;
mod_signal = sqrt(2*Eb/Tb)*cos(2*pi*f1*t + (symbols==1)*pi) + sqrt(2*Eb/Tb)*cos(2*pi*f2*t + (symbols==-1)*pi);

% 添加高斯噪声
for k = 1:length(SNR_dB)
    SNR = 10^(SNR_dB(k)/10);
    noise_power = Eb/SNR;
    noise = sqrt(noise_power/2)*randn(size(mod_signal));
    rx_signal = mod_signal + noise;

    % 解调信号
    demod1 = rx_signal.*cos(2*pi*f1*t);
    demod2 = rx_signal.*cos(2*pi*f2*t);
    s1 = sum(demod1 > 0);
    s2 = sum(demod2 > 0);
    detected_symbols = -1*(s1>s2) + 1*(s1<=s2);

    % 计算误码率
    errors = sum(detected_symbols ~= symbols);
    BER(k) = errors/N;
end

% 绘制误码率曲线
SNR_lin = 10.^(SNR_dB/10);
theory_BER = 0.5*erfc(sqrt(SNR_lin/2));
semilogy(SNR_dB, BER, 'o-', SNR_dB, theory_BER, 'r-');
xlabel('SNR (dB)');
ylabel('BER');
legend('Simulation', 'Theory');

上述代码中，我们首先生成了长度为N的随机二进制数据，然后将其转换为-1和1的符号序列。接着，我们使用2FSK调制将符号序列转换为模拟信号，并添加了高斯噪声。最后，我们使用平方律检测器对解调信号进行解调，并计算误码率。我们还通过绘制误码率曲线将仿真结果与理论误码率曲线进行比较。

请注意，上述代码中的参数设置和函数调用可能需要根据您的具体需求进行调整。