x = np.array([[10,20,90], [9,18,70], [4,12,55],[3,10,50]]) y = np.array([[1], [1], [0], [0]]) 继续给一组x的值,用梯度下降法怎样预测相对于y的值,给出python代码

时间: 2023-11-08 21:04:19 浏览: 121
TXT

Python3X np.load.txt

这是一个二分类问题,可以使用逻辑回归模型进行预测。以下是使用梯度下降法训练逻辑回归模型并进行预测的Python代码: ```python import numpy as np # 定义sigmoid函数 def sigmoid(z): return 1 / (1 + np.exp(-z)) # 定义损失函数 def loss(X, y, w): m = X.shape[0] h = sigmoid(X.dot(w)) J = -1 / m * (y.T.dot(np.log(h)) + (1 - y).T.dot(np.log(1 - h))) return J # 定义梯度函数 def gradient(X, y, w): m = X.shape[0] h = sigmoid(X.dot(w)) grad = 1 / m * X.T.dot(h - y) return grad # 定义梯度下降函数 def gradient_descent(X, y, w_init, lr, num_iter): w = w_init for i in range(num_iter): grad = gradient(X, y, w) w -= lr * grad J = loss(X, y, w) print("Iteration %d | Loss: %f" % (i, J)) return w # 定义预测函数 def predict(X, w): h = sigmoid(X.dot(w)) y_pred = np.round(h) return y_pred # 定义训练数据集和标签 X_train = np.array([[10,20,90], [9,18,70], [4,12,55],[3,10,50]]) y_train = np.array([[1], [1], [0], [0]]) # 添加一列偏置项 X_train = np.hstack((np.ones((X_train.shape[0], 1)), X_train)) # 初始化权重向量 w_init = np.zeros((X_train.shape[1], 1)) # 设置学习率和训练轮数 lr = 0.01 num_iter = 1000 # 训练模型 w = gradient_descent(X_train, y_train, w_init, lr, num_iter) # 预测新数据 X_new = np.array([[1, 5, 30]]) X_new = np.hstack((np.ones((X_new.shape[0], 1)), X_new)) y_pred = predict(X_new, w) print("Predicted label for X_new: %d" % y_pred) ``` 其中,`X_train`是训练数据集,`y_train`是训练标签。在梯度下降函数中,`w_init`是权重向量的初始值,`lr`是学习率,`num_iter`是训练轮数。在预测函数中,`X_new`是新的数据点,`y_pred`是预测的标签。
阅读全文

相关推荐

f = open('G:\jiont\比赛数据2022\charging_data79.csv', encoding='utf-8') data = pd.DataFrame(pd.read_csv(f, encoding='utf-8-sig', low_memory=False)) soc = np.array(data['standard_soc']) # 放电深度DoD current = np.array(data['total_current']) current = [ float(x)/10 for x in current ] all_vol = np.array(data['cell_volt_list']) mileage = np.array(data['mileage']) mileage = [ float(x)/10 for x in mileage ] all_sig_data = cycle_sig(all_vol) all_sig_data = clean_data(all_sig_data) def split_chargedata(chargr_data): a_data = [] all_data = [] for index, m in enumerate(mileage): if index + 1 < len(mileage): if m == mileage[index + 1]: a_data.append(chargr_data[index]) else: a_data.append(chargr_data[index]) all_data.append(a_data) a_data = [] else: all_data.append(a_data) return all_data all_charge_data = split_chargedata(all_sig_data) all_charge_current = split_chargedata(current) all_charge_soc = split_chargedata(soc) dod1 = [] for t in all_charge_soc: dod1.append(t[-1]-t[0]) ind = [] for ind1, t in enumerate(dod1): if t<10: ind.append(ind1) all_charge_data = np.delete(all_charge_data, ind, axis=0) all_charge_current = np.delete(all_charge_current, ind, axis=0) all_charge_soc = np.delete(all_charge_soc, ind, axis=0) ind9 = [5, 13, 25, 35, 47, 55, 81, 84, 86, 88, 89, 92, 94, 101, 111, 115, 116, 126, 157, 162, 167, 174, 180, 198, 200, 216, 237, 245, 261] all_charge_data = np.delete(all_charge_data, ind9, axis=0) all_charge_current = np.delete(all_charge_current, ind9, axis=0) all_charge_soc = np.delete(all_charge_soc, ind9, axis=0)

import numpy as np # 定义参数 n_lags = 31 # 滞后阶数 n_vars = 6 # 变量数量 alpha = 0.05 # 置信水平 # 准备数据 data = np.array([820.95715,819.17877,801.60077,30,26164.9,11351.8], [265.5425,259.05476,257.48619,11.4,12525,4296.5], [696.9681,685.54114,663.32014,47.5,23790.484,8344.8], [4556.1091,440.58995,433.21995,24.6,12931.388,5575.4], [360.08693,353.75386,351.59186,26.9,11944.322,4523], [938.55919,922.25468,894.26468,35.3,27177.893,8287.4], [490.47837,477.35237,385.17474,24.5,14172.1,6650.4], [553.15463,452.35042,425.92277,32.9,16490.17,7795], [740.35759,721.68259,721.68259,15.5,26117.755,7511.7], [1581.99576,1579.50357,1571.23257,65.4,59386.7,15347.2], [1360.91636,1360.20825,1358.11425,66.4,57160.533,8080], [564.06146,560.91611,559.08711,35.2,22361.86,6165.4], [732.17283,727.25063,725.93863,29.7,22177.389,4393.2], [424.12777,424.10579,411.19979,21.6,14691.359,4695.6], [1439.38133,1437.85585,1436.67585,77.3,50123.672,15479], [961.92496,935.21589,931.28189,45.7,28073.9,11273.3], [881.92808,868.65804,832.44504,46.1,27409.15,11224.4], [713.32299,710.75882,707.42682,35.8,24887.111,5164.2], [2657.28891,2599.20299,2515.67859,92,94207.179,19066.4], [420.95033,418.22931,416.80631,25.6,13309.9,7020], [193.92636,193.84936,193.83836,10.9,6133,6139.5], [499.81565,493.73678,485.2468,20.9,13555.897,3412], [951.93942,939.58126,930.049,45.6,27245.608,7752.5], [309.88498,297.05055,295.69055,22.6,11929.038,3903.2], [411.87141,406.63838,389.29638,27.8,12197.085,3834.1], [45.53226,39.24379,55.34631667,7.5,1872.333333,564.3], [532.67282,524.78031,520.89851,24,18041.642,3902], [269.00374,266.96222,211.14422,20.3,7163.069,3515.4], [91.95276,88.77094,85.74583,7.7,1962.8,645.8], [120.60234,116.39872,113.85872,9.8,4227.003,1706.2], [362.98862,350.36495,318.70232,23.7,11615.383,5752.1]) # 计算VAR模型的系数 X = np.zeros((data.shape[0] - n_lags, n_lags * n_vars)) y = np.zeros((data.shape[0] - n_lags, n_vars)) for i in range(n_lags, data.shape[0]): X[i-n_lags, :] = data[i-n_lags:i, :].reshape(1, -1) y[i-n_lags, :] = data[i, :] coefficients = np.linalg.inv(X.T @ X) @ X.T @ y # 计算残差 residuals = y - X @ coefficients # 计算PVAR模型的紧贴矩阵 T = residuals[n_lags:, :] @ residuals[:-n_lags, :].T / (data.shape[0] - n_lags) # 计算PVAR模型的系数 u, s, vh = np.linalg.svd(T) S_inv = np.diag(np.sqrt(s[:n_vars])) @ np.linalg.inv(vh[:n_vars, :]) A = S_inv @ u[:, :n_vars].T @ residuals[n_lags:, :].T # 计算置信区间 t_value = np.abs(np.tinv(alpha/2, data.shape[0]-n_lags-n_vars)) se = np.sqrt((1/(data.shape[0]-n_lags-n_vars)) * (np.sum(residuals[n_lags:, :]**2) / (data.shape[0]-n_lags-n_vars-1))) conf_int = t_value * se print("PVAR模型的系数:\n", A) print("置信区间:[{:.4f}, {:.4f}]".format(A.mean() - conf_int, A.mean() + conf_int))这段代码有什么错误

# -*- coding: utf-8 -*- """ Created on Fri May 26 21:47:54 2023 @author: sleepy student """ import numpy as np jiedian = [] lujing = np.array([[0,1.5,1.6,4,35,44,12,10,24,37,39,33,139,26,16,12,39,49,12,33,28], [1.5,0,3,2.8,33,43,13,12,23,44,37,31,137,25,15,10,37,48,14,31,27], [1.6,3,0,5,33,45,12,10,25,41,36,31,137,25,17,12,37,49,12,31,26], [4,2.8,5,0,34,43,13,13,23,41,38,32,138,26,15,12,38,51,16,30,27], [35,33,33,34,0,43,47,32,33,41,28,7,131,22,25,26,28,50,43,13,13], [44,43,45,43,43,0,54,46,25,73,59,37,108,26,30,34,60,26,57,27,28], [12,13,12,13,47,54,0,17,32,39,36,43,149,36,25,23,42,59,8,41,37], [10,12,10,13,32,46,17,0,32,34,26,31,137,27,25,21,27,58,12,29,22], [24,23,25,23,33,25,32,32,0,64,43,22,120,8,18,20,47,38,41,15,17], [37,44,41,41,41,73,39,34,64,0,29,48,156,53,57,55,29,80,35,46,44], [39,37,36,38,28,59,36,26,43,29,0,24,145,35,44,41,7,72,30,31,28], [33,31,31,32,7,37,43,31,22,48,24,0,122,18,21,22,27,46,39,9,10], [139,137,137,138,131,108,149,137,120,156,145,122,0,117,126,131,145,100,148,115,117], [26,25,25,26,22,26,36,27,8,53,35,18,117,0,11,13,41,35,33,10,9], [16,15,17,15,25,30,25,25,18,57,44,21,126,11,0,11,41,35,33,10,9], [12,10,12,12,26,34,23,21,20,55,41,22,131,13,11,0,18,40,25,21,18], [39,37,37,38,28,60,42,27,47,29,7,27,145,41,41,18,0,69,33,34,33], [49,48,49,51,50,26,59,58,38,80,72,46,100,35,35,40,69,0,60,38,38], [12,14,12,16,43,57,8,12,41,35,30,39,148,33,33,25,33,60,0,36,34], [33,31,31,30,13,27,41,29,15,46,31,9,115,10,10,21,34,38,36,0,5], [28,27,26,27,13,28,37,22,17,44,28,10,117,9,9,18,33,38,34,5,0]]) for i in range(21): lujing[i][i] = 1e9 x = 0 changdu = [] while len(jiedian)<21: min = 1e9 miny = 0 for j in range(21): if lujing[x][j]<min: min = lujing[x][j] miny = j if miny not in jiedian: if x not in jiedian: jiedian.append(x) jiedian.append(miny) changdu.append(lujing[x][miny]) x = miny else: lujing[x][miny] = 1e9 将以上代码转为matlab代码

最新推荐

recommend-type

python 对任意数据和曲线进行拟合并求出函数表达式的三种解决方案

x = np.array([40, 45, 50, 55, 60, 65, 70, 75, 80, 85, 90, 95, 100, 105, 110, 115, 120, 125, 130, 135]) y = ... # 拟合高斯分布 popt, pcov = curve_fit(func, x, y) # 绘制结果 plt.plot(x, y, 's', label=...
recommend-type

基于 .NET 5 + Ant Design Vue 的 Admin Fx.zip

基于 .NET 5 + Ant Design Vue 的 Admin FxColder.Admin.AntdVueWeb后台快速开发框架,.NET5+Ant Design Vue版本代码(GitHub)https://github.com/Coldairarrow/Colder.Admin.AntdVue文档(GitHub)https://github.com/Coldairarrow/Colder.Admin.AntdVue/wiki代码(码云镜像)https ://gitee.com/Coldairarrow/Colder.Admin.AntdVue文档(码云镜像)https://gitee.com/Coldairarrow/Colder.Admin.AntdVue/wikis在线预览地址http://coldairarrow.gitee.io/colder.amin.antdvue.preview.web/ (账号/密码Admin 123456)
recommend-type

Angular实现MarcHayek简历展示应用教程

资源摘要信息:"MarcHayek-CV:我的简历的Angular应用" Angular 应用是一个基于Angular框架开发的前端应用程序。Angular是一个由谷歌(Google)维护和开发的开源前端框架,它使用TypeScript作为主要编程语言,并且是单页面应用程序(SPA)的优秀解决方案。该应用不仅展示了Marc Hayek的个人简历,而且还介绍了如何在本地环境中设置和配置该Angular项目。 知识点详细说明: 1. Angular 应用程序设置: - Angular 应用程序通常依赖于Node.js运行环境,因此首先需要全局安装Node.js包管理器npm。 - 在本案例中,通过npm安装了两个开发工具:bower和gulp。bower是一个前端包管理器,用于管理项目依赖,而gulp则是一个自动化构建工具,用于处理如压缩、编译、单元测试等任务。 2. 本地环境安装步骤: - 安装命令`npm install -g bower`和`npm install --global gulp`用来全局安装这两个工具。 - 使用git命令克隆远程仓库到本地服务器。支持使用SSH方式(`***:marc-hayek/MarcHayek-CV.git`)和HTTPS方式(需要替换为具体用户名,如`git clone ***`)。 3. 配置流程: - 在server文件夹中的config.json文件里,需要添加用户的电子邮件和密码,以便该应用能够通过内置的联系功能发送信息给Marc Hayek。 - 如果想要在本地服务器上运行该应用程序,则需要根据不同的环境配置(开发环境或生产环境)修改config.json文件中的“baseURL”选项。具体而言,开发环境下通常设置为“../build”,生产环境下设置为“../bin”。 4. 使用的技术栈: - JavaScript:虽然没有直接提到,但是由于Angular框架主要是用JavaScript来编写的,因此这是必须理解的核心技术之一。 - TypeScript:Angular使用TypeScript作为开发语言,它是JavaScript的一个超集,添加了静态类型检查等功能。 - Node.js和npm:用于运行JavaScript代码以及管理JavaScript项目的依赖。 - Git:版本控制系统,用于代码的版本管理及协作开发。 5. 关于项目结构: - 该应用的项目文件夹结构可能遵循Angular CLI的典型结构,包含了如下目录:app(存放应用组件)、assets(存放静态资源如图片、样式表等)、environments(存放环境配置文件)、server(存放服务器配置文件如上文的config.json)等。 6. 开发和构建流程: - 开发时,可能会使用Angular CLI来快速生成组件、服务等,并利用热重载等特性进行实时开发。 - 构建应用时,通过gulp等构建工具可以进行代码压缩、ES6转译、单元测试等自动化任务,以确保代码的质量和性能优化。 7. 部署: - 项目最终需要部署到服务器上,配置文件中的“baseURL”选项指明了服务器上的资源基础路径。 8. 关于Git仓库: - 压缩包子文件的名称为MarcHayek-CV-master,表明这是一个使用Git版本控制的仓库,且存在一个名为master的分支,这通常是项目的主分支。 以上知识点围绕Angular应用“MarcHayek-CV:我的简历”的创建、配置、开发、构建及部署流程进行了详细说明,涉及了前端开发中常见的工具、技术及工作流。
recommend-type

管理建模和仿真的文件

管理Boualem Benatallah引用此版本:布阿利姆·贝纳塔拉。管理建模和仿真。约瑟夫-傅立叶大学-格勒诺布尔第一大学,1996年。法语。NNT:电话:00345357HAL ID:电话:00345357https://theses.hal.science/tel-003453572008年12月9日提交HAL是一个多学科的开放存取档案馆,用于存放和传播科学研究论文,无论它们是否被公开。论文可以来自法国或国外的教学和研究机构,也可以来自公共或私人研究中心。L’archive ouverte pluridisciplinaire
recommend-type

深入剖析:内存溢出背后的原因、预防及应急策略(专家版)

![深入剖析:内存溢出背后的原因、预防及应急策略(专家版)](https://d8it4huxumps7.cloudfront.net/uploads/images/65e82a01a4196_dangling_pointer_in_c_2.jpg?d=2000x2000) 参考资源链接:[Net 内存溢出(System.OutOfMemoryException)的常见情况和处理方式总结](https://wenku.csdn.net/doc/6412b784be7fbd1778d4a95f?spm=1055.2635.3001.10343) # 1. 内存溢出的概念及影响 内存溢出,又称
recommend-type

Java中如何对年月日时分秒的日期字符串作如下处理:如何日期分钟介于两个相连的半点之间,就将分钟数调整为前半点

在Java中,你可以使用`java.time`包中的类来处理日期和时间,包括格式化和调整。下面是一个示例,展示了如何根据给定的日期字符串(假设格式为"yyyy-MM-dd HH:mm:ss")进行这样的处理: ```java import java.text.SimpleDateFormat; import java.time.LocalDateTime; import java.time.ZoneId; import java.time.ZonedDateTime; public class Main { public static void main(String[] args
recommend-type

Crossbow Spot最新更新 - 获取Chrome扩展新闻

资源摘要信息:"Crossbow Spot - Latest News Update-crx插件" 该信息是关于一款特定的Google Chrome浏览器扩展程序,名为"Crossbow Spot - Latest News Update"。此插件的目的是帮助用户第一时间获取最新的Crossbow Spot相关信息,它作为一个RSS阅读器,自动聚合并展示Crossbow Spot的最新新闻内容。 从描述中可以提取以下关键知识点: 1. 功能概述: - 扩展程序能让用户领先一步了解Crossbow Spot的最新消息,提供实时更新。 - 它支持自动更新功能,用户不必手动点击即可刷新获取最新资讯。 - 用户界面设计灵活,具有美观的新闻小部件,使得信息的展现既实用又吸引人。 2. 用户体验: - 桌面通知功能,通过Chrome的新通知中心托盘进行实时推送,确保用户不会错过任何重要新闻。 - 提供一个便捷的方式来保持与Crossbow Spot最新动态的同步。 3. 语言支持: - 该插件目前仅支持英语,但开发者已经计划在未来的版本中添加对其他语言的支持。 4. 技术实现: - 此扩展程序是基于RSS Feed实现的,即从Crossbow Spot的RSS源中提取最新新闻。 - 扩展程序利用了Chrome的通知API,以及RSS Feed处理机制来实现新闻的即时推送和展示。 5. 版权与免责声明: - 所有的新闻内容都是通过RSS Feed聚合而来,扩展程序本身不提供原创内容。 - 用户在使用插件时应遵守相关的版权和隐私政策。 6. 安装与使用: - 用户需要从Chrome网上应用店下载.crx格式的插件文件,即Crossbow_Spot_-_Latest_News_Update.crx。 - 安装后,插件会自动运行,并且用户可以对其进行配置以满足个人偏好。 从以上信息可以看出,该扩展程序为那些对Crossbow Spot感兴趣或需要密切跟进其更新的用户提供了一个便捷的解决方案,通过集成RSS源和Chrome通知机制,使得信息获取变得更加高效和及时。这对于需要实时更新信息的用户而言,具有一定的实用价值。同时,插件的未来发展计划中包括了多语言支持,这将使得更多的用户能够使用并从中受益。
recommend-type

"互动学习:行动中的多样性与论文攻读经历"

多样性她- 事实上SCI NCES你的时间表ECOLEDO C Tora SC和NCESPOUR l’Ingén学习互动,互动学习以行动为中心的强化学习学会互动,互动学习,以行动为中心的强化学习计算机科学博士论文于2021年9月28日在Villeneuve d'Asq公开支持马修·瑟林评审团主席法布里斯·勒菲弗尔阿维尼翁大学教授论文指导奥利维尔·皮耶昆谷歌研究教授:智囊团论文联合主任菲利普·普雷教授,大学。里尔/CRISTAL/因里亚报告员奥利维耶·西格德索邦大学报告员卢多维奇·德诺耶教授,Facebook /索邦大学审查员越南圣迈IMT Atlantic高级讲师邀请弗洛里安·斯特鲁布博士,Deepmind对于那些及时看到自己错误的人...3谢谢你首先,我要感谢我的两位博士生导师Olivier和Philippe。奥利维尔,"站在巨人的肩膀上"这句话对你来说完全有意义了。从科学上讲,你知道在这篇论文的(许多)错误中,你是我可以依
recommend-type

【Java内存管理终极指南】:一次性解决内存溢出、泄漏和性能瓶颈

![【Java内存管理终极指南】:一次性解决内存溢出、泄漏和性能瓶颈](https://community.cloudera.com/t5/image/serverpage/image-id/31614iEBC942A7C6D4A6A1/image-size/large?v=v2&px=999) 参考资源链接:[Net 内存溢出(System.OutOfMemoryException)的常见情况和处理方式总结](https://wenku.csdn.net/doc/6412b784be7fbd1778d4a95f?spm=1055.2635.3001.10343) # 1. Java内存模型
recommend-type

c 语言return用法

在C语言中,`return`关键字用于结束函数的执行并返回一个值给函数调用者(如果函数声明了返回类型)。它的基本语法如下: ```c return_type function_name(parameters) { // 函数体内的代码 if (条件) { return value; // 可选的,直接返回一个特定值 } else { // 可能的计算后返回 result = some_computation(); return result; } } ``` 当`return`被执行时,控制权会立即从当前函数转移