ValueError: ข้อผิดพลาดของโดเมนคณิตศาสตร์

ฉันเป็นเพียงการทดสอบตัวอย่างจากวิธีเชิงตัวเลขในงานวิศวกรรมกับงูหลาม

from numpy import zeros, array
from math import sin, log
from newtonRaphson2 import *

def f(x):
    f = zeros(len(x))
    f[0] = sin(x[0]) + x[1]**2 + log(x[2]) - 7.0
    f[1] = 3.0*x[0] + 2.0**x[1] - x[2]**3 + 1.0
    f[2] = x[0] + x[1] + x[2] -5.0
    return f

x = array([1.0, 1.0, 1.0])
print newtonRaphson2(f,x)

เมื่อฉันรันมันแสดงข้อผิดพลาดต่อไปนี้:

File "example NR2method.py", line 8, in f
    f[0] = sin(x[0]) + x[1]**2 + log(x[2]) - 7.0
ValueError: math domain error

ฉันได้ จำกัด มันให้แคบลงในบันทึกเมื่อฉันลบบันทึกและเพิ่มฟังก์ชันอื่นมันใช้งานได้ ฉันคิดว่ามันเป็นเพราะการรบกวนบางอย่างกับฐานฉันไม่สามารถเข้าใจได้ว่าเป็นอย่างไร ใครช่วยแนะนำวิธีแก้ไข

รหัสของคุณกำลังทำlogตัวเลขที่น้อยกว่าหรือเท่ากับศูนย์ นั่นไม่ได้กำหนดไว้ในทางคณิตศาสตร์ดังนั้นlogฟังก์ชันของ Python จึงมีข้อยกเว้น นี่คือตัวอย่าง:

>>> from math import log
>>> log(-1)
Traceback (most recent call last):
  File "<pyshell#59>", line 1, in <module>
    log(-1)
ValueError: math domain error

โดยไม่ทราบว่าnewtonRaphson2ฟังก์ชันของคุณทำหน้าที่อะไรฉันไม่แน่ใจว่าฉันสามารถเดาได้ว่าx[2]ค่าที่ไม่ถูกต้องมาจากไหน แต่หวังว่าสิ่งนี้จะนำคุณไปสู่การติดตามที่ถูกต้อง

คุณกำลังพยายามหาลอการิทึมของสิ่งที่ไม่เป็นบวก

ลอการิทึมหาฐานหลังจากได้รับตัวเลขและกำลังเพิ่มขึ้นเป็น log(0)หมายความว่าบางสิ่งบางอย่างยกกำลังของมี2 0เลขชี้กำลังไม่สามารถทำให้เกิด0* ได้ซึ่งหมายความว่าlog(0)ไม่มีคำตอบจึงโยนmath domain error

* หมายเหตุ: 0^0สามารถส่งผลได้0แต่ยังสามารถส่งผลใน1เวลาเดียวกันได้ ปัญหานี้เป็นที่ถกเถียงกันอย่างหนัก

math.log1pนอกจากนี้คุณยังอาจจะใช้

ตามเอกสารอย่างเป็นทางการ :

math.log1p (x)

ส่งคืนลอการิทึมธรรมชาติของ 1 + x (ฐาน e) ผลลัพธ์คำนวณด้วยวิธีที่แม่นยำสำหรับ x ใกล้ศูนย์

คุณสามารถแปลงกลับเป็นค่าเดิมโดยใช้math.expm1ซึ่งส่งคืนeยกกำลัง x ลบ 1

คุณได้รับข้อผิดพลาดโดเมนคณิตศาสตร์จากสาเหตุใดสาเหตุหนึ่ง: คุณกำลังพยายามใช้ตัวเลขติดลบในฟังก์ชันบันทึกหรือค่าศูนย์