小编And*_*uri的帖子

在谷歌关于矩阵求逆算法时,我发现有几种方法(和意见!)关于如何在代码中执行此操作.我想知道哪种方法最快,或者性能最好的方法,并试图找到答案我什么都没发现.

我知道在某些情况下可以计算伪逆(使用SVD,cholevsky,...),我实际上在我的代码中使用了一些,并且我知道有几次逆不存在等等.很容易找到特定问题的具体答案,但不是对这个巨大的(巨大的!)问题的一般直觉,即矩阵求逆.

所以我的问题是:

对于小型矩阵,哪种方法的性能最佳？精确度？大矩阵怎么样？

我个人的情况是6x6(EDIT:symetric)矩阵必须被反转数千次(是的,是的,具有不同的值)并且我需要高精度,但是肯定速度会非常方便.

请注意,我不是在寻找代码,我会根据自己的情况自行编写任何答案,但我认为这是许多程序员想要了解的问题.

由于这段代码,我使用openCv进行了delaunay三角测量: 示例代码 (在partiluclar draw_subdiv中).然而,当我想显示三角测量时,我得到不属于三角测量的网格和线.这条线是由于三角测量算法开始其工作考虑在"无限"处发布的三角形.

你能解释一下我如何只将网格划入凸包(没有这条线)吗？

显示功能:

 void draw_subdiv(Mat &img, Subdiv2D& subdiv, Scalar delaunay_color)
{

  vector<Vec6f> triangleList;
  subdiv.getTriangleList(triangleList);
  vector<Point> pt(3);

  for(size_t i = 0; i < triangleList.size(); ++i)
    {
      Vec6f t = triangleList[i];

      pt[0] = Point(cvRound(t[0]), cvRound(t[1]));
      pt[1] = Point(cvRound(t[2]), cvRound(t[3]));
      pt[2] = Point(cvRound(t[4]), cvRound(t[5]));

      line(img, pt[0], pt[1], delaunay_color, 1);
      line(img, pt[1], pt[2], delaunay_color, 1);
      line(img, pt[2], pt[0], delaunay_color, 1);
    }
}

主功能 :

Mat image = imread(argv[1], 1);

 ..... ....
 //creat delaunay                                                                                                                                 
 Scalar delaunay_color(255, 255, 255), point_color(0,0,255);
 Rect rect(0,0,image.cols, image.rows); …

正如标题所说,如何在Matlab中向矩阵添加尾随单例维度？

A=ones(3,3,1);给出3x3矩阵,同时A=ones(1,3,3);给出1x3x3矩阵

添加到具体问题:

我有一个应用程序,我有N MxM矩阵,我需要堆叠它们,结果成为一个MxMxN矩阵.但是,N可以是1,如果是,我需要矩阵为MxMx1.

注意:我知道这对Matlab脚本没有多大意义,正如Loren所说,在Matlab中,在非单例之后存在无限的"单例"维度.但是,这是在mex环境中mxGetNumberOfDimensions使用的地方.

问题的简要说明:我使用Newton Raphson算法在多项式中进行根发现,但在某些情况下不起作用.为什么？

我从"c ++中的数字配方"中得到了一个N​​ewton Raphson混合算法,如果New-Raph没有正确收敛(具有低导数值或者如果收敛速度不快),则将其一分为二.

我用几个多项式检查了算法并且它有效.现在我在我的软件内部进行测试,并且我总是遇到一个特定多项式的错误.我的问题是,我不知道为什么这个多项式只是没有得到结果,而其他人做的很多.因为我想改进任何多项式的算法,需要知道哪一个是没有收敛的原因,所以我可以正确对待它.

下面我将发布我可以提供的有关算法和我遇到错误的多项式的所有信息.

多项式:

f(t)= t^4 + 0,557257315256597*t^3 - 3,68254086033178*t^2 +
+ 0,139389107255627*t + 1,75823776590795

它的第一个衍生物:

 f'(t)= 4*t^3 + 1.671771945769790*t^2 - 7.365081720663563*t + 0.139389107255627

情节:

根(由Matlab提供):

  -2.133112008595826          1.371976341295347          0.883715461977390 
  -0.679837109933505

算法:

double rtsafe(double* coeffs, int degree, double x1, double x2,double xacc,double xacc2)
    {
    int j;
    double df,dx,dxold,f,fh,fl;
    double temp,xh,xl,rts;
    double* dcoeffs=dvector(0,degree);
    for(int i=0;i<=degree;i++)
        dcoeffs[i]=0.0;
    PolyDeriv(coeffs,dcoeffs,degree);
    evalPoly(x1,coeffs,degree,&fl);
    evalPoly(x2,coeffs,degree,&fh);
    evalPoly(x2,dcoeffs,degree-1,&df);
if ((fl > 0.0 && fh > 0.0) || (fl < 0.0 && fh < 0.0))
    nrerror("Root must …

我试图为三角测量过程获得一个3x4相机矩阵,但calibrateCamera()只返回3x3和4x1矩阵.

如何从这些矩阵中获取3x4？

提前致谢!!

我通过imfilter命令将图像与LOG过滤器相关联.

z=imfilter(I,fspecial('log',11,1.5)); 
figure
imshow(z);

图片:

假设我想要一个需要做很多事情的CUDA内核,但是有些圆顶参数对所有内核都是不变的.这个参数作为输入传递给主程序,因此它们不能在a中定义#DEFINE.

内核将运行多次(大约65K)并且需要这些参数(以及一些其他输入)来进行数学运算.

我的问题是:将这些常量传递给内核的最快(或者最优雅)方式是什么？

常量是2或3个元素长度float*或int*数组.它们将在这些中约5~10个.

玩具示例:2个常数const1和const2

__global__ void kernelToyExample(int inputdata, ?????){
        value=inputdata*const1[0]+const2[1]/const1[2];
}

好点吗

__global__ void kernelToyExample(int inputdata, float* const1, float* const2){
        value=inputdata*const1[0]+const2[1]/const1[2];
}

要么

__global__ void kernelToyExample(int inputdata, float const1x, float const1y, float const1z, float const2x, float const2y){
        value=inputdata*const1x+const2y/const1z;
}

或者可能在一些全局只读内存中声明它们并让内核从那里读取？如果是这样,L1,L2,全球？哪一个？

有没有更好的方法我不知道？

在特斯拉K40上运行.

我有一些由MATLAB给C++提供的大数组(因此我需要按原样使用它们)需要转换和置换(行 - 市长,列市长问题).

该数组imgaux是double类型的大小size_proj[0]*size_proj[1]*size_proj[2],需要进行浮动,更改值的某些位置.最小的例子如下:

#include <time.h>
#include <stdlib.h>  
   int main(void){
      int size_proj[3];
    size_proj[0] = 512;
    size_proj[1] = 512;
    size_proj[2] = 360;
    size_t num_byte_double = size_proj[0] * size_proj[1] * size_proj[2] * sizeof(double);
    size_t num_byte_float = size_proj[0] * size_proj[1] * size_proj[2] * sizeof(float);

    double *imgaux = (double*) malloc(num_byte_double);
    float* img = (float*) malloc(num_byte_float);


    clock_t begin, end;
    double time_spent;

    begin = clock();
    for (int k = 0; k < size_proj[0]; k++)
        for (int i = 0; i <size_proj[1]; i++) …

我有这个矩阵:

a = [1 2 2 1; 1 1 2 2]

%   1     2     2     1 
%   1     1     2     2

我想找到所有的1并将它们归零.

[~, a_i] = find(a == 1);
a(a_i) = 0

%   0     2     2     1 
%   0     0     2     2

为什么第一排还有1？

我在开发代码时使用MATLABs git支持,经常提交,推送和所有标准的源代码控制.

但是,我只使用MATLABs用户界面,基本上可以通过右键单击文件夹并浏览菜单直到找到正确的选项(见下图).

有没有办法制作MATLABs命令窗口运行git命令,而不是每次都需要通过菜单导航？