function [alpha,theta]=channel_estimatorb(quad_train_seq, inphase_train_seq, orig_train, model, channel)%   [alpha,theta]=channel_estimatorb(quad_train_seq, inphase_train_seq, orig_train, model, channel)%%   Output:%   alpha		        - Alpha in channelmodel alpha = exp(j*theta)
%   theta               - Theta in dito%%   Input:%   quad_train_seq	    - Quadrature part of training-sequence as seen in receiver.
%   inphase_train_seq   - Inphase dito	 %   orig_train  	    - Original training-sequence. (Binary!) (on Antenna 1!)
%							   (not using training info from antenna 2.)%	 model			    - If 0: No Alamouti coding (only sending on antenna 1)
%								if 1: Alamouti coding is used for transmitting antennas.
%	 channel	        - 1 if referring to channel 0 (or 2) in Alamouti's paper.
%								2 if referring to channel 1 (or 3)%   %   Short Theoretical Background for the Function:%%   Only one channel is estimated!!! 
%	 
%	Estimates the one-tap channel. %   (Referring to Appendix A in George J鰊grens master of sciense%   thesis for mathematical theories)%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%     %%% Function part of simulation for Space-Time%%% coding project, group Grey-2001.%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%   Author: Stefan Uppg錼d%   Date: 2001-03-22%   Version: 1.1%   Revision (Name & Date):
%		SU	2001-03-24 Adding variable: model and modifying code for that%						  Also: changed inputargument "sent_training_seq" be binary.
%						  (not BPSK:ed)%       SU  2001-03-27 Added if for calculating theta.%       SU  2001-04-08 Put the mean-expressions on 1/abs(w) and angle(w) instead of on r*conj(r).%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
sent_train_seq = bpsk(orig_train,model);r = complex(inphase_train_seq, quad_train_seq);Sigma_rr = (r .* (conj(r)));Sigma_rd = (r .* sent_train_seq);w = Sigma_rd ./ Sigma_rr;alpha = mean(1./abs(w));        %%% The corresponding values for the channel h = 1/w*.Sigma_rr = mean(r .* (conj(r)));Sigma_rd = mean(r .* sent_train_seq);w = Sigma_rd / Sigma_rr;theta = angle(w);