亚洲欧美第一页_禁久久精品乱码_粉嫩av一区二区三区免费野_久草精品视频

? 歡迎來到蟲蟲下載站！ | ?? 資源下載 ?? 資源專輯 ?? 關(guān)于我們
? 蟲蟲下載站

?? stbc_qpsk_ber.m

?? stbc編碼qpsk調(diào)制的誤碼率,需要的可以下載了學(xué)習(xí)一下
?? M
?? 第 1 頁 / 共 3 頁
字號: 

%**************************************************************************


%QPSK modulation
   
   % Angle [pi/4 3*pi/4 -3*pi/4 -pi/4] corresponds to 
   % Gray code vector [00 10 11 01], respectively.
   table=exp(j*[-3/4*pi 3/4*pi 1/4*pi -1/4*pi]);  % generates QPSK symbols
   table=table([0 1 3 2]+1); % Gray code mapping pattern for QPSK symbols
   inp=reshape(bits_in(i,:),2,full_len/2);
   mod_symbols=table([2 1]*inp+1);  % maps transmitted bits into QPSK symbols
   P=4;% 4 constellation points
   %%%%%%%%%%%%%%%%%%%%%%%%%%%
   temp1=mod_symbols;
   s=table;
    temp=[temp; temp1];
    temp1=0;
end

%ready to transmit symbols of length 'K'
X=temp.';

fr_length=length(X);
 
 % block coding-Alamouti
 x0=X(:,1);% required to verify a 1x1 system
 x1=X;
 x2(:,1)=-conj(X(:,2));
 x2(:,2)=conj(X(:,1));


 % form the channel matrix
 for n=1:N
  %   if channel_model=='AWGN    '
   %      Hr(n,:,:)=ones(fr_length,N);
    % else
         Hr(n,:,:)=(randn(fr_length,N)+j*randn(fr_length,N))/sqrt(2);
   %  end
 end
      
for n=1:M
%transmission matrix   
 H=reshape(Hr(n,:,:),fr_length,N);
 Habs(:,n)=sum(abs(H).^2,2);

%received signal per receiver antenna   
 r1(:,n)=sum(H.*x1,2)/sqrt(N)+sqrt(sigma)*(randn(fr_length,1)+j*randn(fr_length,1));
 r2(:,n)=sum(H.*x2,2)/sqrt(N)+sqrt(sigma)*(randn(fr_length,1)+j*randn(fr_length,1));
 
% demodulate the received signals
 z1(:,n)=r1(:,n).*conj(H(:,1))+conj(r2(:,n)).*H(:,2);
 z2(:,n)=r1(:,n).*conj(H(:,2))-conj(r2(:,n)).*H(:,1);
end

%uncoded(1,1)
r01=H(:,1).*x0+sqrt(sigma)*(randn(fr_length,1)+j*randn(fr_length,1));


%form estimates
for m=1:P
  d01(:,m)=abs(r01-H(:,1)*s(m)).^2; % uncoded signal
%coded signals  
  d1(:,m)=abs(sum(z1,2)-s(m)).^2+(-1+sum(Habs,2))*abs(s(m))^2;
  d2(:,m)=abs(sum(z2,2)-s(m)).^2+(-1+sum(Habs,2))*abs(s(m))^2;
end
     
% determine the minimum of estimates      

%decision for detecting uncoded
      [y0,i0]=min((d01),[],2);
      s0d=s(i0).';
      clear d01 
      
%decision for detecting s1     
      [y1,i1]=min(d1,[],2);
      s1d=s(i1).';
      clear d1
%decision for detecting s2
      [y2,i2]=min(d2,[],2);
      s2d=s(i2).';
      clear d2
% form received symbols    
    Xd=[s1d s2d];

%determine symbol errors   
 error_un(packet_count)=sum(X(:,1)~=s0d);% for uncoded
 temp1=X>0;
 temp2=Xd>0;
 error(packet_count)=sum(sum(temp1~=temp2));% for coded
       end % end of FOR loop for "packet_count"
%calculate FER, SER and BER for current idx

%for uncoded signal
  SER_uncoded(idx)=sum(error_un)/(Num*K);
  BER_uncoded(idx)=SER_uncoded(idx)/BIT;
  FER_uncoded(idx)=SER_uncoded(idx)*K;

%for coded signal
  SER(idx)=sum(error)/(Num*K);
  BER(idx)=SER(idx)/BIT;
  FER(idx)=SER(idx)*K;
  
%increment idx  
  idx=idx + 1; 
end % end of FOR loop for SNR
%semilogy(EbNo,SER,'*')
%hold on
semilogy(EbNo,BER)
%legend('SER','BER')
clc
clear
hold on
%SNR upto 20 dBs
EbNo=[0:2:20];

%N, M: number of transmit and receive antennas

%N=no_tx_antennas;
%M=no_rx_antennas;
N=3;
M=1;

%initialize packet_count
idx=1;
K=1000;  %符號數(shù)
Num=100;
%channel_model='AWGN ';%信道模型是高斯信道
model='Rayleigh';%信道模型是Rayleigh
for SNR=EbNo
    sigma=0.5/(10^(SNR/10));
       for packet_count=1:Num
% we are interested in transmitting 'K' SYMBOLS not bits. Hence, K*2 for QPSK
% etc.
                data=randint(K*2,4);
                BIT=2;
        
tx_bits=data.';
temp1=[];
temp=[];
for i=1:4 %N is actually 3, but we require 4 tx antenna coding to be transmitted from 3 antennas.
   % [temp1 s P]=tx_modulate(tx_bits(i,:),modulation);
   %%%%%%%%%%%%%%%%%%%%%%%%%%
     %%%%%%%%%%%%%%%
   bits_in(i,:)=tx_bits(i,:);
   full_len = length(bits_in(i,:));
    % Angle [pi/4 3*pi/4 -3*pi/4 -pi/4] corresponds to 
   % Gray code vector [00 10 11 01], respectively.
   table=exp(j*[-3/4*pi 3/4*pi 1/4*pi -1/4*pi]);  % generates QPSK symbols
   table=table([0 1 3 2]+1); % Gray code mapping pattern for QPSK symbols
   inp=reshape(bits_in(i,:),2,full_len/2);
   mod_symbols=table([2 1]*inp+1);  % maps transmitted bits into QPSK symbols
   P=4;% 4 constellation points
   temp1=mod_symbols;
   s=table;
   %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    temp=[temp; temp1];
    temp1=0;
end

%ready to transmit symbols of length 'K'
X=temp.';
fr_length=length(X);
 
%block signals in the lth time slot-Block coding for G3 STBC 
 x0=X(:,1);% required to verify a 1x1 system
 x1=X(:,1:3);
 x2(:,1)=-X(:,2); x2(:,2)=X(:,1); x2(:,3)=-X(:,4);
 x3(:,1)=-X(:,3); x3(:,2)=X(:,4); x3(:,3)=X(:,1);
 x4(:,1)=-X(:,4); x4(:,2)=-X(:,3); x4(:,3)=X(:,2);
 x5=conj(x1);
 x6=conj(x2);
 x7=conj(x3);
 x8=conj(x4);
 
 % form the channel matrix
 for n=1:N
    % if channel_model=='AWGN    '
     %    Hr(n,:,:)=ones(fr_length,N);
    % else
         Hr(n,:,:)=(randn(fr_length,N)+j*randn(fr_length,N))/sqrt(2);
     %end
 end
        

for n=1:M
%transmission matrix   
H=reshape(Hr(n,:,:),fr_length,N);
Habs(:,n)=sum(abs(H).^2,2);

 %received signal per receiver antenna   
 r1(:,n)=sum(H.*x1,2)/sqrt(N)+sqrt(sigma)*(randn(fr_length,1)+j*randn(fr_length,1));
 r2(:,n)=sum(H.*x2,2)/sqrt(N)+sqrt(sigma)*(randn(fr_length,1)+j*randn(fr_length,1));
 
 r3(:,n)=sum(H.*x3,2)/sqrt(N)+sqrt(sigma)*(randn(fr_length,1)+j*randn(fr_length,1));
 r4(:,n)=sum(H.*x4,2)/sqrt(N)+sqrt(sigma)*(randn(fr_length,1)+j*randn(fr_length,1));
 
 r5(:,n)=sum(H.*x5,2)/sqrt(N)+sqrt(sigma)*(randn(fr_length,1)+j*randn(fr_length,1));
 r6(:,n)=sum(H.*x6,2)/sqrt(N)+sqrt(sigma)*(randn(fr_length,1)+j*randn(fr_length,1));
 
 r7(:,n)=sum(H.*x7,2)/sqrt(N)+sqrt(sigma)*(randn(fr_length,1)+j*randn(fr_length,1));
 r8(:,n)=sum(H.*x8,2)/sqrt(N)+sqrt(sigma)*(randn(fr_length,1)+j*randn(fr_length,1));

  % demodulate the received signals
  z1_1(:,n)=r1(:,n).*conj(H(:,1))+r2(:,n).*conj(H(:,2))+r3(:,n).*conj(H(:,3));
  z1_2(:,n)=conj(r5(:,n)).*H(:,1)+conj(r6(:,n)).*H(:,2)+conj(r7(:,n)).*H(:,3);
  z1(:,n)=z1_1(:,n)+z1_2(:,n);
  
  z2_1(:,n)=r1(:,n).*conj(H(:,2))-r2(:,n).*conj(H(:,1))+r4(:,n).*conj(H(:,3));
  z2_2(:,n)=conj(r5(:,n)).*H(:,2)-conj(r6(:,n)).*H(:,1)+conj(r8(:,n)).*H(:,3);
  z2(:,n)=z2_1(:,n)+z2_2(:,n);
  
  z3_1(:,n)=r1(:,n).*conj(H(:,3))-r3(:,n).*conj(H(:,1))-r4(:,n).*conj(H(:,2));
  z3_2(:,n)=conj(r5(:,n)).*H(:,3)-conj(r7(:,n)).*H(:,1)-conj(r8(:,n)).*H(:,2);
  z3(:,n)=z3_1(:,n)+z3_2(:,n);
  
  z4_1(:,n)=-r2(:,n).*conj(H(:,3))+r3(:,n).*conj(H(:,2))-r4(:,n).*conj(H(:,1));
  z4_2(:,n)=-conj(r6(:,n)).*H(:,3)+conj(r7(:,n)).*H(:,2)-conj(r8(:,n)).*H(:,1);
  z4(:,n)=z4_1(:,n)+z4_2(:,n);
end

%uncoded(1,1)
r01=H(:,1).*x0+sqrt(sigma)*(randn(fr_length,1)+j*randn(fr_length,1));

  
%form estimates
for m=1:P
   d01(:,m)=abs(r01-H(:,1)*s(m)).^2; %uncoded signal
%coded signals       
   d1(:,m)=abs(sum(z1,2)-s(m)).^2+(-1+sum(Habs,2))*abs(s(m))^2;
   d2(:,m)=abs(sum(z2,2)-s(m)).^2+(-1+sum(Habs,2))*abs(s(m))^2;
   d3(:,m)=abs(sum(z3,2)-s(m)).^2+(-1+sum(Habs,2))*abs(s(m))^2;
   d4(:,m)=abs(sum(z4,2)-s(m)).^2+(-1+sum(Habs,2))*abs(s(m))^2;
end
% determine the minimum of estimates      

%decision for detecting uncoded
      [y0,i0]=min((d01),[],2);
      s0d=s(i0).';
      
      clear d01 
      
%decision for detecting s1     
      [y1,i1]=min(d1,[],2);
      s1d=s(i1).';
      
      clear d1
%decision for detecting s2     
      [y2,i2]=min(d2,[],2);
      s2d=s(i2).';
      
      clear d2
%decision for detecting s3     
      [y3,i3]=min(d3,[],2);
      s3d=s(i3).';
      
      clear d3
%decision for detecting s4     
      [y4,i4]=min(d4,[],2);
      s4d=s(i4).';
      
      clear d4

% form received symbols    
    Xd=[s1d s2d s3d s4d];
%determine symbol errors   
error_un(packet_count)=sum(X(:,1)~=s0d);% for uncoded
temp1=X>0;
temp2=Xd>0;
error(packet_count)=sum(sum(temp1~=temp2));% for coded
       end% end of FOR loop for "packet_count"
%calculate FER, SER and BER for current idx

%for uncoded signal
  SER_uncoded(idx)=sum(error_un)/(Num*K);
  BER_uncoded(idx)=SER_uncoded(idx)/BIT;
  FER_uncoded(idx)=SER_uncoded(idx)*K;

%for coded signal
  SER(idx)=sum(error)/(Num*K);
  BER(idx)=SER(idx)/BIT;
  FER(idx)=SER(idx)*K;
  
%increment idx  
  idx=idx + 1;
  
end% end of FOR loop for SNR
%semilogy(EbNo,SER,'g')
%hold on
semilogy(EbNo,BER,'o')
%legend('SER','BER')
clc
clear
hold on
%SNR upto 20 dBs
EbNo=[0:2:20];

%N, M: number of transmit and receive antennas

N=3;
M=4;

%initialize packet_count
idx=1;

 Num=100;
K=1000;

for SNR=EbNo
    sigma=0.5/(10^(SNR/10));
       for packet_count=1:Num
% we are interested in transmitting 'K' SYMBOLS not bits. Hence, K*2 for QPSK
% etc.
                data=randint(K*2,4);
                BIT=2;
           
        
tx_bits=data.';
temp1=[];
temp=[];
for i=1:4 %N is actually 3, but we require 4 tx antenna coding to be transmitted from 3 antennas.
  %  [temp1 s P]=tx_modulate(tx_bits(i,:),modulation);
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
       bits_in(i,:)=tx_bits(i,:);
   full_len = length(bits_in(i,:));
   %QPSK modulation
   
   % Angle [pi/4 3*pi/4 -3*pi/4 -pi/4] corresponds to 
   % Gray code vector [00 10 11 01], respectively.
   table=exp(j*[-3/4*pi 3/4*pi 1/4*pi -1/4*pi]);  % generates QPSK symbols
   table=table([0 1 3 2]+1); % Gray code mapping pattern for QPSK symbols
   inp=reshape(bits_in(i,:),2,full_len/2);
   mod_symbols=table([2 1]*inp+1);  % maps transmitted bits into QPSK symbols
   P=4;% 4 constellation points
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
      temp1=mod_symbols;
   s=table;
    temp=[temp; temp1];
    temp1=0;
end

%ready to transmit symbols of length 'K'
X=temp.';
fr_length=length(X);
 
%block signals in the lth time slot-Block coding for G3 STBC 
 x0=X(:,1);
 x1=X(:,1:3);
 x2(:,1)=-X(:,2); x2(:,2)=X(:,1); x2(:,3)=-X(:,4);
 x3(:,1)=-X(:,3); x3(:,2)=X(:,4); x3(:,3)=X(:,1);
 x4(:,1)=-X(:,4); x4(:,2)=-X(:,3); x4(:,3)=X(:,2);
 x5=conj(x1);
 x6=conj(x2);
 x7=conj(x3);
 x8=conj(x4);
 
 % form the channel matrix
 for n=1:N
    % if channel_model=='AWGN    '
      %   Hr(n,:,:)=ones(fr_length,N);
     %else
         Hr(n,:,:)=(randn(fr_length,N)+j*randn(fr_length,N))/sqrt(2);
    % end
 end
        

for n=1:N
%transmission matrix   
H=reshape(Hr(n,:,:),fr_length,N);
Habs(:,n)=sum(abs(H).^2,2);

 %received signal per receiver antenna   
 r1(:,n)=sum(H.*x1,2)/sqrt(N)+sqrt(sigma)*(randn(fr_length,1)+j*randn(fr_length,1));
 r2(:,n)=sum(H.*x2,2)/sqrt(N)+sqrt(sigma)*(randn(fr_length,1)+j*randn(fr_length,1));
 
 r3(:,n)=sum(H.*x3,2)/sqrt(N)+sqrt(sigma)*(randn(fr_length,1)+j*randn(fr_length,1));
 r4(:,n)=sum(H.*x4,2)/sqrt(N)+sqrt(sigma)*(randn(fr_length,1)+j*randn(fr_length,1));
 
 r5(:,n)=sum(H.*x5,2)/sqrt(N)+sqrt(sigma)*(randn(fr_length,1)+j*randn(fr_length,1));
 r6(:,n)=sum(H.*x6,2)/sqrt(N)+sqrt(sigma)*(randn(fr_length,1)+j*randn(fr_length,1));
 
 r7(:,n)=sum(H.*x7,2)/sqrt(N)+sqrt(sigma)*(randn(fr_length,1)+j*randn(fr_length,1));
 r8(:,n)=sum(H.*x8,2)/sqrt(N)+sqrt(sigma)*(randn(fr_length,1)+j*randn(fr_length,1));

  % demodulate the received signals
  z1_1(:,n)=r1(:,n).*conj(H(:,1))+r2(:,n).*conj(H(:,2))+r3(:,n).*conj(H(:,3));
  z1_2(:,n)=conj(r5(:,n)).*H(:,1)+conj(r6(:,n)).*H(:,2)+conj(r7(:,n)).*H(:,3);
  z1(:,n)=z1_1(:,n)+z1_2(:,n);
  
  z2_1(:,n)=r1(:,n).*conj(H(:,2))-r2(:,n).*conj(H(:,1))+r4(:,n).*conj(H(:,3));
  z2_2(:,n)=conj(r5(:,n)).*H(:,2)-conj(r6(:,n)).*H(:,1)+conj(r8(:,n)).*H(:,3);
  z2(:,n)=z2_1(:,n)+z2_2(:,n);
  
  z3_1(:,n)=r1(:,n).*conj(H(:,3))-r3(:,n).*conj(H(:,1))-r4(:,n).*conj(H(:,2));
  z3_2(:,n)=conj(r5(:,n)).*H(:,3)-conj(r7(:,n)).*H(:,1)-conj(r8(:,n)).*H(:,2);
  z3(:,n)=z3_1(:,n)+z3_2(:,n);
  
  z4_1(:,n)=-r2(:,n).*conj(H(:,3))+r3(:,n).*conj(H(:,2))-r4(:,n).*conj(H(:,1));
  z4_2(:,n)=-conj(r6(:,n)).*H(:,3)+conj(r7(:,n)).*H(:,2)-conj(r8(:,n)).*H(:,1);
  z4(:,n)=z4_1(:,n)+z4_2(:,n);
  
  
end

%uncoded(1,1)
r01=H(:,1).*x0+sqrt(sigma)*(randn(fr_length,1)+j*randn(fr_length,1));

  
%form estimates
for m=1:P
   d01(:,m)=abs(r01-H(:,1)*s(m)).^2; %uncoded signal
%coded signals       
   d1(:,m)=abs(sum(z1,2)-s(m)).^2+(-1+sum(Habs,2))*abs(s(m))^2;
   d2(:,m)=abs(sum(z2,2)-s(m)).^2+(-1+sum(Habs,2))*abs(s(m))^2;
   d3(:,m)=abs(sum(z3,2)-s(m)).^2+(-1+sum(Habs,2))*abs(s(m))^2;
   d4(:,m)=abs(sum(z4,2)-s(m)).^2+(-1+sum(Habs,2))*abs(s(m))^2;

?? 快捷鍵說明

復(fù)制代碼 Ctrl + C
搜索代碼 Ctrl + F
全屏模式 F11
切換主題 Ctrl + Shift + D
顯示快捷鍵 ?
增大字號 Ctrl + =
減小字號 Ctrl + -
亚洲欧美第一页_禁久久精品乱码_粉嫩av一区二区三区免费野_久草精品视频
欧美国产综合一区二区| 国产乱子伦视频一区二区三区| 日韩国产欧美三级| 精品在线播放午夜| a4yy欧美一区二区三区| 欧美日韩综合在线| 久久网站热最新地址| 自拍av一区二区三区| 男男视频亚洲欧美| 99热99精品| 91精品国产综合久久福利| 国产精品色呦呦| 日本色综合中文字幕| 成人高清视频免费观看| 777a∨成人精品桃花网| 中文字幕在线一区免费| 日日嗨av一区二区三区四区| 成人ar影院免费观看视频| 欧美一区在线视频| 亚洲色图另类专区| 激情丁香综合五月| 欧美日韩一区二区在线视频| 久久亚区不卡日本| 性做久久久久久免费观看| 成人av资源站| 欧美成人激情免费网| 亚洲精品欧美激情| 国产一区二区三区高清播放| 欧美日本不卡视频| 一区在线观看视频| 国产精品影音先锋| 91精品国产色综合久久不卡蜜臀| 亚洲欧美一区二区视频| 久久国产免费看| 欧美军同video69gay| 亚洲欧美激情一区二区| 粉嫩av一区二区三区粉嫩 | 国产精品激情偷乱一区二区∴| 午夜精品视频一区| 一本久久综合亚洲鲁鲁五月天| 国产亚洲一区二区三区| 日韩成人精品在线| 欧美日韩国产区一| 曰韩精品一区二区| 99久久婷婷国产综合精品| 久久久噜噜噜久久人人看| 日韩**一区毛片| 欧美日韩亚洲综合在线 | 欧美做爰猛烈大尺度电影无法无天| 久久久久久久久久电影| 日本午夜精品视频在线观看 | 日本伊人色综合网| 欧美网站一区二区| 亚洲美女屁股眼交3| 99国产精品国产精品久久| 国产清纯白嫩初高生在线观看91 | 久久亚洲一区二区三区四区| 美女任你摸久久| 7777精品伊人久久久大香线蕉经典版下载 | 日本成人在线电影网| 在线免费观看日韩欧美| 亚洲你懂的在线视频| 97精品电影院| 亚洲精选视频在线| 日本高清视频一区二区| 亚洲另类色综合网站| 色综合欧美在线视频区| 亚洲精品国产精华液| 欧美主播一区二区三区美女| 亚洲欧美日韩国产综合| 色综合中文字幕国产| 亚洲免费观看高清在线观看| 91小视频免费观看| 亚洲一区二区高清| 欧美日韩夫妻久久| 日本亚洲电影天堂| 精品免费视频.| 九九在线精品视频| 久久久天堂av| 成人毛片视频在线观看| 亚洲丝袜精品丝袜在线| 色婷婷亚洲综合| 午夜精品123| 日韩欧美一区在线| 奇米影视一区二区三区| 亚洲精品在线电影| 成人理论电影网| 一区二区在线观看视频在线观看| 色婷婷久久99综合精品jk白丝| 亚洲午夜久久久久久久久电影网| 欧美日韩精品一区视频| 日韩不卡手机在线v区| 午夜精品久久久久久久蜜桃app| 91精品国产一区二区三区香蕉| 极品尤物av久久免费看| 久久久91精品国产一区二区精品| youjizz国产精品| 亚洲午夜激情网页| 欧美一区二区在线观看| 国产一区二区精品久久91| 国产精品嫩草久久久久| 欧美在线视频全部完| 蜜桃久久久久久| 欧美国产精品一区二区三区| 91在线视频免费观看| 午夜电影一区二区三区| 久久色.com| 91在线观看下载| 日韩成人精品在线观看| 国产三级欧美三级| 一本大道久久a久久综合婷婷| 日韩中文字幕麻豆| 亚洲国产经典视频| 欧美精三区欧美精三区| 国产成人免费av在线| 亚洲午夜私人影院| 久久品道一品道久久精品| 色婷婷av一区二区| 韩国女主播一区| 一区二区理论电影在线观看| 欧美mv日韩mv| 在线区一区二视频| 国产一区二区三区不卡在线观看| 亚洲男女一区二区三区| 日韩精品一区二区三区视频 | 一区二区三区四区在线| 欧美mv日韩mv国产网站app| 色悠久久久久综合欧美99| 青青草国产精品亚洲专区无| 中文字幕一区二区三区精华液| 91精品国产综合久久久蜜臀粉嫩 | 亚洲电影视频在线| 国产婷婷一区二区| 777久久久精品| 91在线你懂得| 国内精品久久久久影院一蜜桃| 一区二区三区精品视频| 国产欧美日本一区二区三区| 欧美久久一二三四区| 中文字幕av一区 二区| 欧美一级久久久| 欧美亚洲精品一区| 不卡在线观看av| 久久99精品国产| 亚洲成av人影院| 成人欧美一区二区三区白人 | 99久久婷婷国产综合精品电影| 麻豆成人av在线| 五月综合激情网| 自拍偷拍亚洲激情| 国产欧美日韩在线看| 欧美成人性战久久| 欧美午夜免费电影| 99久久久无码国产精品| 国产精品1区2区3区在线观看| 日韩国产欧美在线播放| 亚洲青青青在线视频| 中文字幕成人av| 久久综合久久综合九色| 91精品国产综合久久精品性色 | 亚洲国产精品久久久男人的天堂| 国产精品私人自拍| 日本一区二区三区久久久久久久久不 | 国产亚洲短视频| 精品久久久久一区二区国产| 欧美视频一区二区三区| 波多野结衣中文字幕一区二区三区| 久久97超碰色| 免费成人深夜小野草| 日韩av电影免费观看高清完整版| 亚洲黄色小视频| 亚洲激情中文1区| 亚洲欧美日韩在线| 亚洲少妇30p| 亚洲少妇屁股交4| 综合色中文字幕| 成人欧美一区二区三区白人| 国产精品美女久久久久久2018| 中国av一区二区三区| 国产精品乱人伦| 日韩码欧中文字| 一区二区三区中文字幕精品精品| 一区二区三区四区在线播放 | 一区二区三区在线观看动漫| 日本亚洲欧美天堂免费| 精品卡一卡二卡三卡四在线| 制服丝袜日韩国产| 日韩三区在线观看| 99久久伊人久久99| av综合在线播放| 成人动漫一区二区三区| 国产99精品视频| 一区二区三区在线免费视频| 欧美日韩卡一卡二| 欧美亚洲日本国产| 欧美综合亚洲图片综合区| 色呦呦网站一区| 97久久久精品综合88久久| 久久国产生活片100| 国产精品乱码妇女bbbb|