fs=20000; dt=1/fs; t=0:dt:10; fc=1701.4; fd=11.4; fa=11; g=square(2*pi*fd*t,50); gs=cumsum(g)*dt; rfsk=cos(2*pi*fc*t+2*pi*fa*gs);

这段代码是用 MATLAB 编写的，主要实现了一个频率偏移键控（FSK）调制过程，最终输出一个以正弦波为载波的 FSK 调制信号 rfsk。

具体来说，代码中的 fs=20000 表示采样率为 20000Hz，dt=1/fs 表示采样间隔为 1/20000秒，t=0:dt:10 表示时间范围为 0到10秒，fc=1701.4 表示载波频率为 1701.4Hz，fd=11.4 表示频率偏移为 11.4Hz，fa=11 表示调制信号的幅度为 11。

接下来，使用 square 函数生成一个周期为 1/fd 的方波调制信号 g，并使用 cumsum 函数对 g 进行积分操作，得到一个累积方波调制信号 gs。

最后，将 gs 与 fc 和 fa 带入余弦函数，得到一个以正弦波为载波的 FSK 调制信号 rfsk。