Determinanta matrice u c++?
(1 korsinik/a gleda/ju temu) (1) Gost

Determinanta matrice u c++?


15.09.2011 | 11:23
Lijep pozdrav svima.

Odmah na početku isprika ako već postoji ovakva tema ili ako ne pripada ovdje. Naletjela sam na ovaj forum upisavši u google: determinanta matrice c++. Vidjela sam jednu raspravu iz 2009. slične tematike pa se nadam da ću bar ovdje dobiti odgovor.

Problem: programiranje u c++.

Ja sam početnik, potpuni. Imali smo jedan semestar programiranja u c++ i radili smo sitne programe tipa množenje, računanje kvadratnih jednadžbi, traženje simetričnosti matrice. Na kraju semestra profesor je zatražio da mu napravimo programe kao preduvjet za polaganje ispita.

Nekima je podijelio probleme sa šahovskom pločom, nekima igre riječi (palindromi), i sl.
Meni je zadao determinantu matrice jer to nisam znala na ispitu iz matematike kada smo pisali.

Ugl. sada znam riješiti determinantu kvadratne matrice na papiru, ali prevesti to u c++ mi se čini kao, hm... ne znam.


Zadatak ide ovako: "Unesi n€N, te kvadratnu matricu A. Izračunaj njenu determinantu."

Imam neku skriptu u kojoj piše kako treba započeti programe i sl. Čak ima i jedan dio o matricama. Trebalo mi je dosta da shvatim kako se na početku programa piše a[10][10], a profesor mi je napisao u mailu da ide float a[10][10], a ne int.
Mislim, koliko god se ja to trudila shvatiti, a valjda sam imbecil... šta ću, ne ide mi.

Sjedim za računalom 5-6 sati, piljim u monitor i pokušavam natjerati program da mi pomnoži elemente dijagonale matrice (iako još nisam sredila matricu da ima oblik trokutne matrice). Nikako.

I sada, jedina nada su mi ljudi s interneta, ali ovih dana su svi nekud nestali. Od 4-5 foruma na kojima sam poslala upit nitko ništa.
Nema mi smisla platiti nekome da mi izradi program jer nakon programa slijedi usmeni i profesor pita neka pitanja po kojima vidi jesmo li ili nismo savladali... kak se to zove.. gradivo? O.o

Hm... ispit je u ponedjeljak. Da, imala sam cijelo ljeto. Nisam pogledala ranije taj zadatak. Mislila sam da nam je zadao nešto što se riješi u dvije - tri for petlje (kao što smo imali na satu).o.O

Dobila sam dva primjera programa na nekim web stranicama, ali ljudi su iskusniji programeri koji rade sa pokazivačima (**).
Naš profesor zna da nemamo pojma o pokazivačima pa mi nema smisla nositi mu takav kod, što mi je uostalom i rekao: "za pokazivače bi trebali sati i sati predavanja i objašnjavanja".

No, da skratim... ako je itko uopće čitao sve ovo... O.o


____________________
Determinantu matrice treba riješiti pomoću Gaussa.
Treba ju svesti na trokutastu.
Dakle, treba prvo naći max vrijednost, zamijeniti redove tako da oni koji sadrže max vrijednost (uz funkciju fabs) dođu na mjesto prvog retka. Zatim treba elementima lijevo od dijagonale pribrojiti prvi redak pomnožen sa brojem koji u brojniku sadrži element reda kojem pribrajamo, a u nazivniku broj elementa prvog reda kojem pridajemo suprotni predznak od broja iz onog drugog reda, ako već nema. Onoliko puta koliko su redovi mijenjali mjesta toliko puta se cijeloj matrici mijenja predznak. Na kraju se množe elementi dijagonale.

Komplicirano.

Ugl. ja imam ovaj kod kao najbliže rješenju:
____________________
#include<conio.h>
#include<stdio.h>

int main ()
{
int a[10][10],maxr[10],maxs[10],minr[10],mins[10],i,j,n,k,pom[10],P,b,c,f;
printf ("upisi red matrice: ");scanf("%d",&n);
printf("upisi elemente matrice:\n");
for(i=1;i<=n;i++)
for(j=1;j<=n;j++) {
printf("a[%d][%d]: ",i,j); scanf("%d",&a[i][j]); }
for(i=1;i<=n;i++) {
minr[i]=a[i][1]; maxr[i]=a[i][1];
mins[i]=a[1][i]; maxs[i]=a[1][i];
for (j=1;j<=n;j++) {
if(a[i][j]<minr[i])minr[i]=a[i][j];
if(a[i][j]>maxr[i])maxr[i]=a[i][j];
if(a[j][i]<mins[i])mins[i]=a[j][i];
if(a[j][i]>maxs[i])maxs[i]=a[j][i]; } }
for(i=1;i<=n;i++) {
for(j=1;j<=n;j++)
printf("%d ",a[i][j]);
printf("\n"); }
for(i=1;i<=n;i++)
printf("\n najveci u redu %d je %d, a najmanji %d ",i,maxr[i],minr[i]);
for(i=1;i<=n;i++)
printf("\n najveci u stupcu %d je %d, a najmanji %d ",i,maxs[i],mins[i]);

for (k=a[1][1]; k<=n; k++);
P=1;
P=P*f;
f=k++;
printf("\nUmnozak elemenata dijagonale je:%d", k);
getch ();
}
________________________________________________

Ovaj zadnji for nije ništa doli moj očajnički pokušaj da pomnožim elemente dijagonale matrice.
Naravno, min i neke stvari namjeravam izbaciti.
Čudno mi je što program prepoznaje apsolutnu vrijednost iako nigdje nema funkcije fabs.
_________________________________________________

Od vas očekujem, ako imate volje, vremena i strpljenja da mi objasnite gdje griješim, što bi trebalo dalje. Mislim, nije mi problem pomučiti se i razmišljati, ali eto ne znam ni kuda krenuti.
Trebalo bi ubaciti pokoji potprogram jer će inače čitav kod izgledati vrlo nepregledno, ali greške samo niču (to few/many arguments i sl.)
__________________________________________________

Hvala svima!
16.09.2011 | 00:08
Pozdrav!

Imali smo cijelu trakavicu upravo iniciranu ovom temom. Unutar nje ima konkretno jedan StajpZg-ov post koji će ti pomoći.

www.jabucnjak.hr/forum/programiranje/583...-c-jeziku.html#59599
16.09.2011 | 11:01
Pogledala sam tu temu. Prekopirala sam ono što mi se učinilo važnim, ali to je i dalje daleko od rješenja. Odnosno, imam sada već neki dio programa, ali ono ključno, svođenje matrice na trokutastu stvarno ne znam kako bih.
Netko je isto napisao da je mijenjanje redaka trivijalno. Ja imam tu neke tri for petlje, ali to nikako ne daje rezultat. Ovaj kod daje max vrijednost, množi elemente dijagonale, ali ne mijenja redove. Zahvaljujem na pomoći.

Evo kod:
/*Unosi se n e N, te kvadratana matrica A reda n.
Program treba izračunati determinantnu matrice A.*/
#include<stdio.h>
#include<math.h>

int main()
{
float A[10][10],max,pom[10];
float P=1;
int i,j,n,k;
printf ("Upisi broj n: ");scanf("%d",&n);
printf("Upisi elemente polja:\n");
/*ucitavanje elemenata polja*/
for(i=1;i<n+1;i++)
for(j=1;j<n+1;j++)
{
printf("A[%d][%d]: ",i,j);
scanf("%f",&A[i][j]);
}
/*ispis u obliku tablice*/
printf("Ucitali ste elemente: \n");
for(i=1;i<n+1;i++)
{
for(j=1;j<n+1;j++)
printf("%f ",A[i][j]);
printf("\n"); /*zbog prelaska u novi red*/
}
{
max=fabs(A[1][1]);//tražim max
for(i=1; i<n+1; i++)
if(A[i][1]>max) max=fabs(A[i][1]);
printf("max=%.2f\n", max);
}

{
for (k=1;k<=n;k++)//Zamjena redova (ne radi)!!
pom[k]=A[i][k];
for (k=1;k<=n;k++)
A[i][k]=A[j][k];
for (k=1;k<=n;k++)
A[j][k]=pom[k];
printf("Prvi element matrice je: %.2f\n", A[1][1]);
}

/*Svođenje matrice na trokutastu*/

for (int i=1;i<=n;i++)//Množenje elemenata dijagonale
{
P *= A[i][i];
}
printf("Umnozak elemenata dijagonale iznosi: %.2f\n", P);
scanf(" ");
return 0;
}
  • User
  • Posjetitelj
16.09.2011 | 11:40
Obavezno razbijaj kod na funkcije, bit će lakše i tebi i nama za pročitat i snać se. Nemam baš vremena da riješim cijeli problem (ipak su ispitni rokovi ), ali ti mogu pomoć oko premještaja redaka.

Problem stvarno je trivijalan. Iskoristi sljedeću funkciju:

void zamijeniRetke(float matrica[][N], int indeks1, int indeks2)
{
    int i,j;
    float temp[N]; 
	
    if (indeks1 == indeks2) {
        return;
    }

    // premještaj redaka
    for (i=0; i < N; i++)	{
        temp[i] = matrica[indeks1][i];				
        matrica[indeks1][i] = matrica[indeks2][i];		
        matrica[indeks2][i] = temp[i];
    }
}


indeks1 i indeks2 su indeksi redaka koje zamjenjuješ, npr. ako želiš zamijenit prvi i treći redak onda je indeks1=0, indeks2=2 ili obratno. N je red matrice.

Ispod include-ova moraš definirat sljedeću konstantu da bi funkcija radila:
#define N <red matrice>


tj. ako ti je matrica npr. reda 4 onda stavljaš
#define N 4


I na kraju samo pozivaš funkciju, npr.
zamijeniRetke(matrica, 1, 2);


Bilo bi lakše s pointerima, al kad ih ne smiješ koristit...
Preporučam ti da svugdje koristiš ovu konstantu N (tj. zamijeniš n sa N) jer ćeš lako moć promijenit red matrice kad zatreba
  • User
  • Posjetitelj
16.09.2011 | 12:06
E da, ako koristiš #define N onda gubiš fleksibilnost pri upisivanju proizvoljne matrice, ali valjda ti to ni nije toliko bitno. Neku matricu kao primjer samo hard-kodiraj tako da je ne moraš svaki put unosit, jer pretpostavljam da ti je jedino bitno riješit da ti Gaussova eliminacija funkcionira kako treba
16.09.2011 | 15:01
A zašto bi morala uopće mijenjati retke da bi je svela na trokutastu?

Jednostavnije bi ti bilo ovo:

promotri element na mjestu [0][0]
     svim elementima drugog retka oduzmi vrijednost V, a V = mat[1][0]/mat[0][0]
     svim elementima trećeg retka oduzmi vrijednost V, a V = mat[2][0]/mat[0][0] 
     i tako redom za sve retke ispod prvog. 
i tako redom za sve elemente na glavnoj dijagonali osim posljednjeg. 


Malo ove konstante zamijeni varijablama i,j da bude općenito i time si je svela na trokutastu. Kasnije samo pomnoži elemente na glavnoj dijagonali et voila.
17.09.2011 | 11:11
Hvala ljudi.
Ipak, profesor traži sve standardno. Znači nema cout i cin, a mislim da bi mu bilo čudno da upotrijebim i ovaj define i sl., jer ni to nije koristio u primjerima s predavanja. Napisao mi je na komad papira kako treba izgledati zamjena redaka i mislim da tu ne smijem puno mijenjati od onoga kako je on to zadao jer će biti više nego očito onda da mi je netko drugi riješio zadatak.

Međutim, meni je netko upravo jučer riješio taj zadatak, samo što u tom kodu nema traženja max vrijednosti, ni funkcije fabs.
Nisam toliko napredna da shvatim kako to da program ipak funkcionira iako su izbačene ove dvije stavke, a mislim da će profesor isto skužiti da to nije ono što sam sama trebala napraviti.

Ugl. poslala sam mu jutros taj kod, ali se ne nadam previše oko svega toga.
Još jednom, hvala svima na pomoći. Lijepo je znati da ima ljudi koji žele pomoći nama početnicima.
17.09.2011 | 13:08
Gle, zamjena redaka je od koristi kada LJUDI koriste ovu metodu jer je cilj na elementima glavne dijagonale (osim na zadnjem) dobiti jedinice i onda s tim jedinicama dobiti sve nule u stupcima ispod njih. Tako je lakše, jer kada imaš u jednom stupcu 3 a poviše njega 1 onda ćeš cijeli taj redak s tricom "riješit" tako da od njega oduzmeš trostruku vrijednost retka koji sadrži tu jedinicu. Točno ispod će bit nula a desno šta bude. Dakle zapisala bi tu operaciju (pod pretpostavkom da je redak s jedinicom prvi a s tricom drugi) ovako II-3*I

Obzirom da ti radiš program kojeg izvodi RAČUNALO njemu je svejedno ima li jedinice ili ne, pa tako ako bi imala umjesto trice i jedinice npr. 8 i 5, operacija bi se zapisala kao II - (8/5)*I. Tebi bi to kao čovjeku bilo komplicirano za pratit (stalno neki razlomci) al računalu je to potpuno nebitno. Množio on sa 3 ili sa 8/5 - svejedno.

E sad, ako profa baš inzistira da ti napraviš model koji vjerno preslikava ljudsko određivanje determinante ovom metodom, okej, ali po meni je to ko da radiš u hotelu ko sobarica i svu posteljinu pereš na ruke a pored tebe stoji funkcionalna veš-mašina.
  • User
  • Posjetitelj
17.09.2011 | 21:48
Potpisujem john_st, a i mislim da bi vam profesori trebali bit malo fleksibilniji s obzirom da je ovo ipak doba interneta, foruma, online dokumentacije i sl. Ako si ti želiš olakšat rješenje dodatnim znanjem koje nisi dobila na predavanjima, ne bi smjelo bit razloga da to i ne učiniš, samo šta bi profesor eventualno trebao usmeno preispitat to znanje da vidi da stvarno razumiješ to što si iskodirala. Ovakvo ograničavanje ipak ne vodi ničemu.

Ali dobro, ako nisi neki informatički/računarski smjer, onda ćemo mu oprostit što vas nije naučio pointere i sl.
Moderatori: Bertone
  • Stranica:
  • 1

Vikalica™

Zadnja poruka: pred 17 sati, 19 minuta
  • jmustac: poslano
  • Riba: Javi se meni s detaljima.
  • jmustac: kako obrisati profil na ovom forumu?
  • rusty: @marioart - ovisi o hipervizoru, ali ukratko da. Moguće je dobiti emulaciju x86 na ARM bare metal stroju. UTM hypervisor kojeg osobno koristim nudi u mogućnost. Jasno, to zovu emulacijom, ne virtualizacijom ali eto
  • marioart: @rusty .. moze li se ikako dobiti x86 win kroz virtualizaciju?
  • stefanjos: drlovric. pa odgovaranje poruka na aplikacijama koje namu mac app i ako se puni da ne secem po mob, rijesim preko maca i to je to
  • rusty: Virtualizacija bi bio moj odabir, e sada da li je besplatno ili ne, ovisi o tome kako ćeš tu kopiju licencirati
  • VanjusOS: ima li nekakav besplatni način za imati Windows na Macbook sa M čipovima?
  • drlovric: Za sta ti konkretno mirroring sluzi? Ja mu nisam nasao primjenu :/
  • stefanjos: affinty na popustu 50%. ja sam si uzeo novu licensu. do sada sam bio na v1
  • stefanjos: tako da, meni je to top feature
  • stefanjos: ja imam australski app store i iphone mirroring koristim svaki dan
  • drlovric: To ces probati jednom i nikad vise. Nema razloga za biti razocaran :)
  • Zdravac: - NIJE :/
  • Zdravac: Jel proradio iPhone Mirroring? :)
  • Yonkis: iOS 18.1.1 / macOS Sequoia 15.1.1 [link]
  • kupus: pretpostavljam da to ide meni, znaci fino radi.
  • kupus: hahaha
  • drlovric: Znas kako lijepo klepece. Mozes i mijenjati zvuk. Merak pravi ????????
  • kupus: svejedno hvala smayoo
  • smayoo: Ma nemoj me zajebavat... :D
  • drlovric: Prekrasna aplikacija ako zelite zvuk fancy mehanicke tipkovnice :) [link]
  • drlovric: I ja to ocekujem. Na proljece cu u USA i iskreno se nadam da ce izaci do tada. Da donesem makar dva komada :)
  • Borisone: Osim ako nisi na ti s Tim Cookom, najbliže nagađanju kada je: [link]
  • VanjusOS: zna li netko, hoće li izaći uskoro Air sa M4 čipom?
  • smayoo: Pokušao sam, ali sam odustao kad me zavrtio u krug s registracijom i sign-inom na Epic games. Najprije kaže da može iskoristiti moj AppleID, kad napravi login preko AppleID i kažem da ne želim pokazati svoj AppleID onda me opet svejedno traži login. Tako da... sorry, ne bih.
  • kupus: hvala
  • kupus: moze li netko probat ovaj besplatan virtualan masterpiece na m macovima? [link]
  • ping: @drlovric imaš besplatne konvertere na netu pdf to excel ili konvertiraj u jpeg pa potrazi online free konverter ocr to excel
  • drlovric: nemam Acrobat....smo PDF Exprt i Wondershare.... Raspadne se na oba
  • drlovric: Daj mi mail na poruku pls, pa da ti posaljem....a kakav Pentagon, prosla su ta vremena, glupost neka :)
  • Zdravac: negdje ju uploadaj, pošalji mi link pa ak hoćeš, da vidimo kak će ispasti (osim ako nisu neke secret Pentagon pi*darije :)
  • Zdravac: @drlovric - "export to excel" ?! Adobe acrobat
  • drlovric: Imam fakturu u PDF, htio bih od nje napraviti Excel file. Sta da radim? :)
  • mariolino1971: Pozdrav. Pokušavam poslati oglas no stalno dobivam poruku 413 Request Entity Too Large. U čemu je problem?
  • kupus: UK korisnici vs Apple [link]
  • Yonkis: Jeste li znali da se iPhone sam reboota ako nije otključan u 72 sata? [link]
  • Yonkis: jeste čitali komentare na [link]
  • Ender: takodje zainteresovan za neki IPTV, predlog moze u PM.
  • cariblanco: Ako hočeš ja ti svog "provajdera" kojeg koristim 4-5 godina, pošaljem na pp ?
  • cariblanco: Kako veli Zdravac dobar iptv za 10 E i imaš sve i gdjedaš od kud hočeš. Meni je to jako važno pošto puno putujem, trebam samo dobar net i mirna bosna
  • drlovric: Moj Telemach nema Arene vec Sport Klubove....ali postoji kabelska u zgradi koja ima Arene, preko koaksijalnog kabela.... Pa sam se tehnloski vratio deteljece unatrag. Ali radi evo :)
  • Zdravac: Ali to je paket ex yu kanala, nema samo sport vec sve
  • Zdravac: Kvalitetni, koji mogu preporučiti 10 eur/mj
  • Zdravac: Ima, iptv
  • kupus: hometv to go ili kako već
  • kupus: znači ide i eronet aplikacija
  • kupus: u bih prava ima arenasport ako se ne varam
  • drlovric: Jel ima neki nacin da se HNL Supersport liga gleda online? Neki app? Streaming? Iptv?
  • cariblanco: Prošlo iz petog pokušaja...bit će da je Riba na suhom ;)

Za vikanje moraš biti prijavljen.

Prijava

Prisutni jabučari

Novo na Jabučnjaku

Teme

Poruke

Oglasi

Komentari

Anketa

Koji Mac koristite?

Page Speed 1.72 Seconds

Provided by iJoomla SEO