是这个

好的，那么您可以在代码的最后加入以下内容来计算 `fsk_filt` 和 `st1` 的误码率：

% 计算误码率
error_fsk_filt = sum(abs(fsk_filt-data))/N; % 计算 FSK_filt 的误码率
error_st1 = sum(abs(st1-data))/N; % 计算 ST1 的误码率

% 显示误码率
disp(['FSK_filt 的误码率：', num2str(error_fsk_filt)]);
disp(['ST1 的误码率：', num2str(error_st1)]);

这段代码会计算 `fsk_filt` 和 `st1` 的误码率，并在屏幕上显示结果。完整的代码如下：

% 生成随机比特序列
N = 10000; % 比特序列长度
data = randi([0 1], N, 1); % 随机生成 0 和 1 的比特序列

% FSK 调制
fc1 = 1000; % 第一个载频频率
fc2 = 2000; % 第二个载频频率
fs = 8000; % 采样率
T = 1/fs; % 采样时间间隔
t = 0:T:(N/fs)-T; % 时间向量
s = fskmod(data, fc1, fc2, fs); % FSK 调制

% 信道加噪声
EbNo = 10; % 信噪比
snr = EbNo + 10*log10(log2(numel(unique(data)))); % 转换为信噪比
noisy_signal = awgn(s, snr, 'measured'); % 加入高斯白噪声

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

% ST1 解调
[b, a] = butter(6, 2*fc1/fs, 'low'); % 低通滤波器设计
st1 = filtfilt(b, a, noisy_signal); % 进行滤波
st1 = st1(1:T*fs:end); % 降采样

% 计算误码率
error_fsk_filt = sum(abs(fsk_filt-data))/N; % 计算 FSK_filt 的误码率
error_st1 = sum(abs(st1-data))/N; % 计算 ST1 的误码率

% 显示误码率
disp(['FSK_filt 的误码率：', num2str(error_fsk_filt)]);
disp(['ST1 的误码率：', num2str(error_st1)]);

希望对您有帮助！