function y = channel(sig2,Mt,Mr,x,H,N);
%function y = channel(sig2,Mt,Mr,x,H,N)
%
%信道傳輸模型
%
%輸入： 
%   sig2----噪聲
%   Mt----發射天線數
%   Mr----接收天線數
%   x ----復數形式的輸入信號（對于MIMO，輸入時一個矩陣，每列都是某一時刻天線輸出值）
%   H----頻率選擇性信道
%   N----OFDM幀中的符號數
%
%輸出： 
%   y ----信道輸出（也是 MIMO矩陣，與x矩陣類似）

%產生噪聲序列（每行代表不同天線，列代表時間）噪聲是空間時間獨立的高斯百噪聲
noise = sqrt(sig2)*1/sqrt(2)*(randn(Mt*N,1)+j*randn(Mt*N,1));
y = H*x +noise;