CUDA技术 CUDA框架

昭棠笔记 2023-01-25

中华民谣孙浩-打碟机教程

2022年5月1日发

(作者:火影忍者主题曲邂逅)

Journal of Computer Applications 

ISSN 1001—9081 2012 09.01 

计算机应用,2012,32(9):2466—2469 CODEN JYIIDU http://www.joca.cn 

文章编号:1001—9081(2012)09—2466—04 doi:10.3724/SP.J.1087.2012.02466 

基于CUDA的大规模体行为实时仿真并行实现及优化 

贺毅辉blog营销。,叶晨tieba baidu com,刘志忠www baidu com ,彭 伟 

(解放军理工大学指挥自动化学院泰安网站建设,南京210007) 

( 通信作者hyh.nj@163.eom) 

摘要:体仿真中个体从环境中查相关对象时会导致较高的时间复杂度。要使大规模体能够实时仿真seo每日一贴,必 

须降低模型运算的时间复杂度或者提高计算平台的能力。通过对Biods模型为典型案例进行研究,提出一种基于统一 

计算架构(CUDA)的大规模体行为实时仿真并行实现及优化的方法盛大推广。实现中将个体与GPU逻辑线程一一对应推一把论坛,通过 

将仿真环境离散化来提高相关个体查的效率,通过并行化基数排序法将个体信息组织成具有空间局部性的数组,提高 

图形处理器(GPU)内存带宽的利用率博客发。通过实验验证了该方法将仿真个体的数量提升到CPU方法的约7.3倍。 

关键词:大规模体行为;统一计算架构;并行计算;实时仿真 

中图分类号:TP391.9 文献标志码:A 

Parallel simulation and optimization of CUDA--based real--time huge crowd behavior 

HE Yi.hui YE Chen LIU Zhi.zhongcms是什么,PENG Wei 

(Institute of Command Automation,PLA University ofScience and Technology,?Canjing Jiangsu 210007,China) 

Abstract:That the individuals search relevant objects from the environment may cause high time complexity during the 

crowd simulation.If the crowd should be simulated in real—time.the time complexity of the model needs reducing and the 

computing capability of the simulation platform needs enhancing.In this paper,the Biods model was studied as a typical case 

and a solution of how to parallelize and optimize the real-time huge crowd simulation based on Compute Uniifed Device 

Architecture(CUDA)was presented.Each individual was correspondent to a l0舀cal Graphic Processing Unit(GPU)thread. 

By discretizing simulation environment,the efifciency of searching the relevant individuals was improved.The individual 

information was organized into an array with the spatila locality by parallel radix sort in order to improve the utilization of the 

GPU memo ̄bandwidth.The experiment verifies the solution presented here has improved number of simulation individuals up 

to about 7.3 times as CPU solution. 

Key words:huge crowd behavior;Compute Unified Device Architecture(CUDA);parlalel computing;real—time 

simulation 

0 引言 

理器平台仿真10000个个体的逃散场景达到了45 fps的帧速 

度;文献[6]通过CPU与图形处理器(Graphic Processing Unit, 

体行为是人类及其他生物的一种典型社会性行为,在 

GPU)的混合计算,能够以60 fps的速度仿真1 600个Biods个 

虚拟现实技术中有重要的研究价值和广泛的应用空间网络营销好做吗。研究 

体;文献[7—9]通过完全使用GPU进行体运动计算,避免 

体行为,一个侧重点是根据社会学等理论建立体行为模 

了CPU与GPU间频繁的数据交互,进一步提高了仿真的效 

型:如Biods模型以及考虑社会学意义的人际关系的人类 

率门户网站排名。 

体行为模型等…百度快照。其他类似的模型通过引入心理学效果、社 

大规模体行为仿真主要针对基于人工生命的什么是响应式网站,采用自 

会力等来建立。这些模型重点在于刻画体行为危机公关案例,力求 

底向上的综合方法荆州网,用于研究高层的自组织过程和行为涌现营销文案策划。 

行为的逻辑正确性,为该领域奠定了基础。用计算机仿真,特 

别是实时可视化渲染仿真体行为是该领域研究的另一个侧 

以Biods模型为例,该模型通过对齐、聚集、分离三条规则软推网,使 

重点当阳网。上述模型在建立时又名潍坊站长网,并未考虑实时仿真的需求。以 

个体在达到一定的数量时涌现出聚行为。该模型具备以下 

Biods模型为例,个体对周围环境的感知时都要通过查询所有 

特点网站seo优化软件,即所有的个体遵循相同的行为方式百度电话,区别在于速度、位 

其他个体的信息来辨别出与其相关的个体百度推广托管,时间复杂度为 

置等个体信息自然排名,符合单指令多数据流(Simple Instruction 

O(n )。如果体达到一定规模后app营销,仿真的每一步计算都将导 

Multiple Data,SIMD)特征烟台百度公司,满足用GPU并行计算的条件。继 

致较大的时间耗费,无法满足实时仿真的要求。 

Biods模型之后,出现了基于该模型的各种改进模型,如在 

为使大规模体行为仿真满足实时化的要求手机网站排名优化软件,主要采用 

Biods模型中加入避碰行为、路径规划行为、分组行为等泸州网站建设。这 

的方法有两类。一类是对模型进行优化,降低模型的时间复 

些模型从本质上说南宁seo,虽然个体有了更多的行为类别扬州网站建设,但同类个 

杂度,如文献[2—4]等开发出的一些利用空间特征的算法百度账户优化, 

体仍按照相同的行为运动,SIMD的特征依然存在。所以本文 

通过预处理存储辅助信息以提高个体的辨识效率永康网站建设。另一类方 

以研究Biods模型的并行化实时仿真可以作为基于人工生命 

法是提高仿真平台的计算能力网站制作计划,如文献[5]使用分布式多处 的大规模体行为模型并行化的典型案例。 

收稿日期:2012—02—24;修回日期:2012—05-29。 基金项目:江苏省自然科学基金资助项目(BK2010130原创代写,BK2011120)。 

作者简介:贺毅辉(1973一)google搜索解析,男,河北辛集人写好文章,教授,主要研究方向:作战仿真、决策支持;叶晨(1984一)百度网址首页,男百度推广开户,江苏无锡人google排名优化,工程师,硕士研 

究生为什么打开网页很慢,主要研究方向:系统集成与优化; 刘志忠(1980一),男搜索引擎营销策略,江西吉安人,讲师,博士百度关键词价格,主要研究方向:面向服务计算、人工智能; 彭伟 

(1982一)武汉网站优化,男苏州百度公司,湖南常德人,讲师,硕士,主要研究方向:作战仿真、系统集成危机公关原则。 

第9期 贺毅辉等:基于CUDA的大规模体行为实时仿真并行实现及优化 2467 

1 支持CUDA图形设备的体系结构 

统一计算架构(Compute Uniifed Device Architecture刷百度关键词排名, 

CUDA)H淘宝搜索优化。’是英伟达公司开发的具有并行编程模型和指令集 

架构的通用计算架构。支持CUDA的图形设备Device一般 

分成两部分:由采用SIMD结构的流多处理器(Streaming 

Muhiprocessor)构成的流处理器阵列(Streaming Processor 

Array)以及存储器系统(Memory System)。两部分由GPU片 

上的互联网络(Interconnection Network)实现通信。支持 

CUDA的图形设备体系结构如图1所示。 

4带来的另一个好处是在存储个体安全圆覆盖的单元ID与 

单元包含的个体安全圆ID时可以使用对齐的uint4数据结 

构,使半线程束(hatf-warp)在存取全局内存时可以将16次访 

问合并成一次,避免了没有对齐的访存被编译器拆分成16次 

串行访存引起的性能下降网站运营数据分析。 

2.2基于并行基数排序算法的数据组织 

在Biods系统中的个体对其周边环境的搜索时表现为对 

一片d维空间中具有空间局部性 的连续单元区域中的数 

据进行读取。GPU通过纹理缓存,缓存与被访存单元具有空 

间局部关联单元的纹理内存地址,加速对这种类型的数据访 

图1 支持CUDA图形设备体系结构示意图 

在基于CUDA开发的并行程序中,GPU线程是基本的执 

行单元360搜索引擎怎么样。Device存储器系统分为片内和片外两大类:寄存器 

和共享存储器为片内存储器seo帝国,被流多处理器所拥有;全局内 

存、本地内存、纹理内存、常量内存为片外存储器什么叫营销。纹理内存 

和常量内存在片内有对应的缓存,以提高对片外显存的访问 

速度安卓系统优化。全局内存没有缓存,每次对其访问将导致一次实际的 

显存访问深圳网站推广公司,访问延迟依据显存的速度,GDDR3延迟一般为 

400~600U ̄钟周期。 

2 以Biods模型为典型的并行仿真及优化 

对于Biods模型,一种最直观的并行化方法就是将个体 

与逻辑GPU线程一一对应如何设计网站,仿真数据存放在显卡的全局内存 

上,个体位置映射到顶点缓冲区上,由OpenGL直接可视化渲 

染黑龙江网站建设。然而,采用上述方法的并行化的Biods运行效率很低哪个网站建设好,主 

要原因在于:1)相关个体的搜索量很高草根站长工具,且其中存在大量的 

无效的搜索;2)GPU访存中存在引起接合不良的访存方式珠海网站推广。 

2.1仿真空间划分方法 

将仿真空间离散化能有效降低个体对其他相关个体的搜 

索量。然而相关个体搜索量的降低实际是由增加对离散化后 

的子空间维护成本来换取的,维护量和子空间的数量成正比。 

所以在构建离散环境时必须考虑两者间的平衡,以取得更佳 

的性能。使用空间划分后Biods算法的平均时间复杂度为 

O(n log n) 营销公关。 

本文将仿真设定在二维空间中微指数,Biods系统中个体被视为 

具有安全半径r的圆形粒子南宁搜索引擎优化,仿真空间划分为网格,单元大小 

为安全圆的外接正方形单页优化。这种划分方法的灵感来源于粒子系 

统 物理仿真中的宽阶段碰撞检测 12]小江seotaobao,带来的一个好处是 

在个体相互作用时不会出现个体重叠的不合理现象金山推广员。以安全 

圆的中心界定个体所处的单元序号,在这种情况下英文网站建设,一个个体 

的安全圆最多可以出现在4个(2 新闻危机处理,d是空间划分的维度)” 

单元中,同理四平网,一个单元最多包含4个安全圆百度竞价排名费用,如图2所示。 

图2安全圆覆盖单元数与单元包含安全圆数上限 

安全圆覆盖的单元数与单元包含安全圆数上限确定且为 

存柳州360优化。基于上述原因网站如何赚钱,需要构造个体信息的空间局部性特点。 

构造通过线性时间复杂度 的基数排序算法将个体信息组 

织成按所在单元ID的自然数序列排列的数组来完成。 

在本案例中,基数排序所使用的关键字(key)是单元ID。 

对于有 个二进制位的单元ID,按照每组 个二进制位分组, 

其中L<B,组的数量和排序趟数一致aizhan。 

2.2.1基数统计 

在基数统计阶段需要对每个基数出现的次数进行统计seo178, 

则至少需要radicesNum=2 个32位计数器。将计数器放到 

共享内存上,可以提升内存带宽。每个线程都需要自己的计 

数器来保证计数值不被其他线程修改seo实战密码,共需计数器数为: 

blockNum×threadNumPerBlock×radicesNum 

共享内存有大小限制,如Geforee GTX280上每个流多处 

理器仅有16 KB的该类存储,所以必须根据GPU的硬件特 

性,确定 的大小和每个线程块中线程的数量 

(threadNumPerBlock)。 

减少趟数且不导致共享内存溢出的办法就是共享线程间 

的计数器。将每个线程块中的线程分成线程组,每个组R个线 

程文昌网,则会产生的组数为: 

groupNumPerBlock=threadNumPerBlock/R 

每个组内的线程共享相同的radicesNum个计数器smo优化,通过 

串行化修改组内的基数计数器来避免同组线程间的读写冲 

突。对组内的线程进行串行测试时北京百度优化,避免了线程并行引起的冲 

突绝地求生优化,要避免使用if语句打断程序流水。算法伪代码如下: 

d=((threadlD mod R)==0); 

radixCounter+=d¥appropriateIncrement: 

d=((threadlDmodR)==1); 

radixCounter+=d appropriateIncrement: 

d=((threadlD mod R)==R一1); 

radixCounter+=d appropriatelncrement: 

2.2.2基数求和 

在输出排序结果之前,需要计算输出位置的偏移量,通过 

以下步骤来计算:1)对每个基数网站被挂马检测工具,求它的每个基数计数器的 

前缀和;2)对每个基数网站推广软文案例,求它的基数计数器的总和;3)求所有 

radicesNum个基数计数器的前缀和。 

在第3)步中得到的前缀和给出了输出个体位置信息数 

组中具有相同基数的每一部分的偏移量长沙百度推广。第1)步中求出的前 

缀和给出了这些部分中每一个的偏移量,然后每个线程块会 

将其分配到的个体信息根据这些偏移写入到新的个体信息数 

组中网络营销网络推广。在进行第1)步时,每个线程块处理的基数个数为: 

radicesNumPerBlock=radicesNum/blockNum 

每个基数的基数计数器和线程组的基数计数器的数量一 

致,即每个基数拥有的基数计数器数量为: 

groupsNum:blockNum groupsNumPerBlock 

所以每个块都将其32位的groupsNum存放在共享内存 

计算机应用 第32卷 

上,然后它在分配给它的radicesNumPerBlock个基数上迭代。 持25帧的上限个体数量约为7000,与第一个程序相比仿真 

对线程块中的每个线程而言,都要执行如下操作:1)从全局 

存储器中读取基数的groupsNum个计数器到共享内存中; 

2)求groupsNum个计数器的前缀和;3)全局内存中的 

groupsNum个计数器的值被前面所求的前缀和覆盖;4)在最 

后一个计数器被覆盖之前佛山关键词优化,通过把它的值和上一个前缀和相 

加得出groupsNum个计数器值的总和南通网站优化。 

2.2.3 重排序 

个体数量上限提高到原来的约3.3倍。使用本文提出的并行 

化及优化方法开发的优化版GPU仿真程序寇玉磊,因为设定了仿真 

的帧率上限南昌网站优化,在个体数量小于7 168时收录批量查询,基本稳定在100帧左 

右长沙做网站价格,个体数继续增加后,因为优化后算法的平均时间复杂度为 

0(/Z log n)旺道网络营销软件,帧率以接近直线的方式下降怎么建立自己的博客。当个体数量为 

15 360时,仿真的帧率仍可以达到26帧白帽技术,在满足实时仿真帧率 

要求的前提下,个体数量是无优化GPU仿真的约2.2倍,是 

纯CPU仿真的约7.3倍。 

重新排序在每一趟的最后一个阶段。每个块必须计算 

radicesNum个基数计数器的前缀和以获得输出数组的偏 

移量外贸seo。从全局内存中载入第二阶段得出的blockNum个前缀和 

到共享内存中,对每个前缀和的最后一个值再执行并行的前 

缀求和,并使用它来求出每个前缀和的每个值的合适的偏移 

量。每个线程块都重复地执行这里的并行前缀求和。然后统计分析工具,每 

个块从全局内存中载入groupNumPerBlock×radicesNum个基 

I_l』一 l+无优化版GPU{ ̄真程序 

斗 …_l…_j…L 一l十优化版GPU仿真程序 

t1… 东港网,_j—奋一纯cPu版仿真程序 

… 一 

数计数器到共享内存三亚网,并使用前面计算出的偏移量来放 

置它们。最后,每个线程块都按照和第一阶段相同的顺序diemin,从 

全局内存读取其在第一阶段读取的输入个体信息数组的相同 

部分,每个个体信息都要执行下面的操作:1)读取与线程组 

以及个体信息的基数对应的基数计数器的值;2)按照基数计 

斗i

i  ;l l I\ 

- 

l t_ 

一} iIJ { I一一u 

卜…一 一一}¨…~{.一}…u 

 fI I 

l l{ …I¨…r 

…一

}、} _l…__}\心 . . i l .- . 

仿真Biods数 

一F …… r j…f_…~  {  『_l…{r一 … 

● 

数器指示的偏移量网站建设基础知识,把个体信息写到输出个体信息数组中; 

3)增加基数计数器的值。对于每个半线程束的线程,2)和3) 

的写操作与递增操作必须分别被串行化,以避免同一半线程 

图3各版本Biods仿真程序性能比较 

4 结语 

本文提出了一种针对Biods系统的GPU并行化仿真及优 

束的各线程间存在任何共享内存和全局内存风险宁波seo,并可使用 

同一基数来保持个体信息的顺序。 

在每一步仿真中,经过上述B/L趟的操作,最终可以将个 

体信息组织成按照个体所在单元ID的次序排列的数组。这个 

排序过的个体信息数组将映射到纹理内存上,个体也将以该 

化的方法搜索引擎市场,即将个体与GPU逻辑线程一一对应seo实战密码 pdf,通过将仿真 

空间离散化,降低个体对相关个体的搜索量刷百度权重,通过并行基数排 

序算法,将个体信息数组组织成具有空间局部性的数组网站规划,提高 

内存带宽利用率。经过实验验证,本文方法确实对算法本身 

次序映射到逻辑线程上,为个体对其周边环境搜索的计算提 

供加速上海企业推广。上述过程中并行求前缀和的运算通过引入文献[16一 

l7]设计的算法实现。 

有一定的优化陶然士论文,并通过使用通用GPU计算提升了仿真平台的 

计算能力。 

本文方法可以作为基于人工生命的大规模体仿真并行 

化及优化的一个典型。当模型变得更加复杂,特别是个体依 

行为方式区分存在大量的类别时南京seo,如果将个体与GPU逻辑线 

程一一对应软件优化,则需要在核函数中根据个体类别选择对应的行 

为方式,这样会导致大量的分支语句存在于核函数中。由于 

3 实验验证 

使用上述方法并行化的Biods程序(优化版GPU仿真程 

序)采用Visual Studio 2008+NVIDIA CUDA Toolkit 4.0开发, 

由OpenGL可视化渲染网站推广培训。另外有两个对比程序,可视化渲染 

均由OpenGL进行:第一个由CPU计算个体位置和速度百度快照不更新,仿真 

空间采用连续2D空间(纯CPU版仿真程序);第二个仅将第 

个对比程序的CPU计算工作交给GPU,个体信息的存储放 一

GPU的流控制能力较弱,对分支语句不能如CPU那样进行预 

判优化,GPU程序在运行到分支语句时友情链接交换,大量线程会被挂起, 

造成性能瓶颈seo每日一贴。针对上述问题,如果模型中个体的行为方式 

和区域有关,一种解决方案就是将仿真环境栅格化后将栅格 

到全局内存上(无优化版GPU仿真程序)恒牛网。运行这些Biods 

程序的计算机平台主要配置为:2.83 GHz主频英特尔Q9550 

处理器,1 333 MHz主频4 GB双通道DDR3内存普宁网,NVIDIA 

Geforce GTX280显卡如何提高百度权重。实验结果如图3所示百合seo培训。 

与GPU逻辑线程一一对应竞价单页,仿真并行地更新这些栅格,每个 

栅格中的个体由GPU线程串行更新市场推广策划方案。如何针对这种并行方 

式进行优化将是下一步工作的重点上海seo搜索优化。 

参考文献: 

[1] MUSSE S R陆丰网站,THALMANN D.A hierarchical model for real time 

simulation of virtual human crowds[J】.IEEE Transactions on Visu— 

liazation and Computer Graphics自然排名,2001徐州seo顾问,7(2):152—164. 

图3由Matlab分别对1 024 X (n∈[1,16])个数量的 

Biods仿真5 min统计得到平均帧率绘制得到。从图3可看 

出:纯CPU版仿真程序性能更低,仿真的帧率以二次曲线的 

方式下降,在2048个个体时帧率为27,当个体数继续增加时 

无法达到实时可视化渲染要求的25帧宿州网。造成其低性能的主 

要原因是CPU的计算能力相对较弱,个体对周围环境感知算 

[21 TESCHNER M,HEIDELBERGER B.Optimized spatial hashing for 

collision detection of deformable objects【C】//Proceedings of Vi— 

sion百度刷下拉,Modeling and Visualization.Manchen,Germany:[S.n]seo优化, 

2003:47—54. 

法的时间复杂度为O(n )友情链。无优化版GPU仿真程序较第一个 

已经有了明显的性能提升,除得益于GPU对计算密集型作业 

处理能力的优势外英语网站建设,CUDA与OpenGL的互操作避免了第一个 

程序中可视化渲染中从内存到显存数据复制的瓶颈,但由于 

[3] SHAO W,TERZOPOULOS D.Autonomous pedestrians[c]// 

SCA’05:Proceedings of the 2005 ACM SIGGRAPH/Eurographics 

Symposium On Computer animation.New York:ACM Pressseo优化狼雨,2005: 

19—28. 该程序中个体对周围环境感知算法的时间复杂度仍为 

O(/Z ),所以仿真帧率仍是以二次曲线的方式下降,当个体数 

为6 144时帧率已经降到了30帧,从曲线上判断保证帧率维 

[41 REYNOLDS C.Big fast crowds on PS3[C】//Sandbox’06:Pro— 

ceedings of the 2006 ACM SIGGRAPH Symposium on Videogames. 

第9期 贺毅辉等:基于CUDA的大规模体行为实时仿真并行实现及优化 2469 

(上接第2465页) 

加龙岩做网站公司。对于较大的数据集算法的效率更明显细雨算法,其加速比曲线接 ID3算法并行化是一项重要的课题销售手段。本文在ID3算法的基础 

近线性增加通州网站建设公司,这说明基于MapReduee的H—S—Decision算法可 

上提出了基于属性相似度的决策树算法,并将其应用 

用来解决海量数据的挖掘问题28推广。 

MapReduce平台构建了一种并行性能比较好的决策树算法 

5 

[至三=三至至至 三互] H

—Decision。实验结果表明H—S—Decision算法具有较好的 

分类正确率、较低的时间复杂度和较好并行性能。 

4 

丑 

参考文献: 

3 

景 

[1] 魏红宁.基于SPRINT方法的并行决策树分类研究【J】.计算机 

2 

应用,2005深圳百度竞价排名,25(1):39—41. 

l 

[2] 朱敏衡阳网,万剑怡视频优化,王明文.基于MR的并行决策树分类算法的设计 

节点数 

与实现[J].广西师范大学学报:自然科学版黑龙江网站建设,2011,29(1):82— 

图3 H—S—Decision算法加速比曲线 

84. 

[3] 杨宸铸.基于HADOOP的数据挖掘研究【D】.重庆:重庆大学, 

3)可扩展性汕头网站建设。 

2010:45—46. 

为了反映并行算法集的利用率情况,引入算法效率的 

[4] 李应安.基于MapReduce的聚类算法的并行化研究【D】.广州. 

概念。并行算法的效率是算法的加速比与节点数的比值购买域名和空间,即 

中山大学app推广方法,2010:50—51. 

E=S/Nsns是什么意思,其中S为加速比,Ⅳ为节点数。从图4可看出,算法 

[5] 程苗baidu comn,陈华平.基于Hadoop的Web日志挖掘[J】.计算机工程莱芜网站建设, 

效率整体来说呈下降趋势网站运营方案。从图4可看出D 曲线相对D,、 

2011,37(11):37—38. 

D 平稳一些互联网营销公司,这说明当数据集规模越大,H—S—Decision算法 

[6】 雷万云.云计算技术、平台及应用案例[M].北京:清华大学出版 

的效率曲线越平稳网站建设中图片,所以该算法对海量数据具有较好的可扩 

社seo优化工具,2011:22—224. 

展性深圳公牛通讯。 

[7]LIU YANGurl什么意思,LI MAOZHEN网站建设知识,ALHAM N K.HSIM:A MapReduce 

[至 至至至 三互] simulator in enabling cloud computing【EB/OL】.[2012-01—20]. 

http://www.sciencedirect.com/science/artiele/pii/SO167739X1 10 

00884. 

[8] ALHAM N K百度关键词seo,LI MAOZHEN,LIU YANG.A MapReduce-based 

distributed SVM algorithm of automatic image annotation【J].Com— 

puters and Mathematics with Applications,201 1网上推广平台,62(7):2801— 

2811. 

节点数 【9】 毛国君江苏seo,段立娟青海网站建设,王实阳江网站制作公司,等.数据挖掘原理与算法【M】.北京:清华 

图4 H—S—Decision算法效率曲线 

大学出版社济宁网站建设,2005:117—125. 

4 结语 [1O】韩松来,张辉百度举报中心,周华平.基于关联度函数的决策树分类算法【J]. 

计算机应用,2005网站统计分析,25(11):2655—2657. 

决策树算法的并行性是当前解决海量数据分类问题的研 

[11]陆秋兔子优化大师,程小辉.基于属性相似度的决策树算法[J].计算机工程, 

究重点seo优化的网站,而ID3算法在决策树算法中占着重要地位,如何将 2009天津搜索引擎优化,35(6):82—84. 

马来风光-经典动画片歌曲