毕业设计（论文）-计算机仿真中相似性原理的应用.doc

安徽三联学院毕业论文本科毕业论文（设计、创作）题目： 计算机仿真中相似性原理的应用 学生姓名： 学号： 所在系院： 信息与通信技术系 专业： 电子信息工程 入学时间： 2010 年 9 月导师姓名： 职称/学位： 副教授/硕士/讲师/硕士 导师所在单位： 完成时间： 2014 年 5 月安徽三联学院教务 制 计算机仿真中相似性原理的应用摘 要:计算机仿真技术的有关内容，所针对的系统主要是工程中的；仿真就是用模型（物理模型或数学模型）代替实际系；遵循的基本原则是相似原理，即几何相似、环境相似性。本课程用到的计算机系统仿真中相似性原理的应用技术，仿真是指通过系统模型的试验去研究一个已经存在的、或者是正在研究设计中的系统计算机仿真技术，是一门利用计算机软件模拟实际环境；它是以数学理论为基础，以计算机和各种物理为设施；用仿真理论来研究系统。本课程主要通过对建模方法与原理、仿真算法的学习通过构建程序得出数据表格进行分析；得出来实际生活中和物理现象中完全不同的系统却能演变出相同的规律，体现了仿真中相似性原理。关键词：仿真； 相似性； 系统； 建模The computer simulation of similarity theory in the applicationAbstract: The contents of the computer simulation technology, the system is mainly aimed at engineering; The simulation is to use model (physical model and mathematical model) instead of the actual system; Is similar to follow the basic principles of principle of geometric similarity, similarity of environment. Similarity principle using computer simulation in the course of application technology, the simulation is to point to by system model test to research an already existing, or is working on the computer simulation technology in the design of system, is a use of computer software to simulate the actual environment; It is based on the mathematical theory, with computer and various physical facilities; To study the system in the simulation theory. This course mainly through the modeling method and the principle, simulation algorithm of learning through a builder data table is analyzed; Come out in the real life and physical phenomena completely different system can evolve the same rule, reflects the similarity principle in the simulation.Keywords: simulation; similarity; system; modeling2目录第一章 计算机仿真1第1节 引言1第2节 计算机仿真技术的应用与影响1第2章 计算机仿真的方法2第1节 龙格库塔法2第2节 欧拉法3第3节 计算机仿真方法的比较与选择6第3章 龙格库塔法相似性仿真中的运用7第1节 建模：一阶零状态响应方程7第2节 基于龙格库塔法的程序设计8第3节 VC+6.0程序运行10第4节 实际销售产品中相似性仿真原理的应用17总结及心得体会18致谢19参考文献20I安徽三联学院毕业论文第一章 计算机仿真 第1节 引言 仿真技术涉及的发展和计算机技术的发展。在一般情况下，之前的物理模型实验的外观基于计算机的称为“模拟” ，我是连接到一般等相关领域。在未来，计算机，因为在特定的数字计算机，和计算机模拟，使巨大的存储容量和计算能力在一个较高的速度来计算的复杂计算机仿真技术的值是一个重要的问题这样一来，开发得到了蓬勃发展。随着仿真应用的扩展，电脑模拟外延也扩大。的物理模拟，声景观，操作，即，从核子的各个组成部分的角度，系统，飞机，船舶和航海模拟器培训，监控系统，并也有许多仍然使用：现代化的人喜欢实时计算机数学模拟的训练模拟器。广泛。这些模拟都包括在计算机模拟的范围之内。早在明朝的朱棣在南京做了皇帝以后派亲信大臣到北京盖角楼，大臣领了皇帝的谕旨后，很发愁，不知如何盖这九梁十八柱、七十二条脊的角楼?并说：“如果盖不成，皇帝自然要杀我的头，可是在没杀我的头之前，我就先把你们的头砍了。”工头们对这样的工程都没把握。一天，一位木匠，上大街去乘凉，走看走看，听见老远传来一片蝈蝈的叫声和哟喝：买蝈蝈,眼前一亮于是就买一个了。这个木匠提着蝈蝈笼子回到工地。大伙儿一看就吵嚷起来了，大家心里都很烦，你怎幺还买一笼子蝈蝈来？他把蝈蝈笼子的梁、柱、脊细细地数了一遍又一遍后，蹦起来一拍大腿说：“这不正是九梁十八柱、七十二条脊幺？大伙一听都高兴了。 受这个笼子的启发，建成了现在这个角楼。通过以上的例子可以看出我国古代早就应用了相似性原理仿真在生活实际中有广泛的应用。第2节 计算机仿真技术的应用与影响 计算机仿真的方法的产生是和计算机技术的发展与计算机应用紧密相联的。自从世界上第一台计算机的诞生，到四十时代末期第一台模拟式计算机就被用于三自由度飞机系统的仿真。五十年代末期到六十年代。因为航天科技发展的迫切需要，USA科研人员又发明了混合计算机系统，这使得人们能对较繁杂系统的功能进行仿真研究。七十年代以来，随着数字电子计算机运算速度的很大提高和对应的仿真软件的不断完善，数字电子计算机仿真得到很快发展，其应用范围也由各种工程领域扩展到非工程领域应用范围得到了广泛的提高。第2章 计算机仿真的方法计算机的仿真方法有很多种，大体可分为单步法和多步法。单步的方法主要是在用递推公式求解时，递推公式是按照步进式的，当我们知道前1时刻的数值Yn时就能计算出后一时刻的数值Yn+1, 这种方法被我们称为单步法,单步法是一种能够自动启动的算法，欧拉法就是其中单步法的一种。反之，当在求Yn+1时要用到 Yn, Yn-1, Yn-2 ,等许多个值时, 这种方法被我们称为为多步法,多步法不可以自启的，由于在使用多步法时必须要先用其他的方法计算出该时刻前面的值。本章中我们将介绍欧拉法和龙格库塔法这两种常用的计算机的仿真方法。第1节 龙格库塔法所以对于一阶精度的拉格朗日中值定理有：y(i+1)=y(i)+h*K1K1=f(xi,yi)在用点Xi处的斜率近似值的K1与右边端点Xi+1处的斜率K2的的和除于2作为平均斜率K*的近似值，这么我们就会得到二阶精度的改进拉格朗日中值定理：y(i+1)=y(i)+h*( K1+ K2)/2K1=f(xi,yi)K2=f(x(i)+h),y(i)+h*K1)按次类推，如果我们在区间Xi,Xi+1内选择多个预估几个点上的斜率值K1、K2、Km，和用他们的加权平均数作为平均斜率K*的近似值 ，显然我们能构造出具有很高精度的高阶的计算公式。经过数学的推导和求解，就可以得出四阶龙格库塔公式，也就是在工程中应用广泛的经典龙格库塔算法：y(i+1)=y(i)+h*( K1+ 2*K2 +2*K3+ K4)/6K1=f(x(i),y(i)K2=f(x(i)+h/2,y(i)+h*K1/2)K3=f(x(i)+h/2,y(i)+h*K2/2)K4=f(x(i)+h,y(i)+h*K3)通常以上说的龙格-库塔法是对四阶龙格库塔法而说的，我们可以按照二阶、三阶的情形的出我们常用的标准四阶龙格-库塔法公式。 第2节 欧拉法图1是用欧拉法对两种化学元素反应的模拟的例子。程序 输入K1，K2，a(0),b(0),c(0) T, t和N I0 II1 保存a(I, t),b(I, t) c(I, t) Y N I=N? 打印a,b,c 模拟完毕 图1 化学反应例子模拟程序框图 模拟程序使用C语言编写，程序中的初值如下：K1=0.008/g.min K20.002/min, a(0)=100g, b(0)=50g C(0)=0, T=5mins, t=0.1min, N=50其程序清单如下：#include #include float k1,k2;static float A53,B53,C53,delta,t;void strut(int);main() int i; A1=100.0; B1=50.0; C1=0.0; t=0; delta=0.1; k1=0.008; k2=0.002; for(i=1;i53;i+) printf(%2d,i); printf(%10.2f,%10.2f,%10.2f,%10.2fn,t,Ai,Bi,Ci); strut(i); return; void strut(int i) Ai+1=Ai+(k2*Ci-k1*Ai*Bi)*delta; Bi+1=Bi+(k2*Ci-k1*Ai*Bi)*delta; Ci+1=Ci+2.0*(k1*Ai*Bi-k2*Ci)*delta; t=t+delta; 计算的精度很差, 现在几乎没有人在实际工作中使用，但它导出简单， 这能说明构造数值解法一般是计算公式的基本思想, 模拟的程序也很容易懂得，于是人们愿意用它做为构造数值解法的入门例子。其一般解法如下：设给定微分方程(2.1)在区间(tn,tn+1)上求积分,得 y(tn+1)=y(tn)+f(t,y)dt (2.2)如积分间隔足够小，使得tn与tn+1之间的f (t,y)可近似的看成常数f (tn,yn), 就可以用矩形面积近似地代替在该区间上的曲线积分，于是在tn+1时的积分值为 (2.3) 将上式写成以下差分方程形式: （2.4）这就是欧拉公式。它是一个逐渐递推的差分方程，任意一个新的数值解yn+1都是根据前一个数值解以及它的求导f(tn,yn)值求得的。只要给定初始条件y0及步长h就可根据f(t0,y0)算出y1的值，再以y1算出y2，如此逐步算出y3, y4,，直到满足所需计算的范围才停止计算。欧拉法的基本思路是把连续的时间t分割成等间隔的t, 在这些离散的时刻算得函数值，根据这些值在函数图上可得到一条折线，所以欧拉法又叫折线法，其特点是分析方法简单，计算量小，但计算精度低（后面将讨论欧拉法与其它方法的比较）。下图为欧拉折线法的几何意义。如用梯形面积来代替一个小区间的曲线积分，就可以减少小矩形计算的缺点，提高精度，梯形法计算公式为(2.5)上式为隐式公式，因为公式右端含有yn+1, 这是未知的待求量，故梯形法不能自行启动运算，而要依赖于其它算法的帮助，比如说用欧拉公式法求出初值，算出的近似值，然后将其带入微分方程，计算的近似值，最后利用梯形公式反复带进去。如代一次后就认为求得了近似解，于是得到改进的欧拉法，其公式为 预估 （2.6） 校正 （2.7）上式中第一个为预估公式，第二个为校正公式。通常这种方法称为预估矫正法。在校正公式中计算了两点的斜率，再求其平均值，故计算量比欧拉法要大些。 第3节 计算机仿真方法的比较与选择（1） 一般情况下截断误差阶次越高，则解就越精确，可见所以欧拉法精度最低。（2）因为微分方程的数值解法基本思想是：通过某种离散化手续，将微分方程转化差分方程 （代数方程）来求解，而差分方程的求解还会有计算误差，而这类小扰动在传播过程中会不会恶性增长？这就是差分方程的稳定性地问题。因此，这个问题也是我们在选择积分法中必需考虑的问题。综上所述，龙格库塔法相比于欧拉法更精确、稳定、便捷，所以更能得到我们的信睐。第3章 龙格库塔法相似性仿真中的运用第1节 建模一阶零状态响应电路如图2所示，此电路方程如下：图2. 一阶零状态响应电路图R=100K,Us=20V,Uc=100，Uc(0)=0Rc+Uc=Us第2节 基于龙格库塔法的程序设计龙格库塔法的系统程序框图如图3所示， 主 函 数输 入 及 运 行 输 出 初 始 化 微分方程 数值积分 显 示 打印 绘图 右函数 图3 龙格库塔法程序简要框图程序清单如下：主程序：/Adams法积分函数#include stdio.h#include stdlib.h#include iostream.h void adambm(float A50,float B6,float C50,float D6,float *X,int nx,int nu,int nr,float dt,float *Y, float *U,float &time,int &length)static float DX50,WORK450,XS50;int j; if(length1) diffun(A,B,C,D,nx,nu,nr,X,DX,Y,U); length+; for(j=0;jnx;j+) WORKlength-1j=DXj; return; if(length4) ruk4(dt,A,B,C,D,nx,nu,nr,X,Y,time,U);length+; for(j=0;jnx;j+) WORKlength-1j=DXj; return; for(j=0;jnx;j+)XSj=Xj;Xj+=dt*(2.2916666*WORK3j-2.4583333*WORK2j+1.5416666*WORK1j-0.375*WORK0j); WORK0j=WORK1j; WORK1j=WORK2j; WORK2j=WORK3j; diffun(A,B,C,D,nx,nu,nr,X,DX,Y,U);for(j=0;jnx;j+)Xj=XSj+dt*(0.375*DXj+0.7916666*WORK2j-0.2083333*WORK1j+0.0416666*WORK0j); diffun(A,B,C,D,nx,nu,nr,X,DX,Y,U); for(j=0;jnx;j+) WORK3j=DXj; time+=dt;void output(float *Y,int nr,float T) printf(%3.1ft,T);for(int i=0;inr;i+)printf(%7.6ft,Yi);printf(Y); 第3节 VC+6.0程序运行运行上节中的C语言程序，输入步长、打印步长、最大时间分别为0.01、1、50，方程数n、m、v、w都为1，P0为0.1，U1、U2都为1，得出y0、y1。于是整理出表1数据：表1。 仿真结果数据13.6608625.5361736.8552547.7850358.4299268.9021779.2260589.4548899.61605109.73957119.80953129.86584139.90551149.93345159.96312169.96698179.97674189.98462199.98846209.99287219.99428229.99597239.99716249.99800259.99859269.99901279.99930289.99951299.99965309.99976319.99983329.99988339.99991349.99994359.99996369.99997379.99998389.99999399.99999409.99999419.99999421043104410451046104710481049105010图4是表1数据画出的图形。图4. 由表1数据做出电路响应图形当P0=0.265时的仿真结果如表2所示。表2。当P0=0.265时的仿真结果16.4525728.6017539.4488749.7827759.9143869.9662579.9867089.9947699.99793109.99919119.99968129.99987139.99995149.99998159.9999916101710181019102010211022102310241025102610271028102910301031103210331034103510361037103810391040104110421043104410451046104710481049105010图5是表2数据画出的图形。图5. 由表2数据做出电路响应图形当p0=0.43时结果如表3所示表3。当p0=0.43时的仿真结果数据18.0202529.5645139.9042049.9789359.9953669.9989879.9997889.9999599.9999910101110121013101410151016101710181019102010211022102310241025102610271028102910301031103210331034103510361037103810391040104110421043104410451046104710481049105010图6是表3数据画出的图形。图6. 由表3数据绘出的电路响应图当p0=0.595时的仿真结果如表4所示。表4. 当p0=0.595时的仿真结果数据19.0834729.9066639.9905049.9990359.9999069.9999971081091010101110121013101410151016101710181019102010211022102310241025102610271028102910301031103210331034103510361037103810391040104110421043104410451046104710481049105010图7是表4数据画出的图形。 图7. 由表4数据绘出的电路响应图可以看出通过对对电路的分析可在计算机上实现对电路的相似性仿真。第4节 实际销售产品中相似性仿真原理的应用其实不光是我们电路中可以实现仿真生活中也可以应用相似性加以仿真，假设生活实际中销售保安系统价格3000元，又假定某城镇人口为300万，令H=300万 Y=k（H-y）K=0.43 p0状态方程：G0=P0*H-P0*Y0YO为家庭保安系统套数。输入数据 0.01 仿真步长 0.2打印步长 50仿真总时间1 1 1 10.43 01 1 10令初始条件拥有安保系统而人数为零则可以得出程序运行之后的图8如下图8.以上产品销售仿真的结果图形总结及心得体会通过上面的两个完全不同方面的例子，比较电容充放电和生活实际中的保安系统的销售，可以看出它们是两个完全不同的系统，可以用相同的方法加以仿真，这样做的话会对我们的生活产生巨大的意义，可以帮助我们把复杂的问题简单化，抽象的问题具体化，这样很直观，还可以对一样商品的销售做出一定的预期和一些物理现象加以仿真。总之，完全不同的系统可以实现相同方法的去仿真。同时可以看出上面两个完全不同的系统，可以用相同的方法演变出相同的规律体现了本课题的题目思想仿真中相似性原理的应用。通过这次论文的设计我想到了很多问题，无论是现在还是将来都能应用到相似性原理解决很多事就像我们举的例子商品的销售问题，我们可以进行相似性仿真提前了解商品的销售情况，避免一定损失的风险，为商品的提高和创新提供了依据，同时我们也看到了缺点不可能进行完全的相似性仿真因为涉及的方面太广不够精确，有的系统的相似性仿真的设计需要更精确。致谢 本次论文对学生学习的能力有着巨大的提高，作为即将踏出学校门口的学生在学习阶段有很大的帮助，我学到了全面的应用专业知识和基于知识的过程中实践能力，本设计是一个全面的再学习的最后一个环节，我们都非常努力，并且自身得到提高的过程，因为它反映了大学教育的论文综合水平，学校也对毕业论文的事宜相当重视。在大学学习过程中，不仅是我们的社会，毕业论文，社会和学校生活中解决问题的能力和自主学习能力得以提高，为了完成论文，我尽量把毕业论文工作和实际的毕业设计有机地通过理论与实践很好的相结合起来。 它有助于提高我们的技能。经过一段时间的努力，学生和一些老师的帮助， 最后，我完成了这个重要问题。这次论文对我来说，有一个小小的挑战，实际上，不仅学到很多实践知识，大学的知识，丰富了我们