numpy - 如何在 Julia 中找到线性矩阵方程的最小二乘解?
问题描述
我想知道numpy.linalg.lstsq()
Julia 中是否有替代功能。该函数返回线性矩阵方程的最小二乘解。从 NumPy 堆栈访问元素的正确语法?
蟒蛇示例:
import numpy as np
A = np.vstack([x, np.ones(len(x))]).T
m, c = np.linalg.lstsq(A, y, rcond=None)[0]
m, c
(1.0 -0.95) # may vary
输出图表:
Python 参考:https ://numpy.org/doc/stable/reference/generated/numpy.linalg.lstsq.html
解决方案
在这种情况下,您可以\
从 Julia Base 中使用:
julia> A = [rand(10) ones(10)]
10×2 Matrix{Float64}:
0.637746 1.0
0.296172 1.0
0.795938 1.0
0.611058 1.0
0.737017 1.0
0.992014 1.0
0.914031 1.0
0.522682 1.0
0.3607 1.0
0.934141 1.0
julia> y = A * [1, -1] + rand(10) ./ 10
10-element Vector{Float64}:
-0.34049611405598046
-0.6368145973747783
-0.10597203750574954
-0.2950574213524233
-0.19571807260629853
0.020902316863572645
-0.07310077005612584
-0.40758393396440784
-0.6137424837773662
-0.05149257027230776
julia> A \ y
2-element Vector{Float64}:
0.9582937347300398
-0.9216908912065571
推荐阅读
- python - 如何解决 Python 中 PIL 的“必须作为“序列”警告的问题?
- node.js - 是否可以重命名 pm2 正在使用的文件而不必删除它?
- macos - Vagrant 设置错误:执行“VBoxManage”时出错,这是 Vagrant 用于控制 VirtualBox 的 CLI
- leaflet - 当我尝试打开传单地图时,它只是显示为空白
- python - 为什么培训师在教程中培训时不报告评估指标?
- spring-boot - 可以在@JmsListener SpringBoot 中选择性禁用队列消耗吗?
- javascript - 如何在淡入淡出功能完成时创建回调函数?
- python - 如何根据 Python 中的本地日历计算一天中的实际小时数?
- r - 如何预测输出标签而不是概率?
- tinymce - 图像编辑不起作用,高度宽度值在微型编辑器中的编辑器下不起作用