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 5 sati, 12 minuta
  • drlovric: Ma ide to bez problema. Ja sam uvijek azurirao sve strojeve prvi dan. Nisam imao problema, ali generalno nije pametno :)
  • drpongrac23: [link]
  • drpongrac23: Izgleda da izlazi oko 19.00 sati po našem vremenu. Očekuje ze velika gužva prilikom downloada
  • dpasaric: Možda za neki nebitni stroj za početak...
  • dpasaric: Ne znam, ja ne bih žurio. :)
  • drpongrac23: Sutra macOS Tahoe!
  • dpasaric: Osvrt je na naslovnici...
  • dpasaric: Ja bih samo Siri koja radi kako treba.
  • stefanjos: pa nije da ikome treba, ali radije bi to nego iakkav AI image generator haha
  • drlovric: Top je stvar. Razmazi oci i plus dobijes par sat i baterije usput.
  • dpasaric: Kome to uopće treba? :)
  • drlovric: Slazem se. Natjerala ih AI blamaza a ne dobra volja :)
  • stefanjos: 120hz na obicnom iphoneu. Napokon!
  • stefanjos: Drpongrac. Nema uvrede. Autoagnostik sam ali volim svoj polo :)
  • drpongrac23: Bez uvrede, isprika svakome, imao san bekutanera kockici m3, golfa 4 i polo sdi - sve redom šrot!!!!
  • drpongrac23: @stefanjos - fauwe = krampwagen! Hyundai rulez!!!!
  • drpongrac23: Uskoro event! Euforija počinje!!!! Gemišti na hlađenju!!!!
  • stefanjos: @yonkis, a stvarno je tino haha. samo fauve :P
  • Borisone: [link]
  • Borisone: Kao i svake godine u rujnu, nova generacija iPhona, te eventualno AirPodsi, Apple TV, i sl.
  • Yonkis: A što se predstavlja?
  • dpasaric: Danas Appleov Keynote u 19 sati po našem vremenu!
  • Yonkis: tino1 mora pješke jer ga izdao Bayerische Motoren Werke :)
  • stefanjos: ?
  • Yonkis: [link] tužno
  • smayoo: Rijeka je treći grad u Europi po godišnjoj precipitaciji. Prva je Genova u Italiji, a drugi je Nikšić u Crnoj Gori
  • smayoo: Nije istina. London ima više kišnih dana godišnje, ali Rijeka ima veću godišnju precipitaciju (kumulativnu količinu padalina u mm po m2)
  • hal9000: @yonkis • Navodno je.
  • Yonkis: Jel' ovo istina za Rijeku? [link]
  • smayoo: Riječani i ostali zainteresirani, kavica i ćakula ovu subotu (30.08.) ujutro u 10h u gradskoj kavani Renard
  • drlovric: @vanjus: Barem si imao mogucnost izbora. Odaberes banku koja ima Apple Pay i cao :) biH je doslovce prije tjedan-dva dobila Apple Pay, nije uopce bilo dostupno u drzavi.
  • VanjusOS: pa znam da je Apple pay bio dostupan preko PBZ-a već godinama, ali tek od jučer je to moguće i preko HPB.
  • Yonkis: @smayoo - mogli bi subotu ujutro. Ima li još zainteresiranih, gdje je tino1?
  • Yonkis: @VanjusOS - to misliš kao "napokon" tj sarkastično?
  • smayoo: @yonkis - ja mogu u subotu ujutro. U nedjelju popodne-navečer možda također, ali nisam siguran.
  • VanjusOS: ljudi, tehnološka revolucija! od danas Apple Pay i preko HPB-a!!
  • Yonkis: Jel netko za kakvu kavu ovaj vikend? U Ri mislim :)
  • Yonkis: Glede Tita.. u Splitu: [link]
  • smayoo: @drlovric, nema na čemu, drago mi je. :)
  • mbp2009: Titini kriteriji: onaj koji je bio protiv njega. Kao i kriteriji svakog drugog kriminalca i diktatora :)
  • drlovric: @Smayoo: Isli u plavo bijeli. Sasvim pristojna hrana. Nismo uopce razocarani. Hvala :)
  • jura22: Da, tocno tako. Sakupljao po ulici! Zene, djecu, starce. Po redu. Bez kriterija.
  • smayoo: U samom Osijeku, nažalost, po mom iskustvu ništa nije baš trajno na visokom nivou. Ako bih morao preporučiti bar jedno mjesto, to bi bio "Plavo-bijeli". Klasični restoran, kad vidiš kako su konobari odjeveni, izgleda kao da su još osamdesete, ali hrana je konzistentno vrlo dobra. Ali to je nekakav standardni restoranski jelovnik, nema nekih baš tradicionalno slavonskih jela
  • smayoo: Varga u Bilju, pojeli smo zbilja dobar čobanac prije koju godinu. Još je bolji u Sokaku u Đakovu, kod ergele, ali to ti nije blizu Osijeka
  • drlovric: Sutra sam u Osijeku. Molim lokalce za preporuku za dobru slavonsku hranu.
  • Air: U protivnom je imao vanzemaljsku sreću da je baš smaknuo one koji su mu smetali.
  • Air: @jura22 ja sam nekoj logici siguran da je Tito imao kriterije po kojima je to činio.
  • smayoo: Kako misliš - da je nekog dao smaknuti bez kriterija? Skupljao je nasumično ljude po ulici pa ih vodio na smaknuće?
  • Yonkis: Što je kriterij?
  • jura22: “kriterija”

Za vikanje moraš biti prijavljen.

Prijava

Prisutni jabučari

Novo na Jabučnjaku

Teme

Poruke

Oglasi

Komentari

Anketa

Koji Mac koristite?

Page Speed 1.08 Seconds

Provided by iJoomla SEO