作者:蒋迅 随着高速公路建设的飞速发展,我们在日常生活中所见到的高速公路立交桥越来越多,而且式样也是越来越多。立体交叉是车辆汇集、分叉和转向的核心部分。立交桥的设计的好坏直接影响到行车速度和行车安全。 1922年,法国著名的规划思想家、现代建筑运动创始人之一勒·柯布西埃( Le Corbusier )出版了一本《明天的城市及其规划》( The City of To-morrow and its Planning ),在城市规划研究中首先提出多层、高速的公路立体交叉的思想。现在,立体交叉已经从城市发展到了高速公路。高速公路上的互通式立交桥由高速公路的基本路段、立交桥、匝道、交织区、收费口、监控系统和服务设施组成了一个综合体系。 立交桥的设计不仅体现在它的科学性,也体现在它的美观性。美观的立交桥也会让机动车驾驶员感到有所是从和和有所准备,为驾驶更添一份安全因素。本文以立交桥布局设计中的曲线之美为线索,聊聊相关的数学知识,并用 和 这两个数学软件制图,为读者在以后的旅途中增加一些乐趣。 在公路运输领域里,交汇处通常使用立体化和一个或多个匝道(引道)来实现至少一条高速公路上的交通能通过交叉口而不直接穿过任何其他交通车流。在这里,立交桥扮演着重要的角色。最常见的四方向高速公路立交桥有苜蓿叶型、环状型、涡轮型、风车型和环路型等,另外还有它们的一些 混合型 。我们来一一介绍。由于国际上有靠右和靠左行驶两套系统,而且高速公路与铁路和市区公路也有立交,以下我们只考虑靠右行驶道路并只考虑有四个方向的高速汽车公路的立交桥设计。 2. 苜蓿叶立交桥型 最典型的是苜蓿叶型(cloverleaf interchange)。在这里“cloverleaf”,我们指的是“four-leaf clover”这种植物。苜蓿叶型也称为四叶型和幸运草型。典型的苜蓿叶型交汇有两层,这样使得所有原来需要穿越相交道路的转向都由环形匝道来避免,也就是说,让左转车辆行驶约270度的环道后自右侧切向汇入高速公路。这四条环形匝道就形成了苜蓿叶的形状。苜蓿叶型的优点在 于 它只需要一个立交桥,也就是两层交通。因此建设经费较少。但是这样的交叉口占地面积大,路线迂回较长。更严重的是两环间的路段也容易形成交织路段,直行车辆易受转向车辆干扰,影响了高速公路的运载能力。笔者曾经遇到一位老年妇女在直行道上想上环形匝道却又无法上去结果停在了直路中间,结果躲闪不及而追尾。 这种立交桥最早是美国新泽西州Woodbridge的两条道路交叉处。这也是世界上的第一座立交桥。该立交的平均通量为每昼夜达62500辆,高峰小时交通量达6074辆,即每分钟大约可容许100辆汽车通过。苜蓿叶式在全世界各地都很多。比如下面的南京绕城高速和玄武大道立交(图1(b))。 图1.(a) 苜蓿叶型立交桥的布局。(b)南京玄武大道立交 植物学上,“clover”是三叶草。在西方很多国家(如英国、美国)长有四片叶子的三叶草。四叶草是三叶草的稀有变种(图2(a))。据说大约一万至十万株三叶草中才会有一株是四叶的。西方人认为能找到四叶草是幸运的表现,在日本则认为会得到幸福,所以又称“幸运草”。人们对这四片叶子也赋予了含义。有一种说法是:第一片叶子代表希望(hope)、第二片叶子表示信心(faith)、第三片叶子是爱情(love)、而多出来的第四片叶子则是幸运(luck)的象征。 图2. (a) 四叶草。(b) 四叶玫瑰线。(c) r = 5cos( nθ ) 数学上,我们把这样的曲线叫做“四叶玫瑰线”( Quadrifolium )。它是由极坐标方程 生成的。显然这是当 n = 2 时的玫瑰线 r = a cos( nθ ) 。我们可以很容易地将“四叶玫瑰线”的极坐标方程转换成直角坐标方程 ( x 2 + y 2 ) 3 = 4 a 2 x 2 y 3 。所以它是一个几何亏格为零的代数曲线。但如果我们需要计算它所包含的面积的话,那么还是采用极坐标来计算为宜: 当我们考虑曲线的长度时,则需要用到第二类椭圆积分了。在这里我们只给它的近似值: s = 9.86884… a 。有人说它像是中国结,也有道理。 的表达式是 PolarPlot ,{t,0,2Pi}]。建议读者到 网站上去做出 r = 5cos( nθ ) 的图像并让 n 变动起来,看看能得到一些什么图像。用 做动态模拟时的表达式为: Manipulate ,{t,0,2Pi},PlotRange?5],{n,1,10}] 图2(c)是它的效果图。 3. 环状型立交桥 第二种立交桥是环状型(stack interchange)。环状型也称为定向式(directional interchange)。中文的“环状”与英文的“stack interchange”并没有直接的联系。“stack”的意思是堆,叠加的意思。取这个名字是因为环状型多为数层叠加的原因。所以翻译成“多级立交”更为合适。笔者更倾向于称之为定向式,因为它让左转的车辆保持了左转,而不会像上面苜蓿叶型那样通过右转来实现左转。左转和右转车辆都先从最右车道上匝道,然后二者分离,左转车辆到相对象限里汇入到那里的右转车辆所在匝道,然后一起并入住主车道。环状型立交没有苜蓿叶型容易产生车流交织的缺点,也无需做270度的转弯,但其立交桥层数多,一般多为三层,也有四层和五层的例子, 因此造价相对昂贵,也容易产生视觉上的景观冲击。笔者第一次见到这样的高速系统是在休士顿。那时候中国还完全没有高速公路的概念,所以见到这样的大型立交桥时觉得非常震撼。 第一座四层定向型立交桥在美国洛杉矶市,是州际I-10和US101的交汇处。它的第二、第四层为主干线,每层有六个车道;第一、第三层为左转匝道。其最上一层高出地面14.4米,最下层低于地面6.6米。干线设计车速96公里/小时,匝道设计车速55公里/小时,交通量达75000-100000辆/昼夜,耗资约280万美元。现在在中国也有这种立交模式,比如上海延安东路就有这样一座环状型立交桥(图3(b))。 图3.(a) 环状型立交桥的布局。(b) 上海延安东路的环状型立交桥 不同于苜蓿叶型立交桥,我们没有找到一条漂亮的数学曲线来代表这种形状的立交桥。最接近的应该是“内旋轮线”( Hypotrochoid )。给定一个半径为 R 的褂讪的大圆和一个内切于大圆的半径为 r 的小圆,从这个小圆的圆心出发做一个褂讪在小圆的射线,然后在这条射线上取一个点 P ,点 P 可以在小圆之外。点 P 到小圆中心的距离为 d 。当这个小圆沿着大圆的内边滚动时,点 P 的轨迹就叫作“内旋轮线”。这个曲线的参数方程是: 注意虽然我们把滚动的圆称为小圆,其实我们并不假定 r R 。三个参数 r , R 和 d 之间没有任何限制。它们甚至可以是负数。依据它们的取值的不同,我们可以得到许多不同形状的曲线。比如当 d = r = R /2 时,我们就得到一条线段;当 d = 0 时,我们就得到一个圆。下面是一组对应于不同的参数值 ( R , r , d ) 的“内旋轮线”。我们可以感受到这些曲线是多么地不同。 图4. 取不同参数时的内旋轮线。 内旋轮线在 中的一般表达式是 使用上面的表达式以及 中的 Table,ParametricPlot 和 Grid 等指令,我们就得到了上面的图4。 注意当 R = 1, r = 3/4, d = 5/13 时,“内旋轮线”最接近于环状型立交桥。所不同的是,立交桥的四个“叶子”是尖状的,而“内旋轮线”是光滑的。匝道在接入主线时必须是与主线相切地接入。 既然我们不能用一个数学方程来描述环状型,那么我们干脆把四个左转匝道单独出来,然后只看其中一段。其他匝道都可以通过旋转变换来实现。让我们单独出其中的由东向南的一段左转路线来。对 于 这样的路线,最好的数学公式恐怕是贝济埃曲线( Bezier curve )了。如下图所示,贝济埃曲线可以在给定的两个点上按一点的方向连接,而这正好是匝道接入主线时所要求的。 图5.(a) 环状型立交桥由东向南的一段。(b) 贝济埃曲线示意图。(c) Wolfram给出的一个图案。 Wolfram用贝济埃曲线做出了许多漂亮的 花型图案 。但这些图案都不能满足我们这里的要求(图5.(c))。 我们还可以再进一步,将上面的曲线做一个45 o 逆时针转轴。 于 是我们看到的是近似于一条椭圆曲线( elliptic curve )。数学上,椭圆曲线是一个由代数方程 定义的曲线。下面图6(b)是当 a = 1 和 b = 4 时的椭圆曲线,是用 制作的。用 制作也很简单,我们略过。基 于 椭圆曲线数学,人们开发了一种建 立公开密钥加密的算法 。 图6.(a) 做45 o 逆时针转轴后的一段匝道。(b) 当 a = 1 和 b = 4 时的椭圆曲线。 未完待续。
微积分 偏微分D ,D 多重偏微分D 全微分Dt ,Dt 多重全微分Dt 求隐函数的导数 In := Dt == x^2, x] Out = 10 y Dt + Cos Dt == 2 x In := Solve ] Out = {{Dt - (2 x)/(10 y + Cos )}} 积分 积分Integrate 1}] 多重积分Integrate }] 第一重积分在最后面 数值积分 Nintegrate ; 数值积分 Nintegrate ,以x1,x2...为分割求 上的数值积分 解微分方程 (*y 与 y的区别*) In := DSolve == 2 y , y , x] Out = {{y - E^(2 x) C }} In := DSolve == 2 y , y, x] Out = {{y - Function ]}} (*一元二次方程*) In := DSolve + 2 y' + y == 0, y , x] Out = {{y - E^-x C + E^-x x C }} (*二元一次方程*) In := DSolve == -z' , z == -y' }, {y , z }, x] Out = {{y - 1/2 E^-x (1 + E^(2 x)) C - 1/2 E^-x (-1 + E^(2 x)) C , z - -(1/2) E^-x (-1 + E^(2 x)) C + 1/2 E^-x (1 + E^(2 x)) C }} (*带初始条件解微分方程*) DSolve == y , y == 5}, y , x] Out = {{y - 5 E^x}} In := DSolve == y , y' == 5}, y , x] Out = {{y - E^-x (-5 + C + E^(2 x) C )}} 数值解微分方程 s = NDSolve == y ^2, y == 1}, y, {x, 0, 1}] Plot /. s], {x, .01, 1}]
一、 With 函数的用法 With specifies that in expr occurrences of the symbols x, y, ... should be replaced by x0,y0,... . 关于 With 函数,有两点必需注意 : 1 、在 With 函数中出现的局部变量必需赋初值, With 是不合法的。 2 、在 With 函数中出现的局部变量不可在 expr 中再被赋值。 例 1 : In := With ] With::lvws: Variable x in local variable specification {x} requires a value. Out = With ] 例 2 : In := With Set::setraw: Cannot assign to raw object 5. Out = 4 二、 Module 函数的用法 Module specifies that occurrences of the symbols x, y,... in expr should be treated as local. Module defines initial values for x,... . 关于 Module 函数,有一点必需注意: Module 中的局部变量 {x,y,...} 会相应地被重先定义为 {x$nnn,y$nnn,...} (这里 nnn 代表某一整数,随使用次数变化),同时, expr 中出现的所有 x,y,... 都会相应地被 x$nnn,y$nnn,... 所替代。 例: In := Module ] Module ] Module ] t$115 t$116 t$117 三、 Block 函数的用法 Block specifies that expr is to be evaluated with local values for the symbols x, y,... . Block defines initial local values for x,... . 关于 Block 函数,有一点必需注意: 与 Module 一样同样定义了一组 {x,y,...} 的局部变量,但没有 {x,y,...} 被 {x$nnn,y$nnn,...} 替代这一步。因此局部变量如果没有在 Block 中赋初值,则使用全局变量的值代替。 例: In := i=2; In := Block Out := 2+a In := Module Out := a + i$16 四、三个函数的共同点 三个函数都各自定义了一组局部变量,对局部变量的使用不会导致全局变量的值的变化。 例: In := x = 3; With Module Block x Out = 4 Out = 25 Out = 36 Out = 3 五、三个函数的不同点 1 、在 With 函数中出现的局部变量必需赋初值,且在 With 函数中出现的局部变量不可在 expr 中再被赋值,其它两个函数没此限制。 2 、在 Block 函数中,若定义的局部变量没有被赋初值,则使用全局变量的值,而在 Module 函数中,由于局部变量都被重先定义为新的量变,因此不可能使用全局变量。 如上例所述: In := i=2; In := Block Out := 2+a In := Module Out := a + i$16 3 、在运行 With 时,不管 expr 执不执行, x 、 y 的初值 x0 、 y0 都得首先代入到 expr 中,而 Module 和 Block 是首先执行 expr ,然后再代初值。 例: In := {With ], Module ], Block ]} Out = {Hold , Hold , Hold } 4 、在速度方面, With 的运行速度最快,其次是 Block ,最后是 Module 。 例: In := Do , {10^5}] // Timing Do , {10^5}] // Timing Do , {10^5}] // Timing Out = {0.078, Null} Out = {0.094, Null} Out = {0.25, Null}
科技数学计算 Wolfram Mathematica 7.0 for Students 软件下载 .Wolfram.Mathematica.7.0.for.Students.Win32.rar (452.4MB) .Wolfram.Mathematica.7.0.for.Students.Win32.Keygen.rar (245.1KB) 内容简介 Wolfram Research发布了Mathematica的第7个版本,它如今能做得已经远超过符号计算。Mathematica 7.0最值得关注的功能是引入了并行计算,可以充分发挥多核处理器的计算能力。 其它的新特性包括与Microsoft Word公式编辑器简单的粘贴整合, 数学对象的即时3D模型,可能是有史以来最昂贵的Photoshop克隆——内置图像处理和分析功能等。 Mathematica提供了基因组、化学、气象、天文学、金融、测地学数据的完整支持,令其成为科学研究不可或缺的工具。 在在美国德克萨斯州举行的超级计算大会(SC08)上,Wolfram Research宣布新版通用计算软件Mathematica 7将支持NVIDIA CUDA GPU并行计算技术。 据称,融入CUDA技术的Mathematica 7可以在数学运算、建模、模拟和视觉计算等方面获得10倍甚至100倍的性能提升,而且不需要用户另外学习或者编写C语言代码。 配合NVIDIA同时发布的Tesla个人超级计算机,Mathematica 7可以在桌面上执行复杂的数据计算任务,免除了编写原生C程序或等待公共集群的麻烦,为研究人员节省大量的工作时间。 支持CUDA加速的新版Mathematica 7将在2009年第一季度提供,而刚刚发布的普通版Mathematica 7已经面世,在12个应用领域增加了500多个新功能,包括图像处理、并行高性能计算、矢量场视觉化、布尔数学体系运算、统计模式分析等等。 Mathematica可以说是世界上最强大的通用计算系统,自1988年首次发布以来已经拥有数百万用户,对如何在科技和其它领运用计算机有着深刻的影响,并成为许多机构的标准工具,比如财富50强的所有公司、美国gov15个主要机构、全球最大的50所大学等等。 Mathematica 7支持Windows 2000/XP/Vista、Mac OS X、Linux x86、Solaris、UltraSPARC/x86等平台,建议零售价2495美元(北美)或3120美元(亚洲),还提供1095/1315美元的教育版和140美元的学生版。 Following closely on the dramatic reinvention of Mathematica in 2007, Mathematica 7 continues the momentum of innovation to deliver an array of new capabilities, greatly extending the state of the art in many areas, and bringing a dozen major new application fields into the integrated framework of Mathematica.