做app和做网站相同和区别,这么制作自己的网站,网站建设实践论文,签证网站建设文章目录 0 赛题思路1 问题要求2 假设约定3 符号约定4 建立模型5 模型求解6 实现代码 0 赛题思路
#xff08;赛题出来以后第一时间在CSDN分享#xff09;
https://blog.csdn.net/dc_sinor?typeblog
1 问题要求
现在已知一个教室长为15米#xff0c;宽为12米#xff0… 文章目录 0 赛题思路1 问题要求2 假设约定3 符号约定4 建立模型5 模型求解6 实现代码 0 赛题思路
赛题出来以后第一时间在CSDN分享
https://blog.csdn.net/dc_sinor?typeblog
1 问题要求
现在已知一个教室长为15米宽为12米在距离地面高2.5米的位置均 匀的安放4个光源分别为1、2、3、4各个光源的光照强度均为一个单位如下图 要求
1如何计算教室内任意一点的光照强度光源对目标点的光照强度与该光源到目标点距离的平方成反比与该光源的强度成正比.2画出距离地面1米处各个点的光照强度与位置横纵坐标之间的函数关系曲面图试同时给出一个近似的函数关系式.3假设离地面1米高正是学生桌面的高度如何设计这四个点光源的位置才能使学生对光照的平均满意度达到最高4若将题目中的点光源换成线光源以上(2)、(3)问的结果又如何
对于1、2问假设横向纵向墙壁与光源、光源与光源、光源与墙壁之间的距离是相等的.
2 假设约定
1 光不会通过窗、门等外涉也不考虑光在空气中的消耗即光照强度和不变2 室内不受外界光源影响3 教室高度为2.5米4 不考虑光的反射5 线光源发光是均匀的.
3 符号约定 4 建立模型 5 模型求解 6 实现代码
matlab 实现代码 建议最好用python去实现图会好看一些而且国内当前趋势会逐渐淘汰matlab目前有些学校已经无法使用matlab了
clear
clc
max0;min4;
for i0:0.1:3for j0.1:0.1:4s0;x18i,y15-jx28i,y210jx34-i,y310jx44-i,y45-j for x0:0.1:12for y0:0.1:15for z0:0.1:2.5if x1~x y1~y x2~x y2~y x3~x y3~y x4~x y4~y ss1./((x1-x).^2(y1-y).^2(2.5-z).^2)1./((x2-x).^2(y2-y).^2(2.5-z).^2)1./((x3-x).^2(y3-y).^2(2.5-z).^2)1./((x4-x).^2(y4-y).^2(2.5-z).^2);endendendendk4./s;l0;z1;for x0:0.1:12for y0:0.1:15llk.*(1./((x1-x).^2(y1-y).^2(2.5-z).^2)1/((x2-x).^2(y2-y).^2(2.5-z).^2)1./((x3-x).^2(y3-y).^2(2.5-z).^2)1./((x4-x).^2(y4-y).^2(2.5-z).^2));endendif lmaxmaxl;x11x1;y11y1;x12x2;y12y2;x13x3;y13y3;x14x4;y14y4;endpl./(120.*150);Q0;for x0:0.1:12for y0:0.1:15QQ(k.*(1./((x1-x).^2(y1-y).^2(2.5-z).^2)1./((x2-x).^2(y2-y).^2(2.5-z).^2)1./((x3-x).^2(y3-y).^2(2.5-z).^2)1./((x4-x).^2(y4-y).^2(2.5-z).^2))-p).^2.^(1./2);endendif minQminQ;x21x1;y21y1;x22x2;y22y2;x23x3;y23y3;x24x4;y24y4;endend
end
disp([最大值,x11,num2str(x11), ,y11,num2str(y11), ,x12,num2str(x12), ,y12,num2str(y12), ,x13,num2str(x13), ,y13,num2str(y13), ,x14,num2str(x14), ,y14,num2str(y14)])
disp([最平均,x21,num2str(x21), ,y21,num2str(y21), ,x22,num2str(x22), ,y22,num2str(y22), ,x23,num2str(x23), ,y23,num2str(y23), ,x24,num2str(x24), ,y24,num2str(y24)])
附录二
clear
clc
max0;min4;li4;
for i0:0.1:3for j0.1:0.1:4s0;x18i,y15-jx28i,y210jx34-i,y310jx44-i,y45-j for x0:0.1:12for y0:0.1:15for z0:0.1:2.5if x1~x y1~y x2~x y2~y x3~x y3~y x4~x y4~y ss1./((x1-x).^2(y1-y).^2(2.5-z).^2)1./((x2-x).^2(y2-y).^2(2.5-z).^2)1./((x3-x).^2(y3-y).^2(2.5-z).^2)1./((x4-x).^2(y4-y).^2(2.5-z).^2);endendendendk4./s;l0;z1;e0for x0:0.1:12for y0:0.1:15llk.*(1./((x1-x).^2(y1-y).^2(2.5-z).^2)1/((x2-x).^2(y2-y).^2(2.5-z).^2)1./((x3-x).^2(y3-y).^2(2.5-z).^2)1./((x4-x).^2(y4-y).^2(2.5-z).^2));rk.*(1./((x1-x).^2(y1-y).^2(2.5-z).^2)1/((x2-x).^2(y2-y).^2(2.5-z).^2)1./((x3-x).^2(y3-y).^2(2.5-z).^2)1./((x4-x).^2(y4-y).^2(2.5-z).^2));ee(r-6*10^(-32))^2;endendS(l-0.1278)^2eif SliliSx11x1,y11y1, x12x2,y12y2, x13x3,y13y3, x14x4,y14y4,en4en4
en4
disp([x11,num2str(x11), ,y11,num2str(y11), ,x12,num2str(x12), ,y12,num2str(y12), ,x13,num2str(x13), ,y13,num2str(y13), ,x14,num2str(x14), ,y14,num2str(y14)])
li
