Julia -对从另一个函数返回的矩阵求幂

文章/答案/技术大牛

发布

1回答

julia

我有返回矩阵m1和m2的函数f1和f2，这是使用LinearAlgebra包中的Diagonal, Tridiagonal, SymTridiagonal计算得出的。在一个新的函数f3中，我尝试计算 j = m1 - m2*im m3 = exp(j) 但除非我使用j=Matrix(m1-m2*im)，否则这会给出计算上的Method error，即no matching我的问题是，如何才能以最优的

浏览 19提问于2020-12-02得票数 3

回答已采纳

1回答

函数和曲面图Julia

function、plot、julia、surface

我想在Julia中使用surface绘制一个函数。我设法绘制出所需的函数： x = 0:0.1:4f(x,y) = x^0.2 * y^0.8 surface(x, y, f, camera=(10,30),linealpha=0.3, fc=:heat) 然而，我希望f(*)是一个适当的函数，在这个函数上我也可以进行优化(例如经济学中的效用最大化)。这是我的尝试： function Utility(x1,

浏览 21提问于2020-10-23得票数 1

回答已采纳

4回答

使用Numpy进行矩阵乘法

python、numpy、matrix-multiplication

假设我有一个亲和力矩阵A和一个对角矩阵D，我如何使用nympy在Python中计算拉普拉斯矩阵？这条路对吗？谢谢。

浏览 2提问于2010-08-27得票数 5

回答已采纳

2回答

有没有办法用C创建一个NumPy矩阵？

python、numpy、python-c-api

我阅读了我能找到的关于NumPy C应用程序接口的文档，但仍然不能找出是否有可能用C应用程序接口构造矩阵对象-而不是二维数组。该函数用于处理数学矩阵，如果用户调用矩阵乘法时忘记将该值从数组转换为矩阵(乘法和求幂是矩阵子类唯一的区别)，我不希望出现奇怪的结果。

浏览 4提问于2009-02-21得票数 3

回答已采纳

1回答

在perl中对矩阵求幂

perl、matrix、exponent

我正在寻找一个使用perl对矩阵求幂的函数。我搜索了CPAN，但似乎找不到任何合适的东西。有没有人知道一个合适的包，或者简单的函数？例如，对于以下矩阵($m)： [0.13385,-0.70885,0.16502,0.40998],

浏览 1提问于2012-02-24得票数 2

2回答

boost c++子矩阵和向量类型的逐个元素运算

c++、boost、vector、matrix、linear-algebra

我想在boost矩阵和向量类型上执行按元素的函数，例如，取每个元素的对数，对每个元素求幂，应用特殊函数，如gamma和digamma等(类似于matlab对应用于矩阵和向量的这些函数的处理)。我想为每个想要的函数写一个辅助函数就足够了，但是这看起来很浪费。同样，boost wiki为提供了一些代码，但这似乎相当复杂。建议使用v

浏览 5提问于2009-04-21得票数 7

回答已采纳

1回答

Numpy将两个数组追加在一起

python、arrays、numpy、dimension

我正在尝试创建一个函数，它对2-D矩阵求幂，并将结果保存在3D数组中，其中第一维是指数。这一点很重要，因为我求幂的矩阵的行代表了图上不同顶点的信息。例如，如果我们有A，A^2，A^3，每个都是shape (50,50)，我想要一个矩阵D= (3,50,50)，这样我就可以使用D:,1，：来检索关于节点1的所有信息，并能够对其进行矩阵乘法。我<em

浏览 27提问于2019-10-10得票数 0

1回答

如何在Python中以最快的方式计算矩阵指数？

python、numpy、matrix、scipy、exponential

我想用Python计算矩阵指数。我找到了一种方法("scipy.linalg.expm()")，但它需要很长时间(例如，5000×5000矩阵需要30秒)。matrix_exponential = scipy.linalg.expm(matrix) 在Python中有没有更快的计算矩阵指数的方法？非常提前感谢您。

浏览 103提问于2019-06-13得票数 1

2回答

求非常大N的Fibonacci的第N项

java、fibonacci、largenumber

我已经使用了迭代方法和递归方法，但在这两种情况下，我都没有得到fibo(10^6)更快的结果。复杂度为O(logN)。

浏览 0提问于2013-02-05得票数 1

1回答

svds对某些矩阵不起作用-错误结果

matlab、linear-algebra

下面是我的测试函数： [U,S,V]diff_mat = y - U*S*V';这里有两个非常相似的部分:一个是求随机矩阵的奇异值分解，另一个是求随机稀疏矩阵

浏览 0提问于2013-04-06得票数 1

回答已采纳

1回答

如何对巨型数字求幂和求模

python-3.x、pow

我想把一个数字求幂，然后将它与另一个数字相乘，最后对其进行mod。 Python的pow函数接受3个参数，不能用于我的用例场景，所以我正在尝试寻找快速的替代方案。示例：我知道我能做到 9^10%2 作为 pow(9, 10, 2) 但我做不到 (8*9^10)%2 具有相同的功能我已经尝试过了(数字仅供参考，实数大得多) pow(9, 10)*8 % 2 #tooks too long pow(9, 10, 2)*8 # ret

浏览 16提问于2019-04-08得票数 1

回答已采纳

1回答

示例1-7来自C编程语言(版本2)

我在理解C编程语言给出的示例1.7时遇到了困难。本书将下面的程序描述为这样的程序：“由于C没有像Fortran中的**那样的求幂运算符，让我们通过编写一个函数的幂(m，n)来说明函数定义的机制。也就是说，幂(2，5)的值是32。这个函数不是一个实用的求幂例程，因为它只处理小整数的</e

浏览 16提问于2017-08-14得票数 0

回答已采纳

1回答

具有A是稀疏矩阵且b是朱莉娅中一个稠密向量的线性系统Ax=b的求解

julia、linear-algebra、sparse-matrix

我有一个矩阵A，它是一个以CSC格式存储的稀疏矩阵和一个密集向量b。并且我想解决Ax = b的x系统。我该怎么做呢？有人建议另一个方法吗？使用Julia v1.7.1。另外，矩阵A不是三对角的。请查看的问题以获得错误信息。

浏览 6提问于2022-01-10得票数 1

1回答

更好的算法，使用vcat()或hcat()在一个循环在朱莉娅？

arrays、algorithm、matrix、julia

尽管所有的子矩阵都是以相同的方式创建的，但我可能需要在循环之外构造这个初始矩阵a。如果使用另一个变量来引用有用的行和列，则可能会浪费一些内存分配。如果n有整数平方根，则foo(n)按sqrt(n) * 20矩阵返回

浏览 2提问于2017-01-10得票数 3

回答已采纳

1回答

使用fft对图像进行滤波

image、fft、convolution

我有滤波器矩阵(卷积核等)，但如果我想使用fft，我需要将滤波器和图像矩阵扩展到2的平方矩阵的幂。首先，我用零展开它(所以图片和滤波器的矩阵应该在2的幂矩阵的左上角)，我得到了这个结果然后，我尝试了在不使用fft的情况下应用卷积时对图片矩阵所做的操作，填充图片矩阵的</em

浏览 3提问于2013-01-15得票数 0

2回答

Clojure中的任意精度求幂

math、clojure

有没有办法在Clojure中执行任意精度的求幂？我已经尝试过来自clojure.math.numeric- the的Math/pow和expt函数，但它们都只能返回有限的精度。例如：=> 1.4142135623730951 如何获取更多的数字？

浏览 0提问于2012-12-19得票数 10

回答已采纳

1回答

@NLConstraint和矢量化约束跳转/Julia

optimization、julia-jump

返回的错误为： LoadError: exp is not defined for type Variable.所以一个简单的解决方案是在@NLconstraint函数之外简单地做所有的求幂和求和，这样输入就是一个标量。我想做的另一件事是泰勒级数展开，但这也带来了一个问题，因为它进入了幂大于2的@NLconstraint域，然后我遇到了同样的矢量化问题。所以我感觉卡住了，我觉得如果JuMP允许像@constraint

浏览 0提问于2018-06-16得票数 0

2回答

具有动态规划的所有对最短路径

algorithm、minimum-spanning-tree

全,考虑了加权邻接矩阵本身的乘法，但在这种情况下，我们用加法代替矩阵乘法运算，用极小化代替加法运算。注意，加权邻接矩阵本身的乘积返回一个矩阵，该矩阵包含任意一对节点之间长度为2的最短路径。由此得出的结论是，A到n的幂包含所有最短路径。k≤n {dik(

浏览 7提问于2011-12-05得票数 3

回答已采纳

10回答

Python中的递归函数

python、recursion

1 return fibonacci(number-1) + fibonacci(number-2) 这根据斐波那契级数的(n-1) + (n-2)函数是有意义的。Python如何执行包含另一个递归的递归，而不是在相同的代码行中？'finobacci(number-1)‘是否完成所有的递归，直到它达到'1’，然后它用'fibonacci(number-2)‘做同样的事情并将它们相加？为了便于比较，下面的递

浏览 2提问于2012-12-05得票数 8

回答已采纳

2回答

Julia中不同于转置的Singleton向量

vector、julia、matrix-multiplication

如果我定义A = [1]，我会得到A不等于A'，因为它们属于不同类型： julia> A=[1] 1 1×1 LinearAlgebra.Adjoint{Int64,Array{Int64,1}}:false 如果我定义另一个向量B = [1, 2, 3]，并尝试用A'和A做乘积，我会得到以下输出： B=[1,2,3这似乎是一个*运算符签名的<

浏览 14提问于2020-04-08得票数 3

回答已采纳

点击加载更多