最新数学建模课后习题答案.docVIP

精品文档实验报告 姓名：和家慧 专业:通信工程 学号：20121060248 周一下午78节 方程及方程组的求解 实验一： 一 实验目的：学会初步使用方程模型，掌握非线性方程的求解方法，方程组的求解 函数直接求解法等。TLAB方法，MA二 问题：路灯照明问题。在一条20m宽的道路两侧，分别安装了一只2kw和一只3kw的路灯，它们离地面的高度分别为5m和6m。在漆黑的夜晚，当两只路灯开启时 （1）两只路灯连线的路面上最暗的点和最亮的点在哪里？ （2）如果3kw的路灯的高度可以在3m到9m之间变化，如何路面上最暗点的亮度最大？ （3）如果两只路灯的高度均可以在3m到9m之间变化，结果又如何？ 三 数学模型 解： 根据题意，建立如图模型 y P1=2kw P2=3kw S=20m P2 照度计算公式： ?sinpP1 I?kR1 R2 2rh2 h1 Q （k为照度系数，可取为 1； α2 α P为路灯的功率） x X O S （1）设Q(x,0)点为两盏路灯连线上的任意一点，则两盏路灯在Q点的照度分别为 ??sinpnpsi2112I?kI?k 21 22RR21h2?1ns?i22xRh?? 1R11 1 精品文档． 精品文档h2?222nsi?)x??(sR?h 222R2 点的照度： Q 18hPh10P2112????I(x) 3223322322))((25?x)x((h?)20?x(36?s((h?(?x))21 I(x)的最大值和最小值，所以应先求出函数的极值点要求最暗点和最亮点，即为求函数?3Phx3Ph(s?x)?30x54(20?x)2211???)I(x?2522222555))h(?(s?x))(20h(?x)?x?(25?x)(3621 算法与编程 0)?I(x 利用MATLAB求得x的值时 代码：s=solve((-30*x)/((25+x^2)^(5/2))+(54*(20-x))/((36+(20-x)^2)^(5/2))); s1=vpa(s,8); s1 计算结果 运行结果： s1 = 19 9.338299136 8.538304309-11i .2848997038e-1 8.538304309+11i 因为x=0，选取出有效的x值后，利用MATLAB求出对应的I(x)的值，如下表： x 0 0.028489970 9.3382991 19.976695 20 I(x) 0 0 0 0 0 综上，x=9.33m时，为最暗点；x=19.97m时，为最亮点。 （2）路灯2的高度可以变化时，Q点的照度为关于x和h的二元函数： 2 精品文档． 精品文档hPhh103P22121??I(x,h)?? 233223222322)x)?(h(?(s?x))(h?x)20?(h25(?x)221 I(x,h）同理，求出函数）的极值即为最暗点和最亮点 与（122h3P?IP2220??? h?522223)?x)h?(s?x))(h?(s(222算法与编程 x：matlab