我是python中的新手,我需要一些帮助。
我有一个python脚本来查找圆心的X坐标。如果r1很小,则xA1和xA2是复数(表示错误消息:ValueError:math domain error)。我的问题是,如何将此错误消息用于while语句?我必须复制while语句下面的整个代码吗?是否可以链接到第一个等式xs2 = r3 * math.sin(beta2r)?
感谢您的回答。
罗马
import math
r1=0;
r2=110.5;
r3=212.5;
beta2d=35;
beta1dpoz=45;
beta2r=beta2d*math.pi/180
xs2=r3*math.sin(beta2r)
ys2=r1-r3*math.cos(beta2r)
kruh1=math.pow(xs2,2)+math.pow(ys2,2)
kruh2=math.pow(r1,2)+math.pow(r3,2)-2*r1*r3*math.cos(beta2r)
xA1=1.0/2*((xs2*(1+((math.pow(r1,2)-math.pow(r2,2))/(math.pow(xs2,2)+math.pow(ys2,2)))))+(ys2*(math.sqrt(((2*(math.pow(r1,2)+math.pow(r2,2)))/(math.pow(xs2,2)+math.pow(ys2,2)))-(math.pow(((math.pow(r1,2)-math.pow(r2,2))/(math.pow(xs2,2)+math.pow(ys2,2))),2))-1))))
xA2=1.0/2*((xs2*(1+((math.pow(r1,2)-math.pow(r2,2))/(math.pow(xs2,2)+math.pow(ys2,2)))))-(ys2*(math.sqrt(((2*(math.pow(r1,2)+math.pow(r2,2)))/(math.pow(xs2,2)+math.pow(ys2,2)))-(math.pow(((math.pow(r1,2)-math.pow(r2,2))/(math.pow(xs2,2)+math.pow(ys2,2))),2))-1))))
xAmat=[xA1,xA2]
xA=max(xAmat)
while ('ValueError: math domain error'):
r1=r1+0.1
.
.
.
xA=
答案 0 :(得分:1)
这将增加r1
,直到xA1和xA2有效。
import math
r1=0;
r2=110.5; #vstupny polomer
r3=212.5; #vystupny polomer
beta2d=35; #vystupny uhel
beta1dpoz=45; #vstupny uhel
xA1, xA2 = None, None
#Keep looping till xA1 and xA2 are valid
while xA1 == None and xA2 == None:
beta2r=beta2d*math.pi/180
xs2=r3*math.sin(beta2r)
ys2=r1-r3*math.cos(beta2r)
kruh1=math.pow(xs2,2)+math.pow(ys2,2)
kruh2=math.pow(r1,2)+math.pow(r3,2)-2*r1*r3*math.cos(beta2r)
try:
xA1=1.0/2*((xs2*(1+((math.pow(r1,2)-math.pow(r2,2))/(math.pow(xs2,2)+math.pow(ys2,2)))))+(ys2*(math.sqrt(((2*(math.pow(r1,2)+math.pow(r2,2)))/(math.pow(xs2,2)+math.pow(ys2,2)))-(math.pow(((math.pow(r1,2)-math.pow(r2,2))/(math.pow(xs2,2)+math.pow(ys2,2))),2))-1))))
xA2=1.0/2*((xs2*(1+((math.pow(r1,2)-math.pow(r2,2))/(math.pow(xs2,2)+math.pow(ys2,2)))))-(ys2*(math.sqrt(((2*(math.pow(r1,2)+math.pow(r2,2)))/(math.pow(xs2,2)+math.pow(ys2,2)))-(math.pow(((math.pow(r1,2)-math.pow(r2,2))/(math.pow(xs2,2)+math.pow(ys2,2))),2))-1))))
except ValueError:
r1 += 0.1
xAmat=[xA1,xA2]
xA=max(xAmat)
如果任何其他计算(beta2r, ys2, xs2, kruh1, kruh1
)可能产生ValueError,则将try
语句向上移动以包含这些语句。
7.4. The try statement