数学模型中的反问题逆问题(36页).doc

-数学模型中的反问题逆问题-第 35 页数学模型中的反问题向下运动向上运动风筝数学模型竟赛中有很多涉及反问题。如2010国赛中A题和2011年美赛中A题都涉及反问题。顾名思义，反问题是相对于正问题而言的。正问题的定义为：按着自然顺序来研究事物的演化过程或分布形态，起着由因推果的作用。自然顺序的定义为：不受任何限制和约定俗成的顺序，一般地都认为他们是自然而然的，无须多加解释的。在一般地语境下，认为这些顺序都是是前提条件的。如时间顺序、空间顺序、因果顺序，等等。纯粹的自然顺序的例子是第一，第二，第三这种升序；或者反过来的倒序；约定俗成的例子是上北下南左西右东。反问题的定义为：根据事物的演化结果，由可观测的现象来探求事物的内部规律或所受的外部影响，由表及里，索隐探秘，起着倒果求因的作用。可以看出，正、反两方面都是科学研究的重要内容。但相对正问题，反问题求解难大，计算量大。许多人知道求解问题的思路，但由于选用计算方法不适当，在几天内求不出计算结果，失去获奖机会。尽管一些经典反问题的研究可以追溯很早，反问题这一学科的兴起却是近几十年来的事情。在科学研究中经常要通过间接观测来探求位于不可达、不可触之处的物质的变化规律；生产中经常要根据特定的功能对产品进行设计，或按照某种目的对流程进行控制。这些都可以提出为某种形式的反问题。可见，反问题的产生是科学研究不断深化和工程技术迅猛发展的结果，而计算技术的革命又为它提供了重要的物质基础。 现在，反问题的研究已经遍及现代化生产、生活、研究的各个领域。简单的概括不足以说明问题，我们下面具体介绍一些常见的反问题类型，希望大家能够对它有一个概括的了解.第一节反问题的例子例1 物体下落距离L与时间T，正问题是：已知物体的高度，测量下落时间，即(x). 反问题是：已知物体下落时间，求物体的高度,即(t)。2之前，能用时钟测量物体下落时间，但反过来，给定下落时间，测量物体高度比较难。对于没有读中学的人，能完成时钟测量物体下落时间的试验。但给他物体下落时间，测量物体的下落高度是不容易的事情。例2 年龄与身高。正问题是，根据年龄T，每周岁测身高H，得到身高H与年龄T的关系(T). 反问题是：已知身高H，求年龄T，即求关系(T). 例3速度V与轨道形状(x)，其摩擦系数为，z为高度，初始速度为V0，末速度为(). 正问题是，已知轨道轨道形状(x)。对大学生，正问题能求出来，但反问题有些难。例4 热传导问题（2013美赛A题）设点P到边界的距离为x, 传热系数为a, 温度(). 正问题是：已知距离x, 传热系数a的数值, 求温度T。如用一维传热公式：反问题是：已知温度，求传热系数为a,例5 光电板问题(2012A题)设屋顶面积为D，光电板长为L，宽为H。正问题是：已知D，L，H，求在屋顶上铺设光电板最大数量N。反问题是：已知光电板总铺设面积D*，光电板长L，宽H，求屋顶面积。由上面几个例子，可以在数学上定义正问题为(x)，定义域为D，值域为V。反问题为(y). 由高等数学可知，若函数f(x)在D上是单调的，则反函数g(y)存在且唯一。相对正问题而言，反问题计算量大，选用适当的计算方法是成功求解反问题的关键。因而要求在求反问题之前，要求学生掌握基本的计算方法。第二节计算方法在数学建模中，求解方程的根是经常遇到的。常用求根方法有迭代法，二分法，牛顿法，极小值法，一维寻查法，格子法。迭代法设函数f(x)(x)有一根x*, 则f(x*)=0, 或x*(x*)=0; 或x*(x*); 定义求根的迭代公式为：定理: 若导数g的绝对值小于1, 即|L<1, 则迭代收敛。证：由于x*(x*)，则1*()(x*)() (*)有| 1*|< *|<L2| 1*|<1| x0*|因为L<1, 则极限>0, 故àx*. 证毕。例 求f(x)*x的零点。解：这里g(x)*x, g(x)=2x, 则当<0.5时，(x)|<1, 即<0.5时，迭代公式. 12收敛。取x0=0.1, 计算得X1022=10-2X212=(10-2)2=10-4最后求得àx*=0. 实际上，我们知道0为*x的解，但它还有一解1; 由于|22*12, 则用上面迭代公式(x)*x求不出解1. 它需要构造另一种迭代公式.1()=容易验证当1时，|<1. 取x0=2, 计算得X10=2X21=(2)=2X33=(2)=2最后求得àx*=1.由上面例子可知，对同一函数f(x)，它的不同零点对应的迭代公式不同。在高等数学里，我们已学习下面定理。定理：设f(a)f(b)<0, f(x)在区间上连续可导，则至少有一个()中的点x*，使f(x*)=0.i(b00)à0, àx*. 例. 求f(x)=12 在区间0.5, 2上的零点。2=0.75, f(b0)(2)=1-223, 有f(a0)f(b0)=0.75*(-3)<0, 故在0.5,2上f(x)有一零点x*. 取x1=(a00)/2=(0.5+2)/2=1.25, 有f(x1)(1.25)0.5625, f(a0)f(x1)=0.75*(-0.5625)<0, 则在区间a01中有零点x*, 故取a10=0.5, b11=1.25, x2=(a11)/2=0.875, 计算得f(x2)(0.875)=1-(0.875)2=0.2343, F(a1)f(x2)=0.75*0.2343>0, 则零点x*在区间x21=0.875,1.25中，故取a2=0.875, b2=1.25. 如此计算下去， 当<=0.3时，求得3=1.0625。它为x*的近似值。二分法的计算步骤为：1)输入a00, 误差限；2)若f(a0)f(b0)>0, 无根，停止计算。否则转下一步；3)取x1=(a00)/2, 若f(a0)f(x1)<0, 取a1011; 否则取a11, b10;4)若b11<, 输出近似根x*=(a11)/2; 否则1à0, b1àb0, 转第三步。二分法能用图形来说明，其示意图见图2.1, 图中给出了点a00123, 它们根据二分法计算。由图可知，当二分次数增加时，中间点相互靠近，收敛于零点x*.图2.1 二分法示意图2.1.3 极小值法定型：若x*为f(x)的零点，则它为F(x)2(x)的极小值点。证：由于F(x)非负，F(x*)2(x*)=00=0,  则x*为F(x)的一个极小值点。我们容易得：定理：若F(x) 2(x)(x*)=0, 则f(x*)=0.可见，f(x)的零点计算问题能化为极小值计算问题。它常用一维寻查法求解。一维寻查法比较简单，它的计算步骤为1)输入初始点d0, 步长h, 误差；2)计算函数值F(d0)(d0)(d0);3)若F(d0)<F(d0)(d0); 转第6步。4)若F(d0)>F(d0), 取d10; 否则取d10;5)令d1àd0, hà2h, 转第2步。6)取0, 0, 用二分法求极值点。二分法求极值点的原理与求根原理类似。由下面定理给出：定理：设F(x)在上连续，且0F(x),若c为中的点，且 F(c)<F(a)(b), 则F(x)在上存在极小值点x*. i(b00)à0. 可以证明，àx*为极小值点。由上面讨论可知，求极小值点分为两步，先求极点所在的区间, 然后用二分法逐步缩小区间，求出极小值点。其计算过程可以用图2.2说明。图中函数F只有一个极小值点。 给定初值d0和步长h, 求出d0为最小值，取2h, 计算得d0+2h也为最小，再取4h, 计算得d0+2h也为最. 图2.2 极小值示意图例. 用极小值法求函数f(x)=1*x的零点，x0=1.4, 0.1.解. 令F(x)= f2(x*)=(12)2先用一维寻查法求含有根的区间. 计算F(x0)(1.3)=0.4761; F(x0)(1.4)=0.9216; F(x0)(1.5)=1.5612; 比较3个数值，x1=1.3时0.4761最小。将步长放大2倍，取0.2, 计算F(x1)(1.1)=0.0441; F(x1)(1.3)=0.4761; F(x1)(1.5)=1.5612; 比较3个数值，x2=1.1时0.0441最小。再将步长放大2倍，取0.4, 计算F(x2)(0.7)=0.216; F(x2)(1.1)=0.0441; F(x2)(1.5)=1.5612; 比较3个数值，x3=1.1时0.0441最小。因而取0.7, 1.5. 再用二分法求极值点。取a0=0.70=1.20=0.95, 0.25; 中点x0=(a00)/2=0.825; y0=(c00)/2=1.075; 计算得F(0.7)=0.216; F(0.825)=0.102; F(0.95)=0.0095; F(1.075)=0.02421; F(1.2)=0.1936; 当z1=0.95时，函数F(0.95)=0.0095最小。则取a1=0.825, c1=0.95, b1=1.075; x1=0.8875; y1=1.0125; 计算得F(0.825)=0.102; F(0.8875)=0.04509, F(0.95)=0.0095; F(1.0125)=6.3284, F(1.075)=0.02421; 给定误差=0.1时，若(1)/4<, 输出近似根x*=1.0125.对于高维问题，格子法是求极值点的常用方法。它的思想与二分法类似，基本原理为，给定非负的高维函数(X), 初始点X0, 步长h, 将每个坐标分量加上h和减去h, 求最小值y0F(X0)()(), 和对应的坐标点X1, 若X10，取hà2, 步长减半，否则取hà2h, 步长加倍，再将X1的每个坐标分量加上h和减去h, 求最小值点X2，如此下去，直到步长h<为止。最后为近似最小值点。例.求方程式组的极小解：()-1=0; y()2=0;解：令F()=()-12+y()22取初值点X0=(0,0), 步长0.8; 计算得：F(-0.8,0)=7.24; F(0.8,0)=4.04 F(0,0)=5.0; F(00.8)=3.00; F(0,0.8)=7.355;可知X1=(00.8)为最小值点，取加倍步长1.6, 计算中心得：F(-1.6,0.8)=20.9; F(1.6,0.8)=4.923; F(0,0.8)=3.0; F(02.4)=83.6; F(0,0.8)=7.32;则X21=(00.8)为最小值点，取减半步长0.8,继续计算，最后求得近似极小值点(1.31751.5675).满足误差=0.01.2.1.5 多项式拟合 在反问题计算中，多项式拟合是常用的方法，其基本原理是：给定测量数据(), 1, 2, , 求一个多项式. 012x2+将数据代入得上式可写为矩阵表达式：这里：两边乘以X的转置有故有编程：0,1,2,30,9,35,90(中表示Y的转置矩阵) 1x(i,1)=1 1:3 x(1)()*t(i)(X*X)*X*Y当1时，为线性函数，可以由上式求出具体表达式：式中E(X)为X的平均值，E(Y)为Y的平均值。D(X)为X的方差。上式与最小二乘法得到的结果相同。例. 已知数据(0,0), (1,1), (2,4),(3,8)， 求一元回归函数?解. 我们求得4E(X)=(0+1+2+3)/4=7/4; E(Y)=(0+1+4+8)/4=13/4E()=(0*0+1*1+2*4+3*8)/4=33/4计算得： 1.9524; 0.1667; 则一元回归为：2.1.6 数值积分在数学模型竟赛中，能求出分析解的积分太少，大多只能用数值方法离散计算。设h为步长，0<x1<积分用复化中点公式计算例.计算f(x)2在区间-1,1上的积分。解.取步长1,得区间-1,0,0,1, 中点z00.5, z1=0.5; 积分为实际精确值为0.6667, 相差0.1667; 若取更小的步长，误差将变小。 对于二维积分，其积分区域为D，将区域D划分为若干四边形，第k个四边形中点坐标为, 面积为, 则用复化中点计算例. 计算f()2 2在区域0,2;0,2上的积分。解.取步长1,得子区域0,1;0,1,0,1;1,2, 1,2;0,1,1,2;1,2; 中点为P1=(0.5,0.5); P2=(0.5,1.5)3=(1.5,0.5); P4=(1.5,1.5); 子区域面积1.0; 积分为(P1)S1(P2)S2(P3)S3(P4)S4 = f(P1) (P2)(P3)(P4)=0.52+0.52+0.52+1.52+1.52+0.52+1.52+1.52=10实际精确值为10.667, 相差0.667, 相对误差为0.667/10.667=0.0625. 在可接受范围之内。当然，取小的步长，能缩小误差值。2.1.7 微分方程数值解设初值问题 . y ' (); y(x0)0取步长为h, 点x10, x21, x32, 导数 y(x1)=(y10), 则有. y10(x00). y21(x11)例. 给定初值问题, y=1, y(0)=1; 0.1; 求y(0.2)?解. y10(x00)=1+ y21(x11)=1+第三节反问题数学模型求解例1 物体下落距离L与时间T，正问题是：已知物体的高度，测量下落时间，即(x). 反问题是：已知物体下落时间，求物体的下落距离。解设距离L与时间有多项式关系：23选取取不同距离，测量下落时间，1,2,; 构造函数Q(): Q()=232取()T, 用最小二乘法，求出上式的极小值这里(L12,)T，(X1234)1=(1,1,1) ', X2=(T12,) ', X3=(T1222,2) ', X4=(T1323,3) '.例如测量得一组数据为表3.1 物体下落时间T与距离LTLTLTL219143由上面数据计算得：0.00126; 0.044694.72950.0007893;232.例2 年龄与身高。正问题是，根据年龄T，每周岁测身高H，得到身高H与年龄T的关系(T). 反问题是：已知身高H，求年龄T，即求关系(T). 解设年龄T与身高H的关系为根据下面数据表3.2 年龄T与身高HTHTH1399111248101203571112946612138575131476841415679315165810216174求得30, 9; 即有30+9T例3 速度V与轨道形状(x)，其摩擦系数为，z为高度，初始速度为V0，给定末速度为(),.求轨道形状(x)。解设y0(x0)1(x1),在区间x01,物体速度近似取v0, 由能量守恒定理，设M为物体质量，V1(1)为1 处的速度，有 初始(动+势)能=结束(动+势)能+磨擦力作功其数学表达式为：22 1(a1)S1 =(x10)S1式中a1为线段y0y1与x轴的夹角，为物体在x01段移动时磨擦力作的功，S1为线段的长度。为离心力，定义为：R为曲率半径。取y=(y10), y”=(y2-2y10)2.将区间划分为x01,x12, 0, h为步长，()为对应函数值，当已知第1点的速度1，则第i点的速度为2 12(1)()(1)式中为线段1与x轴的夹角，为物体在x 1段移动时磨擦力作的功，为线段的长度。将区间划分为x011,2,. 当给定函数(x)在点时值, 用上面的公式能计算出末速度V(). 设y12,为y在世点值组成的向量，与末速度V有关系(Y),给定, 得, 则计算出数据(), 1,2, 令 23只要求出常数, 对速度, 我们就可以求出对应. 然后用插值方法求出曲线(x).考虑一个具体的例子，设(x)定义在区间0,1,且f(1)=0.5;末速度4, 初始速度V0=10, 摩擦系数为=0.01,单位取为焦尔，米，牛顿。将0,1划分为两个区间0,.5,.5,1, 由于y0=02=0.5, 则只有一个未知量y1. 改变y1 值，求得V值。表3.3 速度与坐标y1的关系y1v1(0.5)v2(1.0)曲率半径R设速度2与坐标y1的关系 . y12我们求得7.6919; 4.2762; 0.59305; 即有2取2=4, 得y1=0.0759.第四节 具有反问题的数模竟赛题单板滑雪问题A：单板滑雪请设计一个单板滑雪场（现为“半管”或“U型池”）的形状，以便能使熟练的单板滑雪选手最大限度地产生垂直腾空。“垂直腾空“是超出“半管”边缘以上的最大的垂直距离。制定形状时要优化其他可能的要求，如：在空中产生最大的身体扭转。在制定一个“实用”的场地时可能需要权衡哪些因素？解单板滑雪场“U型池”见图4.1, 可以使单板爱好者们，从一面墙到另一面之间移动，跳跃并做花样动作。墙面形状对运动速度有影响，运动轨迹是三维的，为了简化，只考虑二维截面。由于对称性，只考虑一面墙，由于只有墙的形状可以改变，平底部分不能改变，因而墙面形状部分。二维截面的墙面见图4.2, 图中L为墙面宽度，h为高度，h为运动员离开墙面后上升高度。G为运动员身体重量，为离心力。设墙面曲线为(x), 运动员进入初始速度为V0, 末速度为。当f(x)已知时，根据力学原理，能计算出末速度,这是正问题。若给定末速度，求曲线f(x), 这是反问题。可知，反问题比正问题难得多。当曲线f(x)未知时，它要满足条件f(0)=0, f(L), f(0)=0, 1(L)=0. 即运动员沿水平方向进入弧曲线，沿垂直方向离开弧曲线。A1无磨擦当不考虑磨擦损失时，由能量守恒原理，有式中m为运动员质量，, 有H的计算式：故垂直腾空高度. 它说明当不考虑磨擦时，垂直腾空高度与曲线形状无关。A2有磨擦当计及磨擦时，设运动员从点0运动到点1时，速度为v1，由能量守恒原理，有式中为磨擦系数，为作用在曲线法向上分力，为曲线段长度,且有 2 2为切线角，R为曲率半径：取速度0, 由上面公式计算得速度v1. 用类似的方法，将区间0划分为x0=012,. 依次计算出速度v12,. 再改娈曲线值y12,求出对应速度，根据本题要求，要计算出最大速度对应的f(x)，则从所有计算结果中，选出最大的对应值y12,。再用曲线拟合的方法计算f(x). 在实际计算中，=0.1,初始f(x)为圆弧计算结果见表4.1, 由表可知，经过优化后，垂直腾空高度h增加近20%. (:)（） ()（） ()（） ()(m) ()(m) (%)1718192021A3空气阻力运动员速度高速运动时，将产生空气阻力：K为空气阻力系数，A为迎风面积，气体密度，v运动速度。则有能量关系为损失动能、磨擦力作功、势能、阻力作功。然后根据上面方法求出末速度,和垂直腾空高度h。A4扭转要求运动员腾空后，在空中产生最大的身体扭转。产生的方法是：运动员用脚作用于U池。见图4.3。设E为运动员离开U型池作的功，E1为旋转能量，为旋转角速度， E21为向上增加能量，E4为运动到池边，速度为的能量，E542为离开U型池向上运动的能量则有 则在空中停留时间T为式中A为运动员总转动角度。将h代入上式得求得A2最大值时的值再计算出V1，V2，V，得角度为计算结果见表4.2， 由表可知，身体扭转角增加15%。 () (J) () ()() (%)13310020202290131332002051235915 7: 第五节 具有反问题的数模竟赛题储油罐例5.1；2010年国赛A题 储油罐的变位识别与罐容表标定一 内容通常加油站都有若干个储存燃油的地下储油罐，并且一般都有与之配套的“油位计量管理系统”，采用流量计和油位计来测量进/出油量与罐内油位高度等数据，通过预先标定的罐容表（即罐内油位高度与储油量的对应关系）进行实时计算，对于实际储油罐(见图5.1)，试建立罐体变位后标定罐容表的数学模型，即罐内储油量与油位高度及变位参数（纵向倾斜角度a和横向偏转角度b ）之间的一般关系。请利用罐体变位后在进/出油过程中的实际检测数据，根据你们所建立的数学模型确定变位参数，并给出罐体变位后油位高度间隔为10的罐容表标定值。进一步利用实际检测数据来分析检验你们模型的正确性与方法的可靠性。二 逆问题数学模型 本文只考虑储油罐的纵倾和横倾变位，认为储油罐整体是钢性的。设储油罐内液体体积V与油浮显示油位高度h之间的函数关系为式中为纵倾角，为横倾角。再设为高度时的液体体积测量值, ，定义理论计算值与测量值的误差函数： (2.1)由于环境条件，测量值对应的变位倾角是未知的，但可以从上面公式求极值得到: (2.2)当V为倾角的线性函数时： (2.3)上式中，,为已知函数，则有的表达式： (2.4)利用最小二乘法，得极值解为： (2.5)式中： (2.6)对于一般情形，当倾角较小时，取V的线性项有较高的精确度，也能作为迭代初值。由上面讨论可知，当已知罐内液体体积的测量值和理论计算值，变位倾角可以由极值计算。下一节将讨论液体体积理论计算方法。三罐内液体体积理论计算设为储油罐贴体坐标系，平面为储油罐截面积最大的截面所在的平面，z轴垂直于平面。再设液面所在的平面P的法向量为()，它由倾角决定，对应平面方程为 (3.1)式中点(x000)为油浮子在液面上的点。本文假设油浮子始终在液面上。平面将储油罐截为上表面, 和下表面。设为平面0在储油罐内部区域，则罐内液体体积V为： (3.2)根据平面P与油罐上下表面之间的关系，式中F由下面公式表示出： (3.3)由于罐内液体体积V为二重积分，当油罐上下表面用简单函数表示时，我们能求出积分表达式。当罐面比较复杂时，可以用数值方法计算。在下面实例中，可将平面以步长分成N个小区域，将每一个区域对应的油柱近似为直方体，然后对所有N个直方体进行离散求和。四 卧式圆筒球形封头型油罐4.1 上下表面和油面卧式圆筒球形封头型油罐的空间直角坐标系见图1, 图中取为贴体坐标系，中截面在平面上，取为上节中二重积分区域, L为储油罐主体（圆柱体）的长度，r为储油罐主体的截面园半径，q为两端球冠体的高。图1 油罐的空间直角坐标系的建立设油浮子始终作为油面上的点，坐标为 ，这里为油浮显示油位高度，储油罐纵倾角为，横倾角为。4.2 计算结果根据前面定义的误差函数, 利用2010全国大学生数学建模竞赛A题()所给出的实际测量的前302个数据，采用变步长遍历搜索得误差函数最小值对应的油罐纵倾角为，横倾角为. 利用部分测量数据计算参数，另一部分测量数据检验计算的参数是科学研究中常用的方法8。检验变位参数正确性运用求得的与的值得到体积与的对应关系，对A题所给出的实际测量剩余的后299个数据进行检验。则模型所得数据与实际数据的相对误差，如下图10所示。图2 模型计算数据与实验数据相对误差由图中相对误差图可知，模型的计算数据与实验数据误差不超过5%，或计算值与实测值吻合精度高于95%. 实际计算得标准偏差为0.0348, 平均绝对误差小于0.01,可以认为所建立的模型是正确的。4.4 罐容表标定值油罐变位后必需重新标定，根据本文计算公式，计算的标定值见表1, 油位高度间隔取为10。表1 变位后油位高度间隔为10的罐容表标定值油面高度油的容量油面高度油的容量油面高度油的容量 注：表格中油面高度的单位为米（m），油的容量的单位为升（L）五 结论根据油面与油罐表面之间的相对关系，选用油罐贴体坐标系，求出油量与油面高度及变位倾角之间的计算关系式，然后利用测量数据，构造逆问题求解倾角的误差函数。求出了极小值对应的倾角值。最后对罐容表进行了重新标定。文中油量为二重积分值，有多种数值方法能给出满足实际要求精度的计算值，且与用解析积分表达式计算值接近。当油罐表面函数比较简单时，可以求出积分表达式，或部分区域上的积分表达式。但当积分表达式较繁杂时，编写计算程序难度提高。因而，用数值方法计算积分是一个可行的方法。文中方法不仅适用于卧式圆筒球形封头型油罐，也能用于其它类型油罐逆问题计算。第六节 2012B题 太阳能小屋的设计1、问题的重述太阳能小屋， 就是在小屋的顶面及外墙铺设光伏电池，通过光伏电池吸收太阳辐射能量转化为直流电， 在经过逆变器转化为 220 的通用交流电源供居民生活使用，多余的存入电网。光伏电池对太阳的吸收与电池的种类，辐射强度，阳光直射角度等诸多条件有关，而不同种类的光伏电池的价格也相差甚远。根据附件中的数据，在考虑年发电总量又考虑铺设成本的情况下对小屋电池排布进行设计， 计算出小屋光伏电池 35 年寿命期内的发电总量、经济效益（当前民用电价按 0.5 元 计算）及投资的回收年限，并选出合适的逆配器与串并联接线方式。问题 1：根据山西省大同市的气象数据，仅考虑贴附安装方式，选定光伏电池组件，对小屋的部分外表面进行铺设，并根据电池组件分组数量和容量，选配相应的逆变器的容量和数量。问题 2：电池板的朝向与倾角均会影响到光伏电池的工作效率，请选择架空方式安装光伏电池，重新考虑问题 1。问题 3：根据给出的小屋建筑要求，请为大同市重新设计一个小屋，要求画出小屋的外形图， 并对所设计小屋的外表面优化铺设光伏电池，给出铺设及分组连接方式，选配逆变器，计算相应结果。问题一的分析我们在设计电池排布的时候既要考虑太阳辐射的能量，又要考虑铺设成本。然而题中已知每千瓦时的电价，我们可以将设备 35 年寿命内转换的电量转化为收益，然后减去初始成本得到总收益，作为目标函数简化求解。在对每一块墙面上电池板的排布，由于电池板的尺寸，大小和形状皆不经相同， 这样一个排布上的优化问题就涉及到排样算法，十分复杂而且只能不适合于本题的目标函数。 于是我们考虑将矩形拼排的问题转化为长度问题。经过观察发现大部分电池板宽度都在一米左右浮动，于是我们将墙面进行切割，结合实际可用空间形状切割为宽度一米左右的若干个条状得到总长度， 然后要求所有电池板的长度和的数值小于该总长度即可。为保证光伏组件正常工作， 考虑到光伏电源串并联的要求：只允许相同型号的光伏组件进行串联。 多个光伏组件串联后可以再进行并联，并联的光伏组件端电压相差不应超过 10%。最后连接，形成组件阵列。最后考虑到逆变器的额定工作电压（V）范围和功率容量（W）等参数进行分组设计。问题一的求解：首先我们要从房屋的六个表面中选出几个铺设电池板的墙面，而且不考虑四周其他建筑的遮光效应，如果一个面被选定则该面整面都可以用来铺设电池板。我们将六个面编号。图6-1太阳能小屋六个面单位时间接收到的辐射等于每个面在水平面、东、南、西、北五个面投影面积乘以该方向的辐射强度。首先不考虑门窗， 根据题目图中数据，六个面在五个方向的投影面积为：表 1 墙面投影面积表（单位：平方米）根据题中给出的太阳辐射强度，计算全年每个面接收到的能量。在铺设同等面积电池板时，求得全年能接受到的太阳辐射总能量墙面最多，而墙面最少，为了保障铺设电池板的经济效益，我们优先考虑按照、顺序安排电池板。 然后我们从墙面开始，由于实际电池块大部分小于1.114，我们先对墙面进行划分将区域分为 6 个条状长方形， 得到几个长方形长的总和，使用的所有的电池矩形的长的数量和应该小于这六个长方形长的数量和。最后结果为：附件A 编程技术基本上机操作1)打开;2)>>3)在编辑窗中输入程序；4）点击> 5）输入文件名w11, 保存6）屏幕上出现结果3;上机实习程序二；计算 1+2100; 1:10;上机实习程序三计算2+4+620; 如s>10, 输出“优秀”，否则输出“下一次优秀”0; 1:10; 2*k;(s>10)('优秀') (s<=10)('下一次优秀')上机实习程序4计算矩阵乘积：;1 2 34 5 6123*B上机实习程序5计算逆矩阵2 11 2A1(A)线性规划上机实习6在中，用函数只解决下面线性规划问题例1, 求下列线性规划问题写为程序为-3;1;1;1 -2 1;4 -1 -2;113;2 0 -110;0;0;x计算结果为：x = 4 1 9 -2线性规划上机实习7，0-1规划例2: 求下列0-1规划问题解： 程序为：9; i1=1:211-1; i2=1:222-1; i3=1:233-1; f11+3*x2-2*x3; (f1>4); f21+4*x23; (f1>4); f3=2*x23; (f3>6); f412; (f4>1); 3*x12+4*x3; (<z)(1)1(2)2(3)3; 3 21y(4)程序中为进行下一个循环，而不计算后面的语句。程序二：31;4;1 3 -2 1 4 1 0 2 1 1 1 0;4;4;6;1;x，()上面为0-1整数线性规划函数，求解的为目标函数的最小值, 所以要在f前面加个负号。u(1)=1.3175(2)1.5550.0125;0; 1:2; 1:1(j)(j)*i; (1)(2); 1(k,1)(k,2); a(k,3)=(*()-1)2+(y*()2)2; av0=59.8110.01;y1=0.250=02=0.50.5; x0=01=0.52=1; 1:101=0.025*k; 1=(y10)2=(y2-2*y10); 2/(112)1.5; (x10)2+(y10)2)0.5;0.5*v02*y1*(g*0.502*)(w<0);v1(2*w); 1=(y21)2=(y2-2*y10); 2/(112)1.5; (x21)2+(y21)2)0.5;0.5*v02*(y21)*(g*0.512*)(w<0);v2(2*w); v(k,1)1(k,2)1(k,3)2(k,4);v 1:9(i,3)(i,1)=1(i,2)(i,3)*;'*x; b1'*v(:,1)(b);b2*b1