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

关于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 复现论文;matlab代码编写
    • ¥30 echarts 3d地图怎么实现一进来页面散点数据和卡片一起轮播
    • ¥15 数字图像的降噪滤波增强
    • ¥15 心碎了,为啥我的神经网络训练的时候第二个批次反向传播会报错呀,第一个批次都没有问题
    • ¥15 MSR2680-XS路由器频繁卡顿问题
    • ¥15 VB6可以成功读取的文件,用C#读不了
    • ¥15 如何使用micpyhon解析Modbus RTU返回指定站号的湿度值,并确保正确?
    • ¥15 C++ 句柄后台鼠标拖动如何实现
    • ¥15 有人会SIRIUS 5.8.0这个软件吗