generated C/C++ code from `pinv`-command gives different result than the command `pinv` itself

11 次查看（过去 30 天）
Victor Lam
Victor Lam 2018-11-14
评论： Mike Hosea 2018-11-20
I am using Matlab Coder for my Arduino. Inside my Matlab function I use the Matlab-command `pinv`. This Matlab-function is then converted into C/C++ code (let's call it `pinvC`) using the Code Generator from Matlab.
However, I noticed the following. The result of Matlab-command `pinv` inside Matlab function is different to the result of `pinvC`. What is causing this? How can I solve it?

请先登录，再回答此问题。

回答（1 个）

James Tursa
James Tursa 2018-11-14
How much different? The MATLAB code is calling into BLAS and LAPACK library routines to accomplish this. But the Coder is not using that same library code to generate its results, hence the difference.
  2 个评论
Victor Lam
Victor Lam 2018-11-14
The pseudo-inverses are a little different. But after some multiplications the small errors accumulate and the difference is about a factor of 5.
So I tried to compute the pinv beforehand and store it in .mat file and load it inside the function using `coder.load` but then some variables are blowing up to order of 1e3, while they should be below 10.
Could the C/C++ `pinv` version be different, for numerically stability reasons?
Mike Hosea
Mike Hosea 2018-11-20
James is correct. Different stability properties induced by PINV per se does not seem likely, but different PINV implementations will issue different numerical perturbations, which might be magnified by how you use PINV.
You might try an experiment perturbing the output of PINV in MATLAB by some units in the last place to explore what your algorithm does with it in MATLAB, i.e. with no code generation involved. Something like:
X = pinv(A);
X = X + randi([-10,10],size(X)).*eps(X)

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 MATLAB Coder 的更多信息

标签

产品


版本

R2018b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by