机器学习 - 从 Octave 代码在 Python 中实现梯度下降

如何解决机器学习 - 从 Octave 代码在 Python 中实现梯度下降

我正在尝试从头开始在 Python 中实现一个梯度下降函数,我已经在 GNU Octave 中实现并工作了该函数。不幸的是我被卡住了。我摆弄了一段时间并检查了 NumPy 文档,但到目前为止还没有运气。 我知道诸如 scikit-learn 之类的库,但是我的目的是学习从头开始编写这样的函数。也许我的做法是错误的。 您将在下面找到重现错误所需的所有代码。 预先感谢您的帮助。

实际结果:测试失败并出现错误 ->“ValueError:matmul:输入操作数 0 没有足够的维度(有 0,gufunc 核心带有签名 (n?,k),(k,m?)->(n?,m?) 需要 1)"

预期结果:和值为 [5.2148,-0.5733] 的数组

Octave 中的函数 gradientDescent():

function [theta,J_history] = gradientDescent(X,y,theta,alpha,num_iters)
        m = length(y); % number of training examples
        J_history = zeros(num_iters,1);

        for iter = 1:num_iters
          theta = theta - (alpha/m)*X'*(X*theta-y);
          J_history(iter) = computeCost(X,theta);
        end

python 中的函数gradient_descent():

from numpy import zeros

def compute_cost(X,theta):
    m = len(y)
    ans = (X.T @ theta).T - y
    J = (ans @ ans.T) / (2 * m)
    return J[0,0]


def gradient_descent(X,num_iters):
    m = len(y)
    J_history = zeros((num_iters,1),dtype=int)
    for iter in range(num_iters):
        theta = theta - (alpha / m) @ X.T @ (X @ theta - y)
        J_history[iter] = compute_cost(X,theta)
    return theta


测试文件:test_ml_utils.py

import unittest
import numpy as np
from ml.ml_utils import compute_cost,gradient_descent

class TestGradientDescent(unittest.TestCase):
    # TODO: implement tests for Gradient Descent function
    # [theta J_hist] = gradientDescent([1 5; 1 2; 1 4; 1 5],[1 6 4 2]',[0 0]',0.01,1000);
    def test_gradient_descent_00(self):
        X = np.array([[1,5],[1,2],4],5]])
        y = np.array([1,6,4,2])
        theta = np.zeros(2)
        alpha = 0.01
        num_iter = 1000
        r_theta = np.array([5.2148,-0.5733])
        result = gradient_descent(X,num_iter)
        self.assertEqual((round(result,4),r_theta),'Result is wrong!')


if __name__ == '__main__':
    unittest.main()

解决方法

Python 中的 __matmul__ 运算符 @- 绑定得更紧密。这意味着您正在尝试使用操作数 (alpha / m)(标量)和 X.T(实际上是矩阵)进行矩阵乘法。见operator precedence

在 Octave 代码中,(alpha - m) * X' 进行标量乘法,而不是矩阵,因此如果您希望在 Python 中具有相同的行为,请使用 * 而不是 @。这似乎是因为如果一个操作数是标量,Octave 会重载 * 运算符以执行标量乘法,但如果两个操作数都是矩阵,则执行矩阵乘法。

,

添加到 Adam 的答案(对于您遇到的错误是正确的)。

然而,我想更一般地补充一点,如果没有某种提示(无论是以编程方式还是以注释的形式)不同变量采用的维度,这段代码对读者来说毫无意义。

例如,代码中提示 y 可能是二维的,而您正在使用 len 获取其大小。作为这可能如何无声地失败的示例,请考虑:

>>> y = numpy.array([[1,2,3,4,5]])
>>> len( y )
1

而大概你想要

>>> numpy.shape( y )
(1,5)

>>> numpy.size( y )
5

我在您的单元测试中注意到您传递的是 1 级向量而不是 2 级向量,因此结果 y 是 1D 而不是 2D,但由于广播而使用 2D 的 X 进行操作。因此,尽管隐含的逻辑,您的代码仍然可以工作,但是如果没有明确说明这些事情,这是一个等待发生的运行时错误。

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。

相关推荐


依赖报错 idea导入项目后依赖报错,解决方案:https://blog.csdn.net/weixin_42420249/article/details/81191861 依赖版本报错:更换其他版本 无法下载依赖可参考:https://blog.csdn.net/weixin_42628809/a
错误1:代码生成器依赖和mybatis依赖冲突 启动项目时报错如下 2021-12-03 13:33:33.927 ERROR 7228 [ main] o.s.b.d.LoggingFailureAnalysisReporter : *************************** APPL
错误1:gradle项目控制台输出为乱码 # 解决方案:https://blog.csdn.net/weixin_43501566/article/details/112482302 # 在gradle-wrapper.properties 添加以下内容 org.gradle.jvmargs=-Df
错误还原:在查询的过程中,传入的workType为0时,该条件不起作用 <select id="xxx"> SELECT di.id, di.name, di.work_type, di.updated... <where> <if test=&qu
报错如下,gcc版本太低 ^ server.c:5346:31: 错误:‘struct redisServer’没有名为‘server_cpulist’的成员 redisSetCpuAffinity(server.server_cpulist); ^ server.c: 在函数‘hasActiveC
解决方案1 1、改项目中.idea/workspace.xml配置文件,增加dynamic.classpath参数 2、搜索PropertiesComponent,添加如下 <property name="dynamic.classpath" value="tru
删除根组件app.vue中的默认代码后报错:Module Error (from ./node_modules/eslint-loader/index.js): 解决方案:关闭ESlint代码检测,在项目根目录创建vue.config.js,在文件中添加 module.exports = { lin
查看spark默认的python版本 [root@master day27]# pyspark /home/software/spark-2.3.4-bin-hadoop2.7/conf/spark-env.sh: line 2: /usr/local/hadoop/bin/hadoop: No s
使用本地python环境可以成功执行 import pandas as pd import matplotlib.pyplot as plt # 设置字体 plt.rcParams['font.sans-serif'] = ['SimHei'] # 能正确显示负号 p
错误1:Request method ‘DELETE‘ not supported 错误还原:controller层有一个接口,访问该接口时报错:Request method ‘DELETE‘ not supported 错误原因:没有接收到前端传入的参数,修改为如下 参考 错误2:cannot r
错误1:启动docker镜像时报错:Error response from daemon: driver failed programming external connectivity on endpoint quirky_allen 解决方法:重启docker -> systemctl r
错误1:private field ‘xxx‘ is never assigned 按Altʾnter快捷键,选择第2项 参考:https://blog.csdn.net/shi_hong_fei_hei/article/details/88814070 错误2:启动时报错,不能找到主启动类 #
报错如下,通过源不能下载,最后警告pip需升级版本 Requirement already satisfied: pip in c:\users\ychen\appdata\local\programs\python\python310\lib\site-packages (22.0.4) Coll
错误1:maven打包报错 错误还原:使用maven打包项目时报错如下 [ERROR] Failed to execute goal org.apache.maven.plugins:maven-resources-plugin:3.2.0:resources (default-resources)
错误1:服务调用时报错 服务消费者模块assess通过openFeign调用服务提供者模块hires 如下为服务提供者模块hires的控制层接口 @RestController @RequestMapping("/hires") public class FeignControl
错误1:运行项目后报如下错误 解决方案 报错2:Failed to execute goal org.apache.maven.plugins:maven-compiler-plugin:3.8.1:compile (default-compile) on project sb 解决方案:在pom.
参考 错误原因 过滤器或拦截器在生效时,redisTemplate还没有注入 解决方案:在注入容器时就生效 @Component //项目运行时就注入Spring容器 public class RedisBean { @Resource private RedisTemplate<String
使用vite构建项目报错 C:\Users\ychen\work>npm init @vitejs/app @vitejs/create-app is deprecated, use npm init vite instead C:\Users\ychen\AppData\Local\npm-