在下面的代码的最后写计算出fsk_filt和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); % 差分后进行中值滤波

% 误码率计算
error = sum(abs(fsk_filt-data))/N; % 计算 FSK 误码率
disp(['FSK 误码率：', num2str(error)]); % 显示 FSK 误码率

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

% 误码率计算
error = sum(abs(st1-data(1:length(st1))))/length(st1); % 计算 ST1 误码率
disp(['ST1 误码率：', num2str(error)]); % 显示 ST1 误码率

这段代码首先生成了一个随机比特序列 `data`，然后进行 FSK 调制得到调制信号 `s`，加入高斯白噪声后得到有噪声的信号 `noisy_signal`。接着进行 FSK 解调，差分后再进行中值滤波得到 `fsk_filt`，计算 FSK 误码率并显示。然后进行 ST1 解调，滤波后进行降采样得到 `st1`，计算 ST1 误码率并显示。