Matrix Solver
 x   TIP: Přetáhni ikonu na hlavní panel pro připnutí webu

Matrix SolverMatrix Solver

 
Hledat
Vybavení pro Laser Game
Spuštěn Filmový magazín
Laser Game Brno
Laser Game Ostrava

Matrix Solver

Google       Google       22. 4. 2009       16 496×

Následující mikro-článek není vlastně článkem, ale jakýmsi uvedením zdrojového kódu, který má za úkol upravovat matice zadané uživatelem. Možná si řeknete, že ve světě, kde si lze z torrentů zdarma stáhnout SW typu Matlab, je to zbytečné, ale sami uznáte, že program vytvořený pro radost někdy potěší víc než ten "poctivě" ukradený. Uvedený algoritmus navíc narozdíl od uvedeného software nevypisuje pouze výsledek, ale také postup.

Reklama
Reklama

Program Matrix Solver umožňuje zadat vstupní matici příkazem new [a11, a12, a13...; ... aij]. Příklad:

new [1,2,3;3,4,5;5,6]

Přičemž nevyplněné položky v řádku budou nahrazeny nulami.

Jako oddělovače položek na řádku jsou použity čárky, ale akceptovány jsou mezery (tak jako v Matlabu). Jako oddělovače řádků používáme středníky (opět převzato z Matlabu). Po zadání matice je tato matice vypsána a lze s ní dále pracovat.

V tuto chvíli je implementována jediná transformace, a to úprava na horní trojúhelníkový tvar (tedy s nulami pod hlavní diagonálou). Algoritmus vychází z hledání nejmenšího společného násobku (pomocí Eulerovy metody hledání největšího společného dělitele). Algoritmus by bylo vhodné dále doplnit o přehazování řádku, například srovnat podle počtu nul v řádku a dále například přehodit sloupce tak, aby byly nuly převážně v pravé straně matice. Budete-li mít vytvořeny nějaké zajímavé úpravy, určitě se o ně podělte. Opět bych je publikoval na tomto serveru a kdo ví, časem si třeba vytvoříme pěkný SW balík :).

Přikládám kompletní zdrojový kód tak, aby si kdokoliv, kdo by měl nějakou připomínku, mohl zdrojový kód upravit. Projekt je vytvořen ve Visual Studiu 2008 (CS Express). Prosím nepište, že uvedený zdrojový kód je k ničemu, i když si to třeba myslet můžete. Začátečníkům se může hodit. Podaří-li se vám zdrojový kód nějak zajímavě upravit, případně přidat nějakou další transformaci, určitě si to nenechte pro sebe a pošlete zdrojový kód.

Zdrojový kód hlavní části:

using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;

namespace ConsoleApplication1
{
    class Program
    {
        static void Main(string[] args)
        {
            string command;
            int[,] MATRIX=null;
            
            Console.Clear();
            Console.WriteLine("Program pro úpravu matic");
            Console.WriteLine("Autor: Čtenáři serveru Programujte.com \n\n");
            do
            {
                Console.Write(">> ");
                command = Console.ReadLine();                

                command = command.ToLower();
                // pokud příkaz začíná slovem NEW, pak zřejmě půjde o zadání nové matice
                if (command.StartsWith("new"))
                {
                    MATRIX = ParseNewMatrix(command);
                    continue;
                }
                // výpočet matice s nulami pod hlavní diagonálou
                if (command == "diagonal")
                {
                    MATRIX = Diagonal(MATRIX);
                    continue;
                }
                // smaže obrazovku
                if (command == "cls")
                {
                    Console.Clear();
                    continue;
                }
                Console.WriteLine();
            }
            while (command != "exit");
        }

        // zpracuje příkaz NEW pro vložení nové matice
        static int[,] ParseNewMatrix(string cmd)
        {
            char[] col_separator = { ' ', ',' };
            StringSplitOptions col_split_opt = StringSplitOptions.RemoveEmptyEntries;

            // najde začátek matice
            int zacatek_matice = cmd.IndexOf('[');
            int konec_matice = cmd.IndexOf(']');
            // uloží si matici do samostatného řetězce
            string matrix = cmd.Substring(zacatek_matice+1, konec_matice - zacatek_matice -1);
            // rozloží matici do řádků
            string[] radky = matrix.Split(';');
            // zjistí počet řádků
            int pocet_radku = radky.Count();
            // zjistí maximální počet prvků v řádku
            int pocet_prvku = 0;
            for (int i = 0; i < pocet_radku; i++)
            {
                string[] prvky = radky[i].Split(col_separator, col_split_opt);
                if (pocet_prvku < prvky.Count()) pocet_prvku = prvky.Count();
            }
            // vytvoří matici pro uložení výsledků
            int[,] LOCAL_MATRIX = new int[pocet_radku, pocet_prvku];
            // projde všechny řádky
            for (int y = 0; y < pocet_radku; y++)
            {
                // v každém řádku rozdělí položky na jednotlivé buňky
                string[] prvky = radky[y].Split(col_separator, col_split_opt);
                // každou buňku parsuje...
                for (int x = 0; x < pocet_prvku; x++)
                {
                    // pokud byla buňka zadána, vloží ji do matice
                    if (x < prvky.Count())
                    {
                        // pokud je prvek parsovatelný, zapíše se, pokud není, pak se zapíše nula
                        if (!int.TryParse(prvky[x], out LOCAL_MATRIX[y, x]))
                            LOCAL_MATRIX[y, x] = 0;
                    }
                    // pokud nebyla, nahradí ji nulou
                    else
                        LOCAL_MATRIX[y, x] = 0;
                }
            }
            // zobrazí zadanou matici
            ShowMatrix(LOCAL_MATRIX);
            return LOCAL_MATRIX;
        }
        
        // zobrazí matici
        static void ShowMatrix(int[,] LOCAL_MATRIX)
        {
            // zjistí rozměry matic
            int y_max = LOCAL_MATRIX.GetLength(0);
            int x_max = LOCAL_MATRIX.GetLength(1);
            // vypíše matici
            Console.WriteLine("Matice:");
            for (int y = 0; y < y_max; y++)
            {
                for (int x = 0; x < x_max; x++)
                    Console.Write("{0,8:D}",LOCAL_MATRIX[y,x]);
                Console.WriteLine();
            }
            Console.WriteLine();
        }        

        // přetransformuje matici na matici s prvky na hlavní diagonále
        static int[,] Diagonal(int[,] LOCAL_MATRIX)
        {
            // pokud není co zpracovávat, vyhodí chybové hlášení a vyskočí
            if (LOCAL_MATRIX == null)
            {
                Console.WriteLine("Matice je prázdná!!!\n");
                return null;
            }
            // zjistí rozměry matic
            int y_max = LOCAL_MATRIX.GetLength(0);
            int x_max = LOCAL_MATRIX.GetLength(1);
            int[,] TMP_MATRIX = LOCAL_MATRIX;

            int[,] nasobitele = new int[y_max,2];
            // pro celou šířku diagonály provede výpočet...
            
            for (int iii = 0; iii < x_max-1; iii++)
            {
                bool provedena_operace = false;
                // na všech řádcích pod diagonálou
                for (int i = 1+iii; i < y_max; i++)
                {
                    // pokud je položka nulová, sama o sobě přeskočí...
                    if (LOCAL_MATRIX[i, iii] == 0) continue;
                    // vypočte násobitele řádků
                    nasobitele[i, 0] = MyMath.NejmensiSpolecnyNasobek(LOCAL_MATRIX[iii, iii], LOCAL_MATRIX[i, iii]) / LOCAL_MATRIX[iii, iii];
                    nasobitele[i, 1] = -MyMath.NejmensiSpolecnyNasobek(LOCAL_MATRIX[iii, iii], LOCAL_MATRIX[i, iii]) / LOCAL_MATRIX[i, iii];
                    // vypíše postup úprav
                    Console.WriteLine("Postup: {0:D}xřádek {1:D} + {2:D}xřádek {3:D}", nasobitele[i, 0], iii+1,nasobitele[i, 1], i + 1);
                    // propočítá všechny položky v řádku
                    for (int ii = 0; ii < x_max; ii++)
                        LOCAL_MATRIX[i, ii] = nasobitele[i, 0] * LOCAL_MATRIX[iii, ii] + nasobitele[i, 1] * LOCAL_MATRIX[i, ii];
                    provedena_operace = true;
                }
                if (provedena_operace) ShowMatrix(LOCAL_MATRIX);
            }
            return LOCAL_MATRIX;
        }
    }
}

Zdrojový kód matematické knihovny:

using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;

namespace ConsoleApplication1
{
    class MyMath
    {
        // spočítá největší společný dělitel celých čísel
        public static int NejvetsiSpolecnyDelitel(int a, int b)
        {
            int zbytek = 0;
            do
            {
                zbytek = a % b;
                a = b;
                if (zbytek != 0) b = zbytek;
            } while (zbytek != 0);
            return b;
        }

        // spočítá nejmenší společný násobek
        public static int NejmensiSpolecnyNasobek(int a, int b)
        {
            return (a * Math.Abs(b)) / NejvetsiSpolecnyDelitel(a, b);
        }
    }
}
Odkaz na zdrojový kód.

×Odeslání článku na tvůj Kindle

Zadej svůj Kindle e-mail a my ti pošleme článek na tvůj Kindle.
Musíš mít povolený příjem obsahu do svého Kindle z naší e-mailové adresy kindle@programujte.com.

E-mailová adresa (např. novak@kindle.com):

TIP: Pokud chceš dostávat naše články každé ráno do svého Kindle, koukni do sekce Články do Kindle.

Hlasování bylo ukončeno    
0 hlasů
Google
(fotka) Ondřej KarasAutor se věnuje vývoji řídících a komunikačních systémů a tvorbě podpůrných aplikací.
Web    

Nové články

Obrázek ke článku Konference: Moderní informační systémy podporují automatizaci

Konference: Moderní informační systémy podporují automatizaci

Současná situace v šíření onemocnění Covid-19 klade na řadu firem nové nároky a mnohé z nich jsou nyní více než kdy jindy závislé na nejmodernějších informačních technologiích. Proto i v oblasti podnikových informačních systémů vidíme rostoucí důraz na automatizaci nebo na důslednou integraci. Také o těchto trendech se bude mluvit na konferenci Firemní informační systémy, která se koná 24.9.2020 v pražském Kongresovém centru Vavruška na Karlově náměstí.

Reklama
Reklama
Obrázek ke článku Nebezpečí ukrytá v USB: z nuly na škvarek za pět sekund

Nebezpečí ukrytá v USB: z nuly na škvarek za pět sekund

Za cenu šesti dolarů lze celkem bez obtíží koupit nový, líbivě vyhlížející flash disk. Přidaná hodnota, které se vám spolu s ním dostane, už tak moc líbivá není. To, co se před pár sekundami tvářilo jako externí disk, se po připojení k počítači změní v důmyslné elektrické křeslo, které vaše zařízení v onen příslovečný škvarek promění za pár sekund. Cílovou skupinou pro koupi takových zařízení by mohli být záškodníci, kteří by tímto způsobem osnovali pomstu třeba vůči záletnému partnerovi. 

Obrázek ke článku Znalosti, dovednosti i prestižní titul MBA: Jde to i moderně a online

Znalosti, dovednosti i prestižní titul MBA: Jde to i moderně a online

Snad nikdy není špatná příležitost na investici do hodnotného vzdělání. Obzvlášť v případě, že absolvent dovede teoretické poznatky přetavit v praktické dovednosti, využitelné při řešení problémů i v komunikaci. Právě na to se specializuje studijní program MBA Řízení informačních technologií, vyučovaný na Business Institutu.

Obrázek ke článku Coding Bootcamp Praha: Obor IT krize nepoznamenala, žádaní jsou weboví vývojáři

Coding Bootcamp Praha: Obor IT krize nepoznamenala, žádaní jsou weboví vývojáři

Pandemie Covid-19 otřásla trhem práce v základech. Dopady krize pocítilo celkově až 45 % zaměstnanců. Není divu, že čím dál větší jistotu přináší obor IT. Ten zůstal krizí téměř nepoznamenán a při nutnosti začít dělat věci na dálku se ještě více ukázalo, jak moc mnohé firmy kvalitní IT potřebují. Do IT nyní přicházejí začátečníci, kteří v něm vidí lukrativní budoucnost a jistotu, ale i freelanceři a zaměstnanci z oborů zasažených krizí

Reklama autora

Hostujeme u Českého hostingu       ISSN 1801-1586       ⇡ Nahoru Webtea.cz logo © 20032020 Programujte.com
Zasadilo a pěstuje Webtea.cz, šéfredaktor Lukáš Churý