matlab仿真的空时分组码，两根发射天线一根接收天线

%Simulation program to realize DSTBC transmission system on rayleigh channel clear all; %********************** preparation part *************************** Nt=2; %Number of TX Antennas Nr=1; %Number of RX Antennas Num_loop=100000; %Number of data per simulation level=1; %modulation level average=1; %取几次的平均 x1=1/sqrt(Nt); %两参考调制信号 x2=1/sqrt(Nt); X_REF=[x1 x2; -conj(x2) conj(x1)]; X_ref=X_REF; k=1; for ebn0=-2:16 sgma = sqrt((1/Nt)/(2*level*(10^(ebn0/10)))); N_ref = sgma * (randn(Nr,2) + j*randn(Nr,2)); sig_err=0; for jj=1:average for i=1:Num_loop %Starting the simulation %***************************Data generation************************ Source_data=randint(1,Nt,2^level); %****************************modulation**************************** Modul_data=pskmod(Source_data,2^level)/sqrt(Nt); %*****************************coefficient************************** R=Modul_data*X_REF'; %***************************next transmit signal****************** Xt=R*X_ref; x_ref=[Xt(1) Xt(2); -conj(Xt(2)) conj(Xt(1))]; %******************channel***************************************** h=sqrt(0.5)*(randn(Nr,Nt)+j*randn(Nr,Nt)); H=[h(1) conj(h(2));h(2) -conj(h(1))]; N = sgma * (randn(Nr,2) + j*randn(Nr,2)); %******************receive signal********************************** r1=X_ref(1,:)*H + [N_ref(1) conj(N_ref(2))]; r2=X_ref(2,:)*H + [N_ref(2) -conj(N_ref(1))]; %%参考的接收矢量 r3=x_ref(1,:)*H + [N(1) conj(N(2))]; %********************de coefficient******************************** RR=[r3*r1' r3*r2']; %********************de coefficient******************************** Demodul_data=RR*X_REF; %******************decoded signal********************************** Decoding=pskdemod(Demodul_data*sqrt(Nt),2^level); X_ref=x_ref; N_ref=N; %*******************number of err_signal**************************** [number,ratio]=symerr(Source_data,Decoding); sig_err=sig_err+number; end end SER_DSTBC_BPSK_2Tx1Rx(k)=sig_err/(Num_loop*Nt*average) k=k+1; end %********************** Output result *************************** ebn0=-2:16; semilogy(ebn0,SER_DSTBC_BPSK_2Tx1Rx,'-om'); %******************** end of file *********************************