错误使用 fskmod M must be a positive power of 2.

这个错误通常是由于使用了非 2 的幂次的调制阶数导致的。在使用 `fskmod` 函数时，必须将 `M` 设置为 2 的幂次方（例如 2，4，8，16 等），以确保调制阶数是整数。如果 `M` 不是 2 的幂次方，则会出现类似于“M must be a positive power of 2”这样的错误。

要解决这个问题，您可以将 `M` 设置为最接近您所需的调制阶数的 2 的幂次方。例如，如果您需要使用 10 阶的 FSK 调制，则可以将 `M` 设置为 16。这样可以确保调制阶数是整数，并且可以避免出现错误。

以下是一个示例代码，展示了如何将 `M` 设置为最接近所需调制阶数的 2 的幂次方：

% 生成随机比特序列
N = 1000; % 比特序列长度
M = 10; % FSK 调制阶数
data = randi([0 M-1], N, 1); % 随机生成 0 到 M-1 的整数比特序列

% 将 M 设置为最接近的 2 的幂次方
M_pow2 = 2^nextpow2(M);

% FSK 调制
fc1 = 1000; % 第一个载频频率
fc2 = 2000; % 第二个载频频率
fs = 8000; % 采样率
T = 1/fs; % 采样时间间隔
t = 0:T:(N/fs)-T; % 时间向量
s = fskmod(data, M_pow2, fc1, fc2, fs); % FSK 调制，使用最接近的 2 的幂次方

% 显示调制结果
subplot(2,1,1);
plot(t, s);
title('FSK 调制信号');

% FSK 解调
f1 = fskdemod(s, M_pow2, fc1, fs); % 解调得到第一个载频信号
f2 = fskdemod(s, M_pow2, fc2, fs); % 解调得到第二个载频信号
fsk_filt = medfilt1(f1-f2, 5); % 差分后进行中值滤波

% 显示解调结果
subplot(2,1,2);
plot(t, fsk_filt);
title('FSK 解调信号');

在这个示例代码中，我将 `M` 设置为 10，但由于 10 不是 2 的幂次方，我使用 `nextpow2` 函数来找到最接近的 2 的幂次方，即 16。然后，我使用 16 阶的 FSK 调制，并将得到的信号进行解调和滤波处理。