主题:[讨论]难道此题无人能解吗?
ygzzxn
[专家分:0] 发布于 2009-06-07 08:04:00
求解n个方程组:
r1+1/(1/r2+1/r3+1/r4+…+1/rn)=a1 (1)
r2+1/(1/r1+1/r3+1/r4+…+1/rn)=a2 (2)
r3+1/(1/r1+1/r2+1/r4+…+1/rn)=a3 (3)
……………………….
rn+1/(1/r1+1/r2+1/r3+…+1/rn-1)=an (n)
其中,n=3,4,5……..50;r1,r2, r3,......rn未知,a1, a2, a3,……an已知。r,a为大于零的实数。
n个方程组对还是初学者的本人而言实在太难,因此先求解三个方程组,先易后难,试图从中找到规律,再解决n个方程组的问题。试着用穷举的方法写了一段代码,仍无法解出,看来代码有问题。不知道有无高手能求出n个方程组或者3 个,4个,5个方程组的解。
三个方程组:
r1+1/(1/r2+1/r3)=a1 (1)
r2+1/(1/r1+1/r3)=a2 (2)
r3+1/(1/r1+1/r2)=a3 (3)
附,求三个方程组的代码:
input “a1,a2,a3=”;a1,a2,a3
for r1=0.1 to a1 step 0.1
for r2=0.1 to a2 step 0.1
for r3=0.1 to a3 step 0.1
if r1+1/(1/r2+1/r3)=a1 and r2+1/(1/r1+1/r3)=a2 and r3+1/(1/r1+1/r2)=a3 then
print “r1=”;r1,”r2=”;r2,”r3=”;r3
exit for
end if
next r3
next r2
next r1
end
四个方程组:
r1+1/(1/r2+1/r3+1/r4)=a1 (1)
r2+1/(1/r1+1/r3+1/r4)=a2 (2)
r3+1/(1/r1+1/r2+1/r4)=a3 (3)
r4+1/(1/r1+1/r2+1/r3)=a4 (4)
五个方程组:
r1+1/(1/r2+1/r3+1/r4+1/r5)=a1 (1)
r2+1/(1/r1+1/r3+1/r4+1/r5)=a2 (2)
r3+1/(1/r1+1/r2+1/r4+1/r5)=a3 (3)
r4+1/(1/r1+1/r2+1/r3+1/r5)=a4 (4)
r5+1/(1/r1+1/r2+1/r3+1/r4)=a5 (5)
回复列表 (共17个回复)
沙发
ygzzxn [专家分:0] 发布于 2009-06-07 08:10:00
如果各位网友觉得求解n个方程组有难度,可尝试解一下3 个,4个,5个方程组。
下面是本人试着用穷举的方法写了一段代码,仍无法解出,看来代码有问题。不知道有无高手能找出下面代码的问题或者求出3 个,4个,5个方程组的解。
三个方程组:
r1+1/(1/r2+1/r3)=a1 (1)
r2+1/(1/r1+1/r3)=a2 (2)
r3+1/(1/r1+1/r2)=a3 (3)
附,求三个方程组的代码:
input “a1,a2,a3=”;a1,a2,a3
for r1=0.1 to a1 step 0.1
for r2=0.1 to a2 step 0.1
for r3=0.1 to a3 step 0.1
if r1+1/(1/r2+1/r3)=a1 and r2+1/(1/r1+1/r3)=a2 and r3+1/(1/r1+1/r2)=a3 then
print “r1=”;r1,”r2=”;r2,”r3=”;r3
exit for
end if
next r3
next r2
next r1
end
四个方程组:
r1+1/(1/r2+1/r3+1/r4)=a1 (1)
r2+1/(1/r1+1/r3+1/r4)=a2 (2)
r3+1/(1/r1+1/r2+1/r4)=a3 (3)
r4+1/(1/r1+1/r2+1/r3)=a4 (4)
五个方程组:
r1+1/(1/r2+1/r3+1/r4+1/r5)=a1 (1)
r2+1/(1/r1+1/r3+1/r4+1/r5)=a2 (2)
r3+1/(1/r1+1/r2+1/r4+1/r5)=a3 (3)
r4+1/(1/r1+1/r2+1/r3+1/r5)=a4 (4)
r5+1/(1/r1+1/r2+1/r3+1/r4)=a5 (5)
板凳
moz [专家分:37620] 发布于 2009-06-09 01:21:00
先从前两行方程出发分析:
r1+1/(1/r2+1/r3.....1/rn)=a1
把r1移到等号右边得:
1/(1/r2+1/r3.....1/rn)=a1-r1
等号两边倒数得:
1/r2+1/r3+......+rn=1/(a1-r1) 以下称之为A式
第二行方程式同理得:
1/r1+1/r3+......+rn=1/(a2-r2) 以下称之为B式
A式 - B式 得:
1/r2 - 1/r1 = 1/(a1-r1) - 1/(a2-r2)
把减数对调得:
1/r2 + 1/(a2-r2) = 1/r1 + 1/(a1-r1)
各自通分得:
(a2-r2+r2)/(r2*(a2-r2)) = (a1-r1+r1)/(r1*(a1-r1))
得:
a2/(r2*(a2-r2))=a1/(r1*(a1-r1))
等号两边倒数得:
r2*r2/a2-r2 = r1*r1/a1-r1
两边各自是一道一元二次方程式。
根据同样的道理,得:
r1*r1/a1-r1 = r2*r2/a2-r2 = r3*r3/a3-r3 = r4*r4/a4-r4 = r5*r5/a5-r5 = .......
随便给r1套一个数,得出某个结果C
再根据C,和每一个一元二次方程式,计算出下一个rn未知数的(大于零的正值才有效)
一元二次方程有固定的求值函数,单独做一个子函数就可以了。
后面的就简单得多了,我就不说了。
另:题目其实第一天就看见了,邮件反而是多天后才看见的,因为本人脑袋大不如前,已经不太愿意花脑汁去思考问题了。请楼主见谅。
3 楼
ygzzxn [专家分:0] 发布于 2009-06-09 15:03:00
a2/(r2*(a2-r2))=a1/(r1*(a1-r1))
等号两边倒数得:
r2*r2/a2-r2 = r1*r1/a1-r1
MOZ版主,以上两步是不是有误?
4 楼
moz [专家分:37620] 发布于 2009-06-09 20:51:00
有错吗?我没看出来。
应该没错吧,
虽然我已经有十多年没接触过代数方程了,
不过应该还没有退化的。
5 楼
天天和和 [专家分:1420] 发布于 2009-06-11 23:11:00
[quote]有错吗?我没看出来。
应该没错吧,
虽然我已经有十多年没接触过代数方程了,
不过应该还没有退化的。[/quote]
这好像有些不负责任说~
6 楼
moz [专家分:37620] 发布于 2009-06-11 23:17:00
是吧?我也觉得是。太不负责任了。
7 楼
ygzzxn [专家分:0] 发布于 2009-06-12 08:23:00
版主:
你好,以下两步能否再详细推算一下,本人数学知识实在太肤浅了,一直没有整明白!见笑了!
a2/(r2*(a2-r2))=a1/(r1*(a1-r1))
等号两边倒数得:
r2*r2/a2-r2 = r1*r1/a1-r1
8 楼
moz [专家分:37620] 发布于 2009-06-12 08:49:00
原式: a2 / ( r2 * ( a2 - r2 )) = a1 / ( r1* ( a1 - r1 ))
倒数: ( r2 * ( a2 - r2 )) / a2 = ( r1 * ( a1 - r1 )) / a1
去括号: ( r2 * a2 - r2 * r2 ) / a2 = ( r1 * a1 - r1 * r1 )) / a1
去括号: r2 * a2 / a2 - r2 * r2 / a2 = r1 * a1 / a1 - r1 * r1 / a1
约分: r2 - r2 * r2 / a2 = r1 - r1 * r1 / a1
两边同时乘以-1
- r2 + r2 * r2 / a2 = - r1 + r1 * r1 / a1
加法交换率 r2 * r2 / a2 - r2 = r1 * r1 / a1 - r1
我有点老眼昏花,眼花缭乱,看得不太清楚,也许哪个地方看错了也不一定,你要提出来告诉我。
9 楼
ygzzxn [专家分:0] 发布于 2009-06-12 14:30:00
版主很是敬业啊!在此谢过!加分!
可是要得出答案,还有一定难度,看还有那位高手能写出代码?
10 楼
moz [专家分:37620] 发布于 2009-06-12 16:43:00
嗯,要代码是吧?
嗯,抱歉。
我来回复