小白求解答

clear all
close all
snr_in_dB=-4:2:8;
SNR_in_DB=-4:0.1:8;
for i=1:length(SNR_in_DB)
   SNR=10^(SNR_in_DB(i)/10);
   t_e_r(i)=Qfunct(sqrt(2*SNR));
end;
for i=1:length(snr_in_dB),     
   [class1_BER,class2_BER]=tchfsPe_AWGN(snr_in_dB(i));
end;
semilogy(snr_in_dB,class1_BER,'b-',snr_in_dB,class2_BER,'r*-');
title('高斯噪声下TCH/FS信道误码特性');
grid on;
xlabel('信噪比 Eb/No(dB)');
ylabel('误码率 Pe');
legend('class1','class2');

clear all
close all
snr_in_dB=-4:2:8;
SNR_in_DB=-4:0.1:8;
for i=1:length(SNR_in_DB)
   SNR=10^(SNR_in_DB(i)/10);
   t_e_r(i)=Qfunct(sqrt(2*SNR));
end;
for i=1:length(snr_in_dB),     
   [class1_BER,class2_BER]=tchfsPe_AWGN(snr_in_dB(i));
end;
semilogy(snr_in_dB,class1_BER,'b-',snr_in_dB,class2_BER,'r*-');
title('高斯噪声下TCH/FS信道误码特性');
grid on;
xlabel('信噪比 Eb/No(dB)');
ylabel('误码率 Pe');
legend('class1','class2');

function code=GSMconv_EN(msg)
TCH_FS_constraintlength=5;%约束矩阵
TCH_FS_Generator=[23 33];
trellis=poly2trellis(TCH_FS_constraintlength, TCH_FS_Generator);%生成网格
REG_Num=log2(trellis.numStates);
[code1,fstate1]=convenc(msg,trellis);
[code2,fstate2]=convenc(zeros(REG_Num,1),trellis,fstate1);%结尾比特
code=[code1 ;code2];

Error using vertcat
Dimensions of matrices being concatenated are not consistent.
Error in GSMconv_EN (line 8)
code=[code1 ;code2];
Error in tchfsPe_AWGN (line 19)
    inEncode=GSMconv_EN(exEncode);