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 7 sati, 17 minuta
  • smayoo: @mbp2009, imaš točno tu temu u forumu od prije par mjeseci
  • smayoo: Natrijeve baterije postoje u tzv. REDOX izvedbi već 30 godina, prvenstveno za stacionarne sustave, ali nisu bile pogodne za mobilne jer su bitno manje gustoće pohrane od litijevih, a i veće mase, jer su tekućinske. Velika im je prednost to što su ekološki neutralne, jednostavne za proizvesti i požarno nerizične.
  • m@xym: @Zdravac, OD otkrića nečega u laboratoriju do dolaska na tržište (ako se to upće dogodi) prođe između 5 i 10 godina
  • drlovric: O da. Jako veliko dan. Ljudi nisu svjesni bitnosti ali vrijeme ce polazati.
  • Zdravac: @drlovric, jesu, baš je lijepo danas pratiti HRT4 i live prijenos! :)
  • Zdravac: jel senzacionalizam ili ima nešto u ovoj tehnologiji da bi uskoro moglo i u praksu?
  • Zdravac: [link]
  • drlovric: Stigli Rafali ❤️????
  • mbp2009: Ekipa, treba mi preporuka za neki usb-c dongle za M1. Potreban je ulaz za SD karticu, HDMI ili DP (TB2), barem 2 Usb-a i ako je moguće jedan usb-c
  • smayoo: Funkcionalna zamjena je push notifikacija koja preko 5G ide po generičkom protokolu, samo aplikacija na mobitelu mora to podržavati
  • smayoo: Upozorili su me ljudi koji rade u toj branši, radi nekih telemetrijskih rješenja gdje nadzorni sustav SMSom šalje obavijesti korisniku
  • ZeljkoB: smayoo, odakle ti ova informacija za SMS?
  • smayoo: SMS poruke su phase-out. Znači, neki ih operateri više ne nude, niti omogućavaju njihovo slanje. Određen je neki krajnji rok, nakon kojeg više nijedan operater u EU nije dužan osigurati propagiranje SMS poruka. Probaj ručno promijeniti operatera (pretpostavljam da si u roamingu)
  • Matija klc: Ekipa, pomoc. Nalazim se u Rumunjskoj. Primam poruke normalno, no ne mogu poslati poruku. Pricam o obicnoj. Sve drugo radi najnormalnije. Sta napraviti?
  • MoschuS: Ma svugdje ali sve neka mala mjesta ili sela po Hr ili Slo
  • JOHN: Jesu to oni iz Imotski?
  • MoschuS: Ako netko kupuje rabljeno racunalo pripazite se prevaranata na Njuskalu i Bolhi. Otkrio sam ih par, slicni oglasi i kada ih kontaktirate salju iste/slicne racune i serijske brojeve kao potvrdu i onda odjednom obrisu oglas, naprave novi account i pokusavaju dalje uloviti nekoga.
  • Zdravac: :) Na kraju će ispasti da nas ima dosta koji vozimo motore!
  • rkrosl: prodao triumpha, kupio triumpha
  • JOHN: Honda CB650R
  • FiFi: Skoro sam ga jednom zamijenio za Bandita 1200 , ali sam ipak odustao.
  • Zdravac: ...ili Suzuki Bandit :)
  • Zdravac: E, taj motor je nekad bio moj san, da nije uletila yamaha, sigurno bih toga sad vozio :)
  • FiFi: Honda CB900F2 Hornet, 2002. Evo radim servis i nadam se registraciji prvi put nakon tri godine stajanje (djeca ti promijene raspored vožnje)
  • Riba: Problem je sto "sezona" u Irskoj trtaje tjedan dana. :-D Moram vidjeti ima li kakav nacin da ga registriram u Hr.
  • Zdravac: @ZeljkoB Odličan motor, takvi motori su baš po mom guštu! Lijep i ugodan za voziti.
  • Zdravac: @Riba Pa što čekaš, daj, zamjeni ulje, filtere i svječice, napuni bateriju, i lagano pranje. To će biti dovoljno, siguran sam! :)
  • Zdravac: Ajd, super, nisam jedini! Sad ide sezona, i gušt! Polako plan za moto trip na more...
  • ZeljkoB: Honda NC750X, ali sve manje vremena imam za voznje.....
  • Riba: Imam ja BMW 650CS, mozda ga reaktiviram nakon 10 godina stajanja u garazi. A imam i Vespu iz 1959, ta malo duze stoji. :)
  • Zdravac: Evo, ja ću prvi: YAMAHA XJ900 S Diversion, classic iz 1998., ide ko švicarski sat, čista mehanika, pouzdano do bola! :)
  • Zdravac: Jel imamo ovdje kolega ili kolegica, da voze motor? Ako da, koji?
  • zoranowsky: @air, u biti isti je pricnip, samo si sam sebi target :)
  • Air: Polako, polaku, kuda žurba, daj mi prvo link na to kako pročitati vlastite misli :)
  • smayoo: Za one koji ne znaju čitati misli - [link] :D
  • Borisone: Za one koji ne znaju čitati misli, koja je konfiguracija?
  • VanjusOS: pozdrav, treba li netko Mac mini M2, nekorišteno?
  • kupus: Ima takvih naravno. Iz mog skromnog iskustva takvi su u manjini.
  • smayoo: većina online trgovaca odbija slati na poste restante pa čak i na poštanski pretinac
  • kupus: zasto na osobnu adresu*
  • kupus: Zasto na adresu? Naruci na adresu postanskog ureda, tamo ce te cekat
  • drlovric: Imamo li koga iz Zupanje na forumu da bih mogao naruciti paket na njegovu adresu? :)
  • ecvis17: ok, tnx
  • smayoo: Piši u forum
  • ecvis17: pozdrav svima, imam stari MBP koji maksimalno moze voziti BigSur ... malo nisam vise u toku, gdje se danas mogu pronaci Combo installeri? Tnx
  • cariblanco: Čekao je da Hajduk osvoji prvenstvo i na kraju je odustao ;)
  • Yonkis: Od sada ćemo mjeriti vrijeme pre-JM and post-JM time
  • JOHN: Ako je vjerovao, neka mu dragi Bog pomogne i uzme sebi. Ako nije, opet nema veze……sve u svemu mislim da je živio dobar život!
  • smayoo: Ljudi... ode nam Joža Manolić... Nije šala.
  • accom: ne nadjem...

Za vikanje moraš biti prijavljen.

Prijava

Prisutni jabučari

Air, ammadeus, damir1, Ender, Riba, Anonimci (601)

Novo na Jabučnjaku

Teme

Poruke

Anketa

Kupujete li profesionalni Mac?

Čekam novi modularni Mac Pro - 48.5%
Novac nije problem, kupujem iMac Pro - 0.7%
Kupujem Valjak, baš je lijep i tih! - 0%
Kupujem polovni Mac Pro tower - 11.8%
Nadogradit ću postojeći Mac Pro tower - 2.9%
Običan iMac 27" mi je dovoljan za posao - 5.9%
Skromnih sam potreba, Mac mini je zakon! - 7.4%
Radim na terenu, mora biti MacBook Pro - 3.7%
Ne diram ništa, stari Mac služi me odlično - 10.3%
Kupujem PC kantu i prelazim na Windowse! - 8.8%

Ukupno glasova: 136
Anketa je završena dana: 08 Svi 2018 - 12:17
Page Speed 1.21 Seconds

Provided by iJoomla SEO