sttc.asv

nv_encd.m 卷积编码程序 viterbi.m 卷积译码程序 其它的是viterbi.m中用到的子函数 程序来自《现代通信系统－使用matlab》英文版 已经调通！并加上了注释

clc;
clear;
SNR=[0:0.1:20];
for i=1:210
input=randsrc(1,1300,[0 1]);
g=[0 2 1 0;2 2 0 1];
k0=2;
output1=cnv_encd(g,k0,input);
output1=pskmod(output1,4,pi/4);
outuptabs=abs(output1);
E0=sum(outputabs)/651;
N0=(E0/(10^(0.1*SNR(i))))*randn(2,651);
h=mimo_channel(2,2,0.5*10^(-6)*[0:650]);
for i=1:651
    channel_output(:,i)=h(:,:,i)*output1(:,i);
end
channel_output=channel_output+N0;
[decoder_output,survivor_state,cumulated_metric]=viterbi(g,k0,h,channel_output);
error_num=sum(mod(decoder_output+input,2));
errorbit=error_num/1300;
erro(i)=errorbit;
end