Svenska ElektronikForumet

Tips om att återvända, i alla fall från mitt håll, är inte för att undvika någon känd bugg utan mer för att inte få icke-deterministiskt beteende efter att den snurrat i typ tre timmar och malloc börjar få problem pga fragmentering eller nått sånt.

Det här problemet kanske kan bero på någon vektor/matris som borde normaliseras mellan varven? Min erfarenhet av att mangla matriser om och om igen är att om man sparar om någon matris med lite förändringar varje gång så kan små fel växa (ju sämre precision desto fortare går det)...

MATLAB kan väl sätta precision lätt? Hur bra sköter sig MATLAB-koden med typ digits(2)?

Du menar att det kan vara malloc som bråkar med mig? Trots att jag sitter på en PC?

Hur tror du det skulle vara om jag skapar matriserna bara en gång och sedan återanvänder jag dem? Alltså jag allokerar matriserna som jag har behov utav, bara en gång.

Nej, jag tror inte att du nödvändigtvis har problem pga det utan påpekade det bara för att det kan leda till huvudvärk om man har otur (fel som kan var svåra att felsöka för att de dyker upp sporadiskt, inte allomfattande ofrånkomliga fel på något vis).

float har ju 7, kanske 7.5 signifikanta siffror, double har 15. Så precisionen borde blivit bättre med double.
Och du är såklart medveten om att octave avrundar talen när den skriver ut dem - annars testa "format long" eller nåt sånt.
Minnesallokering har inget med precisionen att göra. Att återanvända minnet har flera lobbat för tidigare i tråden.
Så det går bra. Men det blir lite andra trix då som att säkerställa att inte två matriser använder samma minne samtidigt.

Okej! Det blev lite bättre!

Här har jag räknat
\(H = USV^T\)
och sedan
\(H - USV^T\)

Floats

Markering_043.png

Double

Markering_044.png

Notera att Octave får e-18 i precision och Octave har andra algoritmer för att lösa SVD. Jag löser SVD igenom QR-faktorisering. Det är en seg metod, men den är enkel och pedagogisk.

Jag ska testa med Armadillo C++ och jämföra resultaten. STM32 kan ju använda sig av C++, men jag förstår inte varför STM32 har C++ som standard om C++ finns tillgänglig. Orsaken varför C används mest inom inbyggda system har väll med en kostnadsfråga, eller hur?

Enligt tillverkarna av Armadillo så använder Armadillo OpenBLAS som matematikberäknare, och OpenBLAS finns för ARM.

Tror ni det skulle vara några problem att kompilera C++ kod för STM32?

Vilken kostnad ser du?

Ingen. Jag har bara hört att C används mer inom inbyggda system då kompilatorer finns mer tillgängligt för C än C++. Det är ju ett jobb det också att skapa en kompilator.

Jag tycker C++ har fått ett rykte om sig att vara olämpligt för inbäddade system men jag skulle nog säga att de som håller med ryktet har missat designfilosofin i C++: "you pay for what you use".

Mao man måste ha koll på vilka språkfinesser som kan skapa onödiga objekt eller oväntade kopieringar av objekt mm, samt förstå att vissa funktioner i standardbiblioteket kan dra med en massa annat som gör att koden blir onödigt stor osv.

Själv har jag testat att skriva samma funktionalitet i C och C++ (med klasser, virtuella funktioner mm) på Arduino Atmega328p och i de fallen har storleksskillnaden, både i flash och ram, varit obetydlig. Att göra tvärtom, dvs skriva objektorienterat i C går också men är mycket mer plågsamt och pratigt. Många structar och funktionspekare blir det.

Jag har också hört detta rykte, men jag har aldrig hört någon motivering. Det kanske var så förut, att C++ var olämpligt för inbyggda system?

Jag tror det är just risken att råka skapa lite objekt dynamiskt lite här och där som gjort att C++ inte är så poppis bland uC folket, ungefär av samma anledning som att malloc inte är så poppis.
Och en hel del, så här har vi alltid gjort

Om jag tolkar IAR hemsida rätt klarar deras kompilatorer både C och C++ , och de supportar en hel hög plattformar.. (Allt man kan tänka sig utom PIC typ)
https://www.iar.com/iar-embedded-workbench/ klicka på "choose an architecture"

Vad tror du om http://arma.sourceforge.net/docs.html?

AndLi skrev:Jag tror det är just risken att råka skapa lite objekt dynamiskt lite här och där som gjort att C++ inte är så poppis bland uC folket, ungefär av samma anledning som att malloc inte är så poppis.
Och en hel del, så här har vi alltid gjort

Visst är det så. Det är en liten del att folk skriver c++ som att det är ett desktop-system och kommer till slutsatsen att c++ inte fungerar för embedded, och en jättestor del "ja men (random namn) sa ju att det är dåligt.."

Jag gick för några år sen en kurs om embedded-programmering. Jag frågade om vi fick använda c++ istället för c. Läraren skrattade bara och sa "Ja om du tror att det går så gör det. Du ändrar dig nog ganska snabbt."
När jag lämnat in några uppgifter sa han "ja det var ju kul att du försökte men c++ används inte i industrin för den kompilerade koden blir minst 2x större"
Jag skrev en rapport där jag visade att olika c++ vs c-saker kompilerade till exakt samma instruktioner. Hans enda svar var "nej experterna håller inte med". När jag försökte få honom att klargöra vilka dessa "experterna" var blev det knäpptyst.

@Al
Jag ser inga problem som gör att armadillo inte skulle fungera, förutsatt att det finns de funktioner du behöver. Just det har jag ingen aning om.

Även C++ klarade inte ens att få det samma som GNU Octave.

Armadillo med C++

Kod: Markera allt

#include <iostream>
#include <armadillo>
#include "hankel.h"

using namespace std;
using namespace arma;

int main() {

	/*
	 * Declare vector and matrix
	 */
	mat u(1, 72);
	mat y(1, 72);
	mat toe(72, 72);
	mat triup(72, 72);
	mat inverse(72, 72);
	mat g(1, 72); // markov parameters
	vec g_v; // Impulse vector
	mat H1(72, 72); // Hankel 1 of markov parameters
	mat H2(72, 72); // Hankel 2 of markov parameters
	mat H1_half(36,36);

	/*
	 * G(s) = 4/(2s^2 + s + 5) - Model
	 */

	double input[72] = { 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
			5, 5, 5, 5, 5, 5, 5, 5, 5, 5,

			5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
			5, 5, 5, 5, 5,

			5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5 };

	double output[72] = { 0.00000, 3.47145, 6.41826, 4.35654, 2.53855, 3.76202,
			4.88287, 4.15747, 3.46685, 3.89656, 4.32183,

			4.06753, 3.80580, 3.95616, 4.11713, 4.02833, 3.92938, 3.98176,
			4.04256, 4.01169, 3.97436, 3.99253,

			4.01544, 4.00476, 3.99070, 3.99697, 4.00559, 4.00192, 3.99664,
			3.99879, 4.00202, 4.00077, 3.99878,

			3.99952, 4.00073, 4.00030, 3.99956, 3.99981, 4.00026, 4.00012,
			3.99984, 3.99993, 4.00009, 4.00005,

			3.99994, 3.99997, 4.00003, 4.00002, 3.99998, 3.99999, 4.00001,
			4.00001, 3.99999, 4.00000, 4.00000,

			4.00000, 4.00000, 4.00000, 4.00000, 4.00000, 4.00000, 4.00000,
			4.00000, 4.00000, 4.00000, 4.00000,

			4.00000, 4.00000, 4.00000, 4.00000, 4.00000, 4.00000 };

	// Insert
	for (int i = 0; i < 72; i++) {
		u(0, i) = input[i];
		y(0, i) = output[i];
	}

	// Toeplitz
	toe = toeplitz(u);

	// Triangular upper
	triup = trimatu(toe);

	// inverse
	inverse = inv(triup);

	// impulse
	g = y * inverse;

	// Turn it to a vector
	g_v = trans(g);

	// hankel
	H1 = hankel(g_v, 1);
	H2 = hankel(g_v, 2);

	// cut
	H1_half = H1( span(0,35), span(0,35) );


	mat U;
	vec s;
	mat V;

	svd(U,s,V,H1);


	s.print("s: ");


	return 0;
}

Resultatet blev:

s:
1.4958e+00
1.0602e+00
2.1347e-05
1.3254e-05
1.2733e-05
1.2129e-05
1.0989e-05
7.4545e-06
6.7454e-06
6.3479e-06
5.8772e-06
5.8022e-06
5.5698e-06
5.4533e-06
5.0749e-06
4.2402e-06
4.1427e-06
4.0063e-06
3.9360e-06
3.7566e-06
3.7451e-06
3.6993e-06
3.5652e-06
3.4892e-06
3.4648e-06
3.4544e-06
2.6481e-06
2.5976e-06
2.2883e-06
1.8049e-06
1.7767e-06
1.6989e-06
1.4169e-06
1.2371e-06
1.1915e-06
1.1413e-06
9.0625e-07
8.3946e-07
7.3417e-07
7.0425e-07
5.9058e-07
5.7882e-07
4.7199e-07
4.3554e-07
2.4071e-07
9.7150e-08
7.4907e-08
5.6327e-08
5.6258e-08
3.2675e-08
2.9747e-08
3.5307e-10
3.0517e-10
1.1414e-16
1.1414e-16
1.1414e-16
1.1414e-16
1.1414e-16
1.1414e-16
1.1414e-16
1.1414e-16
1.1414e-16
1.1414e-16
1.1414e-16
1.1414e-16
1.1414e-16
1.1414e-16
1.1414e-16
1.1414e-16
1.1414e-16
1.1414e-16
1.1414e-16

Mitt C-bibliotek:

Kod: Markera allt


#include <stdio.h>
#include <stdlib.h>
#include "LinearAlgebra/declareFunctions.h"

int main() {

	/*
	 * G(s) = 4/(2s^2 + s + 5) - Model
	 */

	double input[72] = {  5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,

			   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,

			   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5,   5};


	double output[72] = {0.00000,   3.47145,   6.41826,   4.35654,   2.53855,   3.76202,   4.88287,  4.15747,   3.46685,   3.89656,   4.32183,

			   4.06753,   3.80580,   3.95616,   4.11713,   4.02833,   3.92938,   3.98176,   4.04256,   4.01169,   3.97436,   3.99253,

			   4.01544,   4.00476,   3.99070,   3.99697,   4.00559,   4.00192,   3.99664,   3.99879,   4.00202,   4.00077,   3.99878,

			   3.99952,   4.00073,   4.00030,   3.99956,   3.99981,   4.00026,   4.00012,   3.99984,   3.99993,   4.00009,   4.00005,

			   3.99994,   3.99997,   4.00003,   4.00002,   3.99998,   3.99999,   4.00001,   4.00001,   3.99999,   4.00000,   4.00000,

			   4.00000,   4.00000,   4.00000,   4.00000,   4.00000,   4.00000,   4.00000,   4.00000,   4.00000,   4.00000,   4.00000,

			   4.00000,   4.00000,   4.00000,   4.00000,   4.00000,   4.00000};


	// Create toeplitz matrix
	matrix* toe = toeplitz(input, 72);

	// Create upper triangular matrix of the toeplitz matrix
	matrix* tru = triu(toe, 0);

	// Find the inverse of tru
	matrix* iv = inv(tru);

	// Create vector the horizon - Important! Else, we cannot find the markov-parameters g
    matrix* Y = create(output, 1, 72);

	// Multiply Y with the iv matrix - Find the markov-parameters g
    matrix* G = mul(Y, iv, false);

    // Detete some mats
    freeMatrix(toe);
    freeMatrix(tru);
    freeMatrix(iv);
    freeMatrix(Y);

    /*--------------------------------------*/

    // turn vector G into a normal vector because the function hankel only want 1D vector
    double g[72];
    for(int i = 0; i < 72; i++)
    	g[i] = *(G->data + i);


    // Delete G
    freeMatrix(G);

    // Create hankel matrix H0 and H1
    matrix* H0 = hankel(g, 72, 1);
    matrix* H1 = hankel(g, 72, 2);

    // Cut H1 and H2 to the half - Remember indexing from zero here!
    matrix* H0_half = cut(H0, 0, 35, 0, 35); // 36x36 from 72x72

    // Remove H0
    freeMatrix(H0);


    // Measure the size

    /*
    int n;
    int m;
    size(H0_half, &n, &m);
    printf("The H0_half have the size %dx%d\n\n,", n, m);
    */


    // Do SVD on
    matrix* u =  initMatrix(36,36); // We know the size from the size() command above
    matrix* s =  initMatrix(36,36);
    matrix* v =  initMatrix(36,36);
    svd(H0_half, u, s, v);
    //freeMatrix(H0_half);

    matrix* p = mul(s, tran(v), false);
    matrix* o = mul(u, p, false);
    //printMatrix(sub(H0_half, o));
    printMatrix(s);


    matrix* E = initMatrix(36,36);
    int n = E->row;
    int m = E->column;
    double* e_ptr = E->data;
    double* s_ptr = s->data;
    for(int i = 0; i < n; i++){
    	for(int j = 0; j < m; j++){
    		if(j == i){
    			*((e_ptr + i*n + j)) = 1/sqrt(*(s_ptr + i*n + j));
    		}
    	}
    }
    //printMatrix(u);



	return EXIT_SUCCESS;
}

Kod: Markera allt

 1.495784917789025004 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000
 0.000000000000000000 1.060244148475468995 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000
 0.000000000000000000 0.000000000000000000 0.000019400761171795 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000
 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000011552529611520 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000
 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000011392477463621 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000
 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000009352345453395 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000
 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000007272746558755 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000
 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000006357671418403 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000
 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000005818190973607 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000
 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000005489075329517 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000
 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000004932159190776 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000
 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000004886374883577 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000
 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000004313931603235 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000
 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000004306302153600 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000
 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000003840821399643 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000
 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000003728868524627 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000
 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000003585514083560 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000
 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000003529845886977 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000
 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000003466027087633 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000
 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000003192936832562 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000
 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000003057194448596 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000
 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000002666912744386 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000
 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000002635144423442 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000
 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000002441116481867 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000
 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000002331233014617 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000
 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000002243085195337 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000
 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000001359374862849 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000
 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000867825790150 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000
 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000796750618645 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000
 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000709803530064 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000
 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000624556496871 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000
 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000499296385836 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000
 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000403271341945 0.000000000000000000 0.000000000000000000 0.000000000000000000
 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000301750776794 0.000000000000000000 0.000000000000000000
 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000095266142240 0.000000000000000000
 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000000000000000 0.000000049205780117

Och GNU Octave

Kod: Markera allt

u = linspace(5,5, 72);
t = linspace(0, 71, 72);
G = tf(4, [2, 1, 5]);
y = lsim(G, u, t); close all;

g = y*inv(triu(toeplitz(u)));

% Create the hankel matrix
k = 1
H = hankel(g)(1:length(g)/2,1+k:length(g)/2+k);

[u,s,v] = svdsim(H);
s

Kod: Markera allt

s =

Diagonal Matrix

 Columns 1 through 8:

   1.49578691063427e+00                      0                      0                      0                      0                      0                      0                      0
                      0   1.06024607224369e+00                      0                      0                      0                      0                      0                      0
                      0                      0   1.33255591690263e-15                      0                      0                      0                      0                      0
                      0                      0                      0   1.31787315854601e-15                      0                      0                      0                      0
                      0                      0                      0                      0   1.06167277714223e-15                      0                      0                      0
                      0                      0                      0                      0                      0   1.01367984757753e-15                      0                      0
                      0                      0                      0                      0                      0                      0   9.72307420314964e-16                      0
                      0                      0                      0                      0                      0                      0                      0   8.40989444285943e-16
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0

 Columns 9 through 16:

                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
   7.84581382552389e-16                      0                      0                      0                      0                      0                      0                      0
                      0   7.19899882445797e-16                      0                      0                      0                      0                      0                      0
                      0                      0   6.84337943063624e-16                      0                      0                      0                      0                      0
                      0                      0                      0   6.52662380533485e-16                      0                      0                      0                      0
                      0                      0                      0                      0   5.90031018652976e-16                      0                      0                      0
                      0                      0                      0                      0                      0   5.11850577956612e-16                      0                      0
                      0                      0                      0                      0                      0                      0   4.99214531415125e-16                      0
                      0                      0                      0                      0                      0                      0                      0   4.15463203077949e-16
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0

 Columns 17 through 24:

                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
   3.94340220158252e-16                      0                      0                      0                      0                      0                      0                      0
                      0   3.39904540767023e-16                      0                      0                      0                      0                      0                      0
                      0                      0   2.91838627478874e-16                      0                      0                      0                      0                      0
                      0                      0                      0   2.64709866341420e-16                      0                      0                      0                      0
                      0                      0                      0                      0   2.60750002495710e-16                      0                      0                      0
                      0                      0                      0                      0                      0   2.42707397956568e-16                      0                      0
                      0                      0                      0                      0                      0                      0   2.33268965183877e-16                      0
                      0                      0                      0                      0                      0                      0                      0   2.26449937896559e-16
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0

 Columns 25 through 32:

                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
   1.91419958178306e-16                      0                      0                      0                      0                      0                      0                      0
                      0   1.74228527086487e-16                      0                      0                      0                      0                      0                      0
                      0                      0   1.52749816861709e-16                      0                      0                      0                      0                      0
                      0                      0                      0   1.37245616001238e-16                      0                      0                      0                      0
                      0                      0                      0                      0   1.28865570171581e-16                      0                      0                      0
                      0                      0                      0                      0                      0   7.50766891119711e-17                      0                      0
                      0                      0                      0                      0                      0                      0   6.94131769376700e-17                      0
                      0                      0                      0                      0                      0                      0                      0   6.42272417019620e-17
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0
                      0                      0                      0                      0                      0                      0                      0                      0

 Columns 33 through 36:

                      0                      0                      0                      0
                      0                      0                      0                      0
                      0                      0                      0                      0
                      0                      0                      0                      0
                      0                      0                      0                      0
                      0                      0                      0                      0
                      0                      0                      0                      0
                      0                      0                      0                      0
                      0                      0                      0                      0
                      0                      0                      0                      0
                      0                      0                      0                      0
                      0                      0                      0                      0
                      0                      0                      0                      0
                      0                      0                      0                      0
                      0                      0                      0                      0
                      0                      0                      0                      0
                      0                      0                      0                      0
                      0                      0                      0                      0
                      0                      0                      0                      0
                      0                      0                      0                      0
                      0                      0                      0                      0
                      0                      0                      0                      0
                      0                      0                      0                      0
                      0                      0                      0                      0
                      0                      0                      0                      0
                      0                      0                      0                      0
                      0                      0                      0                      0
                      0                      0                      0                      0
                      0                      0                      0                      0
                      0                      0                      0                      0
                      0                      0                      0                      0
                      0                      0                      0                      0
   4.84239988362713e-17                      0                      0                      0
                      0   2.39198845321293e-17                      0                      0
                      0                      0   2.02231370859494e-17                      0
                      0                      0                      0   6.68357906467677e-18

>>

Det blir samma värden när jag gör min SVD-algoritm i GNU Octave också, dvs SVD som grundar sig på QR-faktorisering.

Så vad kan man dra för slutsats om detta då?
Mina slutsatser är att:
1. Jag behöver inte använda Armadillo.
2. Jag kan använda mitt C bibliotek
3. GNU Octave är bättre på beräkningen
4. Armadillo's SVD och min QR-faktoriserade SVD i C ger samma resultat, men Armadillos SVD är något snabbare och kräver mindre iterationer.
5. Float verkar inte heller göra någon större skillnad i beräkningen heller.

Då kommer jag skriva om min C-kod så att man initialiserar matriserna bara en gång under körningen.

Svenska ElektronikForumet

Matrisberäkningar med för STM32?

Re: Matrisberäkningar med för STM32?

Re: Matrisberäkningar med för STM32?

Re: Matrisberäkningar med för STM32?

Re: Matrisberäkningar med för STM32?

Re: Matrisberäkningar med för STM32?

Re: Matrisberäkningar med för STM32?

Re: Matrisberäkningar med för STM32?

Re: Matrisberäkningar med för STM32?

Re: Matrisberäkningar med för STM32?

Re: Matrisberäkningar med för STM32?

Re: Matrisberäkningar med för STM32?

Re: Matrisberäkningar med för STM32?

Re: Matrisberäkningar med för STM32?

Re: Matrisberäkningar med för STM32?

Re: Matrisberäkningar med för STM32?