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

? 歡迎來到蟲蟲下載站！ | ?? 資源下載 ?? 資源專輯 ?? 關于我們
? 蟲蟲下載站

?? sunview.cpp

?? MFC編寫的一個FIR的程序
?? CPP
字號: 
// sunView.cpp : implementation of the CSunView class
//

#include "stdafx.h"
#include "sun.h"

#include "sunDoc.h"
#include "sunView.h"
#include "complex.h"
#include <math.h>
#define PI_2    2*(double)3.14159265358979323846
#define PI      (double)3.14159265358979323846
#ifdef _DEBUG
#define new DEBUG_NEW
#undef THIS_FILE
static char THIS_FILE[] = __FILE__;
#endif

/////////////////////////////////////////////////////////////////////////////
// CSunView

IMPLEMENT_DYNCREATE(CSunView, CView)

BEGIN_MESSAGE_MAP(CSunView, CView)
	//{{AFX_MSG_MAP(CSunView)
	ON_COMMAND(ID_CHOSE, OnChose)
	ON_COMMAND(ID_DOFIR, OnDofir)
	//}}AFX_MSG_MAP
	// Standard printing commands
	ON_COMMAND(ID_FILE_PRINT, CView::OnFilePrint)
	ON_COMMAND(ID_FILE_PRINT_DIRECT, CView::OnFilePrint)
	ON_COMMAND(ID_FILE_PRINT_PREVIEW, CView::OnFilePrintPreview)
END_MESSAGE_MAP()

/////////////////////////////////////////////////////////////////////////////
// CSunView construction/destruction

CSunView::CSunView()
{
	// TODO: add construction code here

}

CSunView::~CSunView()
{
}

BOOL CSunView::PreCreateWindow(CREATESTRUCT& cs)
{
	// TODO: Modify the Window class or styles here by modifying
	//  the CREATESTRUCT cs

	return CView::PreCreateWindow(cs);
}

/////////////////////////////////////////////////////////////////////////////
// CSunView drawing

void CSunView::OnDraw(CDC* pDC)
{
	CSunDoc* pDoc = GetDocument();
	ASSERT_VALID(pDoc);
	// TODO: add draw code for native data here
	
if(cholg.flag==1)
	{
   
  
   int count;
   count=pow(2,nn);
   int dd=count;
   complex * x;
   double * ms;
   double *phase;
	
   
  


   
   x=new complex[count];
   
   ms=new double[count];
   
   phase=new double[count];

   if(cholg.m_xuan==0)load(x,nn);
   if(cholg.m_xuan==1)sanload(x,nn);
   if(cholg.m_xuan==2)fangload(x,nn);



   



   
   
	CString str;
	int i; 
    
	// 創建畫筆對象
	CPen* pPenRed=new CPen;
	// 紅色畫筆
	pPenRed->CreatePen(PS_SOLID,1,RGB(255,0,0));

	// 創建畫筆對象
	CPen* pPenBlue=new CPen;
    // 藍色畫筆
	pPenBlue->CreatePen(PS_SOLID,1,RGB(0,0, 255));




	
		//原序列垂直軸
		pDC->MoveTo(50,10);
		pDC->LineTo(50,230);

	    //原序列水平軸
	    pDC->MoveTo(50,110);
	    pDC->LineTo(350,110);

     	//繪制原序列X軸箭頭
	    pDC->MoveTo(345,105);
	    pDC->LineTo(350,110);
	    pDC->LineTo(345,115);
	
	    //繪制原序列Y軸箭頭
	    pDC->MoveTo(50,10);
	    pDC->LineTo(45,15);
	    pDC->MoveTo(50,10);
	    pDC->LineTo(55,15);

		str.Format("原序列");
    	pDC->TextOut(70,10,str);

		//FFT序列幅頻垂直軸
		pDC->MoveTo(400,10);
		pDC->LineTo(400,180);

	    //FFT序列幅頻水平軸
	    pDC->MoveTo(400,180);
	    pDC->LineTo(700,180);

     	//繪制FFT序列幅頻X軸箭頭
	    pDC->MoveTo(695,175);
	    pDC->LineTo(700,180);
	    pDC->LineTo(695,185);
	
	    //繪制FFT序列幅頻Y軸箭頭
	    pDC->MoveTo(400,10);
	    pDC->LineTo(395,15);
	    pDC->MoveTo(400,10);
	    pDC->LineTo(405,15);

		str.Format("原序列經FFT所得的幅頻特性");
    	pDC->TextOut(600,10,str);



		//FFT序列相頻垂直軸
		pDC->MoveTo(400,250);
		pDC->LineTo(400,450);

	    //FFT序列相頻水平軸
	    pDC->MoveTo(400,350);
	    pDC->LineTo(700,350);

     	//繪制FFT序列相頻X軸箭頭
	    pDC->MoveTo(695,345);
	    pDC->LineTo(700,350);
	    pDC->LineTo(695,355);
	
	    //繪制FFT序列相頻Y軸箭頭
	    pDC->MoveTo(400,250);
	    pDC->LineTo(395,255);
	    pDC->MoveTo(400,250);
	    pDC->LineTo(405,255);

		str.Format("原序列經FFT所得的相頻特性");
    	pDC->TextOut(420,250,str);
		
		
		
		
		//IFFT序列垂直軸
		pDC->MoveTo(50,250);
		pDC->LineTo(50,450);

	    //IFFT水平軸
	    pDC->MoveTo(50,350);
	    pDC->LineTo(350,350);




       //繪制IFFT序列X軸箭頭
	    pDC->MoveTo(345,345);
	    pDC->LineTo(350,350);
	    pDC->LineTo(345,355);

	
	    //繪制IFFT序列Y軸箭頭
	    pDC->MoveTo(50,250);
	    pDC->LineTo(45,255);
	    pDC->MoveTo(50,250);
	    pDC->LineTo(55,255);

		str.Format("原序列經FFT和IFFT所得序列");
    pDC->TextOut(70,250,str);





		// 選中當前藍色畫筆，并保存以前的畫筆
	    CGdiObject* pOldPen = pDC->SelectObject(pPenBlue);

		for(i=0;i<count;i++)
			{
				pDC->MoveTo(i+50,110);
        	    pDC->LineTo(i+50,110-int(x[i].ShowReal()));
			
			}


     fft(x,nn);

        
		for(i=0;i<count;i++){
			ms[i]=x[i].abs();
		}
       pDC->SelectObject(pPenRed);
	   
		    

			for(i=0;i<count;i++)
			{
			   pDC->MoveTo(i+400,180);
                pDC->LineTo(i+400,180-0.1*int(ms[i])); 

			}

		
		for(i=0;i<count;i++){
			phase[i]=x[i].arg();
		}

		for(i=0;i<count;i++)
			{
				pDC->MoveTo(i+400,350);
                pDC->LineTo(i+400,350-int(50*phase[i])); 
		}
       
		
		ifft(x,nn);
pDC->SelectObject(pPenBlue);

        for(i=0;i<count;i++)
			{
				pDC->MoveTo(i+50,350);
        	    pDC->LineTo(i+50,350-int(x[i].ShowReal())/dd);
			
			}




delete[]ms;
delete[]x;
delete[]phase;


	}





	if(firolg.flag==1){

	CString str;
	int i; 
    
	// 創建畫筆對象
	CPen* pPenRed=new CPen;
	// 紅色畫筆
	pPenRed->CreatePen(PS_SOLID,1,RGB(255,0,0));

	// 創建畫筆對象
	CPen* pPenBlue=new CPen;
    // 藍色畫筆
	pPenBlue->CreatePen(PS_SOLID,1,RGB(0,0, 255));




	
		//原序列垂直軸
		pDC->MoveTo(50,10);
		pDC->LineTo(50,230);

	    //原序列水平軸
	    pDC->MoveTo(50,110);
	    pDC->LineTo(350,110);

     	//繪制原序列X軸箭頭
	    pDC->MoveTo(345,105);
	    pDC->LineTo(350,110);
	    pDC->LineTo(345,115);
	
	    //繪制原序列Y軸箭頭
	    pDC->MoveTo(50,10);
	    pDC->LineTo(45,15);
	    pDC->MoveTo(50,10);
	    pDC->LineTo(55,15);

		str.Format("測試序列X(n)");
    	pDC->TextOut(70,10,str);

		//FFT序列幅頻垂直軸
		pDC->MoveTo(400,10);
		pDC->LineTo(400,180);

	    //FFT序列幅頻水平軸
	    pDC->MoveTo(400,180);
	    pDC->LineTo(700,180);

     	//繪制FFT序列幅頻X軸箭頭
	    pDC->MoveTo(695,175);
	    pDC->LineTo(700,180);
	    pDC->LineTo(695,185);
	
	    //繪制FFT序列幅頻Y軸箭頭
	    pDC->MoveTo(400,10);
	    pDC->LineTo(395,15);
	    pDC->MoveTo(400,10);
	    pDC->LineTo(405,15);

		str.Format("窗函數幅頻特性|W(ejw)|");
    	pDC->TextOut(600,10,str);



		//FFT序列相頻垂直軸
		pDC->MoveTo(400,250);
		pDC->LineTo(400,450);

	    //FFT序列相頻水平軸
	    pDC->MoveTo(400,450);
	    pDC->LineTo(700,450);

     	//繪制FFT序列相頻X軸箭頭
	    pDC->MoveTo(695,445);
	    pDC->LineTo(700,450);
	    pDC->LineTo(695,455);
	
	    //繪制FFT序列相頻Y軸箭頭
	    pDC->MoveTo(400,250);
	    pDC->LineTo(395,255);
	    pDC->MoveTo(400,250);
	    pDC->LineTo(405,255);

		str.Format("濾波器衰減特性20lg(|H(ejw)|");
    	pDC->TextOut(420,250,str);
		
		
		
		
		//IFFT序列垂直軸
		pDC->MoveTo(50,250);
		pDC->LineTo(50,450);

	    //IFFT水平軸
	    pDC->MoveTo(50,350);
	    pDC->LineTo(350,350);




       //繪制IFFT序列X軸箭頭
	    pDC->MoveTo(345,345);
	    pDC->LineTo(350,350);
	    pDC->LineTo(345,355);

	
	    //繪制IFFT序列Y軸箭頭
	    pDC->MoveTo(50,250);
	    pDC->LineTo(45,255);
	    pDC->MoveTo(50,250);
	    pDC->LineTo(55,255);

		str.Format("濾波后序列Y(n)");
    pDC->TextOut(70,250,str);





		// 選中當前藍色畫筆，并保存以前的畫筆
	    CGdiObject* pOldPen = pDC->SelectObject(pPenBlue);

	for(i=0;i<256;i++)
			{
				pDC->MoveTo(i+50,110);
        	    pDC->LineTo(i+50,110-int(50*Xn[i]));
			
			}

for(i=0;i<256;i++)
			{
				pDC->MoveTo(i+50,350);
        	    pDC->LineTo(i+50,350-(50*I[i]));
			
			}

pDC->SelectObject(pPenRed);
for(i=0;i<256;i++)
			{
				pDC->MoveTo(i+400,450);
                pDC->LineTo(i+400,450+int(Hd[i])); 
		}
  for(i=0;i<256;i++)
			{
			   pDC->MoveTo(i+400,180);
                pDC->LineTo(i+400,180-int(100*Hw[i])); 

			}     




	}






}

/////////////////////////////////////////////////////////////////////////////
// CSunView printing

BOOL CSunView::OnPreparePrinting(CPrintInfo* pInfo)
{
	// default preparation
	return DoPreparePrinting(pInfo);
}

void CSunView::OnBeginPrinting(CDC* /*pDC*/, CPrintInfo* /*pInfo*/)
{
	// TODO: add extra initialization before printing
}

void CSunView::OnEndPrinting(CDC* /*pDC*/, CPrintInfo* /*pInfo*/)
{
	// TODO: add cleanup after printing
}

/////////////////////////////////////////////////////////////////////////////
// CSunView diagnostics

#ifdef _DEBUG
void CSunView::AssertValid() const
{
	CView::AssertValid();
}

void CSunView::Dump(CDumpContext& dc) const
{
	CView::Dump(dc);
}

CSunDoc* CSunView::GetDocument() // non-debug version is inline
{
	ASSERT(m_pDocument->IsKindOf(RUNTIME_CLASS(CSunDoc)));
	return (CSunDoc*)m_pDocument;
}
#endif //_DEBUG

/////////////////////////////////////////////////////////////////////////////
// CSunView message handlers





void CSunView::fft(complex *x, int m)
{   complex *w;
    float num=0;
    int n=1;
	complex u,temp,tm,wrecur,w_s;
	complex *xi,*xj,*wptr;
	int i,j,k,l,le,windex;
	double arg;
	if(m==0)return;
	n=(int)pow(2,m);
    le=n/2;
	
	w=(complex*)calloc(le,sizeof(complex));
	
	
	arg=PI/le;//計算pi/le
    xj=w;
	
	w_s.copy(cos(arg),-sin(arg));
	wrecur.copy(1,0);
	for (j=0;j<le;j++){
		*xj=wrecur;	
		xj++;		
        wrecur=wrecur*w_s;
		}

//反向調整數據//
	j=0;
	for(i=1;i<n;i++){
	k=n;	
    for(int a=1;a<=m;a++){
		k=k/2;
		if(k<=j) j=j-k;
		else{j=j+k;break;}
		}

	if(i<j){
		xi=x+i;
		xj=x+j;
		temp =*xj;
		*xj=*xi;
		*xi=temp;
	}
	}
	
//fft計算//
   le=1;
   windex=n/2;
   for(l=0;l<m;l++){
	   le=le*2;
	   wptr=w;	
	   for(j=0;j<le/2;j++){
		   u=*wptr;
		   for(i=j;i<n;i=i+le){
			   xi=x+i;
			   xj=xi+le/2;
			   temp=*xi;
			   *xi=*xi+u*(*xj);
			   *xj=temp-u*(*xj);
		   }
		   wptr=wptr+windex;		   
	   }
	   windex=windex/2;
   }
   free(w);

}





void CSunView::load(complex *x,int n)
{int m=pow(2,n);
	for(int i=0;i<m;i++){
 	x[i].copy(70*sin(i*PI/(m/2)),0);
	}

}




void CSunView::sanload(complex *x,int n)
{int mm;
	int m=1,g=0;
	int d=0;
	mm=pow(2,n);

	for(int i=0;i<mm;i++){
		g=g+m;
		if(g==(mm/4)||g==0){
			m=m*(-1);	
		}
    d=20*g/(mm/16);
 	x[i].copy(d,0);
	}


}



void CSunView::fangload(complex *x,int n)
{int m=0,g=1;
int mm;
mm=pow(2,n);
	for(int i=0;i<mm;i++){
 	m++;
	if(m<(mm/2)){
		g=1;	
	}
	if(m>(mm/2)){
		g=-1;
	}

	
	x[i].copy(50*g,0);
	}

}







void CSunView::OnChose() 
{	if(cholg.DoModal()!=IDOK)
	{
		return;
	}
	// TODO: Add your command handler code here
if(cholg.m_jie==0)nn=4;
if(cholg.m_jie==1)nn=5;
if(cholg.m_jie==2)nn=6;
if(cholg.m_jie==3)nn=7;
if(cholg.m_jie==4)nn=8;
if(cholg.m_jie==5)nn=9;
cholg.flag=1;





Invalidate();

}



void CSunView::ifft(complex *x, int m)
{   complex *w;
    float num=0;
    int n=1;
	complex u,temp,tm,wrecur,w_s;
	complex *xi,*xj,*wptr;
	int i,j,k,l,le,windex;
	double arg;
	if(m==0)return;
	n=(int)pow(2,m);
    le=n/2;
	
	w=(complex*)calloc(le,sizeof(complex));
	if(!w){
		exit(1);
	}	
	
	arg=PI/le;//計算pi/le
    xj=w;
	
	w_s.copy(cos(arg),sin(arg));
	wrecur.copy(1,0);
	for (j=0;j<le;j++){
		*xj=wrecur;	
		xj++;		
        wrecur=wrecur*w_s;
		}
	//ifft計算//
   le=n;
   windex=1;
   for(l=0;l<m;l++){
	   le=le/2;
	   wptr=w;	
	   for(j=0;j<le;j++){
		   u=*wptr;
		   for(i=j;i<n;i=i+2*le){
			   xi=x+i;
			   xj=xi+le;
			   temp=*xi;
			   *xi=((*xi+(*xj)));
			   *xj=((temp-(*xj))*u);
			  

		   }
		   wptr=wptr+windex;		   
	   }
	   windex=windex*2;
   }

//反向調整數據//
	j=0;
	for(i=1;i<n;i++){
	k=n;	
    for(int a=1;a<=m;a++){
		k=k/2;
		if(k<=j) j=j-k;
		else{j=j+k;break;}
		}

	if(i<j){
		xi=x+i;
		xj=x+j;
		temp =*xj;
		*xj=*xi;
		*xi=temp;
	}
	}
   free(w);
}

void CSunView::OnDofir() {

	if(firolg.DoModal()!=IDOK)
	{
		return;
	}
	
     UpdateData();
	 int i;
	 m_a=(firolg.m_nn-1)/2;
     Wn=new double[firolg.m_nn];
	 Hd=new double[firolg.m_nn];
Hw=new double[256];
     Hn=new double[256];
     I=new double[256];
	 Xn=new double[256];
     Hd=new double[256];
Ww=new double[256];
	 xy=new complex[256];
	 orig=new complex[256];
     buf=new complex[256];
	 double Wc=PI*firolg.m_Wc;
	 for(i=0;i<256;i++){	  
		 Hn[i]=0;  Xn[i]=0;
	 }
    


	 for(i=0;i<firolg.m_nn;i++){	 //產生窗函數
		 switch(firolg.m_type){
		 case 0:  Wn[i]=1.0;    break;
		 case 1:  Wn[i]=0.54-0.46*cos(PI_2*i/(firolg.m_nn-1));   break;
		 case 2:  Wn[i]=0.5*(1-cos(PI_2*i/(firolg.m_nn-1)));   break;      
		 case 3:  Wn[i]=0.42-0.5*cos(PI_2*i/(firolg.m_nn-1))+0.08*cos(4*PI*i/(firolg.m_nn-1));   break;
		 default:      break;      
		 }

	// TODO: Add your command handler code here
		if(i!=m_a)
		Hd[i]=sin(Wc*(i-m_a))/(i-m_a)/PI;
		else    Hd[i]=Wc/PI;   
		Hn[i]=Wn[i]*Hd[i];     
	}


 for(int j=0;j<256;j++){
     xy[j].copy(Hn[j],0);
	 }

fft(xy,8);


for(i=0;i<256;i++){
     Hw[i]=xy[i].abs();
}

	
double model;
model=Hw[0];
	
for(i=0;i<256;i++)
		{Ww[i]=Hw[i]/model;
		
			}


for(i=0;i<256;i++)
	{ Hd[i]=20*log10(Ww[i]);	}

double low=PI*firolg.m_low;//待測序列
double high=PI*firolg.m_high;
        for(i=0;i<256;i++)
	 	Xn[i]=cos(i*low)+cos(i*high);

 for(i=0;i<256;i++){
     orig[i].copy(Xn[i],0);
	 }




fft(orig,8);

for(i=0;i<256;i++){//由時域卷積定理:時域卷積==>頻域乘積
	buf[i]=orig[i]*xy[i];
}

ifft(buf,8);//用付里葉反變換求得結果序列Y(n)---存入I

for(i=0;i<256;i++){
	I[i]=buf[i].ShowReal()/256;
}
firolg.flag=1;
Invalidate();
}

?? 快捷鍵說明

復制代碼 Ctrl + C
搜索代碼 Ctrl + F
全屏模式 F11
切換主題 Ctrl + Shift + D
顯示快捷鍵 ?
增大字號 Ctrl + =
減小字號 Ctrl + -
亚洲欧美第一页_禁久久精品乱码_粉嫩av一区二区三区免费野_久草精品视频
麻豆91在线看| 7777精品伊人久久久大香线蕉经典版下载 | 国产精一品亚洲二区在线视频| 国产一区二区三区免费看 | 欧美一级搡bbbb搡bbbb| 欧美怡红院视频| 久久天堂av综合合色蜜桃网| 国产精品美女久久久久高潮| 亚洲国产欧美日韩另类综合 | 色呦呦一区二区三区| 91精选在线观看| 成人免费在线观看入口| 青青草精品视频| 一本久久a久久免费精品不卡| 色综合天天综合网国产成人综合天| 北条麻妃一区二区三区| 日韩午夜激情免费电影| 亚洲卡通动漫在线| 国内精品久久久久影院薰衣草| 91麻豆国产香蕉久久精品| 91久久国产综合久久| 久久精品欧美一区二区三区不卡| 亚洲bdsm女犯bdsm网站| 秋霞电影一区二区| 日本精品免费观看高清观看| 日日夜夜精品视频免费| 91浏览器入口在线观看| 欧美精品黑人性xxxx| 一区二区三区资源| 视频在线观看一区二区三区| 视频一区视频二区在线观看| 毛片av一区二区| 国产精品资源站在线| 欧美精品在线视频| 亚洲一区二区在线观看视频| 99免费精品在线| 国产精品美女久久久久久| 国产精品一区二区久久精品爱涩 | 激情综合色播激情啊| 91麻豆精品国产91久久久久久久久| 日韩女优毛片在线| 1024精品合集| 麻豆一区二区三区| 日韩视频免费直播| 视频一区视频二区中文字幕| 欧美色大人视频| 亚洲成人激情综合网| 国产精品综合一区二区| 在线视频国内一区二区| 337p粉嫩大胆噜噜噜噜噜91av | 91国模大尺度私拍在线视频| 亚洲免费大片在线观看| 色乱码一区二区三区88| 一区二区在线看| 在线观看成人免费视频| 亚洲精品国产视频| 国产美女精品在线| 国产精品免费久久| 蜜臂av日日欢夜夜爽一区| 色综合久久综合| 精品国产凹凸成av人导航| 国模冰冰炮一区二区| 欧美色成人综合| 欧美a级一区二区| 91黄色免费看| 日韩综合一区二区| 欧美成人aa大片| 亚洲国产一区二区视频| eeuss国产一区二区三区| 亚洲另类一区二区| 欧美日韩亚洲综合在线 欧美亚洲特黄一级 | 欧美日韩国产免费一区二区| 中文字幕一区在线观看视频| 精品一区二区免费在线观看| 欧美欧美欧美欧美首页| 久99久精品视频免费观看| 国产日韩精品一区二区三区在线| 日本中文字幕不卡| 日韩毛片一二三区| 国产成a人无v码亚洲福利| 亚洲欧美一区二区三区极速播放| 欧美日韩三级在线| 亚洲第一福利一区| 色噜噜久久综合| 久久福利资源站| 日韩视频在线一区二区| 美女任你摸久久| 亚洲欧洲www| 欧美成人一区二区三区在线观看| aaa欧美色吧激情视频| 日韩国产精品大片| 欧美韩国日本不卡| 欧美日韩国产免费一区二区| 国产精品88888| 午夜精品福利一区二区三区av| 久久尤物电影视频在线观看| 色偷偷88欧美精品久久久| 亚洲激情一二三区| 在线亚洲免费视频| 性久久久久久久久久久久| 国产欧美精品一区二区三区四区| 欧美三级日韩在线| 日韩国产精品久久久久久亚洲| 国产精品毛片无遮挡高清| 99久久免费精品| 一区二区三区欧美亚洲| 欧美三级电影一区| 日韩精品电影在线观看| 欧美大片日本大片免费观看| 狠狠v欧美v日韩v亚洲ⅴ| 亚洲午夜三级在线| 日韩免费高清电影| 国产91丝袜在线18| 国产精品久久毛片av大全日韩| 日韩一级片在线播放| 国产一区二区电影| 免费成人结看片| 亚洲第一二三四区| 最新热久久免费视频| 欧美在线制服丝袜| 91亚洲午夜精品久久久久久| 亚洲超丰满肉感bbw| 中文字幕一区二区三区在线不卡| 欧美性大战久久久久久久| 91美女在线视频| 蜜臀va亚洲va欧美va天堂| 日韩综合一区二区| 视频一区中文字幕| 日韩一区日韩二区| 国产精品久久99| 亚洲欧洲在线观看av| 国产精品免费av| 国产欧美日韩久久| 国产日本亚洲高清| 国产日本欧洲亚洲| 国产精品久久久久一区| 91精品国产手机| 成人av网站免费| 日韩 欧美一区二区三区| 日韩黄色片在线观看| 免费看黄色91| 久久国产日韩欧美精品| 精品一区二区三区在线播放视频| 美女国产一区二区| 亚洲国产日韩一级| 国产精品萝li| 一区二区三区欧美日| 亚洲国产综合人成综合网站| 五月婷婷欧美视频| 日本不卡在线视频| 国产乱码精品一区二区三| 成人手机电影网| 精品在线播放午夜| 国产ts人妖一区二区| 91在线播放网址| 欧美视频中文字幕| 日韩欧美中文字幕制服| 91福利小视频| 欧美精品乱码久久久久久按摩| 欧美一级片在线| 欧美日韩免费不卡视频一区二区三区 | 91麻豆精品国产自产在线观看一区 | 欧美一级二级三级蜜桃| 欧美怡红院视频| 欧美一级高清片在线观看| 久久久另类综合| 亚洲伦理在线精品| 日本vs亚洲vs韩国一区三区| 午夜精品久久久久久久久| 一区二区三区中文字幕电影| 中文字幕亚洲精品在线观看| 亚洲成在人线在线播放| 经典三级在线一区| 蜜臀va亚洲va欧美va天堂| 日韩成人免费在线| 日韩激情在线观看| av高清不卡在线| www.欧美精品一二区| 欧美高清视频在线高清观看mv色露露十八| 奇米四色…亚洲| aaa亚洲精品| 91老司机福利 在线| 色婷婷精品大在线视频| 精品三级在线看| 亚洲最大成人网4388xx| 亚洲黄一区二区三区| 国产一区欧美二区| 精品视频一区二区三区免费| 91久久香蕉国产日韩欧美9色| 精品久久久久99| 精品国内片67194| www日韩大片| 亚洲123区在线观看| 日韩av成人高清| 色婷婷亚洲综合| 欧美激情一区在线观看| 国产精品国产三级国产三级人妇| 蜜臀av一级做a爰片久久| 国产乱妇无码大片在线观看| 宅男噜噜噜66一区二区66|