二、数学建模--初等模型.pptVIP

思考题 试用方差最小模型计算表1-4给出的5单位78名代表的合理分配方案. 量纲齐次原则 等式两端的量纲一致 求摆动周期t的表达式 量纲分析~利用量纲齐次原则寻求物理量之间的关系 例：单摆运动 设物理量t,m,l,g之间有关系式。 ?1, ?2, ?3为待定系数，?为无量纲量 (1)的量纲表达式 对比 为什么假设这种形式？ x,y,z量纲单位缩小a,b,c倍 设p=f (x,y,z) 对x,y,z两组测量值x1,y1,z1， x2,y2,z2 p1=f (x1,y1,z1), p2=f (x2,y2,z2) 单摆运动中t，m，l，g的一般表达式: f (t, m, l, g)=0 y1~y4为待定常数，?为无量纲量 基本解 y =(y1,y2,y3,y4)T=(2, 0, -1, 1)T 设一组物理量 a,b1,b2,…,bs关系为a=f (b1,b2,…,bs)，其中b1,b2,…,bl (l≤s)量纲独立，而a, bl+1,…,bs与b1,b2,…,bl量纲相关，则可的以下无量纲量。 a=f (b1,b2,…,bs) 改写 不随测量单位变化而变化 随着测量单位变化而变化 (*) 但（*）式左端不变 ?定理(Buckingham) 设 f (q1,q2,…,qm)=0是与量纲单位无关的物理定律，X1, X2,…Xn是基本量纲，m?n, q1,q2,…,qm的量纲可表示为: 量纲矩阵记作: A=(aij)n×m 若rank A=r 线性齐次方程组 Ay =0 有m-r个基本解，记作 y=(ys1,ys2,???,ysm)T， s=1,2, ???,m-r F(π1,π2, ???,πm-r)=0 与 f (q1,q2,…,qm)=0 等价, F未定 航船速度v，船体尺寸l，浸没面积s，海水密度?，重力加速度g。 量纲分析示例：波浪对航船的阻力 航船阻力f f (q1,q2,…,qm)=0 ?(g, l, ?, v, s, f )=0 [g]=LT-2, [l]=T, [?]=ML-3, [v]=LT-1, [s]=L2, [f]=LMT-2 A=(aij)n×m m=6, n=3 f (q1,q2,…,qm)=0 ?(g, l, ?, v, s, f )=0 rank A=r rank A=3 Ay =0有m-r个基本解 Ay =0有3个基本解 y=(ys1,ys2,???,ysm)T s=1,2, ???,m-r F(?1, ?2, ???, ?m-r)=0 与 f (q1,q2,…,qm)=0 等价 F(?1, ?2, ?3)=0 与 ?(g, l, ?, v, s, f )=0等价 为得到阻力f 的显式表达式 F=0 ?3 =?(?1, ?2) ?未定 量纲分析法的评注 ? 物理量的选取 ?(???)=0中包括哪些物理量是至关重要的 ? 基本量纲的选取 基本量纲个数n;选哪些基本量纲 ? 基本解的构造 有目的地构造Ay=0的基本解 函数F和无量纲量未定 ? 方法的普适性 不需要特定的专业知识 ? 结果的局限性 2.9.2 量纲分析在物理模拟中的应用 例：航船阻力的物理模拟 通过航船模型确定原型船所受阻力 可得原 型船所 受阻力 已知模 型船所 受阻力 f1, s1, v1, l1, ?1, g1 ~原型船的参数 ~模型船的参数(均已知) f, s, v, l, ?, g (f1未知，其他已知) 注意：二者的?相同 g=g1 按一定尺寸比例造模型船，量测f，可算出 f1~物理模型 2.9.3 无量纲化 例：火箭发射 星球表面竖直发射。初速v，星球半径r，表面重力加速度g。 r 0 m2 g v x m1 研究火箭高度x随时间t的变化规律 t=0时,火箭质量m1，星球质量m2 牛顿第二定律，万有引力定律 ——3个独立参数 x=x(t; r, v, g) 用无量纲化方法减少独立参数个数 变量x, t和独立参数r, v, g的量纲 [x]=L, [t]=T, [r]=L, [v]=LT-1, [g]=LT-2 用参数r,v,g的组合,分别构造与x,t具有相同量纲的xc,tc（特征尺度） ——无量纲变量 如 xc=r, tc=r/v x=x(t; r, v, g)将被简化 xc,tc的不同构造 x=x(t; r, v, g)的不同的简化结果 x=x(t; r, v, g) ?为无量纲量 x=x(t; r, v, g) ?为无量纲量 x=x(t; r, v, g) ?为无量纲量 1）2）3） 的共同点 只含有一个参数——无量纲量? 考察无量纲量 重要差别 在1）2）3）中能否忽略以?为因子的项？ 忽略?项 不能忽略?项 忽略?项 不能忽略?项 忽略?项 原问题 火箭