撞南墙者 2022-09-24 18:23 采纳率: 100%
浏览 22
已结题

numpy中矩阵无法相乘

在使用numpy和sklearn自主实现逻辑回归的过程中,矩阵无法相乘

```python
import sys
import pandas as pd
import numpy as np
import matplotlib.pyplot as plt
from sklearn.datasets import load_breast_cancer
data = load_breast_cancer(return_X_y=True)
X = np.array(data[0])
y = np.array(data[1])

def sigmod(x):
    return 1/(1+pow(np.e,(-x)))

def Logistic_Regression(feature_data,target_data,learning_rate,account):
    m = feature_data.shape[0]
    feature_data = np.hstack((np.full((m,1),1),feature_data))
    m,n = feature_data.shape
    para = np.random.uniform(-1,1,n).reshape(n,1) # n*1
    para = np.mat(para)
    feature_data = np.mat(feature_data)    # m*n
    target_data = np.mat(target_data)
    # =========问题代码,两矩阵阵无法相乘????============#
    d = pd.DataFrame(np.array(feature_data))
    print(d.describe())
    print(type(feature_data), type(para))
    print(feature_data.shape, para.shape)
    print(feature_data @ para)

    Error = (-1/m)*np.sum(np.multiply(target_data,np.log(sigmod(feature_data @ para))) +
                    np.multiply(1-target_data,np.log(1-sigmod(feature_data @ para))))
    count = 1
    error_list = [Error]

    while True:
        grad_vector = (1/m) * [feature_data.T @ (sigmod(feature_data @ para) - target_data)]
        para = para - learning_rate * grad_vector
        Error = (-1 / m) * np.sum(np.multiply(target_data, np.log(sigmod(feature_data @ para))) +
                                np.multiply(1 - target_data, np.log(1 - sigmod(feature_data @ para))))
        error_list.append(Error)
        count = count + 1
        if count == account:
            break
    return para,error_list

para,e = Logistic_Regression(X,y,0.01,100)
plt.plot(e)
plt.show()




###### 无显式报错,但程序中止

###### 试过更改库版本,无果
  • 写回答

1条回答 默认 最新

  • ·星辰大海 2022-09-24 19:16
    关注

    np包的矩阵相乘有两种你需要逐个分清,dot和multiply。
    如果是矩阵的相乘的话需要检查两个矩阵的长和宽是否在每一次的迭代中都符合要求,同时也要避免矩阵中出现None值

    本回答被题主选为最佳回答 , 对您是否有帮助呢?
    评论

报告相同问题?

问题事件

  • 系统已结题 1月8日
  • 已采纳回答 12月31日
  • 创建了问题 9月24日

悬赏问题

  • ¥15 R语言Rstudio突然无法启动
  • ¥15 关于#matlab#的问题:提取2个图像的变量作为另外一个图像像元的移动量,计算新的位置创建新的图像并提取第二个图像的变量到新的图像
  • ¥15 改算法,照着压缩包里边,参考其他代码封装的格式 写到main函数里
  • ¥15 用windows做服务的同志有吗
  • ¥60 求一个简单的网页(标签-安全|关键词-上传)
  • ¥35 lstm时间序列共享单车预测,loss值优化,参数优化算法
  • ¥15 Python中的request,如何使用ssr节点,通过代理requests网页。本人在泰国,需要用大陆ip才能玩网页游戏,合法合规。
  • ¥100 为什么这个恒流源电路不能恒流?
  • ¥15 有偿求跨组件数据流路径图
  • ¥15 写一个方法checkPerson,入参实体类Person,出参布尔值