抒情柠檬水 2022-08-15 10:24 采纳率: 0%
浏览 384
已结题

关于inv报错的一些问题

以下是我尝试牛顿迭代法的代码,问题出现在矩阵求逆时给我的报错
弄了好久也不知道怎么回事

报错

img


numpy.core._exceptions._UFuncInputCastingError: Cannot cast ufunc 'inv' input from dtype('O') to dtype('float64') with casting rule 'same_kind'

from sympy import symbols,diff
import numpy as np
def matrixf(f):
    return np.array([diff((f),x),diff((f),y)]).reshape(2,1)
def jcobe(f,arr):
    return np.array([diff(f[0][0],x).evalf(subs={x:arr[0][0],y:arr[1][0]}),diff(f[0][0],y).evalf(subs={x:arr[0][0],y:arr[1][0]}),
diff(f[1][0],x).evalf(subs={x:arr[0][0],y:arr[1][0]}),diff(f[1][0],y).evalf(subs={x:arr[0][0],y:arr[1][0]})]).reshape(2,2)
x,y=symbols('x y')
f=lambda x,y:x*x+y*y-2*x
f=x*x+y*y-2*x
s=matrixf(f)
x0=np.array([0,0]).reshape(2,1)
j=jcobe(s,x0)
x1=x0-np.dot(np.linalg.inv(j),s)
print(x1)
  • 写回答

0条回答 默认 最新

    报告相同问题?

    问题事件

    • 系统已结题 8月23日
    • 修改了问题 8月15日
    • 修改了问题 8月15日
    • 创建了问题 8月15日

    悬赏问题

    • ¥15 如何让企业微信机器人实现消息汇总整合
    • ¥50 关于#ui#的问题:做yolov8的ui界面出现的问题
    • ¥15 如何用Python爬取各高校教师公开的教育和工作经历
    • ¥15 TLE9879QXA40 电机驱动
    • ¥20 对于工程问题的非线性数学模型进行线性化
    • ¥15 Mirare PLUS 进行密钥认证?(详解)
    • ¥15 物体双站RCS和其组成阵列后的双站RCS关系验证
    • ¥20 想用ollama做一个自己的AI数据库
    • ¥15 关于qualoth编辑及缝合服装领子的问题解决方案探寻
    • ¥15 请问怎么才能复现这样的图呀