| 查看: 682 | 回复: 6 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】请教以下为啥出错？
|
|||
|
新手，正在学习使用matlib，输入以下代码 >> syms d; >> light2=@(x,y,d) 75./(4*pi*((x-d).^2+(y-2).^2+9)); >> eqn=@(d) light2(0,0,d)+light2(0,0,10-d)-light2(5,0,d)-light2(5,0,10-d); >> light2y=@(x,y,d) 75./(4*pi*((x-5).^2+(y-d).^2+9)); >> eqny=@(d) light2y(0,0,d)+light2y(0,0,4-d)-light2y(0,2,d)-light2y(0,2,4-d); >> dint=fzero(eqny,[0 2]) ??? Error using ==> fzero The function values at the interval endpoints must differ in sign. 但是： >> dint=fzero(eqn,[0 2]) dint = 1.4410 我实在搞不懂为什么，请帮忙解答，谢谢谢谢！ |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有5人回复
-
药化及相关博士的申请
已经有3人回复
-
一篇MDPI论文改变了学习工作和生活
已经有4人回复
-
一个化合物的合成路线：CAS:367929-02-0 名称：8β-乙烯基雌二醇
已经有4人回复
-
太白金星有点烦
已经有3人回复
-
中国地质大学（北京）博士招生补录，数理学院材料科学与工程专业和材料与化工专业
已经有3人回复
2楼2009-06-01 22:33:00
sitonyl
铁杆木虫 (文坛精英)
- 应助: 33 (小学生)
- 金币: 12935.9
- 红花: 3
- 帖子: 10448
- 在线: 338.8小时
- 虫号: 722261
- 注册: 2009-03-14
- 性别: GG
- 专业: 导航、制导与传感技术
3楼2009-06-02 08:12:10
★
sunxiao(金币+1,VIP+0):鼓励一下 6-4 02:52
sunxiao(金币+1,VIP+0):鼓励一下 6-4 02:52
|
谢谢楼上指教。继续请教： 这个问题的来源是解答在10*4*3的房间内确定4个100w灯泡在天花板上的安装位置使得亮度最佳。 我的想法是假设4个灯泡安装位置对称，所以可以先确认在x方向的最佳位置，然后再确认y轴方向的最佳位置，这样得到的(x,y,0)就是灯泡的安装位置。 对亮度最佳的理解是所有选择中最暗的点亮度最大。在x轴方向随着灯泡位置的不同会出现两个暗点，即(0,0)和(5,0)，这样可以找到一个零解，x=1.4410. 但y轴方向没有零解，就有问题了！ 是不是我的解题思路就是错的？当然，对matlib不熟也是一个原因，不知道如果把x和y同时考虑的方程怎么写。请各位大侠指教。 |
4楼2009-06-02 16:26:09
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
kuhailangyu(金币+2,VIP+0):谢谢回帖交流 6-5 10:35
小木虫(金币+0.5):给个红包，谢谢回帖交流
kuhailangyu(金币+2,VIP+0):谢谢回帖交流 6-5 10:35
|
X = FZERO(FUN,X0), where X0 is a vector of length 2, assumes X0 is an interval where the sign of FUN(X0(1)) differs from the sign of FUN(X0(2)). An error occurs if this is not true. Calling FZERO with an interval guarantees FZERO will return a value near a point where FUN changes sign. 你给的fzero(eqny,[0 2])中，eqny(0)和eqny(2)结果的符号相同。 应该不同才可以。 |
5楼2009-06-05 09:02:55
6楼2009-06-05 09:07:49
★ ★
sunxiao(金币+2,VIP+0):谢谢分享，欢迎常来 6-6 01:16
sunxiao(金币+2,VIP+0):谢谢分享，欢迎常来 6-6 01:16
|
我放弃了先前的先求x，再求y的想法。但前提“四个灯泡位置对称安放”不变。 以下是新写的程序： % define (d,t) as the position of the bulb syms d; syms t; illm=0;tillm=0;pd=0;pt=0; light4=@(x,y,d,t) 100./(4*pi*((x-d).^2+(y-t).^2+9)); [x,y]=meshgrid(0:0.1:10,0:0.1:4); for d=0:0.1:5 for t=0:0.1:2 tillm=min(light4(x,y,d,t)+light4(x,y,10-d,t)+light4(x,y,d,4-t)+light4(x,y,10-d,4-t)); if illm pd=d;pt=t; end end end 这个是个比较笨的办法，基本就是用了穷举的方法算出来的。有没有更好的办法来解决这个问题，请大家指教。 |
7楼2009-06-05 09:53:06