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 1 dan, 2 sati
  • Yonkis: [link]
  • Yonkis: company recently confirmed a special announcement for March 4. Among the expected products are a low-cost MacBook and a new generation of Studio Display.
  • Gjuroo: Pokoj mu duši
  • smayoo: RIP :(
  • Air: RIP, otišla još jedna legenda -Robert Duvall
  • jura22: Navlakusa za novi macOS.
  • imenso: Nije free, ima pretplatu, neke aplikacije i posebnu za iPad
  • Yonkis: Apple Creator Studio - free - [link]
  • drlovric: Bas tako Mario. Ja sam bio u USA i kupio MBP M1 Pro 16" i najmanju potrebu nemam da ga mijenjam. Pogrijesio sam jedino sto nisam uzeo 14" :)
  • marioart: uredjaji su predobri i predugo traju :)
  • Riba: Zalosno je vidjeti kako mac prodaja stagnira.
  • Yonkis: Ili ovako nešto? [link]
  • Yonkis: Gdje staviti ovakve vijestice? [link]
  • gladhr2: [link] ima li netko pametan da mi malo pomogne sa ovim problemom..?
  • dpasaric: Nema indoor mapa kod nas AFAIK.
  • Yonkis: Jeste li znali za bolji anti-reflective sloj na 17ici? [link]
  • Yonkis: Ima li u Rijeci ili igdje u HR ovo implemetirano? [link]
  • Yonkis: Novi AirTag: [link]
  • Yonkis: Bijah u Smart Doctoru, zamjenska je 100€ a originalna je 170€
  • smayoo: @Yonkis, pitao si kome vjerujemo :)
  • marioart: @Yonkis ovisno o aparatu koji imas, nije ti zgorega ni provjeriti cijene u sluzbenom servisu... ja sam u mediamarktu u njemackoj platio zamjenu za 13pro 88eura ... spojis ugodno sa korisnim i u Graz u shoping i servis
  • Yonkis: A ovi u Rijeci, Smart Doctor? Evo linka:[link]
  • Yonkis: Novi firmware ako nekome treba: [link]
  • Yonkis: Nevjerovatno da je cijena skočila na 150€ a znam da je prije 2 godine bila cca 70€. Prestrašno!
  • smayoo: @Yonkis - nažalost, u Rijeci više nikom. :( U Zagrebu je Davor više puta spomenuo da ima dobra iskustva s Agramservisom
  • VanjusOS: pa ti vratim
  • VanjusOS: ako imaš kabl kojeg trebaš isprobati možeš poslati meni da isprobam spojiti Mac pro
  • VanjusOS: i ja tražim već thunderbolt 2 kabel, iz Osijeka sam
  • zoranowsky: Jel ima netko u ZG za posuditi thunderbolt 2 kabel? Imam issue sa minijem M4 i thunderbolt 27 displayem pa trebam testirati neke kombinacije...
  • Yonkis: Kome vi u Rijeci vjerujete za zamjenu baterije na iPhonu jer u iStyleu je 150 € i čeka se par tjedana
  • Yonkis: :(
  • smayoo: [link]
  • smayoo: RIP Scott Adams :(
  • Gjuroo: Za sad to može zvučati i kao Affinity priča. Not great but not that bad.
  • Soffoklo: Za to vrijeme Apple se lagano pretvara u Adobe - [link]
  • drlovric: Sve je bolje od ove blamaze sa Apple AI i cekanja pametnije Siri koja nikako da dodje :)
  • IkoIko: Apple i Google sklopili ugovor o korištenju Gemini AI... Zanimljiv razvoj situacije
  • ZeljkoB: Spigen - [link]
  • Riba: Moza zato da mozes koristiti BT tipkovnicu i mouse i dok je zaklopljen?
  • user: glupo: zatvorim laptop, BT i dalje upaljen, što ne bi trebalo ovo **** se ugasit kad ga poklopim
  • rusty: date hack je samo na par sekundi. vratiš datum, startaš app čisto da se digne i vratiš datum opet. To je to. Čisto da loada profile i postavke
  • smayoo: Date hack nisam htio raditi jer potrga druge stvari. Krenuo sam instalirati staru verziju LogiOptions preko homebrewa, ali riješili su problem prije nego sam dospio to napraviti
  • rusty: Osobno mislim da je 12h predugo za cert problem I nisam mogao doć na sebe kakav su problem napravili. Ono što je još gore, ovo im nije prvi put za taj isti slučaj. Dakle kako ne naučiš nešto iz ovoga. Dodano, date hack rješava problem da se ne mora drugi put čekati 12h
  • Bertone: A jučer nisam koristio MBP nego T14s pa nisam iskusio veselje sa certifikatima :D
  • Bertone: Imam MXMastera već godinama [link] i jedina dobra stavr kod njega je što ga mogu koristiti sa MBP, Tnink Padom i desktop računalom bez pretjeranog petljanja, no da je (pre)skup za ono što nudi, je.
  • smayoo: I, na kraju, ako za neki komad hardvera koji sam u životu kupio mogu reći da NIJE preskupo plaćen, to je sigurno Logitech MX Master 3. Jučer ujutro jesam bio spreman popi... na njih, ali činjenica je da su problem riješili unutar cca 12h otkad se pojavio.
  • smayoo: Osim toga, hardver radi s bazičnim postavkama, ali svi tweakovi i custom podešenja - ne.
  • smayoo: Radi se o softveru koji mora dobiti dozvolu korisnika da zakopa dosta duboko pod haubu radi komunikacije s periferijom.
  • smayoo: Ne mora biti spojeno na internet. Sigurnosni certifikat aplikacije je bio istekao i Apple ga je blokirao, a Logitech nije pazio na datum isteka certifikata i nije izdao update s novim certifikatom na vrijeme
  • Bertone: Ček malo, preskupo plaćeni hardver ne dela na Apple računalu ako računalo nije spojeno na internet (na stranu sad istekli certifikat) ili sam propustio nešto

Za vikanje moraš biti prijavljen.

Prijava

Novo na Jabučnjaku

Teme

Poruke

Komentari

Anketa

Koji Mac koristite?

Page Speed 0.78 Seconds

Provided by iJoomla SEO