使用tensorflow实现矩阵分解方式


Posted in Python onFebruary 07, 2020

采用最小二乘的求逆方法在大部分情况下是低效率的。特别地,当局镇非常大时效率更低。另外一种实现方法是矩阵分解,此方法使用tensorflow内建的Cholesky矩阵分解法。Cholesky矩阵分解法把一个矩阵分解为上三角矩阵和下三角矩阵,L和L'。求解Ax=b,改写成LL'=b。首先求解Ly=b,然后求解L'x=y得到系数矩阵。

1. 导入编程库,初始化计算图,生成数据集。接着获取矩阵A和b。

>>> import matplotlib.pyplot as plt
>>> import numpy as np

>>> import tensorflow as tf

>>> from tensorflow.python.framework import ops
>>> ops.reset_default_graph()

>>> sess=tf.Session()

>>> x_vals=np.linspace(0,10,100)

>>> y_vals=x_vals+np.random.normal(0,1,100)

>>> x_vals_column=np.transpose(np.matrix(x_vals))
>>> ones_column=np.transpose(np.matrix(np.repeat(1,100)))
>>> A=np.column_stack((x_vals_column,ones_column))
>>> b=np.transpose(np.matrix(y_vals))
>>> A_tensor=tf.constant(A)

>>> b_tensor=tf.constant(b)

2. 找到方阵的Cholesky矩阵分解。

注意:tensorflow的cholesky()函数仅仅返回矩阵分解的下三角矩阵,因为上三角矩阵是下三角矩阵的转置矩阵。

>>> tA_A=tf.matmul(tf.transpose(A_tensor),A_tensor)
>>> L=tf.cholesky(tA_A)
>>> tA_b=tf.matmul(tf.transpose(A_tensor),b)
>>> sol1=tf.matrix_solve(L,tA_b)

>>> sol2=tf.matrix_solve(tf.transpose(L),sol1)

3. 抽取系数

>>> solution_eval=sess.run(sol2)
>>> solution_eval
array([[1.01379067],
    [0.02290901]])
>>> slope=solution_eval[0][0]
>>> y_intercept=solution_eval[1][0]
>>> print('slope:'+str(slope))
slope:1.0137906744047482
>>> print('y_intercept:'+str(y_intercept))
y_intercept:0.022909011828880693
>>> best_fit=[]
>>> for i in x_vals:
...  best_fit.append(slope*i+y_intercept)
...
>>> plt.plot(x_vals,y_vals,'o',label='Data')
[<matplotlib.lines.Line2D object at 0x000001E0A58DD9B0>]
>>> plt.plot(x_vals,best_fit,'r-',label='Best fit line',linewidth=3)
[<matplotlib.lines.Line2D object at 0x000001E0A2DFAF98>]
>>> plt.legend(loc='upper left')
<matplotlib.legend.Legend object at 0x000001E0A58F03C8>

>>> plt.show()

使用tensorflow实现矩阵分解方式

以上这篇使用tensorflow实现矩阵分解方式就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持三水点靠木。

Python 相关文章推荐
Python之父谈Python的未来形式
Jul 01 Python
python基础while循环及if判断的实例讲解
Aug 25 Python
python3操作微信itchat实现发送图片
Feb 24 Python
tensorflow构建BP神经网络的方法
Mar 12 Python
Django1.9 加载通过ImageField上传的图片方法
May 25 Python
windows下pycharm安装、创建文件、配置默认模板
Jul 31 Python
Python模拟浏览器上传文件脚本的方法(Multipart/form-data格式)
Oct 22 Python
Pycharm更换python解释器的方法
Oct 29 Python
Python基于百度云文字识别API
Dec 13 Python
python调用百度地图WEB服务API获取地点对应坐标值
Jan 16 Python
TensorFlow2.X使用图片制作简单的数据集训练模型
Apr 08 Python
Django中使用Celery的方法步骤
Dec 07 Python
python匿名函数lambda原理及实例解析
Feb 07 #Python
Python用input输入列表的实例代码
Feb 07 #Python
详解Python 重学requests发起请求的基本方式
Feb 07 #Python
对tensorflow 中tile函数的使用详解
Feb 07 #Python
TensorFlow2.0矩阵与向量的加减乘实例
Feb 07 #Python
Tensorflow分批量读取数据教程
Feb 07 #Python
python统计字符的个数代码实例
Feb 07 #Python
You might like
php array_slice函数的使用以及参数详解
2008/08/30 PHP
php与XML、XSLT、Mysql的结合运用实现代码
2009/11/19 PHP
php中将时间差转换为字符串提示的实现代码
2011/08/08 PHP
php简单判断两个字符串是否相等的方法
2015/07/13 PHP
PHP list() 将数组中的值赋给变量的简单实例
2016/06/13 PHP
PHP注释语法规范与命名规范详解篇
2018/01/21 PHP
PHP7新功能总结
2019/04/14 PHP
JQuery的ajax基础上的超强GridView展示
2009/09/18 Javascript
设置checkbox为只读(readOnly)的两种方式
2013/10/11 Javascript
JS实现往下不断流动网页背景的方法
2015/02/27 Javascript
jquery实现定时自动轮播特效
2015/12/10 Javascript
JS数组去掉重复数据只保留一条的实现代码
2016/08/11 Javascript
JavaScript实现页面定时刷新(定时器,meta)
2016/10/12 Javascript
jQuery动态生成表格及右键菜单功能示例
2017/01/13 Javascript
基于JavaScript实现全选、不选和反选效果
2017/02/15 Javascript
微信小程序 弹框和模态框实现代码
2017/03/10 Javascript
深入理解Node.js中的进程管理
2017/03/13 Javascript
各种选择框jQuery的选中方法(实例讲解)
2017/06/27 jQuery
vue使用微信JS-SDK实现分享功能
2019/08/23 Javascript
Express 配置HTML页面访问的实现
2020/11/01 Javascript
详解Python中的循环语句的用法
2015/04/09 Python
Python之web模板应用
2017/12/26 Python
python实现自动发送邮件
2018/06/20 Python
实践Vim配置python开发环境
2018/07/02 Python
Python面向对象程序设计示例小结
2019/01/30 Python
利用pyecharts读取csv并进行数据统计可视化的实现
2020/04/17 Python
贝斯特韦斯特酒店集团官网:Best Western
2019/01/03 全球购物
简单介绍Object类的功能、常用方法
2013/10/02 面试题
初中班主任经验交流材料
2014/05/16 职场文书
破坏寝室公物检讨书
2014/11/17 职场文书
投标邀请书范本
2015/02/02 职场文书
化验员岗位职责
2015/02/14 职场文书
辩护意见书
2015/06/04 职场文书
创业计划书之淘宝网店
2019/10/08 职场文书
vue自定义右键菜单之全局实现
2022/04/09 Vue.js
Hive HQL支持2种查询语句风格
2022/06/25 数据库