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

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

?? stbc_qpsk_ber.m

?? stbc編碼qpsk調(diào)制的誤碼率,需要的可以下載了學(xué)習(xí)一下
?? M
?? 第 1 頁 / 共 3 頁
字號: 
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,'-^')
%legend('SER','BER')
clc
clear
hold on
%SNR upto 20 dBs
EbNo=[0:2:20];

%N, M: number of transmit and receive antennas

N=4;
M=1;

%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,N);
                BIT=2;
 
tx_bits=data.';
temp1=[];
temp=[];
for i=1:N
    %[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 G4 STBC 
 x0=X(:,1);% required to verify a 1x1 system
 x1=X(:,1:4);
 x2(:,1)=-X(:,2); x2(:,2)=X(:,1); x2(:,3)=-X(:,4);x2(:,4)=X(:,3);
 x3(:,1)=-X(:,3); x3(:,2)=X(:,4); x3(:,3)=X(:,1);x3(:,4)=-X(:,2);
 x4(:,1)=-X(:,4); x4(:,2)=-X(:,3); x4(:,3)=X(:,2);x4(:,4)=X(:,1);
 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))+r4(:,n).*conj(H(:,4));
  z1_2(:,n)=conj(r5(:,n)).*H(:,1)+conj(r6(:,n)).*H(:,2)+conj(r7(:,n)).*H(:,3)+conj(r8(:,n)).*H(:,4);
  z1(:,n)=z1_1(:,n)+z1_2(:,n);
    z2_1(:,n)=r1(:,n).*conj(H(:,2))-r2(:,n).*conj(H(:,1))-r3(:,n).*conj(H(:,4))+r4(:,n).*conj(H(:,3));
  z2_2(:,n)=conj(r5(:,n)).*H(:,2)-conj(r6(:,n)).*H(:,1)-conj(r7(:,n)).*H(:,4)+conj(r8(:,n)).*H(:,3);
  z2(:,n)=z2_1(:,n)+z2_2(:,n);
  
  z3_1(:,n)=r1(:,n).*conj(H(:,3))+r2(:,n).*conj(H(:,4))-r3(:,n).*conj(H(:,1))-r4(:,n).*conj(H(:,2));
  z3_2(:,n)=conj(r5(:,n)).*H(:,3)+conj(r6(:,n)).*H(:,4)-conj(r7(:,n)).*H(:,1)-conj(r8(:,n)).*H(:,2);
  z3(:,n)=z3_1(:,n)+z3_2(:,n);
  
  z4_1(:,n)=r1(:,n).*conj(H(:,4))-r2(:,n).*conj(H(:,3))+r3(:,n).*conj(H(:,2))-r4(:,n).*conj(H(:,1));
  z4_2(:,n)=conj(r5(:,n)).*H(:,4)-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,'-<')
%legend('SER','BER')
clc
clear
hold on
%SNR upto 20 dBs
EbNo=[0:2:20];

%N, M: number of transmit and receive antennas
N=4;
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,N);
                BIT=2;
        
tx_bits=data.';
temp1=[];
temp=[];
for i=1:N 
   % [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 G4 STBC 
 x0=X(:,1);
 x1=X(:,1:4);
 x2(:,1)=-X(:,2); x2(:,2)=X(:,1); x2(:,3)=-X(:,4);x2(:,4)=X(:,3);
 x3(:,1)=-X(:,3); x3(:,2)=X(:,4); x3(:,3)=X(:,1);x3(:,4)=-X(:,2);
 x4(:,1)=-X(:,4); x4(:,2)=-X(:,3); x4(:,3)=X(:,2);x4(:,4)=X(:,1);
 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))+r4(:,n).*conj(H(:,4));
  z1_2(:,n)=conj(r5(:,n)).*H(:,1)+conj(r6(:,n)).*H(:,2)+conj(r7(:,n)).*H(:,3)+conj(r8(:,n)).*H(:,4);
  z1(:,n)=z1_1(:,n)+z1_2(:,n);
  
  z2_1(:,n)=r1(:,n).*conj(H(:,2))-r2(:,n).*conj(H(:,1))-r3(:,n).*conj(H(:,4))+r4(:,n).*conj(H(:,3));
  z2_2(:,n)=conj(r5(:,n)).*H(:,2)-conj(r6(:,n)).*H(:,1)-conj(r7(:,n)).*H(:,4)+conj(r8(:,n)).*H(:,3);
  z2(:,n)=z2_1(:,n)+z2_2(:,n);
  
  z3_1(:,n)=r1(:,n).*conj(H(:,3))+r2(:,n).*conj(H(:,4))-r3(:,n).*conj(H(:,1))-r4(:,n).*conj(H(:,2));
  z3_2(:,n)=conj(r5(:,n)).*H(:,3)+conj(r6(:,n)).*H(:,4)-conj(r7(:,n)).*H(:,1)-conj(r8(:,n)).*H(:,2);
  z3(:,n)=z3_1(:,n)+z3_2(:,n);
  
  z4_1(:,n)=r1(:,n).*conj(H(:,4))-r2(:,n).*conj(H(:,3))+r3(:,n).*conj(H(:,2))-r4(:,n).*conj(H(:,1));
  z4_2(:,n)=conj(r5(:,n)).*H(:,4)-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,'-+')
%legend('SER','BER')
legend('1x2','1x4','2x1','2x2','3x1','3x4','1x4','4x1','4x4')
clc
clear

?? 快捷鍵說明

復(fù)制代碼 Ctrl + C
搜索代碼 Ctrl + F
全屏模式 F11
切換主題 Ctrl + Shift + D
顯示快捷鍵 ?
增大字號 Ctrl + =
減小字號 Ctrl + -
亚洲欧美第一页_禁久久精品乱码_粉嫩av一区二区三区免费野_久草精品视频
亚洲国产日韩av| 国产欧美一区二区精品仙草咪| 日韩在线一区二区| 精品国精品自拍自在线| 色噜噜夜夜夜综合网| 精品亚洲成a人| 亚洲自拍偷拍九九九| 久久久激情视频| 51精品秘密在线观看| 99久久免费精品| 久久99国产精品久久99果冻传媒| 一区二区三区在线视频观看58| 欧美精品一区二区高清在线观看| 欧美性一二三区| 9久草视频在线视频精品| 国内一区二区视频| 日本va欧美va欧美va精品| 有码一区二区三区| 国产精品美日韩| 久久久久久毛片| 欧美一区二区久久| 欧美主播一区二区三区美女| 成人白浆超碰人人人人| 国产一区美女在线| 麻豆久久久久久| 日韩高清国产一区在线| 亚洲国产欧美在线人成| 日韩美女啊v在线免费观看| 国产三级精品三级在线专区| 欧美成人福利视频| 欧美不卡在线视频| 欧美一区二区三区影视| 欧美日韩高清一区二区不卡| 在线免费av一区| 在线观看日韩电影| 色哟哟亚洲精品| 色哟哟日韩精品| 欧美午夜精品一区二区三区| 欧美在线播放高清精品| 色狠狠色狠狠综合| 欧美在线999| 欧美日韩你懂的| 666欧美在线视频| 欧美一级xxx| 欧美变态tickle挠乳网站| 精品国产乱子伦一区| 久久网站最新地址| 国产欧美日韩另类视频免费观看| 中文天堂在线一区| 中文字幕一区二区三区视频| 亚洲欧美日韩一区二区| 亚洲在线视频网站| 日韩精品免费视频人成| 久久99精品国产麻豆婷婷洗澡| 激情六月婷婷久久| 成人免费视频视频| 91久色porny | 日本欧美在线看| 久久99精品久久久久久久久久久久| 久久狠狠亚洲综合| 成人永久看片免费视频天堂| av在线不卡电影| 在线一区二区观看| 欧美一区二区视频网站| 国产亚洲成av人在线观看导航| 中文字幕一区日韩精品欧美| 亚洲香肠在线观看| 久久综合综合久久综合| 国产成人一区二区精品非洲| 97久久超碰国产精品电影| 欧美色老头old∨ideo| 日韩丝袜美女视频| 中文字幕第一页久久| 一区二区三区在线观看动漫| 日韩1区2区日韩1区2区| 国产精品一品二品| 91国模大尺度私拍在线视频| 欧美一区二区高清| 亚洲欧美在线观看| 日本视频中文字幕一区二区三区| 国产精品一区二区视频| 欧美性欧美巨大黑白大战| 欧美精品一区二| 一区二区三区成人在线视频| 久草热8精品视频在线观看| 波多野洁衣一区| 日韩欧美二区三区| 亚洲精品高清在线| 国产精品一区在线| 欧美日韩精品高清| 日本一区二区不卡视频| 亚洲第一二三四区| 高清国产一区二区| 91麻豆精品91久久久久同性| 国产精品入口麻豆九色| 日日欢夜夜爽一区| caoporm超碰国产精品| 日韩视频国产视频| 亚洲女厕所小便bbb| 国产成人av网站| 欧美一卡二卡在线观看| 亚洲影院在线观看| 成人中文字幕合集| 欧美成人vr18sexvr| 丝袜国产日韩另类美女| 97aⅴ精品视频一二三区| 久久免费精品国产久精品久久久久| 亚洲一区免费视频| 色综合久久综合| 欧美国产精品中文字幕| 久久99久久99小草精品免视看| 欧美写真视频网站| 成人精品国产免费网站| 欧美一级片在线看| 午夜精品久久久久久久蜜桃app| 波多野结衣中文字幕一区| 日韩欧美另类在线| 亚洲成av人片在线观看无码| 一本一道久久a久久精品| 欧美国产在线观看| 国产精品系列在线播放| 日韩三级视频在线看| 婷婷开心激情综合| 欧美日韩一区二区三区视频| 亚洲精品国产精华液| 一本到三区不卡视频| 亚洲欧洲成人精品av97| 成人黄色一级视频| 国产精品天美传媒沈樵| 国产成人免费视频| 国产欧美一二三区| 国产jizzjizz一区二区| 国产喂奶挤奶一区二区三区| 国产在线播精品第三| 久久综合色之久久综合| 韩国午夜理伦三级不卡影院| 精品剧情在线观看| 国产一区二区三区综合| 久久久精品影视| 国产成人免费av在线| 欧美国产日本韩| 91在线视频观看| 亚洲精品第1页| 欧美午夜片在线观看| 午夜电影网亚洲视频| 日韩一区二区三区在线视频| 麻豆视频一区二区| 国产午夜亚洲精品不卡| 成人中文字幕电影| 亚洲靠逼com| 欧美日韩电影一区| 久久精品久久精品| 国产亚洲欧美日韩在线一区| 成人国产精品免费| 亚洲精品视频在线观看免费| 欧美三级在线看| 美女视频免费一区| 久久久www成人免费无遮挡大片| 国产成+人+日韩+欧美+亚洲| 亚洲精品写真福利| 在线观看91av| 国产夫妻精品视频| 亚洲精品国久久99热| 91麻豆精品国产91久久久久| 黑人巨大精品欧美一区| 国产精品久久毛片| 欧美日韩国产一级二级| 久久精品国产秦先生| 2024国产精品| 色就色 综合激情| 久久精品久久精品| 亚洲男人的天堂在线aⅴ视频| 欧美日韩国产综合一区二区| 国产精品一区二区三区99| 1024成人网| 欧美xxxxxxxxx| 色婷婷久久久久swag精品 | 欧美二区三区91| 国产精品一区二区在线播放 | 欧美一a一片一级一片| 国产原创一区二区| 一区二区三区日韩精品| 日韩一区二区精品在线观看| 成人深夜视频在线观看| 偷拍一区二区三区四区| 国产人妖乱国产精品人妖| 欧美午夜精品理论片a级按摩| 国产麻豆午夜三级精品| 亚洲一区二区在线免费观看视频 | 久久视频一区二区| 91久久人澡人人添人人爽欧美| 日韩成人一区二区三区在线观看| 国产精品第五页| 欧美电影免费观看高清完整版在线 | 亚洲蜜桃精久久久久久久| 欧美mv日韩mv国产| 欧美性色黄大片手机版| aaa亚洲精品一二三区| 国产一区不卡精品| 日本一区中文字幕|