青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品

posts - 12,  comments - 54,  trackbacks - 0
一般來(lái)說(shuō)，對(duì)均勻間隔采樣的數(shù)據(jù)的功率譜估計(jì)，可以采用DFT/AR/MA/ARMA/MUSIC等方法估計(jì)；
如果由于實(shí)驗(yàn)條件限制，或者偶爾的數(shù)據(jù)丟失，采樣的數(shù)據(jù)間隔時(shí)間／空間不是均勻間隔的，此時(shí)，大致可以有兩種處理方式：
1）用插值的方法將數(shù)據(jù)轉(zhuǎn)換為等間隔的，然后按等間隔采樣方法處理；這種方法最大的缺陷是，在低頻成分處會(huì)出現(xiàn)假的凸起。
2）采用Lomb方法處理；
本文將給出Lomb方法公式推導(dǎo)的結(jié)果，和C++代碼實(shí)現(xiàn)；至于具體推導(dǎo)過(guò)程，請(qǐng)參閱

Jerry D. Scargle, Studies in Astronomical Time Series Analysis. II. Statistical Aspect of Spectral Analysis of Unevenly Spaced Data, the Astrophysical Journal, vol. 263, pp. 835--853

這里不再多說(shuō)。

Lomb方法

假設(shè)數(shù)據(jù)點(diǎn)集為[;(h_i, t_i);]（所有數(shù)學(xué)公式采用latex語(yǔ)法書(shū)寫(xiě)），

定義其數(shù)學(xué)期望和方差分別為

[;h=\frac{1}{N} \sum_{i} h_i \\ \delta^2=

定義歸一化周期圖為

[;P_N(\omega) = \frac{1}{2 \delta ^2} \{ \frac{[\sum_i (h_i - h) \cos \omega(t_i - \tau)]^2 {\sum_i \cos^2 \omega (t_i - \tau)} + \frac{[\sum_i (h_i - h) \sin \omega(t_i - \tau)]^2 {\sum_i \sin^2 \omega (t_i - \tau)} \} ;]

其中[;\tau;]由下式定義

[; \tan 2 \omega \tau =

Lomb方法的C++實(shí)現(xiàn)：

?1?//lomb.hpp
?2?
?3?#ifndef?LOMB_HPP_INCLUDED__
?4?#define?LOMB_HPP_INCLUDED__
?5?
?6?#include?<vector>
?7?using?std::vector;
?8?
?9?//Warning:
10?//??If?you?would?rather?a?self-defined?type?than
11?//??a?builtin?type,?these?methods?MUST?be?offered:
12?//
13?//??operator=(?const?Type&?other?)
14?//??operator=(?const?long?double&?val?)
15?//??operator+(?const?Type&?lhs,?const?Type&?rhs)
16?//??operator-(?const?Type&?lhs,?const?Type&?rhs)
17?//??operator*(?const?Type&?lhs,?const?Type&?rhs)
18?//??operator/(?const?Type&?lhs,?const?Type&?rhs)
19?//??sin(?const?Type&?val?)
20?//??cos(?const?Type&?val?)
21?//
22?template
23?<
24?????typename?Type?=?long?double
25?>
26?class?Lomb
27?{
28?????public:
29?????????vector<Type>?spectrum()?const;
30?
31?????public:
32?????????Lomb(???const?vector<Type>&?_h,
33?????????????????const?vector<Type>&?_t?);
34?????????~Lomb();
35?
36?????private:
37?????????struct?Lomb_Data;
38?????????Lomb_Data*?lomb_data;
39?};
40?
41?#include?"lomb.tcc"
42?
43?#endif?//?LOMB_HPP_INCLUDED__
44

?


//lomb.tcc

#include?"iterator.hpp"

#include?<cassert>
#include?<cmath>

#include?<numeric>
#include?<iterator>
#include?<algorithm>
using?namespace?std;

template
<
????typename?Type
>
struct?Lomb<Type>::Lomb_Data
{
????vector<Type>?h;
????vector<Type>?t;

????Type?f()?const;
????Type?tau(?const?unsigned?int?)?const;
????Type?e()?const;
????Type?d()?const;

};

template
<
????typename?Type
>
Type?Lomb<Type>::Lomb_Data::f()const
{
????const?Type?max?=?*(?max_element(
????????????t.begin(),?t.end()?)?);
????const?Type?min?=?*(?min_element(
????????????t.begin(),?t.end()?)?);
????const?unsigned?int?N?=?h.size();
????const?Type?ans?=?N?/?(?max?-?min?)?;

????return?ans;
}

template
<
????typename?Type
>
Type?Lomb<Type>::Lomb_Data::tau(?const?unsigned?int?n?)const
{

????const?Type?omega?=
????????2.0L?*?3.14159265358979L?*?(?this?->?f()?)?*?n?;
????const?Type?sin_?=
????????sin2_accumulate(?t.begin(),
?????????????????????????t.end(),
?????????????????????????omega
????????????????);
????const?Type?cos_?=
????????cos2_accumulate(?t.begin(),
?????????????????????????t.end(),
?????????????????????????omega
????????????????);
????const?Type?ans?=?atan(?sin_?/?cos_?)?/?(?2.0L?*?omega?);

????return?ans;
}

template
<
????typename?Type
>
Type?Lomb<Type>::Lomb_Data::e()const
{
????const?Type?ans?=
????????mean(
????????????h.begin(),
????????????h.end(),
????????????Type_Keep<Type>()
????????????);

????return?ans;
}

template
<
????typename?Type
>
Type?Lomb<Type>::Lomb_Data::d()const
{
????const?Type?ans?=
????????variance(
????????????h.begin(),
????????????h.end(),
????????????Type_Keep<Type>()
????????????);

????return?ans;
}

template
<
????typename?Type
>
Lomb<Type>::Lomb(?const?vector<Type>&?_h,
????????????const?vector<Type>&?_t
??????????)
{
????assert(?_t.size()?==?_h.size()?);
????lomb_data?=?new?Lomb_Data;
????lomb_data?->?h?=?_h;
????lomb_data?->?t?=?_t;
}

template
<
????typename?Type
>
Lomb<Type>::~Lomb()
{
????delete?lomb_data;
}

template
<
????typename?Type
>
vector<Type>?Lomb<Type>::spectrum()?const
{
????const?unsigned?int?N?=?(*lomb_data).h.size();

????const?Type?f?=?lomb_data?->?f();
????const?Type?e?=?lomb_data?->?e();

????const?Type?d?=?lomb_data?->?d();

????vector<Type>?ans;

????for?(?unsigned?int?i?=?0;?i?<?N;?++i?)
????{
????????const?Type?tau?=?lomb_data?->?tau(?i?);
????????const?Type?omega?=?6.283185307179586476L?*
????????????????????????????????i?*?f;

????????const?Type?h_cos_square?=
????????????h_cos_square_accumulate(
????????????????????????????????????????(*lomb_data).h.begin(),?(*lomb_data).h.end(),
????????????????????????????????????????(*lomb_data).t.begin(),?(*lomb_data).t.end(),
????????????????????????????????????????e,
????????????????????????????????????????tau,
????????????????????????????????????????omega
????????????????????????????????????);

????????const?Type?cos_square?=
????????????cos_square_accumulate(
????????????????????????????????????(*lomb_data).t.begin(),
????????????????????????????????????(*lomb_data).t.end(),
????????????????????????????????????tau,
????????????????????????????????????omega
?????????????????????????????????);

????????const?Type?h_sin_square?=
????????????h_sin_square_accumulate(
????????????????????????????????????????(*lomb_data).h.begin(),?(*lomb_data).h.end(),
????????????????????????????????????????(*lomb_data).t.begin(),?(*lomb_data).t.end(),
????????????????????????????????????????e,
????????????????????????????????????????tau,
????????????????????????????????????????omega
????????????????????????????????????);

????????const?Type?sin_square?=
????????????sin_square_accumulate(
????????????????????????????????????(*lomb_data).t.begin(),
????????????????????????????????????(*lomb_data).t.end(),
????????????????????????????????????tau,
????????????????????????????????????omega
?????????????????????????????????);

????????ans.push_back(?sqrt(
????????????????????????????????(
????????????????????????????????????h_cos_square/cos_square?+
????????????????????????????????????h_sin_square/sin_square
????????????????????????????????)?/?(2.0L*d)
????????????????????????????)?/?N
?????????????????????);
????}

????return?ans;
}

?

??1?//iterator.hpp
??2?
??3?#ifndef?ITERATOR_HPP_INCLUDED__
??4?#define?ITERATOR_HPP_INCLUDED__
??5?
??6?//this?struct?is?just?employed
??7?//to?keep?the?Type?information
??8?template
??9?<
?10?????typename?Type
?11?>
?12?struct?Type_Keep;
?13?
?14?
?15?//calculate
?16?//??????\sum_{i}?\sin?2?\omega?t_i
?17?template
?18?<
?19?????typename?InputItor,
?20?????typename?Type
?21?>
?22?Type?sin2_accumulate(
?23?????????????????????????InputItor?begin,
?24?????????????????????????InputItor?end,
?25?????????????????????????Type?omega
?26?????????????????????);
?27?
?28?//calculate
?29?//??????\sum_{i}?\cos?2?\omega?t_i
?30?template
?31?<
?32?????typename?InputItor,
?33?????typename?Type
?34?>
?35?Type?cos2_accumulate(
?36?????????????????????????InputItor?begin,
?37?????????????????????????InputItor?end,
?38?????????????????????????Type?omega
?39?????????????????????);
?40?
?41?//calculate
?42?//??????\frac{1}{N}?\sum_{i}?h_i
?43?template
?44?<
?45?????typename?InputItor,
?46?????typename?Type
?47?>
?48?Type?mean(
?49?????????????InputItor?begin,
?50?????????????InputItor?end,
?51?????????????Type_Keep<Type>
?52??????????);
?53?
?54?//calculate
?55?//??????\frac{1}{N-1}?\sum_{i}(h_i?-?h)^{2}
?56?template
?57?<
?58?????typename?InputItor,
?59?????typename?Type
?60?>
?61?Type?variance(
?62?????????????????InputItor?begin,
?63?????????????????InputItor?end,
?64?????????????????Type_Keep<Type>
?65??????????????);
?66?
?67?//calcuate
?68?//???????[\sum_i?(h_i?-?h)?\sin?\omega?(t_i?-?\tau)]^{2}
?69?template
?70?<
?71?????typename?Type,
?72?????typename?InputItor
?73?>
?74?Type?h_sin_square_accumulate(?????InputItor?h_start,?InputItor?h_end,
?75?????????????????????????????????InputItor?t_start,?InputItor?t_end,
?76?????????????????????????????????Type?h,
?77?????????????????????????????????Type?tau,
?78?????????????????????????????????Type?omega
?79?????????????????????????????);
?80?
?81?//calcuate
?82?//???????[\sum_i?(h_i?-?h)?\cos?\omega?(t_i?-?\tau)]^{2}
?83?template
?84?<
?85?????typename?Type,
?86?????typename?InputItor
?87?>
?88?Type?h_cos_square_accumulate(?????InputItor?h_start,?InputItor?h_end,
?89?????????????????????????????????InputItor?t_start,?InputItor?t_end,
?90?????????????????????????????????Type?h,
?91?????????????????????????????????Type?tau,
?92?????????????????????????????????Type?omega
?93?????????????????????????????);
?94?
?95?//calculate
?96?//??????\sum_{i}?\cos?^{2}?\omega?(t_i?-?\tau)
?97?template
?98?<
?99?????typename?Type,
100?????typename?InputItor
101?>
102?Type?cos_square_accumulate(?????InputItor?t_start,
103?????????????????????????????????InputItor?t_end,
104?????????????????????????????????Type?tau,
105?????????????????????????????????Type?omega
106?????????????????????????????);
107?
108?//calculate
109?//??????\sum_{i}?\cos?^{2}?\omega?(t_i?-?\tau)
110?template
111?<
112?????typename?Type,
113?????typename?InputItor
114?>
115?Type?sin_square_accumulate(?????InputItor?t_start,
116?????????????????????????????????InputItor?t_end,
117?????????????????????????????????Type?tau,
118?????????????????????????????????Type?omega
119?????????????????????????????);
120?
121?#include?"iterator.tcc"
122?
123?#endif?//?ITERATOR_HPP_INCLUDED__
124?

?

//iterator.tcc

#include?<cmath>

template
<
????typename?Type
>
struct?Type_Keep
{
????typedef?Type?T;
};

template
<
????typename?Type
>
struct?Sin
{
????Type?operator()(?const?Type?val?)const
????{
????????return?sin(?val?);
????}
};

template
<
????typename?Type
>
struct?Cos
{
????Type?operator()(?const?Type?val?)const
????{
????????return?cos(?val?);
????}
};

template
<
????typename?InputItor,
????typename?Type,
????typename?Function
>
static?Type?f2_accumulate(
????????????????????????????InputItor?begin,
????????????????????????????InputItor?end,
????????????????????????????Type?omega,
????????????????????????????Function?f
????????????????????????)
{
????Type?ans?=?Type();
????while(?begin?!=?end?)
????{
????????ans?+=?f(?2.0L?*?omega?*?(*begin++)?);
????}
????return?ans;
}

template
<
????typename?InputItor,
????typename?Type
>
Type?sin2_accumulate(
????????InputItor?begin,
????????InputItor?end,
????????Type?omega
????????)
{
????return?f2_accumulate(
????????????????????????????begin,
????????????????????????????end,
????????????????????????????omega,
????????????????????????????Sin<Type>()
?????????????????????????);
}

template
<
????typename?InputItor,
????typename?Type
>
Type?cos2_accumulate(
????????InputItor?begin,
????????InputItor?end,
????????Type?omega
????????)
{
????return?f2_accumulate(
????????????????????????????begin,
????????????????????????????end,
????????????????????????????omega,
????????????????????????????Sin<Type>()
????????????????????????);
}

template
<
????typename?InputItor,
????typename?Type
>
Type?mean(
????????InputItor?begin,
????????InputItor?end,
????????Type_Keep<Type>
????????)
{
????Type?sum?=?Type();
????unsigned?int?cnt?=?0;
????while?(?begin?!=?end?)
????{
????????sum?+=?*begin++;
????????++cnt;
????}

????return?sum?/?cnt;
}

template
<
????typename?InputItor,
????typename?Type
>
static?Type?diff_square_accumulate(
????????????????????InputItor?begin,
????????????????????InputItor?end,
????????????????????Type?mean
????????????????????)
{
????Type?ans?=?Type();
????while(?begin?!=?end?)
????{
????????ans?+=?(?mean?-?*begin?)?*?(?mean?-?*begin?)?;
????????++begin;
????}
????return?ans;
}

template
<
????typename?InputItor
>
static?unsigned?int?difference(
????????????????InputItor?begin,
????????????????InputItor?end
????????????????)
{
????unsigned?int?cnt?=?0;
????while(?begin++?!=?end?)
????????++cnt;

????return?cnt;
}

template
<
????typename?InputItor,
????typename?Type
>
Type?variance(
????????????InputItor?begin,
????????????InputItor?end,
????????????Type_Keep<Type>
????????????)
{
????const?Type?average?=
????????mean(?begin,?end,?Type_Keep<long?double>()?);
????const?Type?power?=
????????diff_square_accumulate(?begin,?end,?average?);
????const?unsigned?int?size?=
????????difference(?begin,?end?);
????const?Type?ans?=?power?/?(size-1);

????return?ans;
}

template
<
????typename?Type,
????typename?InputItor,
????typename?Function
>
static?Type?h_f_square_accumulate(?????InputItor?h_start,?InputItor?h_end,
????????????????????????????????????InputItor?t_start,?InputItor?t_end,
????????????????????????????????????Type?h,
????????????????????????????????????Type?tau,
????????????????????????????????????Type?omega,
????????????????????????????????????Function?f
????????????????????????????)
{
????Type?ans?=?Type();
????while(?(?h_start?!=?h_end?)?&&?(?t_start?!=?t_end?)?)
????{
????????ans?+=?(?*h_start++?-?h?)?*?f(?omega?*?(?*t_start++?-?tau)?);
????}
????return?ans*ans;
}


template
<
????typename?Type,
????typename?InputItor
>
Type?h_sin_square_accumulate(?????InputItor?h_start,?InputItor?h_end,
????????????????????????????????InputItor?t_start,?InputItor?t_end,
????????????????????????????????Type?h,
????????????????????????????????Type?tau,
????????????????????????????????Type?omega
????????????????????????????)
{
????return?h_f_square_accumulate(
????????????????????????????????????h_start,?h_end,
????????????????????????????????????t_start,?t_end,
????????????????????????????????????h,
????????????????????????????????????tau,
????????????????????????????????????omega,
????????????????????????????????????Sin<Type>()
????????????????????????????????);
}

template
<
????typename?Type,
????typename?InputItor
>
Type?h_cos_square_accumulate(?????InputItor?h_start,?InputItor?h_end,
????????????????????????????????InputItor?t_start,?InputItor?t_end,
????????????????????????????????Type?h,
????????????????????????????????Type?tau,
????????????????????????????????Type?omega
????????????????????????????)
{
????return?h_f_square_accumulate(
????????????????????????????????????h_start,?h_end,
????????????????????????????????????t_start,?t_end,
????????????????????????????????????h,
????????????????????????????????????tau,
????????????????????????????????????omega,
????????????????????????????????????Cos<Type>()
????????????????????????????????);
}

template
<
????typename?Type,
????typename?InputItor,
????typename?Function
>
static?Type?f_square_accumulate(?????InputItor?t_start,
????????????????????????????????????InputItor?t_end,
????????????????????????????????????Type?tau,
????????????????????????????????????Type?omega,
????????????????????????????????????Function?f
????????????????????????????)
{
????Type?ans?=?Type();
????while(?t_start?!=?t_end?)
????{
????????const?Type?tmp?=?f(?omega?*?(?*t_start++?-?tau?)?);
????????ans?+=?tmp?*?tmp;
????}
????return?ans;
}

template
<
????typename?Type,
????typename?InputItor
>
Type?cos_square_accumulate(?????InputItor?t_start,
????????????????????????????????InputItor?t_end,
????????????????????????????????Type?tau,
????????????????????????????????Type?omega
????????????????????????????)
{
????return?f_square_accumulate(
????????????????????????????????????t_start,
????????????????????????????????????t_end,
????????????????????????????????????tau,
????????????????????????????????????omega,
????????????????????????????????????Cos<Type>()
????????????????????????????????);
}

template
<
????typename?Type,
????typename?InputItor
>
Type?sin_square_accumulate(?????InputItor?t_start,
????????????????????????????????InputItor?t_end,
????????????????????????????????Type?tau,
????????????????????????????????Type?omega
????????????????????????????)
{
????return?f_square_accumulate(
????????????????????????????????????t_start,
????????????????????????????????????t_end,
????????????????????????????????????tau,
????????????????????????????????????omega,
????????????????????????????????????Sin<Type>()
????????????????????????????????);
}

?

實(shí)例

下面這個(gè)例子中，將從文件phi.dat中輸入采樣數(shù)據(jù)，從文件rem.dat中輸入采樣空間間隔，使用Lomb算法處理后，得到的空間頻譜將輸出到文件spectrum.dat中去

?1?#include?"lomb.hpp"
?2?
?3?#include?<iostream>
?4?#include?<fstream>
?5?#include?<vector>
?6?#include?<iterator>
?7?#include?<algorithm>
?8?
?9?using?namespace?std;
10?
11?int?main()
12?{
13?????ifstream?phi_src(?"phi.dat"?);
14?????ifstream?rem_src(?"rem.dat"?);
15?????ofstream?spectrum_dst(?"spectrum.dat"?);
16?
17?????vector<long?double>?h;
18?????vector<long?double>?t;
19?
20?????copy(?istream_iterator<long?double>(phi_src),?istream_iterator<long?double>(),?back_inserter(t)?);
21?????copy(?istream_iterator<long?double>(rem_src),?istream_iterator<long?double>(),?back_inserter(h)?);
22?
23?????Lomb<long?double>*?lomb?=?new?Lomb<long?double>(?h,?t?);
24?
25?????vector<long?double>?v_spectrum?=?lomb?->?spectrum();
26?????copy(?v_spectrum.begin(),?v_spectrum.end(),?ostream_iterator<long?double>(spectrum_dst,?"\n")?);
27?
28?????delete?lomb;
29?????phi_src.close();
30?????rem_src.close();
31?????spectrum_dst.close();
32?
33?????return?0;
34?}
35?

?

備注

Lomb算法的復(fù)雜度是O(n^2)，若采樣數(shù)據(jù)比較長(zhǎng)，可以采用fft方法簡(jiǎn)化復(fù)雜度到O(nlogn)，與大數(shù)乘法中的fft用法一致，此處不再多說(shuō)。




posted on 2009-01-02 21:20 Wang Feng 閱讀(2606) 評(píng)論(3)  編輯 收藏 引用 所屬分類(lèi): Numerical C++

FeedBack:
# re: 非均勻取樣數(shù)據(jù)的功率譜估計(jì)方法
2010-07-09 14:53 | JannieBallard
One knows that men's life seems to be expensive, however some people require money for different things and not every man earns big sums cash. Hence to get fast <a href="http://bestfinance-blog.com/topics/business-loans">business loans</a> and just consolidation loans would be a correct solution.   回復(fù)  更多評(píng)論
  
# re: 非均勻取樣數(shù)據(jù)的功率譜估計(jì)方法
2010-07-18 01:00 | term paper help
There is nothing unnatural If you buy custom essays from the custom research paper service, however this thing would aid you to complete academic results.   回復(fù)  更多評(píng)論
  
# re: 非均勻取樣數(shù)據(jù)的功率譜估計(jì)方法
2010-07-21 11:49 | dissertation writing service
Direct on what you are going to do at the time being and tomorrow, and have loyalty that the extent will advance when it is damn well on tap. Or maybe buy thesis is what you need . Thanks.   回復(fù)  更多評(píng)論
  

<2010年7月>
日一二三四五六
27282930123
45678910
11121314151617
18192021222324
25262728293031
1234567

常用鏈接

留言簿(4)

隨筆分類(lèi)

隨筆檔案

Link List

搜索

  •  

最新評(píng)論

閱讀排行榜

評(píng)論排行榜

青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品

    • <ins id="pjuwb"></ins>
    <blockquote id="pjuwb"><pre id="pjuwb"></pre></blockquote>

    <noscript id="pjuwb"></noscript>
          <sup id="pjuwb"><pre id="pjuwb"></pre></sup>

            <dd id="pjuwb"></dd>
            <abbr id="pjuwb"></abbr>
            
99精品国产高清一区二区|
亚洲激情视频在线观看|
久久精品亚洲一区二区三区浴池
|
99精品视频免费观看|
夜久久久久久|
亚洲精品一二三|
亚洲精品久久久一区二区三区|
久久成人在线|
91久久精品视频|
欧美va亚洲va香蕉在线|
欧美亚洲视频在线看网址|
亚洲国产一区二区三区青草影视|
在线视频亚洲一区|
欧美暴力喷水在线|
免播放器亚洲|
欧美黄污视频|
国产精品日韩在线一区|
亚洲大胆视频|
欧美韩国在线|
欧美国产三区|
欧美在线观看一二区|
国产精品视频区|
国产精品视频不卡|
欧美不卡激情三级在线观看|
韩国视频理论视频久久|
欧美日本韩国一区二区三区|
欧美成人有码|
久久午夜色播影院免费高清|
亚洲国产婷婷综合在线精品|
亚洲美女视频在线免费观看|
亚洲自拍偷拍视频|
尤物yw午夜国产精品视频明星|
国产精品久久久久久超碰
|
欧美搞黄网站|
99国产精品久久久久久久|
能在线观看的日韩av|
亚洲欧美国产精品va在线观看|
亚洲电影毛片|
亚洲伊人网站|
久久精品国产69国产精品亚洲|
午夜精品免费视频|
国产精品av免费在线观看|
欧美日韩免费观看一区|
精品动漫3d一区二区三区|
免费毛片一区二区三区久久久|
国产精品理论片|
在线观看一区二区视频|
久久综合狠狠综合久久综青草|
欧美精品在线视频|
亚洲一区免费网站|
国产精品一区二区久激情瑜伽|
亚洲精品1234|
欧美激情视频一区二区三区不卡|
久久黄色网页|
欧美激情国产精品|
9久草视频在线视频精品|
国产精品高清免费在线观看|
亚洲精品在线视频|
亚洲精选在线观看|
国产欧美一区二区三区在线看蜜臀|
欧美在线一二三四区|
国产精品久久久久一区二区三区|
欧美一区综合|
国产日产欧美a一级在线|
久久综合伊人77777尤物|
亚洲色图自拍|
在线观看三级视频欧美|
国产欧美日韩视频一区二区三区
|
亚洲精品系列|
午夜精品久久久久久久久久久久久|
午夜精品久久99蜜桃的功能介绍|
国产热re99久久6国产精品|
黄色亚洲在线|
欧美成人精品不卡视频在线观看|
欧美高清在线一区二区|
亚洲欧美国产高清|
亚洲美女免费视频|
激情小说亚洲一区|
香蕉av777xxx色综合一区|
中文亚洲欧美|
最新热久久免费视频|
久久久久久穴|
一区免费观看|
一区二区三区国产盗摄|
亚洲视频电影在线|
久久久www|
欧美在线观看视频|
欧美一区在线直播|
性久久久久久久久|
欧美激情精品久久久久久|
亚洲欧美综合网|
国产精品狼人久久影院观看方式|
国产午夜精品视频免费不卡69堂|
免费h精品视频在线播放|
国产精品你懂的|
欧美大片第1页|
久久久久久久高潮|
国产主播一区二区三区四区|
欧美激情中文字幕乱码免费|
欧美日韩国产综合视频在线|
久久亚洲私人国产精品va|
久久久99国产精品免费|
午夜精品久久久|
久久欧美中文字幕|
亚洲欧美综合网|
欧美国产激情|
亚洲视频精选|
久久久高清一区二区三区|
午夜亚洲性色视频|
国内精品久久久|
久久精品在线视频|
狼狼综合久久久久综合网
|
亚洲美女电影在线|
欧美淫片网站|
亚洲午夜激情在线|
韩日视频一区|
亚洲日本成人女熟在线观看|
亚洲午夜电影|
亚洲日韩第九十九页|
乱码第一页成人|
国产精品观看|
亚洲欧美日韩第一区|
国产色视频一区|
国产精品视屏|
麻豆成人精品|
欧美一区亚洲一区|
久久亚洲春色中文字幕|
亚洲免费黄色|
欧美亚洲一区二区在线|
欧美激情女人20p|
亚洲午夜久久久久久久久电影网|
日韩视频免费观看|
亚洲电影自拍|
午夜精品视频网站|
玖玖综合伊人|
亚洲综合精品|
牛牛精品成人免费视频|
亚洲激情视频|
欧美激情国产日韩精品一区18|
亚洲国产精品一区二区三区|
欧美视频二区36p|
亚洲人午夜精品免费|
黄色av日韩|
99精品久久久|
亚洲日本理论电影|
另类亚洲自拍|
午夜激情综合网|
亚洲电影网站|
国产精品乱码|
久久国产精品电影|
久久国产婷婷国产香蕉|
欧美亚洲一级片|
欧美黄污视频|
亚洲一区观看|
亚洲一区二区三区在线观看视频|
亚洲国产精品久久精品怡红院|
亚洲激情在线播放|
国产一区二区欧美|
免费亚洲一区二区|
久久成人18免费网站|
欧美色图首页|
欧美一区二区视频在线|
小黄鸭精品aⅴ导航网站入口|
欧美成人黑人xx视频免费观看|
欧美日韩国产一级片|
亚洲先锋成人|
一区二区三区四区蜜桃|
久久久久久九九九九|
国产精品捆绑调教|
麻豆久久婷婷|
欧美在线播放|
国产精品美女999|
亚洲大片av|
久久久久国产精品午夜一区|
精品91在线|
亚洲精品久久|
国产三级欧美三级|
欧美成人精品激情在线观看|
欧美在线视频导航|
蜜桃久久av|
亚洲视频综合|
亚洲欧美成人网|
国产伦精品一区二区三区高清|
亚洲精品日韩精品|
国产精品伊人日日|
久久精品在线播放|
国产精品免费aⅴ片在线观看|
欧美一区二区三区视频免费播放|
欧美日韩国产成人在线91|
中文日韩欧美|
欧美国产日韩精品|
欧美一区91|
久久三级视频|
欧美成人激情视频免费观看|
亚洲日本aⅴ片在线观看香蕉|
一本大道av伊人久久综合|
国产视频在线一区二区|
亚洲午夜日本在线观看|
亚洲夜间福利|
国产麻豆午夜三级精品|