边际、弹性分析(经济数学建模课件(西安交通大学,戴雪峰)

边际、弹性分析(经济数学建模课件(西安交通大学,戴雪峰)
边际、弹性分析(经济数学建模课件(西安交通大学,戴雪峰)

一、边际分析

边际的概念.

如果一个经济指标y 是另一个经济指标x 的函数)(x f y =,那么当自变量有改变量x ?时,对应有函数的改变量y ?.在经济学中,当自变量在x 处有一个单位改变量时,所对应的函数改变量为该函数所表示的经济指标在x 处的边际量.例如当生产量在x 单位水平时的边际成本,就是在已生产x 单位产品水平上,再多生产一个单位产品时总成本的改变量,或者可以说是再多生产一个单位产品所花费的成本.

设x 的改变量为x ?时,经济变量y 的改变量为y ?=)()(x f x x f -?+,则相应于x ?,y 的平均变化率是

x

x f x x f x y ?-?+=??)()( 由边际的概念,在上式中取1=?x 或1-=?x 就可得到边际量的表达式.但边际概念的定义和计算使我们想到能否用函数)(x f y =的导数作为y 的边际量呢?如果按纯粹的数学概念来讲,似乎行不通,因为导数定义要求自变量增量必须趋向于零,而实际问题中自变量x 的经济意义通常是按计件的产量或销量作为单位的,改变量为小数且趋于零不合乎实际.但我们可以这样考虑,对于现代企业来讲,其产销量的数额和一个单位产品相比是一个很大数目,1个单位常常是其中微不足道的量,可以认为改变一个单位的这种增量是趋近于零的.正是这个缘故,在经济理论研究中,总是用导数

x

x f x x f x f x ?-?+='→?)()(lim

)(0 表示经济变量y 的边际量,即认为)(x f '的经济意义是自变量在x 处有单位改变量时所引起函数y 的改变数量.

1.边际成本 在经济学中,边际成本定义为产量为x 时再增加一个单位产量时所增加的成本.

成本函数的平均变化率为

x

x C x x C x C ?-?+=??)()( 它表示产量由x 变到x +x ?时,成本函数的平均改变量. 当成本函数()C x 可导时,根据导数定义,成本函数在x 处变化率为

x x C x x C x C x ?-?+='→?)()(lim

)(0 在经济上我们认为)(x C '就是边际成本.因此,边际成本)(x C '是成本函数)(x C 关于产量x 的一阶导数.,它近似等于产量为x 时再生产一个单位产品所需增加的成本,即

)()1()()(x C x C x C x C -+=?≈'

在实际问题中企业为了生产要有厂房、机械、设备等固定资产,在短期成本函数中作为

固定成本0C ,它是常数,而生产中使用劳力,原料、材料、水电等方面的投入随产量x 的变化而改变,生产的这部分成本是可变成本,以)(1x C 记,于是成本函数可表示为

)()(10x C C x C +=

此时边际成本为

)()()()(11

0x C x C C x C '='+'=' 由此,边际成本与固定成本无关,它等于边际可变成本.

在实际经济量化分析问题中,经常将产量为x 时的边际成本)(x C '和此时已花费的平均成本

x

x C )(做比较,由两者的意义知道,如果边际成本小于平均成本,则可以再增加产量以降低平均成本,反之如果边际成本大于平均成本,可以考虑削减产量以降低平均成本.由此可知,当边际成本等于平均成本时可使产品的平均成本最低. 2.边际收入和边际利润

在经济学中,边际收入定义为销量为x 时再多销售一个单位产品时所增加的收入. 设收入函数)(x R R =是可导的,收入函数的变化率是

x

x R x x R x R x ?-?+='→?)()(lim

)(0 同边际成本道理一样,我们认为)(x R '就是边际收入.因此,边际收入)(x R '是收入函数)(x R 关于产量x 的一阶导数.,

它近似等于销量为x 时再销售一个单位产品所增加(或减少)的收入.即

)()1()()(x R x R x R x R -+=?≈'

设利润函数为)(x L L =,由于利润函数是收入函数与成本函数之差,即

)()()(x C x R x L -=

则边际利润是

)()()(x C x R x L '-'='

因此,边际利润)(x L '是利润函数)(x L 关于产量x 的一阶导数,它近似等于销量为x 时再销售一个单位产品所增加(或减少)的利润.

在经济学中还经常用到边际效用,边际产量、边际劳动生产率等概念,它和边际成本、边际收入、边际利润的经济解释方法大同小异,在此不再阐述.

下面用具体例子说明边际概念在实际问题中的意义和作用.

例 1 设某企业的产品成本函数和收入函数分别为5

2003000)(2x x x C +

+=和20

350)(2

x x x R +=,其中x 为产量,单位为件,)(x C 和)(x R 的单位为千元,求: (1)边际成本、边际收入、边际利润; (2)产量20=x 时的收入和利润,并求此时的边际收入和边际利润,解释其经济意义. 解 由边际的定义有

(1)边际成本 x x C 52200)(+=' 边际收入 10350)(x x R +=' 边际利润 x x C x R x L 10

3150)()()(-='-'=' (2)当产量为20件时,其收入和利润为

702020

)20(20350)20(2

=+?=R (千元) 6070807020)20()20()20(-=-=-=C R L (千元)

其边际收入与边际利润为

35210

20350)20(=+

='R (千元/件) 144208352)20()20()20(=-='-'='C R L (千元/件) 上面计算说明,在生产20件产品的水平上,再把产品都销售的利润为负值,即发生了亏损,亏损值为60千元;而此时的边际收入较大,即生产一件产品收入为352千元,从而得利润144千元.这样以来,该企业的生产水平由20件变到21件时,就将由亏损60千元的局面转变到盈利8460144=-千元的局面,故应该再增加产量.

二、弹性分析

一个简单引例.设2

x y =,当x 由10变到11时,y 由100变到121.显然,自变量和函数的绝对改变量分别是x ?=1,y ?=21,而它们的相对改变量x x ?和y y ?分别为 x x ?=%10101= y y ?=%21100

21= 这表明,当自变量x 由10变到11的相对变动为10%时,函数y 的相对变动为21%,这时两个相对改变量的比为 1.2%

10%21==

??=x x y y

E 解释E 的意义:x =10时,当x 改变1%时,y 平均改变2.1%,我们称E 为从x =10到x =11

时函数2x y =的平均相对变化率,也称为平均意义下函数2x y =的弹性.

这个大小度量了)(x f 对x 变化反应的强烈程度.特别是在经济学中,定量描述一个经济变量对另一个经济变量变化的反应程度对科学决策至关重要.

如果极限 00000000/)(/)]()([lim /)(/lim

x x x f x f x x f x x x f y x x ?-?+=??→?→? 存在,则称此极限值为函数)(x f y =在点x 0处的点弹性,记为0

x x Ex Ey

=,

=???=→?=x y x f x Ex

Ey

x x x )(lim 0000)()(000x f x f x ' 称)()

(x f x f x Ex Ey '=为函数)(x f y =在区间Ⅰ的点弹性函数,简称弹性函数. 而称

00000/)(/)]()([/)(/x x x f x f x x f x x x f y ?-?+=?? 为函数)(x f y =在以x 0与x 0+x ?为端点的区间上的弧弹性.

弧弹性表达了函数)(x f 当自变量x 从x 0变到x 0+x ?时函数的平均相对变化率,而点弹性正是函数)(x f 在点x 0处的相对变化率.

例2 求指数函数)1,0(≠>=a a a y x 的弹性函数.

解 因为a a y x ln ='

所以

a x a

x a a y x y Ex Ey x x ln ln =?='=. 1. 需求弹性

函数的弹性表达了函数)(x f 在x 处的相对变化率,粗略来说,就是当自变量的值每改变百分之一所引起函数变化的百分数.

需求弹性就是在需求分析中经常用来测定需求对价格反应程度的一个经济指标. 设某商品的市场需求量Q 是价格p 的函数:)(p Q Q =,)(p Q 是可导函数,则称

Q Q

p p Q p Q p Ep EQ '='=)()( 为该商品的需求价格弹性,简称为需求弹性,记为p ε.

可以这样解释p ε的经济意义;当商品的价格为p 时,价格改变1%时需求量变化的百分数.

为什么不使用变化率而要使用这种相对变化率来表达价格改变对需求量的反应呢?由弹性定义看到,弹性与量纲无关,需求弹性与需求量和价格所用的计量单位无关.以对水果的需求为例,在我国将以m 公斤/元来度量,在美国将以n 公斤/美元来度量,这就无法比较两国需求对价格的反应.正因为弹性可不受计量单位的限制,所以在经济活动分析中广泛采用,除需求价格弹性,还有收入价格弹性,成本产量弹性等.

由经济理论知道,一般商品的需求函数为价格的减函数,从而0)(<'p Q ,这说明需求价格弹性p ε一般是负的.由此,当商品的价格上涨(或下跌)1%时,需求量将下跌(或上涨)约%p ε,因此在经济学中,比较商品需求弹性的大小时,是指弹性的绝对值p ε,一般在经济分析中将需求弹性记为p p εε-=. 当1=p ε时,称为单位弹性,此时商品需求量变动的百分比与价格变动的百分比相等;

当1>p ε时,称为高弹性,此时商品需求量变动的百分比高于价格变动的百分比,价格的变动对需求量的影响比较大;当1

在商品经济中,商品经营者关心的是提价(0>?p )或降价(0

设收入函数为R ,则pQ R =,此时边际收入为

Q p Q p R '+=')()1(Q Q

p Q '+

=)1(p Q ε+= (2) 当p ?很小时,有

p Q p p R R p ?+=?'≈?)1()(ε p Q p ?-=)1(ε (3) 由此可知,当1>p ε(高弹性)时,商品降价时(0?R ,即降价可使收入增加,商品提价时(0>?p ),0>?R ,即提价将使总收入减少. 当1

上述分析使我们看到,根据商品需求弹性的不同,应制定不同的价格政策,以使收入快速增长.

例3 设某种产品的需求量Q 与价格p 的关系为 p p Q )4

1(1600)(= (1)求需求弹性;

(2)当产品的价格10=p 时再增加1%,求该产品需求量变化情况.

解 (1)由需求弹性公式

'?????

??='=p p p p Q Q p )41(1600)41(1600ε p p 39.14

1ln -≈= 需求弹性为-1.39p ,说明产品价格p 增加1%时,需求量Q 将减少1.39p %.

(2)当产品价格10=p 时,有

9.131039.1-=?-=p ε

这表示价格10=p 时,价格增加1%,产品需求量将减少13.9%;如果价格降低1%,产品的需求量将增加13.9%.这也表明此商品的需求弹性是高弹性的,适当降价会使销量大增.

例4 已知某企业的产品需求弹性为2.1,如果该企业准备明年降价10%,问这种商品的销量预期会增加多少?总收益预期会增加多少?

题中价格的改变量是相对量,所以所求的销量和总收益的改变也采用相对改变量. 解 由需求函数弹性定义知,当p ?较小时

p

Q Q p dp dQ Q p p ???≈?=ε 即

p p Q Q p ?≈?ε 故当1.2=p ε,1.0-=?p p 时,有 %21)1.0(1.2=-?-≈?Q

Q 因为R =PQ ,由(3)式有

p Q p Q R R p ??-≈?)1(εp p p ?-=)1(ε 当1.2=p ε时,有

%11)1.0()1.21(=-?-≈?R R 可见,明年企业若降价10%,企业销量将增加21%,收入将增加11%.

统计西安交大期末考试试题(含答案)

西安交大统计学考试试卷 一、单项选择题(每小题2 分,共20 分) 1.在企业统计中,下列统计标志中属于数量标志的是(C) A、文化程度 B、职业 C、月工资 D、行业 2.下列属于相对数的综合指标有(B ) A、国民收入 B、人均国民收入 C、国内生产净值 D、设备台数 3.有三个企业的年利润额分别是5000 万元、8000 万元和3900 万元,则这句话中有(B)个变量? A、0 个 B、两个 C、1 个 D、3 个 4.下列变量中属于连续型变量的是(A ) A、身高 B、产品件数 C、企业人数 D、产品品种 5.下列各项中,属于时点指标的有(A ) A、库存额 B、总收入 C、平均收入 D、人均收入 6.典型调查是(B )确定调查单位的 A、随机 B、主观 C、随意 D 盲目 7.总体标准差未知时总体均值的假设检验要用到(A ): A、Z 统计量 B、t 统计量 C、统计量 D、X 统计量 8.把样本总体中全部单位数的集合称为(A ) A、样本 B、小总体 C、样本容量 D、总体容量 9.概率的取值范围是p(D ) A、大于1 B、大于-1 C、小于1 D、在0 与1 之间 10.算术平均数的离差之和等于(A ) A、零 B、1 C、-1 D、2 二、多项选择题(每小题2 分,共10 分。每题全部答对才给分,否则不计分) 1.数据的计量尺度包括(ABCD ): A、定类尺度 B、定序尺度 C、定距尺度 D、定比尺度 E、测量尺度 2.下列属于连续型变量的有(BE ): A、工人人数 B、商品销售额 C、商品库存额 D、商品库存量 E、总产值 3.测量变量离中趋势的指标有(ABE ) A、极差 B、平均差 C、几何平均数 D、众数 E、标准差 4.在工业企业的设备调查中(BDE ) A、工业企业是调查对象 B、工业企业的所有设备是调查对象 C、每台设备是 填报单位D、每台设备是调查单位E、每个工业企业是填报单位 5.下列平均数中,容易受数列中极端值影响的平均数有(ABC ) A、算术平均数 B、调和平均数 C、几何平均数 D、中位数 E、众数 三、判断题(在正确答案后写“对”,在错误答案后写“错”。每小题1 分,共10 分) 1、“性别”是品质标志。(对)

西安交大景思睿、张鸣远版《流体力学》复习资料

流体力学知识要点 第一章 流体及其主要物理性质 1. 流体的连续介质模型 a) 流体的定义:任何微小的剪切力都会导致连续变形的物质 b) 质点:含有足够多分子数,并且具有确定宏观统计特征的分子集合。 c) 连续介质模型:(欧拉)假定组成流体的最小物理实体是流体质点而不是流体分子, 即:流体是由无穷多个、无穷小的、紧密毗邻、连绵不断的流体质点所组成的一种绝无间隙的连续介质。 2. 流体的主要物理性质 a) 流体的密度:表征流体在空间某点质量的密集程度 i. 密度:'lim V V m V ('V 特征体积,此时具有统计平均特性和确定性) ii. 比容:1v b) 压缩性:当作用在一定量流体上的压强增加时,其体积将减小, 用单位压强所引起 的体积变化率表示 i. 压缩性系数b : /b dV V dp ii. 体积弹性模量E :1 /b dp Vdp E dV V dV (Pa) v dp E d (1/)(1/)/V dp Vdp dp dp m dp dV d dV d d m 对气体: (等温 E p ;等熵 E p ,一般 1.4 ) 对液体,无明确比例 可压缩流体和不可压缩流体 液体的体积弹性模量值大,液体平衡和运动的绝大多数问题可以用不 可压缩流体解决。 气体的体积弹性模量值小,气体平衡和运动的大多数问题需要按可压 缩流体来解决。 c) 流体的粘性:是流体抵抗剪切变形或相对运动的一种固有属性,表现为流体内摩擦 i. 粘性内摩擦力产生的原因: 分子间吸引力(内聚力)产生阻力 分子不规则运动的动量交换产生的阻力 ii. 牛顿粘性实验

U U F A F A h h 牛顿内摩擦定律: /U F A h (μ 动力粘性系数,Pa ·s ) du d dy dt (d dt 角变形率) iii. 粘性系数 动力粘性系数 Pa ·s 运动粘性系数 2 /m s iv. 影响粘性的因素 压强:0p p e 正相关 温度:液体温度大粘度小 气体温度大粘度大 v. 理想流体:不具有粘性(对应粘性流体,一切实际流体都具有粘性) vi. 牛顿流体:满足牛顿内摩擦定律的流体(对应非牛顿流体,不满足牛顿内摩擦 定律) 3. 作用在流体上的力 ( 表面力 质量力) a) 表面力:作用在所取的流体分离体表面上的力。即分离体以外的流体通过接触面作 用在分离体上的力(压力,粘性力) b) 质量力:外力场作用在流体质点上的非接触力,在流体质量均匀情况下又称体积力。 质量力与外力场的强度和流体的分布有关,与它周围的微元体积无关。(重力) 4. 理想流体中的压力与方向无关 a) ,,p p x y z (即理想流体中任一点流体静压强的大小与其作用的面在空间的 方位无关,只是该点坐标的函数) b) 流体静压强的方向沿作用面的内法线方向。 第二章 流体静力学 1. 流体静压强及其特性 a) 流体静压强的方向沿作用面的内法线方向。 b) 静止流体中任一点上不论来自何方的静压强均相等,所以在静止流体中流体静压强 是空间坐标的连续函数。 2. 静止流体平衡微分方程式(欧拉平衡微分方程式) a) 1 0g a p dy du dt dy dudt dt d /

西安交通大学计算方法B上机试题

1.计算以下和式:01421181 84858616n n S n n n n ∞ =?? =--- ?++++??∑ ,要求: (1)若保留11个有效数字,给出计算结果,并评价计算的算法; (2)若要保留30个有效数字,则又将如何进行计算。 (1)题目分析 该题是对无穷级数求和,因此在使用matlab 进行累加时需要一个累加的终止条件。这里令?? ? ??+-+-+-+= 681581482184161n n n n a n n ,则 ()()1.016 1 6855844864816114851384128698161 681581482184161148113811282984161111<< ? ??? ????? ??++++++???? ????? ??++++++=??? ????? ??+-+-+-+??? ????? ??+-+-+-+=+++n n n n n n n n n n n n n n n n a a n n n n n n 故近似取其误差为1+≈k a ε,并且有m -1m -111021 21 ?=?=≈+βεk a , (2)算法依据 使用matlab 编程时用digits 函数和vpa 函数来控制位数。 (3)Matlab 运行程序 %%保留11位有效数字 k1=11; s1=0;%用于存储这一步计算值 for n=0:50 a=(1/16^n)*(4/(8*n+1)-2/(8*n+4)-1/(8*n+5)-1/(8*n+6)); n1=n-1; if a<=0.5*10^(1-k1) break end end; for i=0:1:n1 t=(1/16^i)*(4/(8*i+1)-2/(8*i+4)-1/(8*i+5)-1/(8*i+6)); s1=s1+t; end s11=vpa(s1,k1); disp('保留11位有效数字的结果为:');disp(s11); disp('此时n 值为:');disp(n1); %%保留30位有效数字 clear all; k2=30;

西安交大流体力学题与答案

一、 选择题(略) 二、 判断题(略) 三、 简答题 1.等压面是水平面的条件是什么? :①连续介质 ② 同一介质 ③ 单一重力作用下. 2. 同一容器中装两种液体,且21ρρ?,在容器侧壁装了两根测压管。试问:图中所标明的测压管中液面位置对吗?为什么?

C (c) 盛有不同种类溶液的连通器 D C D 水 油 B B (b) 连通器被隔断 A A (a) 连通容器 解:不对,(右测压管液面要低一些,从点压强的大小分析) 3. 图中三种不同情况,试问:A-A 、B-B 、C-C 、D-D 中哪个是等压面?哪个不是等压面?为什么? :( a )A-A 是 (b )B-B 不是 (c )C-C 不是, D-D 是。 四、作图题(略) 五、计算题(解题思路与答案) 1. 已知某点绝对压强为80kN/m 2 ,当地大气压强p a =98kN/m 2 。试将该点绝对压强、相对压强和真空压强用水柱及水银柱表示。 解: 用水柱高表示 (1)该点绝对压强:8.16mH 2o (2)该点相对压强:-1.84mH 2o (3)该点真空压强:1.84mH 2o

用水银柱高表示 (1)该点绝对压强:599.1mm H g (2)该点相对压强:-135.4 mm H g (3)该点真空压强:135.4 mm H g 2. 一封闭水箱自由表面上气体压强p 0=25kN/m 2 ,h 1=5m ,h 2=2m 。求A 、B 两点的静水压强。 解:由压强基本公式gh p p ρ+= 0求解 A p = 7.551 mH 2o (74 kN/m 2 ) B p = 4.551 mH 2o (44.6 kN/m 2 ) 3 如图所示为一复式水银测压计,已知m 3.21=?,m 2.12=?, m 5.23=?,m 4.14=?,m 5.15 =?(改为3.5m)。试求水箱液面上的绝 对压强0p =?

2009西安交通大学高等代数考研真题

西安交通大学2009年攻读硕士学位研究生入学考试试题 科目代码:818 科目名称:高等代数 一 (20分)计算行列式: 000 00 0000 00000n D αβαβαβαβαβαβαβαβ +++=+ + 二 (20分)已知12(0,1,0),(3,2,2)T T αα==-,是线性方程组 1231231 2321341x x x x x x ax bx cx d -+=-??++=??++=? 的两个解,求此方程组的全部解. 三 (20)当t 取什么值时,下面二次型是正定的: 222123123121323(,,)42106f x x x x x x tx x x x x x =+++++ 四(15分)设3阶实对称矩阵A 有特征值1231,1λλλ=-==,A 的属于特征值-1的特征向量1(0,1,1)T ξ=,矩阵32B A A E =-+,其中E 为3阶单位阵(下同),问: (1) 1ξ是否为B 的特征向量?求B 的所有特征值和特征向量; (2) 求矩阵B . 五(15分)设,1200000,,,,00,,,00a c x W a a b c R W y x y z R c b z z ????????????????=∈=∈???????????????????????? (1) 求12W W +; (2) 记12W W W =+,试求空间3W 使得33()M R W W =⊕(其中3()M R 为实数域 上3阶矩阵全体),并说明理由. 六(15分)设向量组12,,,r ααα线性无关,而12,,,,,r αααβγ线性相关.证明:

要么β与γ中至少有一个可被12,,,r ααα线性表出,要么12,,,,r αααβ与12,,,,r αααγ等价. 七(15分)设A 为(1)n n ?+阶常数矩阵,X 为(1)n n +?阶未知数矩阵.试证明矩阵方程AX E =有解的充要条件为()r A n =. 八(10)若12,αα是数域F 上的二维线性空间2()V F 的基,σ和τ是2()V F 上的线性变换,且满足 112212121212,,(),()σαβσαβτααββτααββ==+=+-=- 试证:στ=. 九(10)设A 和B 是两个n 阶实正交矩阵,并且det()det()A B =-.证明 ()r A B n +<. 十(10分)证明A 可与一个对角矩阵相似的充要条件是:对于A 的任意特征值i λ,方程组 2()0i E A X λ-=与()0i E A X λ-= 是同解的,其中11(,,,)n n X x x x =.需要更多试题请https://www.360docs.net/doc/5010240821.html,/exam.taoba -//maths :http 高等代数试题分数分布: 行列式:20分(1); 线性方程组:35分(2); 矩阵:15分(1); 二次型:20分(1); 线性空间:15分(1); 欧几里得空间:10分(1) 线性变换:35分(3)

数学建模之回归分析法

什么是回归分析 回归分析(regression analysis)是确定两种或两种以上变数间相互依赖的定量关系的一种统计分析方法。运用十分广泛,回归分析按照涉及的自变量的多少,可分为一元回归分析和多元回归分析;按照自变量和因变量之间的关系类型,可分为线性回归分析和非线性回归分析。如果在回归分析中,只包括一个自变量和一个因变量,且二者的关系可用一条直线近似表示,这种回归分析称为一元线性回归分析。如果回归分析中包括两个或两个以上的自变量,且因变量和自变量之间是线性关系,则称为多元线性回归分析。 回归分析之一多元线性回归模型案例解析 多元线性回归,主要是研究一个因变量与多个自变量之间的相关关系,跟一元回归原理差不多,区别在于影响因素(自变量)更多些而已,例如:一元线性回归方程为: 毫无疑问,多元线性回归方程应该为: 上图中的x1, x2, xp分别代表“自变量”Xp截止,代表有P个自变量,如果有“N组样本,那么这个多元线性回归,将会组成一个矩阵,如下图所示: 那么,多元线性回归方程矩阵形式为: 其中:代表随机误差,其中随机误差分为:可解释的误差和不可解释的误差,随机误差必须满足以下四个条件,多元线性方程才有意义(一元线性方程也一样) 1:服成正太分布,即指:随机误差必须是服成正太分别的随机变量。 2:无偏性假设,即指:期望值为0 3:同共方差性假设,即指,所有的随机误差变量方差都相等 4:独立性假设,即指:所有的随机误差变量都相互独立,可以用协方差解释。

今天跟大家一起讨论一下,SPSS---多元线性回归的具体操作过程,下面以教程教程数据为例,分析汽车特征与汽车销售量之间的关系。通过分析汽车特征跟汽车销售量的关系,建立拟合多元线性回归模型。数据如下图所示:(数据可以先用excel建立再通过spss打开) 点击“分析”——回归——线性——进入如下图所示的界面:

西安交通大学-流体力学-期末总复习

第一章流体及其主要物理性质 一概念 流体:在任何微小剪切力持续作用下连续变形的物体。 连续介质模型:假定组成流体的最小物质实体是流体质点,流体是由无限多个流体质点组成,质点之间不存在间隙。 适用条件:分子平均自由程远小于流动问题特征尺寸。 不适用条件:稀薄气体,激波层内等。 粘性:流体抵抗剪切变形(相对运动)的一种属性。 流体层间无相对运动时不表现粘性 牛顿平板实验: 两板间流体速度: 剪切力,即: 或: 式中与板间流体的种类、流体的温度、压强有关,成为液体的动力粘性系数,简称粘性系数。 流体做任意层状流动: 牛顿内摩擦定律的数学表达式 式中是角变形率或角变形速度。 凡符合牛顿内摩擦定律的流体称为牛顿流体。 产生粘性机理: 同一种流体的动力粘性系数与流体的温度有很大的关系,而受压强的影响很小。液体与气体的产生粘性机理不一样,液体的粘性主要取决于分子间的距离和分子间的吸引力,故温度升

高粘性下降;气体的粘性主要取决于分子气体热运动所产生的动量交换,故温度升高,其粘性增大。 在国际单位中: μ反应流体真实粘性的大小 运动粘性系数: 物理单位是: 粘性系数等于零的流体称为理想流体或无粘流体。 工程上常用体积弹性模量衡量流体的可压缩性,体积弹性模量定义为: 体积弹性模量的量纲和压强相同,是或。 流体的压缩性越大,则越大,即越小;反之,可压缩性越小,则越 小,即越大。 体积弹性模量又可以表示为: 等温体积弹性模量: 等熵体积弹性模量: 由: 当很小或者很大,由或者二者兼得,则此时流体的密度相对变化量就很小。如果忽 略流体密度的变化,不考虑流体的可压缩性的影响,这种简化的模型称为不可压缩流体,其密度可视为常量;反之,考虑密度为变量或压缩性影响的流体,称为可压缩流体。 不可压缩流体:液体低速流动的气体 可压缩流体:气体水下爆炸和水锤现象的液体

回归分析在数学建模中的应用

摘要 回归分析和方差分析是探究和处理相关关系的两个重要的分支,其中回归分析方法是预测方面最常用的数学方法,它是利用统计数据来确定变量之间的关系,并且依据这种关系来预测未来的发展趋势。本文主要介绍了一元线性回归分析方法和多元线性回归分析方法的一般思想方法和一般步骤,并且用它们来研究和分析我们在生活中常遇到的一些难以用函数形式确定的变量之间的关系。在解决的过程中,建立回归方程,再通过该回归方程进行预测。 关键词:多元线性回归分析;参数估计;F检验

回归分析在数学建模中的应用 Abstract Regression analysis and analysis of variance is the inquiry and processing of the correlation between two important branches, wherein the regression analysis method is the most commonly used mathematical prediction method, it is the use of statistical data to determine the relationship between the variables, and based on this relationship predict future trends. introduces a linear regression analysis and multiple linear regression analysis method general way of thinking and the general steps, and use them to research and analysis that we encounter in our life, are difficult to determine as a function relationship between the variables in the solving process, the regression equation is established by the regression equation to predict. Keywords:Multiple linear regression analysis; parameter estimation;inspection II

数学建模-回归分析-多元回归分析

1、 多元线性回归在回归分析中,如果有两个或两个以上的自变量,就称为 多元回归。事实上,一种现象常常是与多个因素相联系的,由多个自变量的最优组合共同来预测或估计因变量,比只用一个自变量进行预测或估计更有效,更符合实际。 在实际经济问题中,一个变量往往受到多个变量的影响。例如,家庭消费支出,除了受家庭可支配收入的影响外,还受诸如家庭所有的财富、物价水平、金融机构存款利息等多种因素的影响,表现在线性回归模型中的解释变量有多个。这样的模型被称为多元线性回归模型。(multivariable linear regression model ) 多元线性回归模型的一般形式为: 其中k 为解释变量的数目,j β (j=1,2,…,k)称为回归系数(regression coefficient)。上式也被称为总体回归函数的随机表达式。它的非随机表达式为: j β也被称为偏回归系数(partial regression coefficient)。 2、 多元线性回归计算模型 多元性回归模型的参数估计,同一元线性回归方程一样,也是在要求误差平方和(Σe)为最小的前提下,用最小二乘法或最大似然估计法求解参数。 设( 11 x , 12 x ,…, 1p x , 1 y ),…,( 1 n x , 2 n x ,…, np x , n y )是一个样本, 用最大似然估计法估计参数: 达 到最小。

把(4)式化简可得: 引入矩阵: 方程组(5)可以化简得: 可得最大似然估计值:

3、Matlab 多元线性回归的实现 多元线性回归在Matlab 中主要实现方法如下: (1)b=regress(Y, X ) 确定回归系数的点估计值 其中 (2)[b,bint,r,rint,stats]=regress(Y,X,alpha)求回归系数的点估计和区间估计、并检 验回归模型 ①bint 表示回归系数的区间估计. ②r 表示残差 ③rint 表示置信区间 ④stats 表示用于检验回归模型的统计量,有三个数值:相关系数r2、F 值、与F 对应的 概率p 说明:相关系数r2越接近1,说明回归方程越显著;F>F1-alpha(p,n-p-1) 时拒绝H0,F 越大,说明回归方程越显著;与F 对应的概率p<α 时拒绝H0,回归模型成立。 ⑤alpha 表示显著性水平(缺省时为0.05) (3)rcoplot(r,rint) 画出残差及其置信区间

高等代数与解析几何同济答案

高等代数与解析几何同济答案 【篇一:大学所有课程课后答案】 资料打开方法:按住 ctrl键,在你需要的资料上用鼠标左键单击 资料搜索方法:ctrl+f 输入关键词查找你要的资料 【数学】 o o o o o o o o o o o o o o o o o

习题答案 o o o o o o o o o o o o o o o o o o o o 【计算机/网络/信息】 o

o o o 【经济/金融/营销/管理/电子商务】 o o o o o o o o o o o o 【 o o o o o

【篇二:各门课程课后答案】 式]《会计学原理》同步练习题答案 [word格式]《成本会计》习题及答案(自学推荐,23页) [word格式]《成本会计》配套习题集参考答案 [word格式]《实用成本会计》习题答案 [word格式]《会计电算化》教材习题答案(09年) [jpg格式]会计从业《基础会计》课后答案 [word格式]《现代西方经济学(微观经济学)》笔记与课后习题详解(第3版,宋承先) [word格式]《宏观经济学》习题答案(第七版,多恩布什) [word格式]《国际贸易》课后习题答案(海闻 p.林德特王新奎) [pdf格式]《西方经济学》习题答案(第三版,高鸿业)可直接打印 [word格式]《金融工程》课后题答案(郑振龙版) [word格式]《宏观经济学》课后答案(布兰查德版) [jpg格式]《投资学》课后习题答案(英文版,牛逼版) [pdf格式]《投资学》课后习题答案(博迪,第四版) [word格式]《微观经济学》课后答案(高鸿业版) [word格式]《公司理财》课后答案(英文版,第六版) [word格式]《国际经济学》教师手册及课后习题答案(克鲁格曼,第六版)

西安交通大学计算方法B大作业

计算方法上机报告 姓名: 学号: 班级:

目录 题目一------------------------------------------------------------------------------------------ - 4 - 1.1题目内容 ---------------------------------------------------------------------------- - 4 - 1.2算法思想 ---------------------------------------------------------------------------- - 4 - 1.3Matlab源程序----------------------------------------------------------------------- - 5 - 1.4计算结果及总结 ------------------------------------------------------------------- - 5 - 题目二------------------------------------------------------------------------------------------ - 7 - 2.1题目内容 ---------------------------------------------------------------------------- - 7 - 2.2算法思想 ---------------------------------------------------------------------------- - 7 - 2.3 Matlab源程序---------------------------------------------------------------------- - 8 - 2.4计算结果及总结 ------------------------------------------------------------------- - 9 - 题目三----------------------------------------------------------------------------------------- - 11 - 3.1题目内容 --------------------------------------------------------------------------- - 11 - 3.2算法思想 --------------------------------------------------------------------------- - 11 - 3.3Matlab源程序---------------------------------------------------------------------- - 13 - 3.4计算结果及总结 ------------------------------------------------------------------ - 14 - 题目四----------------------------------------------------------------------------------------- - 15 - 4.1题目内容 --------------------------------------------------------------------------- - 15 - 4.2算法思想 --------------------------------------------------------------------------- - 15 - 4.3Matlab源程序---------------------------------------------------------------------- - 15 - 4.4计算结果及总结 ------------------------------------------------------------------ - 16 - 题目五----------------------------------------------------------------------------------------- - 18 -

数学建模多元回归模型

实习报告书 学生姓名: 学号: 学院名称: 专业名称: 实习时间: 2014年 06 月 05 日 第六次实验报告要求 实验目的: 掌握多元线性回归模型的原理,多元线性回归模型的建立、估计、检验及解释变量的增减的方法,以及运用相应的Matlab软件的函数计算。 实验内容: 已知某市粮食年销售量、常住人口、人均收入、肉、蛋、鱼的销售数据,见表1。请选择恰当的解释变量和恰当的模型,建立粮食年销售量的回归模型,并对其进行估计和检验。

表1 某市粮食年销售量、常住人口、人均收入、肉、蛋、鱼的销售数据 年份粮食年销售 量Y/万吨 常住人口 X2/万人 人均收 入X3/ 元 肉销售 量X4/万 吨 蛋销售 量X5/ 万吨 鱼虾销 售量 X6/万吨 197498.45560.20153.20 6.53 1.23 1.89 1975100.70603.11190.009.12 1.30 2.03 1976102.80668.05240.308.10 1.80 2.71 1977133.95715.47301.1210.10 2.09 3.00 1978140.13724.27361.0010.93 2.39 3.29 1979143.11736.13420.0011.85 3.90 5.24 1980146.15748.91491.7612.28 5.13 6.83 1981144.60760.32501.0013.50 5.418.36 1982148.94774.92529.2015.29 6.0910.07

1983158.55785.30552.7218.107.9712.57 1984169.68795.50771.1619.6110.1815.12 1985162.14804.80811.8017.2211.7918.25 1986170.09814.94988.4318.6011.5420.59 1987178.69828.731094.6 523.5311.6823.37 实验要求: 撰写实验报告,参考第10章中牙膏销售量,软件开发人员的薪金两个案例,写出建模过程,包括以下步骤 1.分析影响因变量Y的主要影响因素及经济意义; 影响因变量Y的主要影响因素有常住人口数量,城市中人口越多,需要的粮食数量就越多,粮食的年销售量就会相应增加。粮食销量还和人均收入有关,人均收入增加了,居民所能购买的粮食数量也会相应增加。另外,肉类销量、蛋销售量、鱼虾销售量也会对粮食的销售量有影响,这些销量增加了,也表示居民的饮食结构也在发生变化,生活水平在提高,所以相应的,生活水平提升了,居民也有能力购买更多的粮食。

西交大计算方法上机报告

计算方法(B)实验报告 姓名: 学号: 学院: 专业:

实验一 三对角方程组Tx f =的求解 一、 实验目的 掌握三对角方程组Tx f =求解的方法。 二、 实验内容 求三对角方程组Tx f =的解,其中: 4 -1 -1 4 -1 -1 4 1 -1 4T ????????=?? ?? ???? , 3223f ?? ? ? ?= ? ? ??? 三、 算法组织 设系数矩阵为三对角矩阵 11222333111 b c a b c a b c a b c b n n n n T ---???????? =?????? ?????? 则方程组Tx f =称为三对角方程组。 设矩阵T 非奇异,T 可分解为T=LU ,其中L 为下三角矩阵,U 为单位上三角矩阵,记 1 1 212 313 1 1 1111 ,11n n n n n r l r l r L U l r l μμμμμ---???? ? ? ? ? ? ?== ? ? ? ? ? ? ? ? ? ?? ? ? ? 可先依次求出,L U 中的元素后,令Ux y =,先求解下三角方程组Ly f =得出 y ,再求解上三角方程组Ux y =。 追赶法的算法组织如下: 1.输入三对角矩阵T 和右端向量f ;

2.将Tx f =压缩为四个一维数组{}{}{}{}i i i i a b c d 、、、,{}{}{}i i i a b c 、、是T 的三对角线性方程组的三个对角,{}i d 是右端向量。将分解矩阵压缩为三个一维数组 {}{}{}i i i l r μ、、。 3.对T 做Crout 分解(也可以用Doolittle 分解)导出追赶法的计算步骤如下: 1111,b r c μ== for 2i n = 111, , ,i i i i i i i i i i i i i l a b a r r c y d l y μμ---==-==- end 4.回代求解x /n n n x y μ= for 11i n =- 1()/i i i i i x y c x μ+=- end 5. 停止,输出结果。 四、 MATLAB 程序 MATLAB 程序见附件1. 五、 结果及分析 实验结果为: (1.0000 1.0000 1.0000 1.0000)T x =

数学建模统计模型

数学建模

论文题目: 一个医药公司的新药研究部门为了掌握一种新止痛剂的疗效,设计了一个药物试验,给患有同种疾病的病人使用这种新止痛剂的以下4个剂量中的某一个:2 g,5 g,7 g和10 g,并记录每个病人病痛明显减轻的时间(以分钟计). 为了解新药的疗效与病人性别和血压有什么关系,试验过程中研究人员把病人按性别及血压的低、中、高三档平均分配来进行测试. 通过比较每个病人血压的历史数据,从低到高分成3组,分别记作,和. 实验结束后,公司的记录结果见下表(性别以0表示女,1表示男). 请你为该公司建立一个数学模型,根据病人用药的剂量、性别和血压组别,预测出服药后病痛明显减轻的时间.

一、摘要 在农某医药公司为了掌握一种新止痛药的疗效,设计了一个药物实验,通过观测病人性别、血压和用药剂量与病痛时间的关系,预测服药后病痛明显减轻的时间。我们运用数学统计工具m i n i t a b软件,对用药剂量,性别和血压组别与病痛减轻

时间之间的数据进行深层次地处理并加以讨论概率值P (是否<)和拟合度R-S q的值是否更大(越大,说明模型越好)。 首先,假设用药剂量、性别和血压组别与病痛减轻时间之间具有线性关系,我们建立了模型Ⅰ。对模型Ⅰ用m i n i t a b 软件进行回归分析,结果偏差较大,说明不是单纯的线性关系,然后对不同性别分开讨论,增加血压和用药剂量的交叉项,我们在模型Ⅰ的基础上建立了模型Ⅱ,用m i n i t a b软件进行回归分析后,用药剂量对病痛减轻时间不显着,于是我们有引进了用药剂量的平方项,改进模型Ⅱ建立了模型Ⅲ,用m i n i t a b 软件进行回归分析后,结果合理。最终确定了女性病人服药后病痛减轻时间与用药剂量、性别和血压组别的关系模型: Y=1x 3x 1x 3x 2 1 x 对模型Ⅱ和模型Ⅲ关于男性病人用m i n i t a b软件进行回归分析,结果偏差依然较大,于是改进模型Ⅲ建立了模型Ⅳ,用m i n i t a b软件进行回归分析后,结果合理。最终确定了男性病人服药后病痛减轻时间与用药剂量、性别和血压组别的关系模 型:Y=1x1x 3x 2 1 x关键词止痛剂药剂量性别病痛减轻时 间

西交计算方法A上机大作业

计算方法A 上机大作业 1. 共轭梯度法求解线性方程组 算法原理:由定理3.4.1可知系数矩阵A 是对称正定矩阵的线性方程组Ax=b 的解与求解二次函数1()2 T T f x x Ax b x =-极小点具有等价性,所以可以利用共轭梯度法求解1()2 T T f x x Ax b x = -的极小点来达到求解Ax=b 的目的。 共轭梯度法在形式上具有迭代法的特征,在给定初始值情况下,根据迭代公式: (1)()()k k k k x x d α+=+ 产生的迭代序列(1)(2)(3)x x x ,,,... 在无舍入误差假定下,最多经过n 次迭代,就可求得()f x 的最小值,也就是方程Ax=b 的解。 首先导出最佳步长k α的计算式。 假设迭代点()k x 和搜索方向()k d 已经给定,便可以通过()()()() k k f x d φαα=+的极小化 ()()min ()()k k f x d φαα=+ 来求得,根据多元复合函数的求导法则得: ()()()'()()k k T k f x d d φαα=?+ 令'()0φα=,得到: ()() ()()k T k k k T k r d d Ad α=,其中()()k k r b Ax =- 然后确定搜索方向()k d 。给定初始向量(0)x 后,由于负梯度方向是函数下降最快的方向,故第一次迭代取搜索方向(0) (0)(0)(0)()d r f x b Ax ==-?=-。令 (1)(0)00x x d α=+ 其中(0)(0)0(0)(0) T T r d d Ad α=。第二次迭代时,从(1) x 出发的搜索方向不再取(1)r ,而是选取(1) (1)(0)0d r d β=+,使得(1)d 与(0)d 是关于矩阵A 的共轭向量,由此可 求得参数0β:

西安交通大学入学测试机考《高等数学一(专升本)》模拟题及答案

西安交通大学入学测试机考 专升本高数(一)模拟题1、题目Z1-1(2)() A.A B.B C.C D.D 标准答案:B 2、题目1-1(2)() A.A B.B C.C D.D 标准答案:C 3、题目1-2(2)() A.A B.B C.C D.D 标准答案:A

4、题目1-3(2)() A.A B.B C.C D.D 标准答案:B 5、题目6-1:(2)() A.A B.B C.C D.D 标准答案:B 6、题目1-4(2)() A.A B.B C.C D.D 标准答案:D 7、题目1-5(2)() A.A B.B

C.C D.D 标准答案:C 8、题目1-6(2)() A.A B.B C.C D.D 标准答案:D 9、题目1-7(2)() A.A B.B C.C D.D 标准答案:B 10、题目1-8(2)() A.A B.B C.C D.D 标准答案:A 11、题目1-9(2)()

A.A B.B C.C D.D 标准答案:B 12、题目1-10(2)() A.A B.B C.C D.D 标准答案:B 13、题目6-2:(2)() A.A B.B C.C D.D 标准答案:C 14、题目2-1(2)()

A.A B.B C.C D.D 标准答案:D 15、题目2-2(2)() A.A B.B C.C D.D 标准答案:C 16、题目2-3(2)() A.A B.B C.C D.D 标准答案:D 17、题目6-3:(2)() A.A B.B C.C D.D 标准答案:C

18、题目2-4(2)() A.A B.B C.C D.D 标准答案:B 19、题目6-4:(2)() A.A B.B C.C D.D 标准答案:C 20、题目2-5(2)() A.A B.B C.C D.D 标准答案:A 21、题目2-6(2)()

西安交通大学计算方法A实验报告

实验一 矩阵的分解 一、实验目的 掌握矩阵的分解原理和一般方法,学会利用矩阵分解直接求解线性方程组。 二、实验内容 求矩阵() 2020 =ij A α?的T LDL 分解与Cholesky 分解,其中 ,min(,),ij i i j i j i j α=?=? ≠? 。 三、问题分析 1. Cholesky 分解 Cholesky 分解是针对被分解矩阵为对称正定的情况给出的。 分解步骤如下: 11g =1111/y b g =,1111i i g g α= 2i n = ; DO 2j n = jj g = IF 0jj g < STOP ,JUMP TO (5) DO 1i j n =+ 1 1j ij ik kj k ij jj g g g g α-=??- ? ? ?=∑ ji ij g g = 1 1j i ik k k i jj b g y y g -=??- ? ? ?=∑ END DO END DO

2. T LDL 分解 T LDL 分解是针对Cholesky 分解中的开平方运算进行的改进。 分解步骤如下: 11i i r α=,1111/i i r r r =,11y b = 1i n = DO 2i n = DO j i n = 1 1i ij ij ik kj k r l r α-=??=- ??? ∑ /ji ij ii l r r = 1 1i i i ik k k y b l b -=??=- ??? ∑ END DO END DO 四、matlab 求解 分别写出T LDL 分解和Cholesky 分解的函数程序gaijinsqrt.m 和.cholesky m ,调用格 式如下: 1. [index,x,r]=gaijinsqrt(A,b) 参数说明: A 和b 分别是线性代数方程组Ax =b 的系数矩阵和右端向量;输出x 为解向量。 [index,x,g]=Cholesky(A,b) 参数说明: A 和b 分别是线性代数方程组Ax =b 的系数矩阵和右端向量;输出x 为解向量。 然后写出主程序2homework .m 如下: %生成矩阵A A=zeros(20,20); for i=1:20 for j=1:20 if i~=j if i>j A(i,j)=j; else A(i,j)=i; end

西安交通大学计算方法B大作业资料

计算方法上机报告 姓名: 学号:

班级: 目录 题目一-------------------------------------------------------------------- 4- 1.1题目内容 ---------------------------------------------------------- 4- 1.2算法思想 ----------------------------------------------------------- 4- 1.3Matlab 源程序--------------------------------------------------- 5 - 1.4计算结果及总结 ------------------------------------------------ 5- 题目二---------------------------------------------------------------- 7- 2.1题目内容 ----------------------------------------------------------- 7- 2.2算法思想 ------------------------------------------------------------ 7 2.3 Matlab 源程序 -------------------------------------------------- 8 - 2.4计算结果及总结 ------------------------------------------------ 9- 题目三------------------------------------------------------------------- 11- 3.1题目内容 --------------------------------------------------------- 11- 3.2算法思想 --------------------------------------------------------- 11- 3.3Matlab 源程序------------------------------------------------------- 13- 3.4计算结果及总结 -------------------------------------------------- 14- 题目四------------------------------------------------------------------- 15- 4.1题目内容 --------------------------------------------------------- 15- 4.2算法思想 --------------------------------------------------------- 15- 4.3Matlab 源程序------------------------------------------------------- 15- 4.4计算结果及总结 ----------------------------------------------- 16- 题目五------------------------------------------------------------------- 18- 5.1题目内容 --------------------------------------------------------- 18- 5.2算法思想 ---------------------------------------------------------- 18 5.3 Matlab 源程序 ------------------------------------------------------ 18 5.3.1 非压缩带状对角方程组----------------------------------- 18 - 5.3.2压缩带状对角方程组--------------------------------------- 20- 5.4实验结果及分析 ----------------------------------------------- 22-

相关文档
最新文档