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 2 sati, 3 minuta
  • Performita: Kalimera, Dobro jutro... ovo za Apple Business rolls zvuči kao izvrsna vijest!!!
  • Borisone: Apple Business rolls out to 200 countries today - brace for impact - even Croatia!
  • Yonkis: iOS 26.4.1 and iPadOS 26.4.1 This update has no published CVE entries.
  • stefanjos: Treba mi as u rukavu kada me netko pita gdje da ide haha :D
  • smayoo: Ima ih nekoliko koji se oglašavaju kao iPhone serviseri, nijedan mi nije baš ulio povjerenje, ali nemam neposrednih loših iskustava.
  • stefanjos: gdje nositi iphone na popravak u rijeci kada nema jm elektronike?
  • Yonkis: @smayoo Nema na čemu, evo ja upisah.
  • smayoo: @Yonkis, evo, otključano, hvala na upozorenju!
  • dpasaric: Booo?
  • Yonkis: “ Lažu li nas :) ...nastavak sa vikalice.” - tema zaključana a zadnji post završava sa pitanjem?
  • smayoo: Sretan Uskrs!
  • Yonkis: Sretan Uskrs braćo! I sestre!
  • zoranowsky: Sretan Uskrs!
  • Gjuroo: Sretan Uskrs!
  • Miro Spiro: sretan Uskrs
  • user: da bar je propao kao što su Microsoft i ostali predviđali, hvala što si pokrenuo najgori trend u techu Steve.
  • Yonkis: [link] Who’s in Apple’s top 10? Here’s the full list of the most influential people of all time
  • Yonkis: [link] netko fino složio povijest svih modela iPhonea, jako lijepo.
  • Yonkis: Apple employee #8 Chris Espinosa on working his whole life at one company: [link]
  • Ender: vrlo moguce, ja video to na nasem sajtu 02. apr. sorry.
  • Yonkis: @ender članak je prekopiran odavdje [link] i po tome ispada da je 1.april šala.
  • Ender: AMD kupuje Intel: istorijski preokret koji menja čitavu x86 industriju [link]
  • dpasaric: Super je prošlo! :) Sretno im!
  • miomika: Možda glupo pitanje, dal će se vidjeti putovanje te rakete iz .hr prema mjesecu? zasada nije oblačno
  • Ender: USA lansiranje rakete sa ljudskom posadom ka mesecu
  • Ender: [link]
  • Ender: [link]
  • CoffeePod: @yonkis iStyle takvo ruglo od firme. Pogotovo dućan u Splitu. Osoblje nema pojma, morao sam ih sve školovat :P
  • dpasaric: Večeras tko je u Rijeci u Piki na proslavu 50 godina! :)
  • Borisone: Onaj dan u godini?
  • Yonkis: zato iStyle daje popust od 2.27% ne od 3% ne od 2.5% nego 2.27% Presmiješno.
  • Tito: Kupio sam 3 MBA i 2 iphonea
  • ZoPaj: 50% off all Apple products in all Apple stores for Apple's 50th anniversary... today only! apple.com
  • dpasaric: [link]
  • dpasaric: 50 godina Applea u srijedu u 19.30 u Pikiju!
  • jura22: Nema mjesta za ekspanziju jer je tako Apple odlucio. Sigurno bi se mogla graficka ustekati ili vise njih. Mozda memorija i nema smisla.
  • vjeko1: Realno me ne čudi ovo za Mac Pro jer desktop lineup nije imao previše smisla sa Studio i Pro. S druge strane, arhitektura sa Apple Silicon je takva da nema previše mjesta za expanziju tako da.. RIP
  • Piko: Samo da velim da sam "ukrao" ideju za vikalicu .... i napravio svoju implementaciju na [link] ....
  • Soffoklo: Apple odustaje od Mac Proa - [link]
  • Yonkis: @CoffeePod koristiš li ono punjenje do 80% ili si light user? Moj iPhone 15 je nakon 26 mj izgubio samo 6% a 2 različita SIMa (priv posl). Tako da 1% je sasvim ok.
  • drlovric: @Coffe: Po gruboj racunici moze ti trajati do kraja stoljeca taj telefon. Nista ne azuriraj :)
  • CoffeePod: Ja sqm i dalje na iOS18 i ne planiram jos updejtat. Nakon 9 mjeseci baterija na 16PRO pala je samo za 1% od max kapaciteta
  • ZoPaj: eh, da je bar više tih OS updatea bez novih featuresa:) optimiziranih i očišćenih od bugova pored famoznih sigurnosnih patcheva.. mislim da Snow Leopard nije imao ni jednu jedinu novost? 26.4 mi se čini brz i stabilan...
  • miomika: 26.4 je vani, ali ništa od Geminija. Samo APlle Music i fora s koncertima. Neznam dal je vrijedno i skidati
  • Borisone: [link]
  • Borisone: Jel to Apple MD dostupan od 14.4. u HR? [link]
  • zoranowsky: @rak21, ja nisam, upravo sam dobio info da mi je narudžba stornirana zbog ograničenih količina. Netko je očito zeznuo kalkulaciju...
  • Yonkis: WWDC26 June 8–12, 2026 [link]
  • rak21: Je li netko ulovioacbook pro M5 pro za 2200 u instaru? Baš me zanima jel prošla narudžba ili su stornirali jjer su bili jeftiniji 500€ i više od drugih.
  • smayoo: Ja sam nedavno imao tako raspadnut prikaz oglasa u Firefoxu, ispalo je da je problem bio u Ghosteryju

Za vikanje moraš biti prijavljen.

Prijava

Novo na Jabučnjaku

Teme

Poruke

Komentari

Anketa

Koji Mac koristite?

Page Speed 0.67 Seconds

Provided by iJoomla SEO