克重计算方法
系数纬纱纱支纬密经纱纱支经密克重*??
? ??+=+树脂用量+涂胶用量
坯布系数：25.2 成品系数：24.5 树脂用量：15-20克 涂胶用量：
1、PA 胶：15-25克
2、EVA 胶：20-30克
3、LDPE ：25-30克
4、HDPE ：25-30克
纱支与克重
克重计算相关 平方克重与盎司的关系 克重--面料的克重一般为平方米面料重量的克数，克重是针织面料的一个重要的技术指标，粗纺毛呢通常也把克重作为重要的技术指标。牛仔面料的克重一般用"盎司（OZ）"来表达，即每平方码面料重量的盎司数，如7盎司、12盎司牛仔布等 1 yard 码=3 feet 英尺=0.9144 metre 米 1 ounce 盎司=16 drams 打兰=28.35 grams 克 1盎司/码=28.35克/（0.9144 米*0.9144 米）=33.9克/平方米 11盎司/码=373克/平方米 针织布类的克重与支数识别方法 支数指的是一镑重（454克）的棉纱所具有的 840码的长度的个数。支数越大，纱线越细。如48S比23细。克重指每平方米的面料所具有的重量克数。一般来说，相同品质的纱线，支数越大，克重越大。 坯布的克重与成品克重换算，如以公式计算，与实际可能有很大偏差，主要原因分两方面，胚布设定可用公式计算，只要条件愈接近，结果愈准碓，因织造为物理性改变，由胚布加工为成品布，经过染整多个工序，织物经过物理性及化学性改变，变量增大，期间涉及各种不同类型之机械，工序等等，形成各企业有着不同之规格标准，因此坯布的克重与成品克重换算，每个企业均有其标准。 坯布的克重与成品克重换算 如何建立换算公式？ 建议将企业以往的胚布至成品布之结果收集，再用计算机筛选最多之结果作为标准，再按企业生产情况，评估其实用性，再订下换算公式。 实例参考：（一般布类，一般织染整工艺） 胚重 == 订单成品重 x 胚重系数 布类胚重系数 平纹 92% 双珠 88%
双位衣 89% 单珠 93% 单位衣 89% 拉架平纹 100% 毛巾 106% 双面 88% 拉架罗纹 91% 罗纹 83% 拉架灯蕊 92% 灯蕊 87% ⑴、坯布克重：是指织物单位面积的重量，一般以每平方米的克数表示（g/m2）。商业上最简单计算方法。 例如：棉织物一般在70～250 g/m2范围内； 凡立丁（精梳）毛织物重量在185 g/m2 左右； 厚花呢（粗梳）为280g/m2左右； 薄形丝织物一般在20～100 g/m2之间。 坯布克重也是考核织物物理性能的品德指标之一，这对毛织物来说更为注重。在精梳毛织物、毛型化纤织物的《部颁》标准中规定：平方米重偏轻不能超过允许公差（规定为5%），定为一等品；超过5%而不超过7%，为二等品。 ⑵、按《部颁标准》计算方法（应扣除回潮率因素）： 式中：GK --公定回潮率下的平方米重量（g/m2） GO --试样干重（g） WK --试样的公定回潮率（%） L --试样长度（cm） B --试样宽度（cm）
数值计算方法比较
有限差分方法(FDM:Finite Difference Method)是计算机数值模拟最早采用的方法，至今仍被广泛运用。该方法将求解域划分为差分网格，用有限个网格节点代替连续的求解域。有限差分法以Taylor级数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。有限差分法主要集中在依赖于时间的问题（双曲型和抛物型方程）。有限差分法方面的经典文献有Richtmeyer & Morton的《Difference Methods for Initial-Value Problems》；R. LeVeque《Finite Difference Method for Differential Equations》；《Numerical Methods for C onservation Laws》。 注：差分格式： （1）从格式的精度来划分，有一阶格式、二阶格式和高阶格式。 （2）从差分的空间形式来考虑，可分为中心格式和逆风格式。 （3）考虑时间因子的影响，差分格式还可以分为显格式、隐格式、显隐交替格式等。 目前常见的差分格式，主要是上述几种形式的组合，不同的组合构成不同的差分格式。差分方法主要适用于有结构网格，网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。 构造差分的方法： 构造差分的方法有多种形式，目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式：一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等，其中前两种格式为一阶计算精度，后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合，可以组合成不同的差分计算格式。 有限差分法的不足：由于采用的是直交网格，因此较难适应区域形状的任意性，而且区分不出场函数在区域中的轻重缓急之差异，缺乏统一有效的处理自然边值条件和内边值条件的方法，难以构造高精度(指收敛阶)差分格式，除非允许差分方程联系更多的节点(这又进一步增加处理边值条件韵困难)。另外它还有编制不出通用程序的困难。 有限差分法的优点：该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念 直观，表达简单，精度可选而且在一个时间步内，对于一个给定点来说其相关的空间点只是 与该相邻的几点，而不是全部的空间点。是发展较早且比较成熟的数值方法 广义差分法（有限体积法）（GDM:Generalized Difference Method）：1953年，Mac—Neal 利用积分插值法(也称积分均衡法)建立了三角网格上的差分格 式，这就是以后通称的不规划网格上的差分法．这种方法的几何误差小，特别是给出了处理自然边值条件(及内边值条件)的有效方法，堪称差分法的一大进步。1978年，李荣华利用有限元空间和对偶单元上特征函数的推广——局部Taylor展式的公项，将积分插值法改写成广义Galerkin法形式，从而将不规则网格差分法推广为广义差分法.其基本思路是，将计算区域划分为一系列不重复的控制体积，并使每个网格点周围有
统计学计算题例题及计算分析
计算分析题解答参考 1.1．某厂三个车间一季度生产情况如下： 计算一季度三个车间产量平均计划完成百分比和平均单位产品成本。 解：平均计划完成百分比＝实际产量/计划产量＝733/（198/0.9+315/1.05+220/1.1） ＝101.81% 平均单位产量成本 X=∑xf/∑f=(15*198+10*315+8*220)/733 =10.75(元/件) 1.2．某企业产品的有关资料如下： 试分别计算该企业产品98年、99年的平均单位产品成本。 解：该企业98年平均单位产品成本 x=∑xf/∑f=(25*1500+28*1020+32*980)/3500 =27.83(元/件) 该企业99年平均单位产品成本x=∑xf /∑(m/x)=101060/(24500/25+28560/28+48000/32) =28.87(元/件) 年某月甲、乙两市场三种商品价格、销售量和销售额资料如下： 1.3．1999 解：三种商品在甲市场上的平均价格x=∑xf/∑f=(105*700+120*900+137*1100)/2700 =123.04(元/件) 三种商品在乙市场上的平均价格x=∑m/∑(m/x)=317900/(126000/105+96000/120+95900/137) =117.74(元/件) 2.1．某车间有甲、乙两个生产小组，甲组平均每个工人的日产量为22件，标准差为 3.5件；乙组工人日产量资料：
试比较甲、乙两生产小组中的哪个组的日产量更有代表性？ 解：∵X 甲=22件 σ甲=3.5件 ∴V 甲=σ甲/ X 甲=3.5/22=15.91% 列表计算乙组的数据资料如下: ∵x 乙=∑xf/∑f=(11*10+14*20+17*30+20*40)/100 =17(件) σ乙= √[∑(x-x)2 f]/∑f =√900/100 =3(件) ∴V 乙=σ乙/ x 乙=3/17=17.65% 由于V 甲＜V 乙,故甲生产小组的日产量更有代表性。 2.2．有甲、乙两个品种的粮食作物，经播种实验后得知甲品种的平均产量为998斤，标准差为162.7斤；乙品种实验的资料如下： 试研究两个品种的平均亩产量，确定哪一个品种具有较大稳定性，更有推广价值？ 解：∵x 甲=998斤 σ甲=162.7斤 ∴V 甲=σ甲/ x 甲=162.7/998=16.30% 列表计算乙品种的数据资料如下:
坯布克重计算方法
坯布克重计算方法 坯布规格190T涤塔夫68DX68D20筘*2穿*30梭门副165cm 以上算法只实用与不加捻产品，算出来是米克重哦．当然不加捻产品算克重都可以套用．如果是织造厂算原料的用量，在纬的门副上至少+１０ＣＭ．经不用加. 一．原料用量计算 定义：织一米坯布需经纬原料多少克。 1.长纤类：经用量（g/m）＝总经跟数×(D/9000)×1.1 纬用量（g/m）＝坯布纬密×上机门幅×(D/9000)×1.1 或＝成品纬密×成品门幅×(D/9000)×1.1 注：1.1＝1+10％，10％为织缩率+损耗，一般FDY取1.08，DTY取1.12 例1涤塔夫上机门幅168cm17号筘3穿坯布纬密28根/cm (168×17#×3入)
原料：63DFDY×63DFDY 求：原料用量 解：经用量=168×17×3×63/9000×1.08=64.77g/m 纬用量＝168×28×63/9000×1.08＝35.56g/m 即理论米克重＝64.77+35.56＝100g/m 例2. 75DFDY+100DDTY五枚缎 12＃×5入×16630根/cm 求：原料用量 解：经用量＝12×5×166×75/9000×1.08=89.64g/m 纬用量＝30×166×100/9000×1.12=61.97g/m 2.短纤类：经用量＝0.64984×（经密/经线纱支）×门幅 纬用量＝0.64984×（纬密/纬线纱支）×门幅 注：也可把短纤换算成长纤,用长纤公式来计算，D=5315/S
例1. 70D锦纶FDY×21S棉133×7259”，求：原料用量解：经用量＝133×59×70/9000×1.08＝65.91g/m 纬用量＝0.64984×72/21×59＝131.4g/m 例2.经：70D锦纶DTY+30D氨纶 纬：30S棉+40D氨纶86 ×5859” 求：原料用量及原料所占比例 解：氨纶有个拉缩比，一般30D的氨纶按10-13D计算，40D按15-18D 计算，拉缩比一般可取 3，如40/3＝13.33，按15-18算 经用量锦纶＝70×86×59/9000×1.12＝44.2 经用量氨纶＝86×59×13/9000×1.1＝8.06 纬用量棉＝58×59×5315/30/9000×1.1＝74.1 纬用量氨纶＝58×59×18/9000×1.1＝7.53 总＝44.2+8.06+74.1+7.53＝133.89g/m 氨纶比例＝（8.06+7.53）/133.89＝11.64％
数值分析(计算方法)总结
第一章绪论 误差来源：模型误差、观测误差、截断误差（方法误差）、舍入误差 是的绝对误差，是的误差，为的绝对误差限（或误差限） 为的相对误差，当较小时，令 相对误差绝对值得上限称为相对误差限记为：即： 绝对误差有量纲，而相对误差无量纲 若近似值的绝对误差限为某一位上的半个单位，且该位直到的第一位非零数字共有n位，则称近似值有n位有效数字，或说精确到该位。 例：设x==…那么,则有效数字为1位，即个位上的3，或说精确到个位。 科学计数法：记有n位有效数字，精确到。 由有效数字求相对误差限：设近似值有n位有效数字，则其相对误差限为 由相对误差限求有效数字：设近似值的相对误差限为为则它有n位有效数字 令 1.x+y近似值为和的误差（限）等于误差（限）的 和 2.x-y近似值为 3.xy近似值为 4. 1．避免两相近数相减 2．避免用绝对值很小的数作除数 3．避免大数吃小数 4．尽量减少计算工作量 第二章非线性方程求根 1.逐步搜索法 设f (a) <0, f (b)> 0，有根区间为(a, b)，从x0=a出发，按某个预定步长(例如h=(b-a)/N)
一步一步向右跨，每跨一步进行一次根的搜索，即判别f(x k)=f(a+kh)的符号，若f(x k)>0(而 f(x k-1)<0),则有根区间缩小为[x k-1,x k] (若f(x k)=0，x k即为所求根), 然后从x k-1出发，把搜索步长再缩小，重复上面步骤，直到满足精度：|x k-x k-1|
地方时计算方法及试题精选(DOC)
关于地方时的计算 一．地方时计算的一般步骤： 1．找两地的经度差： (1)如果已知地和要求地同在东经或同在西经，则： 经度差=经度大的度数—经度小的度数 (2)如果已知地和要求地不同是东经或西经，则： 经度差=两经度和（和小于180°时） 或经度差=（180°—两经度和）。（在两经度和大于180°时） 2．把经度差转化为地方时差，即： 地方时差=经度差÷15°/H 3．根据要求地在已知地的东西位置关系，加减地方时差，即：要求点在已知点的东方，加地方时差；如要求点在已知点西方，则减地方时差。 二．东西位置关系的判断： （1）同是东经，度数越大越靠东。即：度数大的在东。 （2）是西经，度数越大越靠西。即：度数大的在西。 （3）一个东经一个西经，如果和小180°，东经在东西经在西；如果和大于180°，则经度差=（360°—和），东经在西，西经在东；如果和等于180，则亦东亦西。 三．应用举例： 1、固定点计算 【例1】两地同在东经或西经 已知：A点120°E,地方时为10：00，求B点60°E的地方时。 分析：因为A、B两点同是东经，所以，A、B两点的经度差=120°－60°=60° 地方时差=60°÷15°/H=4小时 因为A、B两点同是东经，度数越大越靠东，要求B点60°E比A点120°E小，所以，B点在A点的西方，应减地方时差。 所以，B点地方时为10：00—4小时=6：00 【例2】两地分属东西经 A、已知:A点110°E的地方时为10:00,求B点30°W的地方时. 分析:Ａ在东经，Ｂ在西经，110°+30°=140°＜180°，所以经度差=140°，且A点东经在东，B 点西经在西，A、B两点的地方时差=140°÷15°/H=9小时20分，B点在西方， 所以，B点的地方时为10：00—9小时20分=00：40。 B、已知A点100°E的地方时为8：00，求B点90°W的地方时。 分析：A点为东经，B点为西经，100°+90°=190°＞180°， 则A、，B两点的经度差=360°—190°=170°，且A点东经在西，B点西经在东。 所以，A、B两点的地方时差=170°÷15°/H=11小时20分，B点在A点的东方， 所以B点的地方时为8：00+11小时20分=19：20。 C、已知A点100°E的地方8：00，求B点80°W的地方时。 分析：A点为100°E，B点为80°W，则100°+80°=180°，亦东亦西，即：可以说B点在A 点的东方，也可以说B点在A点的西方，A，B两点的地方时差为180÷15/H=12小时。 所以B点的地方时为8：00+12小时=20：00或8：00—12小时，不够减，在日期中借一天24小时来，即24小时+8：00—12小时=20：00。 2、变化点计算 【例1】一架飞机于10月1日17时从我国上海（东八区）飞往美国旧金山（西八区），需飞行14小时。到达目的地时，当地时间是（） A. 10月2日15时 B. 10月2日3时 C. 10月1日15时 D. 10月1日3时
尺寸链计算方法
第十章装配精度与加工精度分析任何机械产品及其零部件的设计，都必须满足使用要求所限定的设计指标，如传动关系、几何结构及承载能力等等。此外，还必须进行几何精度设计。几何精度设计就是在充分考虑产品的装配技术要求与零件加工工艺要求的前提下，合理地确定零件的几何量公差。这样，产品才能获得尽可能高的性能价格比，创造出最佳的经济效益。进行装配精度与加工精度分析以及它们之间关系的分析，可以运用尺寸链原理及计算方法。我国业已发布这方面的国家标准GB5847—86《尺寸链计算方法》，供设计时参考使用。 第一节尺寸链的基本概念 一、有关尺寸链的术语及定义 1．尺寸链 在机器装配或零件加工过程中，由相互连接的尺寸形成的封闭尺寸组，称为尺寸链。尺寸链分为装配尺寸链和工艺尺寸链两种形式。 （a）齿轮部件（b）尺寸链图（c）尺寸链图 图10-1 装配尺寸链示例 图10-1a为某齿轮部件图。齿轮3在位置固定的轴1上回转。按装配技术规范，齿轮左右端面与挡环2和4之间应有间隙。现将此间隙集中于齿轮右端面与挡环4左端面之间，用符号A0表示。装配后，由齿轮3的宽度A1、挡环2的宽度A2、轴上轴肩到轴槽右侧面的距离A3、弹簧卡环5的宽度A4及挡环4的宽度A5、间隙A0依次相互连接，构成封闭尺寸组，形成一个尺寸链。这个尺寸链可表示为图10-1b与图10-1c两种形式。上述尺寸链由不同零件的设计尺寸所形成，称为装配尺寸链。 图10-2a为某轴零件图（局部）。该图上标注轴径B1与键槽深度B2。键槽加工顺序如图10-2b所示：车削轴外圆到尺寸C1，铣键槽深度到尺寸C2，磨削轴外圆到尺寸C3（即图10-2a中的尺寸B1），要求磨削后自然形成尺寸C0（即图10-2a 中的键槽深度尺寸B2）。在这个过程中，加工尺寸C1、C2、C3和完工后尺寸C0构成封闭尺寸组，形成一个尺寸链。该尺寸链由同一零件的几个工艺尺寸构成，称为工艺尺寸链。
数值计算方法大作业
目录 第一章非线性方程求根 (3) 1.1迭代法 (3) 1.2牛顿法 (4) 1.3弦截法 (5) 1.4二分法 (6) 第二章插值 (7) 2.1线性插值 (7) 2.2二次插值 (8) 2.3拉格朗日插值 (9) 2.4分段线性插值 (10) 2.5分段二次插值 (11) 第三章数值积分 (13) 3.1复化矩形积分法 (13) 3.2复化梯形积分法 (14) 3.3辛普森积分法 (15) 3.4变步长梯形积分法 (16) 第四章线性方程组数值法 (17) 4.1约当消去法 (17) 4.2高斯消去法 (18) 4.3三角分解法 (20)
4.4雅可比迭代法 (21) 4.5高斯—赛德尔迭代法 (23) 第五章常积分方程数值法 (25) 5.1显示欧拉公式法 (25) 5.2欧拉公式预测校正法 (26) 5.3改进欧拉公式法 (27) 5.4四阶龙格—库塔法 (28)
数值计算方法 第一章非线性方程求根 1.1迭代法 程序代码： Private Sub Command1_Click() x0 = Val(InputBox("请输入初始值x0")) ep = Val(InputBox(请输入误差限ep)) f = 0 While f = 0 X1 = (Exp(2 * x0) - x0) / 5 If Abs(X1 - x0) < ep Then Print X1 f = 1 Else x0 = X1 End If Wend End Sub 例：求f(x)=e2x-6x=0在x=0.5附近的根（ep=10-10）
1.2牛顿法 程序代码： Private Sub Command1_Click() b = Val(InputBox("请输入被开方数x0")) ep = Val(InputBox(请输入误差限ep)) f = 0 While f = 0 X1 = x0 - (x0 ^ 2 - b) / (2 * b) If Abs(X1 - x0) < ep Then Print X1 f = 1 Else x0 = X1 End If Wend End Sub 例：求56的值。（ep=10-10）
盈亏问题计算公式+例题分析(打印版)
数学运算:盈亏问题计算公式 把若干物体平均分给一定数量得对象，并不就是每次都能正好分完。 如果物体还有剩余，就叫盈; 如果物体不够分，就叫亏。 凡就是研究盈与亏这一类算法得应用题就叫盈亏问题。 盈亏问题得常见题型为给出某物体得两种分配标准与结果，来求物体数量与参与分配得对象数量。由于每次分配都可能出现刚好分完、多余或不足这三种情况，那么就会有多种结果得组合，这里以一道典型得盈亏问题对三种情况得几种组合加以说明。 注意：公司中两次每人分配数得差也就就是大分减小分 一、基础盈亏问题 1、一盈一亏（不够）【一次有余（盈），一次不够（亏）】可用公式：（盈+亏）÷（两次每人分配数得差）=人数。例如，“小朋友分桃子，每人10个少9个，每人8个多7个。问：有多少个小朋友与多少个桃子？” 解：（7+9）÷（10-8）=16÷2=8（个）………………人数 10×8-9=80-9=71（个）………………………桃子 或8×8+7=64+7=71（个）（答略) 测试：如果每人分9 个苹果，就剩下10 个苹果;如果每人分12 个苹果，就少20 个苹果。 2、两次皆盈（余），可用公式：（大盈-小盈）÷（两次每人分配数得差）=人数。 例如，“士兵背子弹作行军训练，每人背45发，多680发；若每人背50发，则还多200发。问：有士兵多少人？有子弹多少发？” 解：（680-200）÷（50-45）=480÷5=96（人） 45×96+680=5000（发）或50×96+200=5000（发）（答略） 测试：如果每人分8 个苹果，就剩下20 个苹果;如果每人分7 个苹果，就剩下30 个苹果。 3、两次皆亏（不够），可用公式：（大亏-小亏）÷（两次每人分配数得差）=人数。 例如，“将一批本子发给学生，每人发10本，差90本；若每人发8本，则仍差8本。有多少学生与多少本本子？”解：（90-8）÷（10-8）=82÷2=41（人）10×41-90=320（本）（答略） 测试：如果每人分11 个苹果，就少10 个苹果;如果每人分13 个苹果，就少30 个苹果。 4、一盈一尽（刚好分完），可用公式：盈÷（两次每人分配数得差）=人数。 测试：如果每人分6 个苹果，就剩下40 个苹果;如果每人分10 个苹果，就刚好分完。 5、一亏一尽（刚好分完），可用公式：亏÷（两次每人分配数得差）=人数。 测试：如果每人分14 个苹果，就少40 个苹果;如果每人分10 个苹果，就刚好分完。 由上面得问题，我们归纳出盈亏问题得公式： 【提示】解决这类问题得关键就是要抓住两次分配时盈亏总量得变化，经过比对后，再来进行计算。 【例题1】某班去划船，如果每只船坐4 人，就会少3 只船;如果每只船坐6 人，还有2 人留在岸边。问有多少个同学? （） A、30 B、31 C、32 D、33 解析：此题答案为C。 设小船有x 只，根据人数不变列方程：4(x+3)=6x+2，解得x=5。 所以有同学6×5+2=32 人。 盈亏问题例题讲解：
三重积分的计算方法与例题
三重积分的计算方法： 三重积分的计算是化为三次积分进行的。其实质是计算一个定积分（一重积分）和一个二重积分。从顺序看： 如果先做定积分?2 1),,(z z dz z y x f ，再做二重积分??D d y x F σ),(，就是“投 影法”，也即“先一后二”。步骤为：找Ω及在xoy 面投影域D 。多D 上一点（x,y ）“穿线”确定z 的积分限，完成了“先一”这一步（定积分）；进而按二重积分的计算步骤计算投影域D 上的二重积分，完成“后二”这一步。σd dz z y x f dv z y x f D z z ??????Ω =2 1]),,([),,( 如果先做二重积分??z D d z y x f σ),,(再做定积分?2 1 )(c c dz z F ，就是“截面 法”，也即“先二后一”。步骤为：确定Ω位于平面21c z c z ==与之间，即],[21c c z ∈，过z 作平行于xoy 面的平面截Ω，截面z D 。区域z D 的边界曲面都是z 的函数。计算区域z D 上的二重积分??z D d z y x f σ),,(，完成 了“先二”这一步（二重积分）；进而计算定积分?2 1 )(c c dz z F ，完成“后 一”这一步。dz d z y x f dv z y x f c c D z ]),,([),,(2 1σ??????Ω = 当被积函数f （z ）仅为z 的函数（与x,y 无关），且z D 的面积)(z σ容易求出时，“截面法”尤为方便。 为了简化积分的计算，还有如何选择适当的坐标系计算的问题。可以按以下几点考虑：将积分区域Ω投影到xoy 面，得投影区域D(平面) （1） D 是X 型或Y 型，可选择直角坐标系计算（当Ω的边界曲
数值分析论文
题目：论数值分析在数学建模中的应用 学院: 机械自动化学院 专业: 机械设计及理论 学号: 学生姓名: 日期: 2011年12月5日
论数值分析在数学建模中的应用 摘要 为了满足科技发展对科学研究和工程技术人员用数学理论解决实际的能力的要求，讨论了数值分析在数学建模中的应用。数值分析不仅应用模型求解的过程中，它对模型的建立也具有较强的指导性。研究数值分析中插值拟合，解线性方程组，数值积分等方法在模型建立、求解以及误差分析中的应用，使数值分析作为一种工具更好的解决实际问题。 关键词 数值分析；数学建模；线性方程组；微分方程 the Application of Numerical Analysis in Methmetical Modeling Han Y u-tao 1 Bai Y ang 2 Tian Lu 2 Liu De-zheng 2 （1 College of Science ，Tianjin University of Commerce ，Tianjin ，300134 2 College of Science ，Tianjin University of Commerce ，Tianjin ，300134) Abstract In order to meet the technological scientific researchers who use mathematical theory to solve practical problems, the use of numerical analysis in mathematical modeling is discussed.Numerical analysis not only solve the model,but also relatively guide the model.Research on some numerical methods in numerical analysis which usually used in mathmetical modeling and error analysis will be a better way to solve practical problems. Key Words Numerical Analysis ；Mathematical Modeling; Linear Equations ；differential equation 1. 引言 数值分析主要介绍现代科学计算中常用的数值计算方法及其基本原理，研究并解决数值问题的近似解，是数学理论与计算机和实际问题的有机结合[1]。随着科学技术的迅速发展，运用数学方法解决科学研究和工程技术领域中的实际问题，已经得到普遍重视。数学建模是数值分析联系实际的桥梁。在数学建模过程中，无论是模型的建立还是模型的求解都要用到数值分析课程中所涉及的算法,如插值方法、最小二乘法、拟合法等，那么如何在数学建模中正确的应用数值分析内容，就成了解决实际问题的关键。 2. 数值分析在模型建立中的应用 在实际中，许多问题所研究的变量都是离散的形式，所建立的模型也是离散的。例如，对经济进行动态的分析时，一般总是根据一些计划的周期期末的指标值判断某经济计划执行的如何。有些实际问题即可建立连续模型，也可建立离散模型，但在研究中，并不能时时刻刻统计它，而是在某些特定时刻获得统计数据。例如，人口普查统计是一个时段的人口增长量，通过这个时段人口数量变化规律建立离散模型来预测未来人口。另一方面，对常见的微分方程、积分方程为了求解，往往需要将连续模型转化成离散模型。将连续模型转化成离散模型，最常用的方法就是建立差分方程。 以非负整数k 表示时间，记k x 为变量x 在时刻k 的取值，则称k k k x x x -=?+1为k x 的一阶差分，称k k k k k x x x x x +-=??=?++1222)(为k x 的二阶差分。类似课求出k x 的n 阶差分k n x ?。由k ，k x ，及k x 的差分给出的方程称为差分方程[2]。例如在研究节食与运动模型时，发现人们往往采取节食与运动方式消耗体内存储的脂肪，引起体重下降，达到减肥目的。通常制定减肥计划以周为时间单位比较方便，所以采用差分方程模型进行讨论。记第k 周末体重为)(k w ，第k 周吸收热量为)(k c ，热量转换系数α，代谢消耗系数β，在不考虑运动情况下体重变化的模型
数值计算方法》试题集及答案
《计算方法》期中复习试题 一、填空题： 1、已知3.1)3(,2.1)2(,0.1)1(===f f f ，则用辛普生（辛卜生）公式计算求得 ?≈3 1 _________ )(dx x f ,用三点式求得≈')1(f 。 答案：2.367，0.25 2、1)3(,2)2(,1)1(==-=f f f ，则过这三点的二次插值多项式中2 x 的系数为 ，拉 格朗日插值多项式为 。 答案：-1， )2)(1(21 )3)(1(2)3)(2(21)(2--------= x x x x x x x L 3、近似值*0.231x =关于真值229.0=x 有( 2 )位有效数字； 4、设)(x f 可微,求方程)(x f x =的牛顿迭代格式是( )； 答案 )(1)(1n n n n n x f x f x x x '--- =+ 5、对1)(3 ++=x x x f ,差商=]3,2,1,0[f ( 1 ),=]4,3,2,1,0[f ( 0 )； 6、计算方法主要研究( 截断 )误差和( 舍入 )误差； 7、用二分法求非线性方程 f (x )=0在区间(a ,b )内的根时，二分n 次后的误差限为 ( 1 2+-n a b )； 8、已知f (1)＝2，f (2)＝3，f (4)＝5.9，则二次Newton 插值多项式中x 2系数为( 0.15 )； 11、 两点式高斯型求积公式?1 d )(x x f ≈( ?++-≈1 )] 321 3()3213([21d )(f f x x f )，代数精度 为( 5 )； 12、 为了使计算 32)1(6 )1(41310-- -+-+ =x x x y 的乘除法次数尽量地少，应将该表达 式改写为 11 ,))64(3(10-= -++=x t t t t y ，为了减少舍入误差，应将表达式1999 2001-
各种利息计算方法例题[]
各种利息计算方法例题 利息计算基本公式：利息=本金×利率×存期=本金×天数×日利率=本金×月数×月利率 税后利息=利息×80% 天数计算=月×30天+另头天数（如4月24日即为144天）利率表示法：%代表年利率，‰代表月利率，万分比代表日利率。 1、活期储蓄存单：按实际存期有一天算一天，大小月要调整。现行日利率为每天0.2元。 例：2006年2月18日存入的活期存单一张，金额为1000元，于06年05月08日支取。问应实付多少利息？ 解：（158-78-1）天×0.1万×0.2元×80%=1.26元 2、定期存款利息计算： A、提前支取按活期存单的计算方法计算。 B、到期支取的利息=本金×年利率×年数 C、过期支取的利息=到期息＋过期息（到期息参照B，过期息参照A）实付利息=应付利息×80% 例：※2006年03月16日存入一年期存款一笔，金额为50000元，于2006年9月3日支取，利率为2.25%,问应付给储户本息多少？ 解：实付息=（273-106+4）天×5万×0.2元×80%=136.80元 本息合计=50000+136.8=50136.80元 ※2001年6月16日存入五年期存款一笔,金额为20000元,利率为2.88%,
于2006年6月16日支取,问应实付多少利息? 解：实付息=20000×2.88%×5年×80%=2304元. ※2003年01年27日存入三年期存款一笔,金额为12000元,利率2.52%,于2006年6月16日支取,问实付利息为多少? 解：到期息=12000×2.52%×3年=907.2元 过期息=(196-57+1)×1.2万×0.2元=33.60元 实付利息=(到期息+过期息)×80%=(907.2+34.08)×0.8=752.64元. 3、利随本清贷款利息计算：方法与活期存单一样，按头际天数有一天算一天。逾期归还的，逾期部分按每天3/万计算。（现行计算方法是按原订利率的50%计算罚息） ※例：某户于2006年2月3日向信用社借款30000元，利率为10.8‰,定于2006年8月10日归还,若贷户于2006年7月3日前来归还贷款时,问应支付多少利息? 解：利息=（213-63+0）天×(10.8‰÷30)×30000元=1620元. ※例：某户于2005年10月11日向信用社借款100000元,利率为9.87‰,定于2006年5月10日到期,贷户于2006年6月15日前来归还贷款,问应支付多少利息? 解：利息=（160+360-311+2）天×100000元×(9.87‰÷30)+(195-160+1)天×100000元×(9.87‰÷30×1.5)=6941.90+1776.60=8718.50元 4、定活两便利息计算：存期不足三个月按活期存款利率计算。三个月以上六个月以下的整个存期按定期三个月的利率打六折计算，六个月以上一年以下的整个存期按定期六个月的利率打六折计算，超过一年的整个存期都按一年期利率
尺寸链典型案例计算分析报告
尺寸链分析报告 工艺过程： 1、橡胶圈由分离机构从直振中拉出到固定位置。 2、视觉拍照，找橡胶套中心位置。 3、机器人理线工位辅助理线，配合机器人夹具将探 头sensor 线理直好插入橡胶圈。 已知条件： 1、橡胶套的内圆公差中心半径公差(理论中心与实际安装中 心的差值)mm A 15.015.010+-=2、探头的外圆半径公差(理论中心与实际安装中心的差 值)mm A 05.005.020+-=3、机器手抓取重复放置精度(理论中心与实际安装中心的差 值)mm A 05.0030+=4、相机本身引导误差mm A 05.00 40+=5、人工示教的容差mm A 2.01.050++=问题描述： 已知安装探头sensor 时机器探头中心与硅胶套中心的偏差， 即半径差值0.5mm ，即（探头能够安装进去橡胶套的最大偏 差值0.5mm 能够安装成功） 求：安装探头sensor 时机器探头中心与硅胶套中心的偏差， 即半径差值0 A 求解：根据题意，增环：1A ，2A ，3A ，4A ，5A ， 减环：无 封闭环：0 A 方法：尺寸链计算步骤及方法（统计法） 1.尺寸链的分析建立如图：
2.计算封闭环的基本尺寸： 封闭环的基本尺寸等于所有增环的基本尺寸和减去所有减环的基本尺寸和。 0=A 3.计算封闭环的公差： 批量生产条件下，组成环与封闭环的实际偏差均服从正态分布，且实际尺寸分布范围与公差带宽度一致。此时，封闭环的公差平方值等于所有组成环公差平方值之和。 4 .0, 16.01.005.005.01.03.0, 022222202 52423222120==++++=++++=T T T T T T T T 公差：公差：公差：4.计算封闭环的中间偏差。 封闭环中间偏差等于所有增环中间偏差之和减去所有减环中间偏差之和。 注：中间偏差等于上下偏差代数和再除以2.2 .0, 15.0025.0025.000, 00543210=?++++=??+?+?+?+?=?中间偏差：中间偏差：中间偏差：5.计算封闭环的极限偏差。 上偏差等于中间偏差加上二分之一公差值；下偏差等于中间偏差减去二分之一公差值。4.00 00000000, 02/4.02.02/)(, 4.02/4.02.02/)(+==-=-?==+=+?=A T A EI T A ES 偏差：下偏差：上偏差：答：满足装配精度要求,最大公差0.4mm,小于理论偏差0.5mm. 例2：感谢百度作者qq1473114691经验分享的方法： (2) (3)
数值分析学习方法
第一章 1霍纳（horner）方法： 输入=c + bn*c bn?1*c b3*c b2*c b1*c an an?1 an?2 ……a2 a1 a0 bn bn?1 bn?2 b2 b1 b0 answer p（x）=b0 该方法用于解决多项式求值问题=anxn+an?1xn?1+an?2xn?2+……+a2x2+a1x+a0 ? 2 注：p为近似值 p（x） 绝对误差： ?|ep?|p?p ?||p?p rp? |p| 相对误差： ?|101?d|p?p rp?? |p|2 有效数字: （d为有效数字，为满足条件的最大整数） 3 big oh(精度的计算)： o(h?)+o(h?)=o(h?); o(hm)+o(hn)=o(hr) [r=min{p,q}]; o(hp)o(hq)=o(hs) [s=q+p]; 第二章 2.1 求解x=g(x)的迭代法用迭代规则 ，可得到序 列值｛｝。设函数g 满足 y 定义在得 。如果对于所有 x ，则函数g 在 ，映射y=g(x)的范围 内有一个不动点; 此外，设 ，存在正常数k<1，使 内，且对于所有x，则函数g 在 内有唯一的不动点p。 ，（ii）k是一个正常数， 。如果对于所有 定理2.3 设有（i）g，g ’（iii ） 如果对于所有x在
这种情况下，p成为排斥不动点，而且迭代显示出局部发散 性。波理 尔 查 . 诺 二 分 法 ( 二 分 法 定) <收敛速度较慢> 试值（位）法：<条件与二分法一样但改为寻求过点(a,f(a))和(b,f(b))的割线l与 x轴的交点(c,0)> 应注意 越来越 小，但可能不趋近于0，所以二分法的终止判别条件不适合于试值法 . f(pk?1) 其中k=1,2,……证明：用 f(pk?1) 牛顿—拉夫森迭代函数：pk?g(pk?1)?pk?1? 泰勒多项式证明 第三章线性方程组的解法对于给定的解线性方程组ax=b a11x1 ? a12x2 ? ? ? a1nxn ? b1 a21x1 ? a22x2 ? ? ? a2nxn ? b2 ? an1x1 ? an2x2 ? ? ? annxn ? bn 一gauss elimination （高斯消元法第一步forward elimination 第二步 substitution 二lu factorization 第一步 a = lu 原方程变为lux=y ； 第二步令ux=y,则ly = b由下三角解出y；第三步 ux=y,又上三角解出x ； 三iterative methods（迭代法） a11x1 ? a12x2 ? ? ? a1nxn ? b1 a21x1 ? a22x2 ? ? ? a2nxn ? b2? ） back 初始值 0,x0,?,x0x1n2 四 jacobi method 1.选择初始值 2.迭代方程为 0,x0,?,x0x1n2 k?1? x1k?1 ? x2
坯布克重计算方法
坯布克重计算方法 坯布规格 190T涤塔夫 68DX68D 20筘*2穿*30梭门副165cm 以上算法只实用与不加捻产品，算出来是米克重哦．当然不加捻产品算克重都可以套用．如果是织造厂算原料的用量，在纬的门副上至少+１０ＣＭ．经不用加. 一．原料用量计算 定义：织一米坯布需经纬原料多少克。 1. 长纤类：经用量（g/m）＝总经跟数×(D/9000)× 纬用量（g/m）＝坯布纬密×上机门幅×(D/9000)× 或＝成品纬密×成品门幅×(D/9000)× 注：＝1+10％，10％为织缩率+损耗，一般FDY取，DTY取例1 涤塔夫上机门幅168cm 17号筘3穿坯布纬密28根/cm (168×17#×3入) 原料：63DFDY×63DFDY 求：原料用量 解：经用量=168×17×3×63/9000×=m
纬用量＝168×28×63/9000×＝m 即理论米克重＝+＝100g/m 例2. 75DFDY+100DDTY 五枚缎 12＃×5入×16630根/cm 求：原料用量 解：经用量＝12×5×166×75/9000×=m 纬用量＝30×166×100/9000×=m 2. 短纤类：经用量＝×（经密/经线纱支）×门幅 纬用量＝×（纬密/纬线纱支）×门幅 注：也可把短纤换算成长纤,用长纤公式来计算，D=5315/S 例1. 70D锦纶FDY×21S棉133×7259”，求：原料用量解：经用量＝133×59×70/9000×＝m 纬用量＝×72/21×59＝m 例2. 经：70D锦纶DTY+30D氨纶
纬：30S棉+40D氨纶86 ×5859” 求：原料用量及原料所占比例 解：氨纶有个拉缩比，一般30D的氨纶按10-13D计算，40D按15-18D 计算，拉缩比一般可取 3，如40/3＝，按15-18算 经用量锦纶＝70×86×59/9000×＝ 经用量氨纶＝86×59×13/9000×＝ 纬用量棉＝58×59×5315/30/9000×＝ 纬用量氨纶＝58×59×18/9000×＝ 总＝+++＝m 氨纶比例＝（+）/＝％ 棉比例＝＝％ 锦纶比例＝＝33％ 3. 加捻类：经用量＝总经根数×D/9000××(1+捻缩率) ×（1+蒸缩率） 纬用量＝坯布纬密×上机门幅×D/9000××(1+捻缩率) ×（1+蒸缩率） 捻缩率
数值分析答案
数值分析模拟试题1 注：计算题取小数点后四位。 注：计算题取小数点后四位。 1. (10分)利用Gauss-Legendre 求积公式 ?-++-≈1 1)7746.0(5556.0)0(8889.0)7746 .0(5556.0)(f f f dx x f 导出求积分 3()f x dx -?的三点高斯型求积公式。 2. (15分)写出求解线性代数方程组 123121322531272 x x x x x x x -+=??-+=-??+=? 的Gauss-Seidel 迭代格式，并分析此格式的敛散性。 3. (15分) 设矩阵21011000201010A ??????=?????， （1）试计算 ||||A ∞。 （2）用Householder 变换阵H 将A 相似约化为上Hessenberg 阵，即HAH 为上 Hessenberg 阵。 4. (10分) 求关于点集{}1,2,3,4的正交多项式 {}012(),(),()x x x ???。
5. (10分)用最小二乘法确定一条经过原点的二次曲线，使之拟合下列数据 1.0 2.0 3.0 4.00.8 1.5 1.8 2.0i i x y ??? 6. (20分)给出数据点： 0 134 19156 i i x y =??=? （1）用012,,x x x 构造二次Lagrange 插值多项式2()L x ，并计算 1.5x =的近似值2(1.5)L 。 （2）用123,,x x x 构造二次Newton 插值多项式2()N x ，并计算 1.5x =的近似值2(1.5)N 。 （3）用事后误差估计方法估计2(1.5)L 、2(1.5)N 的误差。 7．(10分) 设矩阵A 可逆，A δ为A 的误差矩阵，证明：当11A A δ-<时， A A δ+也可逆。 8．(10分)设()f x 四阶连续可导，0,0,1,2.i x x ih i =+=试建立如下数值微分 公式 ''01212()2()()()f x f x f x f x h -+≈ 并推导该公式的截断误差。