Loading [MathJax]/jax/element/mml/optable/GeneralPunctuation.js
Advanced Search
Volume 39 Issue 3
Mar.  2017
Turn off MathJax
Article Contents
Shen Ji-zhong, Shao Zhi-long, Jiang Zheng-ke. Design of Low Voltage Low Power Current-Mode CMOS Circuits Based on Parallel Switches[J]. Journal of Electronics & Information Technology, 2004, 26(8): 1325-1331.
Citation: LIU Cheng, ZHAO Zhan, DU Lidong, FANG Zhen. Method of Measuring Two-dimensional Wind Based on Diametrical Pressure Differences Developed by Flow Around Cylinder[J]. Journal of Electronics & Information Technology, 2017, 39(3): 737-742. doi: 10.11999/JEIT160468

Method of Measuring Two-dimensional Wind Based on Diametrical Pressure Differences Developed by Flow Around Cylinder

doi: 10.11999/JEIT160468
Funds:

Special Fund for Public Welfare (GYHY- 201006044), The National Natural Science Foundation of China (51305423)

  • Received Date: 2016-05-09
  • Rev Recd Date: 2016-09-23
  • Publish Date: 2017-03-19
  • The cylindrical two-dimensional wind sensor employing flow meters is not accurate enough for surface wind measurement. By analyzing the data from wind tunnel experiments, a model is proposed to describe the diametrical pressure differences developed by the flow around a circular cylinder. A method is derived from the model for two-dimensional wind measurement by detecting the diametrical pressure differences. When the proposed method is applied to the data from wind tunnel tests in range of 2~40 m/s, the relative wind speed errors and the wind direction errors are no more than (0.2+0.03 V) m/s and 5 respectively. The proposed method is more accurate than the cylindrical two-dimensional wind sensor using flow meters. Without moving parts, the method is immune to mechanical wear and inertia.
  • 弹道导弹已成为现代战争中的主导武器,而导弹飞行中段往往会伴有诱饵,弹头和诱饵等目标微动特性因其质量分布不同而存在明显差异,因此估计目标的微动参数成为导弹防御系统目标识别重要技术[1]。随着高分辨雷达技术和现代信号处理技术的不断发展,空间目标的微动参数估计已由基于RCS序列的目标微动参数估计逐渐发展到基于微多普勒序列[2,3]、高分辨距离像序列的目标参数估计等技术。结合相位测距等精细化处理技术,利用宽带雷达高分辨序列距离像的参数估计,近年来获得了广泛关注[4,5]

    目标微动是在3维空间进行的,关键散射点瞬时位置蕴含丰富的目标运动几何信息,为3维干涉测量实现目标微动参数估计提供基础[69]。本文考虑分布式相参雷达体制多天线接收系统的空间分集[10],利用接收波程相位差对目标瞬时空间位置进行测量,进而获得目标3维运动几何参数估计。相较于高分辨距离像序列的参数估计,宽带3维干涉测量技术可为微动目标分析在宽带测距的基础上提供干涉测角信息,为目标微动几何参数估计提供更多自由度和更高稳健性。文献[11]利用目标微动在宽带系统下回波信号的调制特征差异来区分各散射点的回波信号,通过干涉处理重构微动目标的3维坐标。文献[12]采用对目标距离-慢时间像进行干涉处理的方法重构出弹道目标的真实3维坐标并有效估计出目标的微动参数与结构参数。考虑到实际应用中,目标会不可避免地偏离电轴,即在斜视情况下[13,14],上述文献算法重构目标3维坐标将发生坐标扭曲,对目标几何运动参数估计应用受限。

    针对以上问题,本文提出了一种斜视校正干涉测量微动几何参数估计方法。首先分别对L型分布式相参雷达体制下各天线接收距离-慢时间回波进行干涉处理,通过求解2元2次非线性方程组并利用坐标变换来校正斜视引起的误差,重构出目标真实的3维运动轨迹;接着利用αβ滤波对重构出的3维运动轨迹进行平滑[15],以达到改善3维轨迹估计的波动和阶跃误差;最后利用平滑后的3维运动轨迹优化求解弹道目标微动几何参数。仿真结果表明本文算法可有效重构目标瞬时3维空间位置,验证了算法具有较好的噪声稳定性。该算法可广泛应用于弹道目标的参数估计、目标识别等实际任务中。

    L型分布式相参雷达与弹道目标几何关系如图1所示,(A,X,Y,Z)为雷达坐标系,(O,x,y,z)为参考坐标系,与雷达坐标系平行。收发一体天线A位于雷达坐标系坐标原点,且与接收天线B, C构成沿X轴和Y轴方向相互垂直的干涉基线,基线长L。目标中心O在雷达坐标系下的坐标为(X,Y,Z),其俯仰角及方位角分别是βα,目标任意散射点K在参考坐标系下坐标为(xk,yk,zk)图2为弹道目标微动模型示意图,H为锥体高,r为底面半径,d为质心距锥顶的距离,自旋轴ws和锥旋轴wc的夹角即为进动角θ。锥体弹道目标通常认为只有3个点起作用,分别是图2中的p1, p2, p3。实际由于遮挡效应,目标的等效散射中心只有锥顶散射点p1和靠近雷达视线的锥底散射点p2

    图  1  雷达与弹道目标几何关系
    图  2  弹道目标微动模型示意图

    对目标各散射点回波进行解线频调处理,并消除包络斜置、残余视频项并距离维脉冲压缩后,天线A接收到的散射点K在快时间频域-慢时间的回波信号表示为

    SAd(f,tm)=σkTpsinc{Tp[f+2μcRΔAk(tm)]}exp[j4πλRΔAk(tm)] (1)

    其中,σk为射点K的散射系数,Tp为脉冲重复周期,μ为调频率,tm为慢时间,λ为波长,c为光速。RΔAk(tm)=RAk(tm)RAref(tm), RAk(tm)表示tm时刻散射点K到天线A的距离,本文选择目标中心到各天线的距离做为参考距离,即RAref(tm)=RAO(tm)。同理对天线B和天线C的接收回波信号进行相同处理,可得到散射点K在快时间频域-慢时间的回波信号SBd(f,tm)SCd(f,tm)。对A, B, C 3天线的接收回波信号进行干涉处理,即可得到tm时刻散射点K处的干涉相位差,其中angle()表示取复数相位

    ΔφAB=angle(SAd(f,tm)SBd(f,tm))=2πλ(RΔAkRΔBk) (2)
    ΔφAC=angle(SAd(f,tm)SCd(f,tm))=2πλ(RΔAkRΔCk) (3)

    实际3维干涉测量技术应用中,目标通常会偏离电轴,即斜视测量。当斜视达到一定程度时,两对天线测出的坐标信息不能直接对应散射点真实坐标值,并且散射点到发射天线的斜距也不能近似为其z轴坐标。忽略斜视问题得到的3维坐标必然会存在误差,还会带来相位模糊问题。下面将对斜视情况下干涉3维测量的误差来源进行分析讨论[14],并提出相应方法校正斜视带来的误差。

    以天线对AB接收信号模型为例,根据图1中目标与天线的几何关系,有

    RAkRBk=2L(X+xk)L2RAk+RBk (4)
    RAORBO=2LXL2RAO+RBO (5)

    进一步可推导出

    ΔRAB=RΔAkRΔBk=(RAkRAO)(RBkRBO)=(RAkRBk)(RAORBO)=2Lxk(RAO+RBO)(RAk+RBk)(RAO+RBO)+2LX(RAO+RBORAkRBk)(RAk+RBk)(RAO+RBO)L2(RAO+RBORAkRBk)(RAk+RBk)(RAO+RBO) (6)

    其中,Xtm时刻目标中心在雷达坐标系的横坐标,xktm时刻散射点K在参考坐标系的横坐标。当远场正视情况下,可近似有RAkRBkRAORBOR0,通过近似后可忽略式(6)的后两项。当目标在斜视情况下,由于X较大导致式(6)中的第2项误差不可忽略,第2项和第3项即为斜视引起的误差相位项。现有方法大多通过利用中心相位补偿误差相位项来解决斜视问题[13],然而中心相位易受噪声的影响,进而影响估计准确性。本文通过求解一个2元2次非线性方程组直接获得散射点x维和y维坐标值,具有较高的运算效率。

    根据上文分析,斜视时由于误差相位项不可忽略,下面介绍通过求解一个2元2次非线性方程组校正斜视引起的误差。

    类似式(6),有

    ΔRAC=2Lxk(RAO+RCO)(RAk+RCk)(RAO+RCO)+2LX(RAO+RCORAkRCk)(RAk+RCk)(RAO+RCO)L2(RAO+RCORAkRCk)(RAk+RCk)(RAO+RCO) (7)

    式(6)和式(7)中可做如下近似:RAO+RCORAO+RBORAk+RCkRAk+RBk2R0。根据目标与天线的几何关系,并做泰勒展开有

    f1=RAO+RBORAkRBk=LykR01R0(2xkX+2ykY+2zkZ+x2k+y2k+z2k)f2=RAO+RCORAkRCk=LxkR01R0(2xkX+2ykY+2zkZ+x2k+y2k+z2k)}   (8)

    将式(6)—式(8)代入式(2)和式(3)中并作近似后可得

    ΔφAB=2πLλR0yk+2YL2λR20Lπf1ΔφAC=2πLλR0xk+2XL2λR20Lπf2} (9)

    式(8)中f1f2中的zkA天线接收信号的距离-慢时间维图像中散射点K所在距离单元坐标表示。则式(9)为关于xkyk2元2次非线性方程组,通过穷举优化搜索求得方程组的解,即为斜视时散射点的x轴和y轴坐标。

    斜视时雷达视线与目标纵轴存在较大夹角,散射点到天线的斜距不能直接近似为z轴坐标,在此进行坐标变换即可完成坐标校正。如图3,对参考坐标系(x,y,z)做两次坐标变换后得到雷达视线坐标系(w,u,v), K为目标上一散射点,O为目标中心,βα分别是俯仰角和方位角。

    图  3  坐标系变换示意图

    首先以z轴为旋转轴,将xoy顺时针旋转π/2αwou,将oxyz坐标系变换到owuz坐标系,再以w轴为旋转轴,将zou顺时针旋转βvou,将owuz坐标系变换到owuv坐标系。通过两次坐标变换将参考坐标系转换成雷达视线坐标系,即

    [wuv]=M1M2[xyz] (10)

    其中M1M2分别为两次坐标变换的旋转矩阵,通过求解式(10)可得

    z=(vysinβy)(1/cosβ) (11)

    根据式(11)求得斜视状态下散射点Kz轴坐标zk。由此,即可得到斜视情况下任意时刻第K个散射点在参考坐标系中的真实坐标。对目标上所有散射点进行上述处理后即可重构出微动目标运动轨迹。

    由于噪声和距离分辨率的影响,上述步骤获得的散射点重构3维运动轨迹曲线会出现波动和阶跃现象,与真实运动轨迹偏差过大,若利用其直接进行微动参数估计,会导致较大的参数估计误差。因此需要对重构3维运动轨迹进行平滑滤波,以减少误差带来的不利影响。微动目标散射中心运动模型不易建立,本文用CV模型进行近似。本文目标的角速度和角加速度变化较小,利用αβ滤波对目标的真实运动信息的影响不大。因此利用αβ滤波对重构3维运动轨迹曲线进行处理,αβ滤波的数学原理在此不再赘述。经αβ滤波后获得经平滑处理后的3维运动轨迹[15]

    由上节重构目标3维运动轨迹可以看出,目标的微动几何参数与3维运动轨迹有直接关系。本节在获取目标真实3维坐标曲线的基础上首先估计进动频率,然后结合重构运动轨迹3维图像的几何关系估计进动角和几何参数。

    对于进动目标而言,锥顶散射点的微动变化曲线为正弦形式,频域为单频信息,频点位置由进动频率确定,可直接对目标进动频率进行估计。

    对上一节估计所得的锥顶散射点坐标曲线xk(tm)做傅里叶变换,F=FFT(xk(tm)), F即代表坐标曲线频谱,则目标进动频率估计值可表示为:ˆfc=max(F)

    利用3维干涉测量技术对弹道目标回波进行处理后,可得到在观测时间内与目标真实空间位置一致的各散射点3维运动轨迹图像,利用该图像可以提取目标进动角及几何参数[12]

    根据图2的几何关系,设向量n为锥顶散射点运动轨迹平面的法向量,KP分别为某一时刻锥顶和锥底散射点所在位置,则向量OK为该时刻自旋轴所在向量,为锥体目标母线长,那么进动角\theta 可表示为

    \theta = \arccos \frac{{{n} \cdot {OK}}}{{\left\| {n} \right\| \cdot \left\| {{OK}} \right\|}} (12)

    半锥角\eta 即为母线与自旋轴的夹角,可表示为

    \eta = \arccos \frac{{{KP} \cdot {KO}}}{{\left\| {{KP}} \right\| \cdot \left\| {{KO}} \right\|}} (13)

    根据几何关系,锥体高度H即为母线长在自旋轴上的投影,可表示为

    H = \left\| {{KP}} \right\| \cdot \cos \eta (14)

    根据勾股定理即可得出底面半径r

    r = \sqrt {{{\left\| {{KP}} \right\|}^2} - {H^2}} (15)

    至此完成斜视测量下的宽带雷达3维干涉测量目标微动参数估计。为了更清楚地说明本文算法,具体流程如图4所示。

    图  4  斜视背景下目标微动参数估计流程图

    本实验所估计弹道目标形状如图2,主要目标参数如表1所示,主要雷达参数如表2所示。

    表  1  实验锥体目标参数
    参数数值
    锥体高度H (m)0.96
    质心到锥顶距离d (m)0.64
    底面半径r (m)0.25
    自旋频率fs (Hz)2
    锥旋频率fz (Hz)3
    进动角θ (°)10
    下载: 导出CSV 
    | 显示表格
    表  2  雷达系统的主要参数
    参数数值
    载频f0 (GHz)10
    带宽B (GHz)1
    脉冲宽度tp (μs)10
    脉冲重复周期prf (Hz)100
    驻留时间T (s)1
    基线长L (m)200
    下载: 导出CSV 
    | 显示表格

    初始时刻目标中心在雷达坐标系下的坐标为(2 km, 2 km, 10 km)。为对比本文所提方法重构3维运动轨迹的精确性,在本文设置的实验条件下,使用参考文献[12]中的估计方法,直接通过干涉测量求解各散射点的x轴和y轴坐标,并根据A天线接收回波距离-慢时间像求得各散射点z轴坐标。理论3维坐标如图5所示,文献[12]干涉测量重构3维坐标如图6所示。

    图  5  理论3维坐标
    图  6  文献[12]干涉测量重构3维坐标

    图6可以看出,由于斜视角的存在,利用文献[12]中所提方法直接通过干涉测量技术重构的3维坐标与真实值有较大偏差。在相同的实验条件下,利用本文算法,在干涉测量的基础上通过求解式(9)和式(11)完成斜视校正,得到的3维重构坐标如图7所示。可以看出,图7中的重构坐标更接近理论值。

    图  7  斜视矫正后3维重构坐标

    由于图7中重构的3维坐标存在波动和阶跃现象,不利于后续微动参数估计。故采用\alpha {\tiny{-}} \beta 滤波对重构出的3维坐标曲线进行平滑,结果如图8所示,与图7对比可以看出平滑处理能够明显改善重构坐标曲线中的波动和阶跃现象。根据图9散射点微动轨迹可以看出平滑处理后的图9(c)较未平滑处理的图9(b)更接近理论值。

    图  8  平滑后3维重构坐标
    图  9  散射点微动轨迹

    结合瞬时3维坐标估计值{F_{\rm e}}({t_m}),瞬时3维坐标理论值{F_{\rm r}}({t_m}),给出式(16)定义估计结果的相似度:

    A = \left( {1 - \frac{{\displaystyle\sum\limits_{{t_m}} {\left| {{F_{\rm r}}({t_m}) - {F_{\rm e}}({t_m})} \right|} }}{{\displaystyle\sum\limits_{{t_m}} {\left| {{F_{\rm r}}({t_m})} \right|} }}} \right) \times 100\% (16)

    实验是在信噪比为20 dB条件下进行,所得实验结果相似度记录在表3中。

    表  3  曲线估计实验结果相似度(%)
    锥顶散射点x锥顶散射点y锥顶散射点z锥底散射点x锥底散射点y锥底散射点z
    文献[12]干涉测量68.2660.4288.7262.9364.9182.73
    斜视校正干涉测量88.5690.4196.5990.8491.2487.28
    下载: 导出CSV 
    | 显示表格

    表3图8中可以看出,相较于文献[12]中所提算法,本文算法能够有效校正斜视带来的重构坐标误差,有较高的曲线估计相似度。

    对锥顶散射点的坐标曲线进行傅里叶变换即可求得目标进动频率{f_{\rm c}} = 3\;{\rm Hz}。在所估计出的3维微动轨迹曲线中取任一时刻锥顶散射点所在位置K和锥底散射点所在位置P,通过求解式(12)—式(15)即可得到目标进动角及几何参数估计值。参数估计结果如表4所示。可以看出,由未平滑处理的重构运动轨迹估计目标参数结果与真实值偏差过大。

    表  4  目标参数估计结果
    进动频率{f_c}\left( {\rm Hz} \right)进动角\theta \left( {^ \circ } \right)半锥角\eta \left( {^ \circ } \right)锥体高度H\left( {\rm m} \right)底面半径r\left( {\rm m} \right)
    真实值310.0014.590.960.25
    平滑处理估计值39.6015.250.950.26
    未平滑处理估计值312.4918.170.870.31
    下载: 导出CSV 
    | 显示表格

    上述实验在SNR为20 dB的实验条件下得出,为验证方法的稳定性,设置SNR=–15~20 dB进行噪声稳健性分析,并以式(17)定义估计相对误差

    Q = \frac{{\left| {a - \bar a} \right|}}{a} \times 100\% (17)

    其中,a\bar a分别代表目标参数真实值与估计值。其他条件不变的情况下,在各信噪比下进行50次蒙特卡罗仿真实验,各参数估计相对误差如图10所示。

    图  10  SNR对参数估计影响

    表4可以看出,各参数估计值与实验模型参数设置值都较为接近,从图10可以看出,SNR较低时,由于\alpha {\tiny{-}}\beta 滤波器的输出受噪声影响较大,进而影响了参数估计精度。在SNR高于5 dB条件下,相对误差最大的参数估计精度也在5%以内,可以说明本文算法有较高稳定性。

    本文以L型分布式相参雷达作用在锥体弹道目标为模型,利用斜视3维干涉测量技术处理回波信号,通过求解一个2元2次非线性方程组和坐标变换重构斜视情况下散射点3维运动轨迹,并在此基础上通过滤波和优化求解实现目标微动几何参数估计。实验结果表明该方法通过斜视3维干涉测量可有效重构散射点3维空间位置,并且具有较高的稳健性。但目前本文仅针对单个光滑锥体目标在进动模式下的参数估计,而实际情况中目标的个数、形状和运动方式会更加复杂,重构运动轨迹和参数的估计也会更加困难,后续工作中可对上述问题进行进一步研究。

  • 孙学金, 王晓蕾, 李浩, 等. 大气探测学[M]. 北京: 气象出版社, 2009: 172-247.
    SUN Xuejin, WANG Xiaolei, LI Hao, et al. Atmospheric Observation[M]. Beijing: China Meteorological Press, 2009: 172247.
    BRUSCHI P, DEI M, and PIOTTO M. A low-power 2-D wind sensor based on integrated flow meters[J]. IEEE Sensors Journal, 2009, 9(12): 16881696. doi: 10.1109/JSEN.2009. 2030652.
    ZHU Yanqing, CHEN Bei, QIN Ming, et al. 2-D micromachined thermal wind sensorsA review[J]. IEEE Internet of Things Journal, 2014, 1(3): 216232. doi: 10.1109/ JIOT.2014.2319296.
    VERESHCHAGINA E, TIGGELAAR R M, SANDERS R G P, et al. Low power micro-calorimetric sensors for analysis of gaseous samples[J]. Sensors and Actuators B: Chemical, 2015, 206: 772787. doi: 10.1016/j.snb.2014.08.077.
    LIU Haobing, LIN Nay, PAN Shanshan, et al. High sensitivity, miniature, full 2-D anemometer based on MEMS hot-film sensors[J]. IEEE Sensors Journal, 2013, 13(5): 19141920. doi: 10.1109/JSEN.2012.2236014.
    KOWALSKI L, JIMENEZ V, DOMINGUEZ-PUMAR M, et al. Low pressure spherical thermal anemometer for space missions[C]. 2013 IEEE Sensors, Baltimore, Maryland, USA, 2013: 14. doi: 10.1109/ICSENS.2013.6688493.
    KOWALSKI L, ATIENZA M T, GORRETA S, et al. Spherical wind sensor for the atmosphere of Mars[J]. IEEE Sensors Journal, 2016, 16(7): 18871897. doi: 10.1109/JSEN. 2015.2509168.
    DONG Ziqiang, CHEN Jingjing, QIN Yukun, et al. Fabrication of a micromachined two-dimensional wind sensor by Au-Au wafer bonding technology[J]. Journal of Microelectromechanical Systems, 2012, 21(2): 467475. doi: 10.1109/JMEMS.2011.2179014.
    ZHU Yanqing, CHEN Bei, GAO Di, et al. A robust and low-power 2-D thermal wind sensor based on a glass-in-silicon reflow process[J]. Microsystem Technologies, 2015, 22(1): 151162. doi: 10.1007/s00542-015-2423-9.
    GEORGIOU D P and KYPROS F M. Fabrication and calibration of a sub-miniature 5-hole probe with embedded pressure sensors for use in extremely confined and complex flow areas in turbomachinery research facilities[J]. Flow Measurement and Instrumentation, 2014, 39: 5463. doi: 10.1016/j.flowmeasinst.2014.07.005.
    PIOTTO M, DEI M, PENNELLI G, et al. A miniaturized 2D solid state anemometer based on thermal flow sensors[J]. Procedia Chemistry, 2009, 1(1): 14631466. doi: 10.1016/j. proche.2009.07.365.
    PIOTTO M, PENNELLI G, and BRUSCHI P. Fabrication and characterization of a directional anemometer based on a single chip MEMS flow sensor[J]. Microelectronic Engineering, 2011, 88(8): 22142217. doi: 10.1016/j.mee.2010.11.009.
    张兆顺, 崔桂香. 流体力学[M]. 北京: 清华大学出版社, 1999: 167168.
    ZHANG Zhaoshun and CUI Guixiang. Hydromechanics[M]. Beijing: Tsinghua University Press, 1999: 167168.
    SCHLICHTING H. Boundary-Layer Theory[M]. 7th Edition, New York: McGraw-Hill, 1979: 1223.
    ZURELL Cory. Aerodynamics of a circular cylinder inclined to airflow and wind-induced vibrations of dry, inclined cables at high wind speeds[D]. [Ph.D. dissertation], University of Ottawa, 2004: 2225.
    CARE I and ARENAS M. On the impact of anemometer size on the velocity field in a closed wind tunnel[J]. Flow Measurement and Instrumentation, 2015, 44: 210. doi: 10.1016/j.flowmeasinst.2014.11.007.
    CHEN Wenli, LI Hui, and HU Hui. An experimental study on a suction flow control method to reduce the unsteadiness of the wind loads acting on a circular cylinder[J]. Experiments in Fluids, 2014, 55(4): 14. doi: 10.1007/s00348-014-1707-7.
    LYSENKO D A, ERTESVG I S, and RIAN K E. Large- eddy simulation of the flow over a circular cylinder at Reynolds number 2104[J]. Flow, Turbulence and Combustion, 2014, 92(3): 673698. doi: 10.1007/s10494-013- 9509-1.
    GIACOMO P. Equation for the determination of the density of moist air (1981)[J]. Metrologia, 1982, 18(3): 3340. doi: 10.1088/0026-1394/18/3/011.
    LIU Cheng, DU Lidong, and ZHAO Zhan. A directional cylindrical anemometer with four sets of differential pressure sensors[J]. Review of Scientific Instruments, 2016, 87(3): 035105. doi: 10.1063/1.4943222.
  • Cited by

    Periodical cited type(2)

    1. 冯存前,李江,黄大荣,胡晓伟,韩立珣. 弹道中段不同平动多目标的平动参数估计方法. 电子与信息学报. 2021(03): 564-571 . 本站查看
    2. 王红敏,宁生科,董涛. 基于红外触发技术的外弹道目标微动参数估计研究. 兵器装备工程学报. 2021(07): 92-97 .

    Other cited types(5)

  • 加载中

Catalog

    通讯作者: 陈斌, bchen63@163.com
    • 1. 

      沈阳化工大学材料科学与工程学院 沈阳 110142

    1. 本站搜索
    2. 百度学术搜索
    3. 万方数据库搜索
    4. CNKI搜索

    Article Metrics

    Article views (1147) PDF downloads(342) Cited by()
    Proportional views
    Related

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return