搜档网
当前位置:搜档网 › .欧拉方程具有伽利略变换的不变性

.欧拉方程具有伽利略变换的不变性

.欧拉方程具有伽利略变换的不变性

3.欧拉方程具有伽利略变换的不变性 在理想流体力学中动力学基本方程是欧拉方程:

或 2v 1(v )v G (1)v 11v v (v)G 2P t P t ρρ?+??=-?????????+?-???=-?? 下面证明欧拉方程在惯性坐标系变换下的协变性:

在方程(1)中G 、ρ、P 、t 是不变量,可直接变换为G /、ρ/、P /、t /

;v 变换为v /+u .其中u 是常矢,故

v u v u v u v v v v u v v '??+'=+'??+'=??'

?'?=?+'?=??)()()()(;)(t t t 再考虑算符▽z y x ??+??+??=k j i

的坐标变换,单位矢i 、j 、k 都是不变量,可用i /、j /、k /代入,y 、z 用y /、z /代入.但

x x t u x ut x x x x x x '

????-=?-?'??=?'?'??=??)1()(, 当算符▽所作用场量为压强P 时,t 与x 可认为是独立坐标,从而P P x x x t '?'=?'

??=??=??,,0 当算符▽作用于场量v 时,t 与x 是相关的,x v dt

dx t x ==??,从而x v u x x '??-=??)1( =?∴?'??'+'-?'='??'-?'='??'+'??'+'??'-x u v u x v u z y x v u x

x x i i k j i )1(-----(2) 将(2)式代入

?'?'='

??-?'?+?'?'='??+'?+'-?'?+?'?'='??'+'-

?'?+'=??+'=??v u v i u v u v i u v u v v x u x u v u x u v u x x )()()()( 欧拉方程最终变换为:

P t '?''-'='?'?'+'?'?ρ1)(G v v v .可见,欧拉方程在x /系中的形式与在x 系中形式完全相同.

欧拉方程在惯性坐标系变换下协变是意料

中的,因为欧拉方程是牛顿运动定律在流体力

学中的表达,而牛顿运动定律对伽利略变换是

协变的,故对欧拉方程自然也协变.

欧拉方程的求解教材

欧拉方程的求解 1.引言 在数学研究领域,我们经常会看到以数学家名字命名的概念、公式、定理等等,让人敬佩跟羡慕.但是,迄今为止,哪位数学家的名字出现得最多呢?他就是数学史上与阿基米德、牛顿、高斯齐名的“四杰”之一,人称“分析学的化身”的盲人数学家欧拉(Leonhard Euler,1707--1783). 几乎在每一个数学领域都可以看到他的名字,譬如我们熟悉的“欧拉线”、“欧拉圆”、“欧拉公式”、“欧拉定理”、“欧拉函数”、“欧拉积分”、“欧拉变换”、“欧拉常数” 欧拉还是许多数学符号的发明者,例如用π表示 圆周率、e 表示自然对数的底、()f x 表示函数、∑表示求和、i 表示虚数单位 以欧拉命名的数学名词有很多,本文主要讲解以欧拉命名的方程即“欧拉方程”. 在文献[1]中,关于欧拉方程的求解通常采用的是变量变换的方法.变量变换法就是将所求的欧拉方程化为常系数齐次线性微分方程,然后再来求解这个常系数齐次线性微分方程的解,亦即求其形如K y x =的解,进而求得欧拉方程的解. 但有些欧拉方程在用变量变换法求解时比较困难.本文在所学的欧拉方程的求解的基础上,对欧拉方程进行了简单的分类,并针对不同阶的欧拉方程的求解给出了不同的定理.最后在每类欧拉方程后面给出了典型的例题加以说明. 2.几类欧拉方程的求解 定义1 形状为 ()1(1)110n n n n n n y a x y a xy a y x ---'++++= (1) 的方程称为欧拉方程. (其中1a ,2a , ,1n a -,n a 为常数)

2.1二阶齐次欧拉方程的求解(求形如K y x =的解) 二阶齐次欧拉方程: 2120x y a xy a y '''++=. (2) (其中1a ,2a 为已知常数) 我们注意到,方程(2)的左边y ''、y '和y 的系数都是幂函数(分别是2x 、1a x 和02a x ) ,且其次依次降低一次.所以根据幂函数求导的性质,我们用幂函数K y x =来尝试,看能否选取适当的常数K ,使得K y x =满足方程(2). 对K y x =求一、二阶导数,并带入方程(2),得 212()0K K K K K x a Kx a x -++= 或 212[(1)]0K K a K a x +-+=, 消去K x ,有 212(1)0K a K a +-+=. (3) 定义2 以K 为未知数的一元二次方程(3)称为二阶齐次欧拉方程(2)的特征方程. 由此可见,只要常数K 满足特征方程(3),则幂函数K y x =就是方程(2)的解. 于是,对于方程(2)的通解,我们有如下结论: 定理1 方程(2)的通解为 (i) 1112ln K K y c x c x x =+, (12K K =是方程(3)的相等的实根) (ii)1212K K x c x y c +=, (12K K ≠是方程(3)的不等的实根) (iii)12cos(ln )sin(ln )x x c x x y c ααββ+=.(1,2K i αβ=±是方程(3)的一对共轭复根) (其中1c 、2c 为任意常数)

微分方程数值解法

《微分方程数值解法》 【摘要】自然界与工程技术中的很多现象,可以归结为微分方程定解问题。其中,常微分方程求解是微分方程的重要基础内容。但是,对于许多的微分方程,往往很难得到甚至不存在精确的解析表达式,这时候,数值解提供了一个很好的解决思路。,针对于此,本文对常微分方程数值解法进行了简单研究,主要讨论了一些常用的数值解法,如欧拉法、改进的欧拉法、Runge —Kutta 方法、Adams 预估校正法以及勒让德谱方法等,通过具体的算例,结合MA TLAB 求解画图,初步给出了一般常微分方程数值解法的求解过程。同时,通过对各种方法的误差分析,让大家对各种方法的特点和适用范围有一个直观的感受。 【关键词】 常微分方程 数值解法 MA TLAB 误差分析 引言 在我国高校,《微分方程数值解法》作为对数学基础知识要求较高且应用非常广泛的一门课程,不仅 在数学专业,其他的理工科专业的本科及研究生教育中开设这门课程.近四十年来,《微分方程数值解法》不论在理论上还是在方法上都获得了很大的发展.同时,由于微分方程是描述物理、化学和生物现象的数学模型基础,且它的一些最新应用已经扩展到经济、金融预测、图像处理及其他领域 在实际应用中,通过相应的微分方程模型解决具体问题,采用数值方法求得方程的近似解,使具体问题迎刃而解。 2 欧拉法和改进的欧拉法 2.1 欧拉法 2.1.1 欧拉法介绍 首先,我们考虑如下的一阶常微分方程初值问题 ???==0 0)() ,('y x y y x f y (2--1) 事实上,对于更复杂的常微分方程组或者高阶常微分方程,只需要将x 看做向量,(2--1)就成了一个一阶常微分方程组,而高阶常微分方程也可以通过降阶化成一个一阶常微分方程组。 欧拉方法是解常微分方程初值问题最简单最古老的一种数值方法,其基本思路就是把(2--1)中的导数项'y 用差商逼近,从而将一个微分方程转化为一个代数方程,以便求解。 设在[]b a ,中取等距节点h ,因为在节点n x 点上,由(2--1)可得:

对于欧拉方程的理解

关于欧拉方程的理解 1755年,瑞士数学家L.欧拉在《流体运动的一般原理》一书中首先提出这个方程。 形如:)(1)1(11)(x f y x p y x p y x n n n n n ='+++--- (1) 的方程称为欧拉方程, 其中n p p p ,,,21 为常数。 欧拉方程的特点是: 方程中各项未知函数导数的阶数与其乘积因子自变量的幂次相同。 现阶段欧拉方程的应用领域很广,现只结合流体力学来探讨我对于欧拉方程的理解。 欧拉方程提出采用了连续介质的概念,把静力学中压力的概念推广到了运动流体中。 流体静力学着重研究流体在外力作用下处于平衡状态的规律及其在工程实际中的应用。 这里所指的静止包括绝对静止和相对静止两种。以地球作为惯性参考坐标系,当流体相对于惯性坐标系静止时,称流体处于绝对静止状态;当流体相对于非惯性参考坐标系静止时,称流体处于相对静止状态。 流体处于静止或相对静止状态,两者都表现不出黏性作用,即切向应力都等于零。所以,流体静力学中所得的结论,无论对实际流体还是理想流体都是适用的。 流体静压强的特性 1静压强的方向—沿作用面的内法线方向 2任一点的流体静压强的大小与作用面的方向无关,只与该点的位置有关

由上图可以推到出流体平衡微分方程式,即欧拉平衡方程 x y z p f x p f y p f z ρρρ??=?????=?????=??? 当流体处于平衡状态时,单位体积质量力在某一轴向上的分力,与压强沿该轴的递增率相平衡。 这里的fx 、fy 、fz 是流体质量力在x 、y 、z 轴上的投影,且质量力中包含以下两项:重力和惯性力。在这里如果假定fx 、fy 、fz 仅仅是重力在三个坐标轴上的投影,那么惯性力在x 、y 、z 轴上的投影分别为:-du/dt ,-dv/dt 和-dw/dt 。于是,上式便可写成 d d d d d d x y z u p f t x v p f t y w p f t z ρρρ????-= ???? ??????-=? ??? ??????-=? ??? ?? 上式整理后可得:

最新薛定谔方程及其解法

关于薛定谔方程 一.定义及重要性 薛定谔方程(Schrdinger equation)是由奥地利物理学家薛定谔提 出的量子力学中的一个基本方程,也是量子力学的一个基本假定, 其正确性只能靠实验来检验。是将物质波的概念和波动方程相结合 建立的二阶偏微分方程,可描述微观粒子的运动,每个微观系统都 有一个相应的薛定谔方程式,通过解方程可得到波函数的具体形式 以及对应的能量,从而了解微观系统的性质。 薛定谔方程是量子力学最基本的方程,亦是量子力学的一个基 本假定,它的正确性只能靠实验来检验。 二.表达式 三.定态方程 ()() 2 2 2 V r E r m η ψψ + ?? -?= ?? ?? 所谓势场,就是粒子在其中会有势能的场,比如电场就是一个带电粒子的势场;所谓定态,就是假设波函数不随时间变化。 其中,E是粒子本身的能量;v(x,y,z)是描述势场的函数,假设不随时间变化。

2 2 22222 z y x ??????++=? 可化为 d 0)(222 =-+ψψv E h m dx 薛定谔方程的解法 一. 初值解法;欧拉法,龙格库塔法 二. 边值解法;差分法,打靶法,有限元法 龙格库塔法(对欧拉法的完善) 给定初值问题 ). ()()((3) ) ,(),()( ,,(2) )(),( 311212 2111021h O t y t y hk y h t f k y t f k k c k c h y y y c c a y b t a y t f dt dy i i i i i i i i =-???????++==++==?????=≤≤=++的局部截断误差使以下数值解法的值及确定常数ββα βα

常微分方程作业欧拉法与改进欧拉法

P77 31.利用改进欧拉方法计算下列初值问题,并画出近似解的草图:dy + =t = t y y ≤ ≤ ,2 ;5.0 0,3 )0( )1(= ,1 ? dt 代码: %改进欧拉法 function Euler(t0,y0,inv,h) n=round(inv(2)-inv(1))/h; t(1)=t0; y(1)=y0; for i=1:n y1(i+1)=y(i)+h*fun(t(i),y(i)); t(i+1)=t(i)+h; y(i+1)=y(i)+1/2*h*(fun(t(i),y(i))+ fun(t(i+1),y1(i+1))) end plot(t,y,'*r') function y=fun(t,y); y=y+1; 调用:Euler(0,3,[0,2],0.5) 得到解析解:hold on; y=dsolve('Dy=y+1','(y(0)=3)','t'); ezplot(y,[0,2]) 图像:

dy y =t - t y ;2.0 t = ≤ )0( 0,5.0 ,4 )2(2= ≤ ? ,2 dt 代码: function Euler1(t0,y0,inv,h) n=round(inv(2)-inv(1))/h; t(1)=t0; y(1)=y0; for i=1:n y1(i+1)=y(i)+h*fun(t(i),y(i)); t(i+1)=t(i)+h; y(i+1)=y(i)+1/2*h*(fun(t(i),y(i))+ fun(t(i+1),y1(i+1))) end plot(t,y,'*r') function y=fun(t,y); y=y^2-4*t; 调用: Euler1(0,0.5,[0,2],0.2) 图像:

微分方程常用的两种数值解法:欧拉方法与龙格—库塔法

四川师范大学本科毕业论文 微分方程常用的两种数值解法:欧拉方法与龙 格—库塔法 学生姓名XXX 院系名称数学与软件科学学院 专业名称信息与计算科学 班级2006级 4 班 学号20060640XX 指导教师Xxx 四川师范大学教务处 二○一○年五月

微分方程常用的两种数值解法:欧拉方法与龙格—库塔法 学生姓名:xxx 指导教师:xx 【内容摘要】微分方程是最有生命力的数学分支,在自然科学的许多领域中,都 会遇到常微分方程的求解问题。当前计算机的发展为常微分方程的应用及理论研究提供了非常有力的工具,利用计算机解微分方程主要使用数值方法,欧拉方法和龙格——库塔方法是求解微分方程最典型常用的数值方法。本文详细研究了这两类数值计算方法的构造过程,分析了它们的优缺点,以及它们的收敛性,相容性,及稳定性。讨论了步长的变化对数值方法的影响和系数不同的同阶龙格—库塔方法的差别。通过编制C程序在计算机上实现这两类方法及对一些典型算例的结果分析比较,能更深切体会它们的功能,优缺点及适用场合,从而在实际应用中能对不同类型和不同要求的常微分方程会选取适当的求解方法。 关键词:显式单步法欧拉(Euler)方法龙格—库塔(Runge—Kutta)方法截断误差收敛性 Two commonly used numerical solution of differential equations:Euler method and Runge - Kutta method Student Name: Xiong Shiying Tutor:Zhang Li 【Abstract】The differential equation is the most vitality branch in mathematics. In many domains of natural science, we can meet the ordinary differential equation solution question. Currently, the development of computer has provided the extremely powerful tool for the ordinary differential equation application and the fundamental research, the computer solving differential equation mainly uses value method. The Euler method and the Runge—Kutta method are the most typical commonly value method to solve the differential equation. This article dissects the structure process of these two kinds of values commonly value method to solve the analyses their good and bad points, to their astringency, the compatibility, and the stability has made the proof. At the same time, the article discuss the length of stride to the numerical method changing influence and the difference of the coefficient different same step Runge—kutta method. Through establishing C program on the computer can realize these two kind of methods, Anglicizing some models of calculate example result can sincerely realize their function, the advantage and disadvantage points and the suitable situation, thus the suitable solution method can be selected to solve the different type and the

伽利略科学研究方法

伽利略(1564—1642)是意大利物理学家、天文学家,近代实验科学的创始人。他的一生完全献给了科学事业,他所取得的巨大成就开创了近代物理的新纪元。正如爱因斯坦所说:“伽利略的发现以及他所应用的科学推理方法是人类思想史上最伟大的成就之一,标志着物理学的开端。”他1632年发表了《关于两个世界体系的对话》一书,1638年又发表了《关于两门新科学的对话》一书。伽利略在长达几十年的科学研究工作期间开创了许多物理学研究方法,对今天的科学研究人员来说,在科技创新方面仍有重要的指导作用。本文试图从这些方面进行一些探究。1伽利略主要的力学研究工作自由落体问题的研究是伽利略力学研究的突破口。当时在力学问题上流行的是亚里士多德的理论。亚里士多德认为落体以匀速下落,其速度的大小与落体的重量成正比。伽利略首先指出了亚里士多德落体观念的逻辑矛盾。他假定一根不太长的绳子,两端分别系着一块石头,这两块石头的重量不同。那么,这两块石头将以什么速度下落呢?按照亚里士多德的观念,它们的重量是大小两块石头重量之和,所以它们下落的速度比任一块石头单独下落的速度都要快。另一方面,也从亚里士多德的观念出发,大石头下落得快,小石头下落得慢,则当两石头串联在一起时,会出现这样的情况:大石头快落在下,小石头慢落在上,在大石头带动下,小石头比单独下落时要快些,而大石头被小石头拖后,使之比单独下落时要慢些。即同是应用亚里士多德的观念,得出的是以上两种自相

矛盾的结果。所以亚里士多德关于落体的理论从逻辑推理上就不攻自破了。还是眼见为实,伽利略知道仅用逻辑推理是不够的,还必须用人们能够观察到的事实来驳斥亚里士多德的落体观念。相传1589年伽利略登上了意大利的比萨斜塔,让10磅重和1磅重的两个球同时下落。塔下的人都看到,这两个重量不同的球几乎是同时落地的。而根据亚里士多德的落体观念,当大球落到地面时,小球只下落到1/10的高度,这显然不符合眼见的事实。做这个实验之后,伽利略想到,有人会说物体下落速度虽然不是同重量成正比,但重物看起来总是比轻物似乎要落得快一些。由于比萨斜塔只有56米高,相对高度而言,球下落得太快了,肉眼不容易看出两者的差距。所以伽利略就想“冲淡引力”,让球落得慢一些。这样就可以比较容易地得到两个重量不同的球究竟是先后落地,还是同时落地的结论。伽利略是通过斜面实验来达到“冲淡引力”的设想的。他在长约8米的木板上,刻着一条光滑的槽,并放置成一斜面,斜面的夹角可以随意调控。他使重量不同的小球在同一高度沿斜面同时滚下。夹角越小,小球滚得就越慢。这就好比冲淡了引力。伽利略发现,重量不同的球在相同的斜面上滚动的速度是相同的,与斜面的夹角的大小无关。当斜面夹角为90度时,小球的滚动就成了自由下落。于是他得出结论:物体自由下落的速度同其重量无关。伽利略为了研究落体运动,不断人为地调整木板与水平面的夹角,观察小球在人为控制下运动,这本身就是一种典型的科学实验(这个实验曾在2002年被英国著名的《物理学世界》杂志的广大读者评为历史上“最美丽”的十大物理实验之一)。伽利略在

实验三 定态薛定谔方程的矩阵解法

实验三 定态薛定谔方程的矩阵解法 一.实验目的 1.掌握定态薛定谔方程的矩阵解法。 2.掌握几种矩阵特征值问题数值解法的原理,会调用相应的子程序求解具体问题。 二.实验内容 1.问题描述 以/2ω/()m ω为长度单位,一维谐振子的哈密顿量为 2 202d H x dx =-+, 其本征值为21n E n =+,本证波函数为 2 /2)()n n x H x ?=-, 其中()n H x 为厄米多项式,满足递推关系 11()2()2()n n n H x xH x nH x +-=-。 用矩阵方法求 2 22d H x x dx =-++ 的本证能量和相应的波函数。 2.问题分析 H E ψψ= 0()|j j j t c ψ?∞ ==>∑ 0||i i j i j i j c E c x Ec ??∞ =+<>=∑ 11|j j j x ???-+>=>>

11||||j j j j x x ????-+<>= <>= 0010010 112111,211,11,1 n n n n n n n n n n n n E x c c x E x c c E x E x c c x E c c -------?????????????????????????=??????????????????????? ? 3.程序编写 子程序及调用方法见《FORTRAN 常用算法程序集(第二版)》第三章 徐士良,P97 4.实验要求 ◆用恰当的算法求解以上实对称三对角矩阵的特征值问题。 ◆取n=8,给出H 的全部特征值和相应的特征向量。 5.实验步骤 ● 启动软件开发环境Microsoft Developer Studio 。 ● 创建新工作区shiyan03。 ● 创建新项目xm3。 ● 创建源程序文件xm3.f90,编辑输入源程序文本。 ● 编译、构建、运行、调试程序。 6.实验结果 程序设计:

欧拉及改进的欧拉法求解常微分方程

生物信息技术0801 徐聪U200812594 #include #include void f1(double *y,double *x,double *yy) { y[0]=2.0; x[0]=0.0; yy[0]=2.0; for(int i=1;i<=9;i++) { x[i]=x[i-1]+0.2; y[i]=y[i-1]+0.2*(y[i-1]-x[i-1]); yy[i]=x[i]+1+exp(x[i]); printf("若x=%f,计算值是%f,真实值是%f,截断误差是%f\n ",x[i],y[i],yy[i],y[i]-yy[i]); } }; void f2(double *y,double *x,double *yy) { y[0]=1.0; x[0]=0.0; yy[0]=1.0; for(int i=1;i<=9;i++) { x[i]=x[i-1]+0.2; y[i]=y[i-1]+0.2*(2*y[i-1]+x[i-1]*x[i-1]); yy[i]=-0.5*(x[i]*x[i]+x[i]+0.5)+1.25*exp(2*x[i]); printf("若x=%f,计算值是%f,真实值是%f,截断误差是%f\n ",x[i],y[i],yy[i],y[i]-yy[i]); } }; void f3(double *y,double *x,double *yy,double *y0) { y[0]=2.0; x[0]=0.0; yy[0]=2.0; for(int i=1;i<=9;i++) { x[i]=x[i-1]+0.2; y0[i]=y[i-1]+0.2*(y[i-1]-x[i-1]); y[i]=y[i-1]+0.1*(y[i-1]-x[i-1]+y0[i-1]-x[i-1]);

研究性学习伽利略的研究艺术

伽利略科学研究方法探究 伽利略(1564—1642)是意大利物理学家、天文学家,近代实验科学的创始人。他的一生完全献给了科学事业,他所取得的巨大成就开创了近代物理的新纪元。正如爱因斯坦所说:“伽利略的发现以及他所应用的科学推理方法是人类思想史上最伟大的成就之一,标志着物理学的开端。”他1632年发表了《关于两个世界体系的对话》一书,1638年又发表了《关于两门新科学的对话》一书。伽利略在长达几十年的科学研究工作期间开创了许多物理学研究方法,对今天的科学研究人员来说,在科技创新方面仍有重要的指导作用。本文试图从这些方面进行一些探究。1伽利略主要的力学研究工作自由落体问题的研究是伽利略力学研究的突破口。当时在力学问题上流行的是亚里士多德的理论。亚里士多德认为落体以匀速下落,其速度的大小与落体的重量成正比。伽利略首先指出了亚里士多德落体观念的逻辑矛盾。他假定一根不太长的绳子,两端分别系着一块石头,这两块石头的重量不同。那么,这两块石头将以什么速度下落呢?按照亚里士多德的观念,它们的重量是大小两块石头重量之和,所以它们下落的速度比任一块石头单独下落的速度都要快。另一方面,也从亚里士多德的观念出发,大石头下落得快,小石头下落得慢,则当两石头串联在一起时,会出现这样的情况:大石头快落在下,小石头慢落在上,在大石头带动下,小石头比单独下落时要快些,而大石头被小石头拖后,使之比单独下落时要慢些。即同是应用亚里士多德的观念,得出的是以上两种自相矛盾的结果。所以亚里士多德关于落体的理论从逻辑推理上就不攻自破了。还是眼见为实,伽利略知道仅用逻辑推理是不够的,还必须用人们能够观察到的事实来驳斥亚里士多德的落体观念。相传1589年伽利略登上了意大利的比萨斜塔,让10磅重和1 磅重的两个球同时下落。塔下的人都看到,这两个重量不同的球几乎是同时落地的。而根据亚里士多德的落体观念,当大球落到地面时,小球只下落到1/10的高度,这显然不符合眼见的事实。做这个实验之后,伽利略想到,有人会说物体下落速度虽然不是同重量成正比,但重物看起来总是比轻物似乎要落得快一些。由于比萨斜塔只有5 6米高,相对高度而言,球下落得太快了,肉眼不容易看出两者的差距。所以伽利略就想“冲淡引力”,让球落得慢一些。这样就可以比较容易地得到两个重量不同的球究竟是先后落地,还是同时落地的结论。伽利略是通过斜面实验来达到“冲淡引力”的设想的。他在长约8米的木板上,刻着一条光滑的槽,并放置成一斜面,斜面的夹角可以随意调控。他使重量不同的小球在同一高度沿斜面同时滚下。夹角越小,小球滚得就越慢。这就好比冲淡了引力。伽利略发现,重量不同的球在相同的斜面上滚动的速度是相同的,与斜面的夹角的大小无关。当斜面夹角为90度时,小球的滚动就成了自由下落。于是他得出结论:物体自由下落的速度同其重量无关。伽利略为了研究落体运动,不断人为地调整木板与水平面的夹角,观察小球在人为控制下运动,这本身就是一种典型的科学实验(这个实验曾在2002年被英国著名的《物理学世界》杂志的广大读者评为历史上“最美丽”的十大物理实验之一)。伽利略在斜塔上不能“冲淡引力”,他在家里通过斜面实验就可以做到这一点,弥补了斜塔观察活动的不足。 伽利略在斜面实验的基础上,利用数学的方法,确定了路程与时间的数量关系为:s正比于t2,这就是时间平方定律。在他的实验记录上,有这样两列数字,不同的下落时间t=1,2,3,4……;物体下落的距离之比s=1:4:9:16……。从这两列数的比例关系,伽利略证明沿斜面下滑的物体正在做匀加速运动,从而也证明了自由落体运动是匀加速直线运动。斜面实验还使伽利略发现了惯性定理。他做了两个斜面,上面都刻有一条光滑的槽,让小球从第一个斜面滚下,再爬上第二个斜面。他发现,当小球在第二个斜

伽利略对自由落体运动的研究

2.6 伽利略对自由落体运动的研究 教材分析 本节内容是让学生了解并学习伽利略研究自由落体运动的科学思维方法和巧妙的实验构思.教材编写的脉络清楚,逻辑推理严谨,文字表述生动、通俗易懂,因此,适合于学生自主学习. 本节是新教材注重过程与方法、情感态度和价值观的一个标志性内容.过去的教学过分注重对知识与技能的掌握,而忽略了对科学精神、科学研究方法的培养.因此,能否通过这节课的学习让学生体会到人类对自然世界的探究思想和方法,感受到一位伟大的科学家的高尚情操,就成为这节课最终的目标.为了更好地落实新课标的精神,该教学策略采用了先让学生收集相关资料,在课堂上经过讨论和发表见解,充实和完善伽利略的研究过程与方法.引导学生一步步体会伽利略严谨的科学态度、不畏强权的探索精神和正确地解决问题的思路,树立正确的科学观念. 教学目标 (一)教学目标: 1、了解伽利略对自由落体运动的研究思路和方法; 2、能够合理设计实验,并将实验数据用图线法处理。 (二)过程和方法: 1、经历伽利略对自由落体运动的研究方法,感悟科学探究的方法; 2、分组进行科学探究活动,完成实验操作; 3、培养学生进行数学推理和图象处理数据的能力。 (三)情感目标: 1、激发了学生学习伽利略敢于向权威挑战,善于观察思考,知难而进的优秀品质; 2、培养学生耐心细致的意志品质,创新思想和互相协作的精神。 教学重点 通过重现重大发现的历史过程,让学生亲临其境探究伽利略对自由落体运动研究的实验,学习其科学思维方法和巧妙的实验构思。 教学难点 1、当无法验证自由落体运动是初速度为零的匀加速直线运动时,如何引导学生巧妙设计斜面实验间接 验证; 2、引导学生在实验过程中怎样进行合理猜想、数学推理、合理外推等重要方法。 教学过程 (一)预习检查、总结疑惑

欧拉方程

泛函的欧拉方程(by zhengpin1390) (二)、泛函的欧拉方程 欧拉方程是泛函极值条件的微分表达式,求解泛函的欧拉方程,即可得到使泛函取极值的驻函数,将变分问题转化为微分问题。 (1)最简单的欧拉方程: 设函数F(x,y,y') 是三个变量的连续函数,且点(x,y)位于有界闭区域B 内,则对形如 的变分,若其满足以下条件: c) 在有界闭区域B内存在某条特定曲线y。(x) ,使泛函取极值,且此曲线具有二阶连续导数。 则函数y。(x) 满足微分方程: 上式即为泛函Q[y]的欧拉方程。 (2)含有自变函数高阶倒数的泛函的欧拉方程 一般来说,对于下述泛函: 在类似条件下,可以得到对应的欧拉方程为: (3)含有多个自变函数的泛函的欧拉方程

对于下述泛函: 其欧拉方程组为: (4)多元函数的泛函及其欧拉方程 此处仅考虑二元函数的情况,对如下所示多元函数的泛函: 其欧拉方程为: 泛函分析 泛函分析是研究拓扑线性空间到拓扑线性空间之间满足各种拓扑和 代数条件的映射的分支学科。它是20世纪30年代形成的。从变分法、微分方程、积分方程、函数论以及量子物理等的研究中发展起来的,它运用几何学、代数学的观点和方法研究分析学的课题,可看作无限维的分析学。 泛函分析的产生 十九世纪以来,数学的发展进入了一个新的阶段。这就是,由于对欧几里得第五公设的研究,引出了非欧几何这门新的学科;对于代数方程求解的一般思考,最后建立并发展了群论;对数学分析的研究又建立了集合论。这些新的理论都为用统一的观点把古典分析的基本概念和方法一般化准备了条件。

本世纪初,瑞典数学家弗列特荷姆和法国数学家阿达玛发表的著作中,出现了把分析学一般化的萌芽。随后,希尔伯特和海令哲来创了“希尔伯特空间”的研究。到了二十年代,在数学界已经逐渐形成了一般分析学,也就是泛函分析的基本概念。 由于分析学中许多新部门的形成,揭示出分析、代数、集合的许多概念和方法常常存在相似的地方。比如,代数方程求根和微分方程求解都可以应用逐次逼近法,并且解的存在和唯一性条件也极其相似。这种相似在积分方程论中表现得就更为突出了。泛函分析的产生正是和这种情况有关,有些乍看起来很不相干的东西,都存在着类似的地方。因此它启发人们从这些类似的东西中探寻一般的真正属于本质的东西。 非欧几何的确立拓广了人们对空间的认知,n维空间几何的产生允许我们把多变函数用几何学的语言解释成多维空间的影响。这样,就显示出了分析和几何之间的相似的地方,同时存在着把分析几何化的一种可能性。这种可能性要求把几何概念进一步推广,以至最后把欧氏空间扩充成无穷维数的空间。 这时候,函数概念被赋予了更为一般的意义,古典分析中的函数概念是指两个数集之间所建立的一种对应关系。现代数学的发展却是要求建立两个任意集合之间的某种对应关系。 这里我们先介绍一下算子的概念。算子也叫算符,在数学上,把无限维空间到无限维空间的变换叫做算子。

第8章 常微分方程数值解法 本章主要内容: 1.欧拉法

第8章 常微分方程数值解法 本章主要内容: 1.欧拉法、改进欧拉法. 2.龙格-库塔法。 3.单步法的收敛性与稳定性。 重点、难点 一、微分方程的数值解法 在工程技术或自然科学中,我们会遇到的许多微分方程的问题,而我们只能对其中具有较简单形式的微分方程才能够求出它们的精确解。对于大量的微分方程问题我们需要考虑求它们的满足一定精度要求的近似解的方法,称为微分方程的数值解法。本章我们主要 讨论常微分方程初值问题?????==00 )() ,(y x y y x f dx dy 的数值解法。 数值解法的基本思想是:在常微分方程初值问题解的存在区间[a,b]内,取n+1个节点a=x 0<x 1<…<x N =b (其中差h n = x n –x n-1称为步长,一般取h 为常数,即等步长),在这些节点上把常微分方程的初值问题离散化为差分方程的相应问题,再求出这些点的上的差分方程值作为相应的微分方程的近似值(满足精度要求)。 二、欧拉法与改进欧拉法 欧拉法与改进欧拉法是用数值积分方法对微分方程进行离散化的一种方法。 将常微分方程),(y x f y ='变为() *+=?++1 1))(,()()(n x n x n n dt t y t f x y x y 1.欧拉法(欧拉折线法) 欧拉法是求解常微分方程初值问题的一种最简单的数值解法。 欧拉法的基本思想:用左矩阵公式计算(*)式右端积分,则得欧拉法的计算公式为:N a b h N n y x hf y y n n n n -= -=+=+)1,...,1,0(),(1 欧拉法局部截断误差 11121 )(2 ++++≤≤''=n n n n n x x y h R ξξ或简记为O (h 2)。

伽利略对落体运动的研究

伽利略对落体运动的研究 秦岭中学新课程高中物理导学案 班级姓名小组上课时间:年月日 课题§2.1伽利略对落体运动的研究(阅读课型) 学前 预习 目标 提示学习 目标1、知道伽利略研究落体运动的艰难过程,理解伽利略对物体运动的研究方法在科学发展和人类进步上的重大意义,学会用实验加科学推理方法研究问题。 2、通过查阅资料,结合小组内合作探究,认识研究自然规律的科学方法。 3、通过合作探究和成果分享,体验学习乐趣。 重点 提示本节课重点是伽利略对落体运动研究的思路及方法。 难点 突破通过老师讲解和小组讨论,分步讨论伽利略研究的科学推理方法。 教法

选用小组学习讨论为主,教师讲解为辅 学法 推荐自主查阅资料,按步探究;课堂小组交流,归纳总结。 学生预习案(自主完成部分) 课本 知识 预习一类 知识 二类 知识 三类 知识 自主 学习 探究 【探究要点】 1、查阅相关资料后回答什么是“佯谬”?“伽利略佯谬”巧妙否定亚里士多德论断的对你有什么启发? 佯谬就是指看上去是一个错误,但实际上是正确的。 通过逻辑推理,让自相矛盾的结果来证明理论的错误之处,是辩论时常常应用的有力武器。

自主 学习 探究 2.在对落体运动进行研究的过程中,伽利略成功地 使用了突出主要因素,忽略和排除次要因素的科学研究 方法。请你谈一谈他是突出了哪种主要因素,又是忽略 和排除了哪些次要因素的。 主要因素:物体质量对速度的影响; 次要因素:空气阻力、物体的形状对速度的影响 3.在研究落体运动时,伽利略为什么要设计那个著 名的“斜面实验”? 为了克服落体运动中物体下落快,不便测量的困难。伽利略对落体运动的研究思路可以如何概括?你从中可 以得到何种启示。 研究思路概括:问题→猜想→数学推理→实验验证 →合理外推→得出结论。 给出了我们以后研究问题的基本方法,给出了进行 科学探究的一般方法和思路。设计一个表格对课本P46 提供的伽利略手稿中记录的一组实验数据进行分析,看 看你能否得出s∝t2的结论。 t时间单位:12距离单位:32130298526824119216002104

欧拉方程的求解

欧拉方程的求解 1、引言 在数学研究领域,我们经常会瞧到以数学家名字命名的概念、公式、定理等等,让人敬佩跟羡慕、但就是,迄今为止,哪位数学家的名字出现得最多呢?她就就是数学史上与阿基米德、牛顿、高斯齐名的“四杰”之一,人称“分析学的化身”的盲人数学家欧拉(Leonhard Euler,1707--1783)、 几乎在每一个数学领域都可以瞧到她的名字,譬如我们熟悉的“欧拉线”、“欧拉圆”、“欧拉公式”、“欧拉定理”、“欧拉函数”、“欧拉积分”、“欧拉变换”、“欧拉常数”L L 欧拉还就是许多数学符号的发明者,例如用π表示圆周率、e 表示自然对数的底、()f x 表示函数、∑表示求与、i 表示虚数单位L L 以欧拉命名的数学名词有很多,本文主要讲解以欧拉命名的方程即“欧拉方程”、 在文献[1]中,关于欧拉方程的求解通常采用的就是变量变换的方法、变量变换法就就是将所求的欧拉方程化为常系数齐次线性微分方程,然后再来求解这个常系数齐次线性微分方程的解,亦即求其形如K y x =的解,进而求得欧拉方程的解、 但有些欧拉方程在用变量变换法求解时比较困难、本文在所学的欧拉方程的求解的基础上,对欧拉方程进行了简单的分类,并针对不同阶的欧拉方程的求解给出了不同的定理、最后在每类欧拉方程后面给出了典型的例题加以说明、 2、几类欧拉方程的求解 定义1 形状为 ()1(1)110n n n n n n y a x y a xy a y x ---'++++=L (1) 的方程称为欧拉方程、 (其中1a ,2a ,L ,1n a -,n a 为常数)

2、1二阶齐次欧拉方程的求解(求形如K y x =的解) 二阶齐次欧拉方程: 2120x y a xy a y '''++=、 (2) (其中1a ,2a 为已知常数) 我们注意到,方程(2)的左边y ''、y '与y 的系数都就是幂函数(分别就是 2x 、1a x 与02a x ),且其次依次降低一次、 所以根据幂函数求导的性质,我们用幂函数K y x =来尝试,瞧能否选取适当的常数K ,使得K y x =满足方程(2)、 对K y x =求一、二阶导数,并带入方程(2),得 212()0K K K K K x a Kx a x -++= 或 212[(1)]0K K a K a x +-+=, 消去K x ,有 212(1)0K a K a +-+=、 (3) 定义 2 以K 为未知数的一元二次方程(3)称为二阶齐次欧拉方程(2)的特征方程、 由此可见,只要常数K 满足特征方程(3),则幂函数K y x =就就是方程(2)的解、 于就是,对于方程(2)的通解,我们有如下结论: 定理1 方程(2)的通解为 (i) 1112ln K K y c x c x x =+, (12K K =就是方程(3)的相等的实根) (ii)1212K K x c x y c +=, (12K K ≠就是方程(3)的不等的实根) (iii)12cos(ln )sin(ln )x x c x x y c ααββ+=、(1,2K i αβ=±就是方程(3)的一对

微分方程数值解

微 分方程数值解及其应用 绪论 自然界中的许多事物的运动和变化规律都可以用微分方程来描述,因此对工程和科学技术中的实际问题的研究中, 常常需要求解微分方程.但往往只有少数较简单和典型的微分方程可求出其解析解,在大多数情况下,只能用近似法求解,数值解法是一类重要的近似方法.本文主要讨论一阶常微分方程的初值问题的数值解法,探讨这些算法在处理来自生活实际问题中的应用,并结合MATLAB 软件,动手编程予以解决. 1 微分方程的初值问题[1] 1.1 预备知识 在对生活实际问题的研究中,通常需要考虑一阶微分方程的初值问题 00(,)()dy f x y dx y x y ?=???=? (1) 这里(),f x y 是矩形区域R :00,x x a y y b -≤-≤上的连续函数. 对初值问题(1)需要考虑以下问题:方程是否一定有解呢?若有解,有多少个解呢?下面给出相关的概念与定理. 定义1 Lipschitz 条件[1][2]:矩形区域R :00,x x a y y b -≤-≤上的连续函数(),f x y 若满足:存在常数0L >,使得不等式()()1212,,f x y f x y L y y -≤-对所有()()12,,,x y x y R ∈都成立,则称(),f x y 在R 上关于y 满足Lipschitz 条件. 定理 1 解的存在唯一性定理[1][3]:设f 在区域()}{,,D x y a x b y R =≤≤∈上连续,关于y 满足Lipschitz 条件,则对任意的[]00,,∈∈x a b y R ,常微分方程初值问题(1)当[],x a b ∈时存在唯一的连续解()y x . 该定理保证若一个函数(),f x y 关于y 满足Lipschitz 条件,它所对应的微分方程的初值问题就有唯一解.在解的存在唯一性得到保证的前提下,自然要考虑方程的求

欧拉方程的求解

欧拉方程的求解 1. 引言 在数学研究领域,我们经常会看到以数学家名字命名的概念、公式、定理等等,让人敬佩跟羡慕. 但是,迄今为止,哪位数学家的名字出现得最多呢?他就是数学史上与阿基米德、牛顿、高斯齐名的“四杰”之一,人称“分析学的化身”的盲人数学家欧拉( Leonhard Euler,1707--1783 ) . 几乎在每一个数学领域都可以看到他的名字,譬如我们熟悉的“欧拉线”、“欧拉圆”、“欧拉公式”、“欧拉定理”、“欧拉函数”、“欧拉积分”、“欧拉变换”、“欧拉常数” L L 欧拉还是许多数学符号的发明者,例如用表示圆周率、e表示自然对数的底、f(x)表示函数、表示求和、i表示虚数单位L L 以欧拉命名的数学名词有很多,本文主要讲解以欧拉命名的方程即“欧拉方程”. 在文献[1] 中,关于欧拉方程的求解通常采用的是变量变换的方法. 变量变换法就是将所求的欧拉方程化为常系数齐次线性微分方程,然后再来求解这个常系数齐次线性微分方程的解,亦即求其形如y x K的解,进而求得欧拉方程的解. 但有些欧拉方程在用变量变换法求解时比较困难. 本文在所学的欧拉方程的求解的基础上,对欧拉方程进行了简单的分类,并针对不同阶的欧拉方程的求解给出了不同的定理. 最后在每类欧拉方程后面给出了典型的例题加以说明. 2. 几类欧拉方程的求解 定义 1 形状为 n (n) n 1 ( n 1) n y(n)a1x n 1y(n 1)L a n 1xy a n y 0 (1) x 的方程称为欧拉方程. (其中a i, a2, L , a ni, a.为常数)

2.1 二阶齐次欧拉方程的求解(求形如 y x K 的解) 二阶齐次欧拉方程: x 2y a i xy a 2y 0. ( 其中 a 1, a 2 为已知常数) 我们注意到,方程(2)的左边y 、y 和y 的系数都是幕函数(分别是x 2 a i x 和a 2X °),且其次依次降低一次.所以根据幕函数求导的性质,我们用幕 函数y x K 来尝试,看能否选取适当的常数 K ,使得y x K 满足方程(2). x K 求一、二阶导数,并带入方程(2),得 由此可见,只要常数K 满足特征方程(3),则幕函数y x K 就是方程(2) 共轭复根) (其中C i 、c 为任意常数) 证明(i )若特征方程(3)有两个相等的实根:? K 2,贝U 2) 消去 x K ,有 (K 2 [K 2 K 2 定义 2 以 K 为未知数的 的特征方程. K)X K (a 1 (a 1 KK a i Kx a 2 x 0 K i)K a 2]x K 0, 1)K a 2 0. 3) 元二次方程( 3)称为二阶齐次欧拉方程( 2) 的解. 于是,对于方程( 2)的通解, 定理 i 方程( 2)的通解为 y c i x Ki 我们有如下结论: (i) c 2X K1 ln X , (K i K 2是方程(3)的相等的实根) (ii) K 1 y c 1X 1 c2X K2 K i K 2是方程(3)的不等的实根) (iii) y c 1 X cos( ln X) c 2X sin( ln X). (K 1,2 i 是方程( 3)的一对

相关主题