数学建模的实验报告分析.docVIP

  1. 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
  2. 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  3. 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  4. 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  5. 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  6. 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  7. 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
问题 路灯照明问题。在一条20m宽的道路两侧,分别安装了一只2kw和一只3kw的路灯,它们离地面的高度分别为5m和6m。在漆黑的夜晚,当两只路灯开启时,两只路灯连线的路面上最暗的点和最亮的点在哪里?如果3kw的路灯的高度可以在3m到9m之间变化,如何路面上最暗点的亮度最大?如果两只路灯的高度均可以在3m到9m之间变化,结果又如何? 数学模型 已知P1为2kw的路灯,P2为3kw的路灯,以地面为X轴,路灯P1为Y轴,建立平面直角坐标系。其中,P1、P2高度分别为h1、h2,水平距离为S=20m。设有一点Q(x,0),P1、P2分别与其相距R1、R2。如下图示。 经查阅资料得,光照强度公式为:,设光照强度k=1。则,两个路灯在Q点的光照强度分别为: 其中: R12=h12+x2 R22=h22+(S-x)2 则Q点的光照强度Ix=I1+I2 分别按照题目中的不同要求,带入不同数值,求导,令导数为零,求得极值,进一步分析对比,求得最值。 算法与编程 当h1=5m,h2=6m时: symptoms x y x=0:0.1:20; y=10./sqrt((25.+x.^2)^3)+18./sqrt((36.+(20-x).^2).^3); plot(x,y) grid on; 在图中的0-20米范围内可得到路灯在路面照明的最亮点和最暗点 对Ix求导: syms x f=10./sqrt((25.+x.^2)^3)+18./sqrt((36.+(20-x).^2).^3) 运用MATLAB求出极值点 s=solve((-30*x)/((25+x^2)^(5/2))+(54*(20-x))/((36+(20-x)^2)^(5/2))); s1=vpa(s,8) s1 = -1 8.5383043+11.615790*i 19.976696 9.3382991 8.5383043-11.615790*i ③根据实际要求,x应为正实数,选择19.9767、9.3383、0.02849三个数值,通过MATLAB计算出相应的I值: syms x I=10/(25+x^2)^(3/2)+18/(36+(20-x)^2)^(3/2); subs(I,x,19.9767) subs(I,x,9.3383) subs(I,x,0.02849) ans = 0.0845 ans = 0.0182 ans = x 0.02849 9.33829 19.9766 I 0.0820 0.0182 0.0845 综上,在19.3米时有最亮点;在9.33米时有最暗点 2.当h1=5m,3mh29m时: 对h2求偏导,并令其为0: ②运用MATLAB求出极值点 solve(3/((h^2+(20-x)^2)^(3/2))-3*(3*h^2)/((h^2+(20-x)^2)^(5/2))=0) ans = 20+2^(1/2)*h 20-2^(1/2)*h 对x求偏导,并令其为0: 通过MATLAB,将步骤②中计算出的关于h2的表达式带入上式,并求出h2的值; solve(-30*(20-2^(1/2)*h)/((25+(20-2^(1/2)*h)^2)^(5/2))+9*h*(20-(20-2^(1/2)*h))/((h^2+(20-(20-2^(1/2)*h))^2)^(5/2))=0) ans = 7.4223928896768612557104509932965 ⑤通过MATLAB,利用已求得的h2,计算得到x,并进一步计算得到I h=7.42239; x=20-2^(1/2)*h I=10/((25+x^2)^(3/2))+(3*h)/((h^2+(20-x)^2)^(3/2)) x = 9.5032 I = 0.0186 3.当h1,h2均在3m-9m之间时: ①同上,通过MATLAB求解下面的方程组: solve(p1/(h1^2+x^2)^(3/2)-3*p1*h1^2/(h1^2+x^2)^(5/2)) solve(3/((h^2+(20-x)^2)^(3/2))-3*(3*h^2)/((h^2+(20-x)^2)^(5/2))=0) ans = 2^(1/2)*h1 -2^(1/2)*h1 ans = 20+2^(1/2)*h 20-2^(1/2)*h ②根据实际,选

文档评论(0)

baobei + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档