回 帖 发 新 帖 刷新版面

主题:[讨论]难道此题无人能解吗?

求解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个回复)

沙发


如果各位网友觉得求解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)

 

板凳

先从前两行方程出发分析:
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 楼


a2/(r2*(a2-r2))=a1/(r1*(a1-r1))
等号两边倒数得:
r2*r2/a2-r2 = r1*r1/a1-r1
MOZ版主,以上两步是不是有误?

4 楼

有错吗?我没看出来。
应该没错吧,
虽然我已经有十多年没接触过代数方程了,
不过应该还没有退化的。

5 楼

[quote]有错吗?我没看出来。
应该没错吧,
虽然我已经有十多年没接触过代数方程了,
不过应该还没有退化的。[/quote]




这好像有些不负责任说~

6 楼

是吧?我也觉得是。太不负责任了。

7 楼

版主:
  你好,以下两步能否再详细推算一下,本人数学知识实在太肤浅了,一直没有整明白!见笑了!
a2/(r2*(a2-r2))=a1/(r1*(a1-r1))
等号两边倒数得:
r2*r2/a2-r2 = r1*r1/a1-r1

8 楼

原式:      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 楼

版主很是敬业啊!在此谢过!加分!
可是要得出答案,还有一定难度,看还有那位高手能写出代码?

10 楼

嗯,要代码是吧?
嗯,抱歉。

我来回复

您尚未登录,请登录后再回复。点此登录或注册