由于某种原因，我不知道，这似乎只通过外部循环一次，有人愿意帮助吗？以下是我到目前为止所拥有的代码：
echelon <- function(mat){
noswaps <- 0
for(c in 1:dim(mat)[2]){
if(sum(abs(mat[,c])) != 0){
pivotcolumn <- mat[,c]
for(r in 1:dim(mat)[1]){
if(mat[r,c] != 0){
pivot <- mat[r,c]
浏览 3提问于2014-03-15得票数 0
我是BLAS/Lapack的新用户，我想知道有没有一个例程可以进行高斯消元甚至高斯-乔丹消元？我在谷歌上搜索并查看了他们的记录，但仍然找不到。
非常感谢你对我的帮助！
浏览 2提问于2011-08-08得票数 5
回答已采纳
我正在学习Steven C. Charpa的书中的数值方法。书中说" Gauss -Siedel比高斯消元法使用更少的内存，因为它不在矩阵中存储"0“值”，然而，书中写的算法处理与高斯消元法相同的矩阵。我不明白为什么Gauss-Siedel使用更少的内存。我在网上搜索了这个问题，人们都说同样的话，但没有人解释是怎么回事。
注:我可以在书中分享算法，如果不会有版权的问题。
浏览 0提问于2014-05-07得票数 1
2回答
矩阵求逆GPU MATLAB
、、、、
我正在寻找一种方法来找到一个矩阵的伪逆，这样它就可以在GPU上完成。SVD/QR很难并行化，并且不被MATLAB的GPU支持，但似乎LU虽然可以并行运行，但也不被MATLAB的GPU支持。我比较了它的性能，它似乎比在单核CPU上运行要慢。
我正在寻找一个我可以使用的伪逆(甚至是方阵的正则逆)。根据Matlab，使用mldivide ()执行适用于GPU的高斯消除。
我尝试使用A\I，但不幸的是它不能在GPU上高效运行。
有没有人能帮我找到一个优化的并行逻辑单元或高斯消除的代码？
我听说过MAGMA包，但它似乎有很多工作要安装和编译，我真的需要这个简单的东西。
C++代码也是受欢迎的。
谢谢，吉
浏览 2提问于2013-06-22得票数 1
回答已采纳
我已经写了一个高斯消除的代码。但它只适用于N>450和N<500，其中N是行数。
对于少于450，它需要更多的时间。
对于500或更大的值，则不显示任何内容。
#include <conio.h>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
#include <omp.h>
int main() {
printf("\nGauss -\n");
doub
浏览 0提问于2020-10-03得票数 0
1回答
见鬼，我在计算只有1或0的二元矩阵的秩时遇到了一些问题。二进制矩阵的秩将基于使用布尔运算XOR进行的行约简。让我们看看XOR操作：
1 xor 1 =0
1 xor 0= 1
0 xor 0= 0
0 xor 1= 1
给定一个二进制矩阵为
A =
1 1 0 0 0 0
1 0 0 0 0 1
0 1 0 0 0 1
我们可以看到第三行等于第一行与第二行的xor。因此，矩阵A的秩只有2，而不是用matlab函数的秩3。我有一种方法可以使用下
浏览 1提问于2014-12-05得票数 1
我将在Fortran中做一个标准的高斯消除过程，我想知道是否已经有了一些现成的解决方案？最好是以函数的形式。我看了很多fortran的数字书，但似乎没有什么好的解决方案。
浏览 2提问于2013-10-22得票数 2
我刚刚打印了大小为n*n的矩阵
代码是(假设已经定义了所有矩阵)：
public class ReadContents {
public static void main(String args[]) {
System.out.println();
System.out.println("The Matrix Is:");
double mat[][]=new double[col][col];
for(k=0;k<col;k++) {
for(p=0;p<col;p
浏览 2提问于2011-05-06得票数 1
我有一个看起来像这样的系统：
z1 = 5*x1 + x2*cos(x3)
z2 = x1*sin(x3) + 3*x2
z3 = 3*x1 - 2*x2
这是一个微分方程系统的变换(只是给出一些背景，我不认为你真的需要知道这一点)。现在我想找出这个变换的逆，它看起来像这样：
x1 = ...(stuff with z1, z2, z3)...
x2 = ...
x3 = ...
现在我想用这个系统手动找到它相对容易，但我必须用7变量7方程系统做同样的事情，这会变得相当单调乏味。
我找不到让Matlab做这件事的方法，有人能帮上忙吗？
浏览 0提问于2016-04-02得票数 3
3回答
我不知道该选择哪个库(适用于windows)：LAPACK++、Armadillo、IT++、Eigen，或者其他库？
我所需要做的就是检查一个大的(大约10,000*10,000)矩阵是否可对角化，如果是的话，得到对角线和可逆矩阵，使得D=(P^(-1))*A*P。这必须尽快完成。我不知道该使用哪个库。
此外，我很高兴了解这些库的优缺点。
浏览 5提问于2013-10-26得票数 2
2回答
如何使用Cholesky分解计算矩阵的逆的行列式。我发现直接计算矩阵的行列式不是一个好主意。那么，有人能提供一些见解吗？
浏览 0提问于2012-05-04得票数 2
回答已采纳
1回答
我正试图解决算法问题。为了方便起见，我将其复制如下：
给出了一个尺寸为N (1 <= N，M <= 1000)的网格G，它只包含1s和0。如果我们选择一个单元格，这将切换相邻单元格(以及该单元格本身)中的值。两个差异于1行或1列的单元格被认为是相邻的(即对角线单元格不相邻)。我们的目标是在我们需要选择的单元格位置找到一个包含1s的网格G‘，以便将G中的所有单元格转换为0(我们不必选择的单元格被标记为0)。给定任何G(在这个问题中)，G‘被保证存在。
注意:网格中没有环绕。
例如，如果G如下所示：
000
100
000
如果我们选择中间单元，G将变成：
010
011
01
浏览 3提问于2017-11-27得票数 1
2回答
我必须做一个DC电路生成器应用程序，比如或。我做了所有的构建逻辑，但我不知道如何分析生成的电路。电路在内存中表示为一个图，其中每个节点都是一个电路元件(灯泡、电阻器..)如果有必要，我想知道通过两边的电流来点亮灯泡。Java中有一个库可以做到这一点吗？或者任何算法？或者有人能给我解释一下如何解决这个问题？
我知道我必须找到电路回路(图回路)并应用欧姆定律，但一些回路会共享一些电路元件，我不知道下一步该做什么。
浏览 2提问于2016-07-25得票数 0
3回答
我想使用最小二乘技术实现多项式曲线拟合，但要使用各种误差函数，即不仅仅是最小二乘。在MATLAB中有什么方法可以做到这一点吗？(我想比较不同误差函数的结果。我还想使用正则化，我需要改变误差函数)。
你能分享一些资源(MATLAB/C++)，这些资源可以帮助你在没有内置函数的情况下实现曲线拟合吗？我只能找到那些使用高斯消除的方法-这与最小二乘拟合相同吗？
浏览 2提问于2017-09-08得票数 1
4回答
我已经为基本操作构建了一个普通的分数类。唯一的问题是，由于有大量的操作(我正在做高斯消元法)，分子或分母都有溢出。
我有一个100的方程式，所以一个100 x 100的矩阵。并且最终结果需要精确到10^-6。我该怎么办？
浏览 0提问于2011-06-28得票数 1
1回答
我使用Gallager方法为n = 20、j = 3和k = 4生成了一个LDPC n = 20。
奇偶校验矩阵的定义如下：
问题是，我不知道如何进一步生成需要codeword的generator matrix。高斯消去法毫无帮助。你能提出一个更好的方法来寻找generator matrix。提前谢谢。
浏览 2提问于2016-06-17得票数 3
1回答
最近，我被介绍到优化线性问题的单纯形算法中。我想我知道它是如何工作的，但我不知道，为什么会这么复杂。
为什么我们必须使用枢轴行/列/元素。如果我正确地看到了这一点，那就是高斯消去，它以特定的顺序消除元素。为什么命令很重要？是数字噪音造成的吗？如果一切都是精确计算的话，那么排序就不重要了，对吧？
在我到目前为止遇到的所有例子中(不是太多，不太公平)，所有的松弛变量最终都被设置为零。为什么我们一开始就引入它们，而不把所有的不等式直接改写成方程呢？
浏览 0提问于2014-03-05得票数 1
回答已采纳
3回答
如何在MATLAB中实现函数lu(A)，使L*U直接为A，同时也得到真实的L矩阵？
当我使用[L,U] = lu(A)时，MATLAB没有给出正确的L矩阵。当我使用[L,U,P] = lu( A )时，我需要实现P*A = L*U，但我只想将L*U相乘以接收A。
浏览 1提问于2016-12-15得票数 7
回答已采纳
2回答
如何反转只有平移和旋转，没有缩放的4x3矩阵？你会用什么来做OpenGL矩阵求逆(只是没有缩放)？
浏览 42提问于2021-09-29得票数 0
1回答
我有一组方程式
1=x⊕y⊕z
1=x⊕y⊕w
0=x⊕w⊕z
1=w⊕y⊕z
我正在尝试实现高斯消元法来解决这个系统，如所述，用异或代替除法、减法和乘法，但它给出了错误的answer..the正确答案是(x，y，z，w)=(0,1,0,0)
我做错了什么？
public static void ComputeCoefficents(byte[,] X, byte[] Y)
{
int I, J, K, K1, N;
N = Y.Length;
for (K = 0; K < N; K++)
{
浏览 2提问于2012-07-14得票数 2
1回答
非矩形和N维数组的真实示例是什么？我知道你可以使用三维游戏，但不确定什么时候你会使用更多，我也从来没有看到过非矩形阵列的例子。
浏览 0提问于2012-07-05得票数 1
Code = { 000, 011 , 101, 110 }
generator matrix = {011, 101}
这里给出了代码矩阵，我应该找到生成矩阵。但我对确定生成矩阵一无所知。有人能解释一下这件事吗？
浏览 0提问于2012-07-02得票数 0
我需要使用4个点来绘制一个图像。我将这4个点存储在一个数组中，这些位置会发生变化。现在，我需要绘制一个反映这些位置的图像。context.drawImage仅接受1个positon。
应该是这样的，4个位置：
x,y----------x,y
| |
| |
| |
| |
| |
x,y----------x,y
浏览 3提问于2012-05-19得票数 2
回答已采纳
1回答
可以使用Javascript在Node.js中进行图像变形吗？我正在尝试找出是否有可能做这样的事情。有没有什么node.js库可以用来实现图像变形？
浏览 8提问于2012-10-05得票数 0
回答已采纳
我想在Java语言中使用Math.ulp(double)方法计算一系列加法、乘法和除法的浮点舍入误差。根据最后一位单元的维基页面，似乎一个浮点计算的误差，例如2+3或2*3将分别是0.5*ulp( 2+3 )或0.5*ulp( 2*3 )，其中2*3和2+3是浮点计算。但是，将这些错误加在一起并不能说明我在最终产品中得到的实际错误。例如，说2+3*4 = 0.5*ulp(2+3*4) + 0.5*ulp(3*4)的最大误差似乎不能解释我得到的实际误差。因此，我很困惑，也许我误解了Math.ulp(double)，或者我需要使用某种相对错误。我不知道。有没有人能给我解释一下，给我举几个浮点和精确
浏览 3提问于2015-12-26得票数 0
3回答
我有以下问题：
给定空间中的n个点，我正在搜索一个通过这些点的超平面。
这类问题最简单的例子是两个点(x_1=0，x_2=0)和(1，-1)，我希望返回1*x_1+1*x_2=0。
我的点将是32位整数的n元组。所需超平面的系数a_i a_1 x_1 + a_2 x_2 + ... =c也必须是32位整数。如果超平面不能以这种方式定义，我想要报告这一点。
我的项目是用c++编写的。
我可能会自己编写代码，但我预计这将是相当多的工作。此外，我的直觉是，这是一个足够普遍的问题，可能会有一个开源库来解决我的问题。有谁知道可以解决我的问题的库吗？
提前感谢！
浏览 1提问于2012-07-22得票数 5
回答已采纳
1回答
我正在尝试获得一个点的笛卡尔坐标，为此我使用了一个简单的函数来转移窗口坐标，类似于(100,200)。我想把它转换成笛卡尔坐标，问题是窗口大小是可变的，所以我真的不知道如何实现它。
编辑：
操作：
我尝试通过执行以下操作来使用guUnPorject
GLdouble modelMatrix[16];
glGetDoublev(GL_MODELVIEW_MATRIX,modelMatrix);
GLdouble projMatrix[16];
glGetDoublev(GL_PROJECTION_MATRIX,projMatrix);
double position[3];
gluUnPr
浏览 2提问于2011-12-04得票数 2
回答已采纳
有谁知道我如何在Rascal中用愚蠢的方式保留关键字吗？
由于GitHub ()上报告的#968号问题，我不能简单地使用单引号。这个问题很可能会得到解决。目前，我正在使用一种语言，在这种语言中，关键字不区分大小写，并试图分析代码，在这些代码中，显然有人在关键字中随机应用camelcase或所有大写字母，然而，我正在寻找解决办法。现在，我只是用三种方式定义每个关键字。由于camelcase的应用中存在错误(比如GroupbyGroup，它应该是camelcase中的GroupByGroup )，所以我不得不用更多不同的方式重新定义关键字。因此，我希望有人知道一个更好的解决办法(或者可以明显地解决
浏览 11提问于2016-06-24得票数 0
回答已采纳
1回答
我在密码学应用程序中遇到了以下问题:我已经给出了一组线性同余
a[1]*x[1]+a[2]*x[2]+a[3]*x[3] == d[1] (mod p)
b[1]*x[1]+b[2]*x[2]+b[3]*x[3] == d[2] (mod p)
c[1]*x[1]+c[2]*x[2]+c[3]*x[3] == d[3] (mod p)
这里，x是未知数，a，b，c，d是给定的
系统很可能是欠确定的，所以我有一个很大的解决方案空间。我需要一个算法，使用伪随机数生成器(或失败)找到该问题的等分布解决方案(即在解空间中等分布)。
据我所知，在我的线性代数课程中，大多数线性方程组的标准算法都不能直接应
浏览 3提问于2011-12-14得票数 2
回答已采纳
3回答
我有一个相当基础的数学问题，但诀窍是我需要在C++中使用它。我正在遵循维基百科上给出的伪代码。这是我的尝试：
createMatrixForAllSolutions(*this);
std::cout << equationMatrix.to_string() << endl;
bool solved = false;
int rows = equationMatrix.getRows();
int cols = equationMatrix.getCols();
int i = 0;
int j = 0;
int maxi = 0;
double current =
浏览 0提问于2011-10-09得票数 1
我不可能是第一个想到这个问题的人，而且肯定有一个没有人设计密码系统来解决这个问题的理由。
让我们把MECC定义为椭圆曲线点的矩阵，MI定义为非负整数的矩阵。
给出了MECC G、MI v和MECC P=v\cdot G。量子计算机如何从P和G中计算出D4？
浏览 0提问于2020-09-28得票数 2
回答已采纳
6回答
在C++中，我想做一些类似的事情：
int n = get_int_from_user();
char* matrix = new char[n][n];
matrix[0][0] = 'c';
//...
matrix[n][n] = 'a';
delete [][] matrix;
但这当然行不通。做类似的事情最好的方法是什么？我已经看到了一些解决方案，但它们看起来相当混乱。
浏览 12提问于2008-11-02得票数 20
回答已采纳
2回答
我正在尝试使用OpenMP并行化一个简单的c++代码，该代码通过高斯消元法求解线性方程组。
我的代码的相关部分是：
#include <iostream>
#include <time.h>
using namespace std;
#define nl "\n"
void LinearSolve(double **& M, double *& V, const int N, bool parallel, int threads){
//...
for (int i=0;i<N;i++){
#p
浏览 0提问于2015-02-26得票数 2
2回答
我有一个数据点的网格，目前我使用双线性插值来查找网格中的缺失点。我被指向了Kriging的方向，也就是最好的线性无偏估计器，但我找不到好的源代码或代数解释。有没有人知道我可以使用的其他插值方法？
--更新@Sam Greenhalgh我考虑过双三次插值，但我使用我找到的代码示例得到的结果似乎不太对劲。
下面是Bicubic的代码示例
注意，我是用C#编写代码的，但我也欢迎来自其他语言的示例。
//array 4
double cubicInterpolate(double[] p, double x)
{
return p[1] + 0.5 * x *
浏览 10提问于2014-11-21得票数 3
除了Lisp方言之外，还有其他基于S表达式或波兰语前缀表示法的语言吗？
它不必是函数式编程语言，也不必以某种方式与Lisp相关。
当我说Lisp方言时，我指的是Arc、AutoLISP、Clojure、Common Lisp、Emacs Lisp、EuLisp、Franz Lisp、Interlisp、ISLISP、LeLisp、Maclisp、MDL、Newlisp、NIL、Picolisp、可移植标准Lisp、球拍、方案、SKILL、Spice Lisp、T、XLISP、Zetalisp (来自页面)。
有许多完全不同的语言使用类似C的语法或类似Python的语法，但是对于S表达式或波兰语
浏览 28提问于2013-06-17得票数 0
2回答
System.Drawing.Drawing2D.Matrix有一个非常方便的构造函数
public Matrix(RectangleF rect, PointF[] plgpts)
基本上，你可以给出一个矩形和一个转换后的矩形的平行四边形，然后你就会得到相应的转换。在我的例子中，这是非常有用的。
不幸的是，我需要双精度的矩阵。有人知道如何从矩形和平行四边形计算System.Windows.Media.Matrix吗？
浏览 0提问于2012-09-12得票数 0
我正在寻找一种有效的方法来实现一个简单的滤波器，它具有一个随时间变化的系数，并由与输入信号长度相同的向量指定。
以下是所需行为的简单实现：
def myfilter(signal, weights):
output = np.empty_like(weights)
val = signal[0]
for i in range(len(signal)):
val += weights[i]*(signal[i] - val)
output[i] = val
return output
weights = np.random.uni
浏览 15提问于2017-01-27得票数 2
回答已采纳
我需要解一组形式为Ax = B的联立方程。我使用了numpy.linalg.solve函数，输入A和B，但是得到了误差LinAlgError:数组的最后2维必须是正方形。我该怎么解决这个问题？
这是我的密码：
A = matrix([[v1x, v2x], [v1y, v2y], [v1z, v2z]])
print A
B = [(p2x-p1x-nmag[0]), (p2y-p1y-nmag[1]), (p2z-p1z-nmag[2])]
print B
x = numpy.linalg.solve(A, B)
矩阵/向量的值是在代码的前面计算的，这很好，但是值是：
A =
(-0.5
浏览 1提问于2018-02-21得票数 13
1回答
我使用gf256库来做伽罗瓦域的数学运算，并且我把它放在一个数值矩阵中。但是当使用它调用np.linalg.inv()时，它会抛出一个错误。
这就是摘要，下面是详细信息：
import numpy as np
from gf256 import GF256 as gf
npgf = np.vectorize(gf)
arr = np.identity(4, np.uint8) * 10
gfarr = npgf(arr)
经过这一切之后，gfarr看起来像这样
array([[GF256(0b00001010), GF256(0b00000000), GF256(0b00000000),
浏览 15提问于2018-08-21得票数 0
2回答
在“黑书”中，数值计算第三版，给出了求解线性方程组的Gauss算法。后面是一个关于计算LU分解的部分，然后使用这个部分来求解线性方程组(参见LUdcmp::solve on第53页)。不幸的是，这本书并没有解释为什么一种方法比另一种方法更好。这两种方法是等同的，还是有理由在特定情况下选择一种方法而另一种方法？
浏览 2提问于2016-01-22得票数 6
回答已采纳
3回答
/* Program to demonstrate gaussian <strong class="highlight">elimination</strong>
on a set of linear simultaneous equations
*/
#include <iostream>
#include <cmath>
#include <vector>
using namespace std;
const double eps = 1.e-15;
/*Preliminary pivoting s
浏览 0提问于2011-03-05得票数 1
回答已采纳
1回答
实体识别和消歧api
、、、、
我正在寻找实体识别和消除歧义的API或库。我发现了拼图算法，但在单个词的层面上做了WSD，所以它不能识别“金融机构”或“利率”作为一个单位，即使这些单位可以在wordnet中找到并具有同义词集(例如。参见“wordnet中的金融机构”。
然而，我找到了BabelNet，它有一个图形化的web界面就可以做到这一点。但是，它在其当前API v.3中没有提供该功能
您是否知道提供此类功能的API或库？我的意思是将“金融机构”或“利率”作为一个单位，并在此基础上进行意义上的消歧！？
提前谢谢。
浏览 17提问于2015-06-03得票数 1
我偶然发现了一个代码，它使用低性能系列for循环来计算所谓的"Crout因式分解“。如果您需要了解更多，我认为这个概念类似于什么是。乍一看，代码是一个for循环，它可以通过omp指令变得并行：
SparseMtrx *Skyline :: factorized()
{
// Returns the receiver in U(transp).D.U Crout factorization form.
...
// this loop is very expensive and only uses single thread:
浏览 5提问于2022-04-07得票数 0
回答已采纳
2回答
简略版:如何才能将连续高斯圈成真正的离散高斯(通常表示为\mathcal{D}_{\mathbb{Z},\alpha q})？其目标是从连续LWE降为真正离散LWE，并将其与从\textsf{GapSVP}_\gamma到连续LWE的约简相结合。
更长的版本:在Reg05中，他们认为噪声的离散高斯(表示\bar{\Psi}_\alpha，有时是\lfloor\mathcal{D}_{\alpha q}\rceil)是一个“奇怪的高斯”：它可以从参数\alpha q的连续高斯中得到，你只需把它旋转到最接近的整数模q。它们还证明了，如果你能求解连续的线性微分方程，你就可以求解\textsf{Gap
浏览 0提问于2021-03-06得票数 7
回答已采纳
3回答
我被困在这里了。我有一个大小为NxN的矩阵，存储在一个双精度数组中。然后我想删除给定的一列，比如说第一列。所以我创建了一个新的双精度数组，大小为NxN-1，并将第一个矩阵中的值复制到第二个矩阵中，当然第一列除外。然后我想将第一个数组设置为第二个数组。我在这里是空白的。
double matrix[N][N]
//fill up the matrix code here...
// remove first column of array
double newMatrix[N][N-1];
for(i = 0; i < N; i++){
for(j = 1; j < N;
浏览 1提问于2013-03-08得票数 2
回答已采纳
1回答
用C++求解一组模方程组
、、、
我正在研究c++中的二次筛分算法。经过高斯消除后，我需要求解一组模块方程，例如：
(1) b + c = 0 mod 2
(2) a + c = 0 mod 2
在这里，符号=的意思是“是一致的”。我正在处理矩阵，如下所示：。如果有人有任何想法，如何实现这样一个函数，以解决这些方程，我将不胜感激。
浏览 1提问于2015-11-27得票数 2
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
腾讯会议
实时音视频活动推荐
腾讯云开发者
