Dotaz: Python vytvoreni matice

[guus]

Příloha:

• r.jpg (476854 bytů)

ahoj mam otazku ohlende vytvoreni matice (seznamu) v pythonu. mam toto zadani viz priloha

prvni cast jsem zvladl naprogramovat-kdyby nekdo potreboval napsal jsem program na vypocet soustavy rovnic Gauss-Jordanovou eliminaci

Nevim si ale rady s tou druhou casti...Tedy jak vytvorit tu matici. Nejake navrhy?

Odpovědi

[guus]

Příloha:

• ori.py (3145 bytů)

tady je jeste ten skript na reseni soustavy Gauss-Jordan eliminaci

[Jan Šimák]

Pokud není vyžadována vlastní implementace matice, tak je tu modul Numpy, který operace s maticemi řeší http://docs.scipy.org/doc/numpy/reference/routines.matlib.html

[guus]

dekuji za odpoved

nejsem jsi zcela jist zda modul numpy vyresi muj druhy ukol...modul ktery uvadis je spise jiz na pocitani s maticemi... ja nejprve potrebuji matici dle zadaneho schematu vytvorit...

ale dekuji jeste nejake dalsi navrhy?

PS prosim kouknete se do zadani aby bylo zrejme o co mi jde

[Jan Šimák]

ja nejprve potrebuji matici dle zadaneho schematu vytvorit...

A k tomu je právě určena třída http://docs.scipy.org/doc/numpy/reference/generated/numpy.matrix.html#numpy.matrix.
Na zadání jsem se díval, proto jsem poslal odkaz na numpi, který umí matici reprezentovat a provádět nad ní matematické operace.

[guus]

asi si nerozumime...nebo ja nerozumim spravne zadani...potrbuji udelat proceduru ktera vztvori matici na zaklade porovnani indexu prvku matice...tedz napriklad jak to chapu ja udelam nejaky cyklus a ten v okamziku kdy bude v bode 0 ve schematu odporu zapise do matice 2, protoze je v rohu prom kdyz prechazim z bodu 0 do bodu 1 pres odpor do matice se zapise -1 ...ted jsem v bode 1 nic se tu nedeje (v zadani klasifikovano jako "jinak") tedy se zapise 0...potom jdu z 1 do 2 pres odpor ted -1 jsem v bode 2 tam se nic nedeje takze se zapise 0...atd...potom se uy nic nedeje takze budou samy 0. Tim jsem dostal prvni radek matice tedy: 2 -1 0 -1 0 0 0 0 ' 0

Matici je mozne take dostat resenim kirchoffovych rovnic ktere jsou taky naznaceny v zadani pro N=5..Bude-li odpor R=0 dostanu primo napeti U...kdyz rovnice upravim dostanu take primo hodnoty matice kterou potrebuji...

Resenim teto matice pomoci G-J eliminacni metody ktrerou jsem naprogramoval dostanu vysledny odpor

Opet dekuji za odkaz. Ale myslim si ze to neni zrovna to co potrebuji. Modul munpy.matrix umoznuje uz primo upravovat vygenerovanou matici. Ale to ja nepotrebuji .Ja potrebuji napsat proceduru ktera vygeneruje jenom tu matici a tu potom poslu na G-J eliminici. Navic by bylo dobre kdyby to bylo udelano pokud mozno jako pure python bez modulu...

Nevim porad jestli si rozumime...

dekuji

[rADOn]

Tim jsem dostal prvni radek matice tedy: 2 -1 0 -1 0 0 0 0 ' 0

Cili mas pole hodnot, a az budes mit radky vsechny tak mas pole poli a to je (gasp) matice. Numpy matice je "jen" optimalnejsi na pamet a ma nejaky dalsi uzitecny funkce.

[guus]

no super...ale ja nevim jak se prave dostat k tem radkum ..k tem polim...nevim jak naprogramovat ten algoritmus pro vytvoreni pole poli tedy matice

[guus]

resp vim jak to funguje teoreticky vim jak se ta matice ma sestavit ale nevim jak to naprogramovat...proste bych potreboval pomoci nebo nakopnout s druhou casti ukolu v zadani...

[guus]

jo a nechapu proc je diskuze oznacena jako vyresena ...ona neni vyresena!!!

[rADOn]

V tom pripade jsi tady spatne. Tohle neni web "napiste za me domaci ukol". Jestli mas nejakej konkretni problem tak ho popis a ukaz kod.

[guus]

cast kodu resp prvni ukol mam vyresen...napsal jsem gauss-jordanovu eliminaci pro reseni soustavy rovnic.,,kod je v priloze ori.py...zasekl jsem se na reseni druheho prikladu...nikde nerikam ze chci vyresit domaci ukol...navic to neni ani domaci ukol ale pouze cviceni pro me...chtel jsem jen slysel nejake navrhy reseni jak matici sestavit...treba bomoci for cyklu nebo nejake podminky...protoze me nic nenapada jak to napsat obratil jsem se sem ... no a prominte ...reseni ve forme odkazu na modul .numpy ktery stejne neresi muj problem ...no..

[rADOn]

Mea culpa, ted jsem znova kouknul na to zadani a vidim ze jsem nejak nepochopil ze ten druhej bod neni soucasti prvniho. Omlouvam se, vyresni jsem zrusil, to druhy zadani nejak nechapu.

Abys na tom nebyl uplne skodnej, tak 2d pole muzes inicializovat snaze:

jednotka = [[0] * velikost for i in range(velikost)]

Doporucuju podivat se na list comprehension, imo by se ten kod dal zkratit tak na petinu

[guus]

nihil factum

no to je prave to co jsem psal...ze si nerozumime co chci...ano je pravdou ze se da kod zkratit ale jak jiz jsem napsal chtel jsem program napsal jak pure (cisty) python tedy bez modulu.

co se tyka inicializace 2D pole i to jsem zvladl...muj problem tkvi v tom ze nevim jak naprogramovat to aby skript prirazoval patricna cisla do matice...tedy aby udelal to ze kdyz napr dob 0 bude v rohu dal do matice 2 prechod z 0 do 1 zapsal -1 atd...proste aby byla splnena podnimka ktera je udana v zadani pro matici Aij...a to pro zadani N... tedy neni mi jasne jak rozhodnout o tom zda je bod (odpor) na rohu na strane uvnitr atd...

dekuji za cas

[guus]

jinak prvni a druhy bod zadani mezi sebou samozrejme souvisi...tim ze jsem napsal prvni bod tedy ze umim resit soustavu rovnic pomoci G-J eliminace mi pomuze s tim ze ve druhem bode sestavim matici proudu a napeti a tuto matici poslu na G-J eliminaci a tim dostanu vysledny odpor...

ja nevim jak to lepe vyjadrit...vse je v zadani

[guus]

a cemu nerozumis na zadani toho druheho prikladu?

[rADOn]

podivas se jestli jsi na krajnich souradnicich?

def sousedi(x, y, xsize, ysize):
  return 4 - (x == 0) + (x == xsize) + (y == 0) + (y == ysize)

V podstate to mas v tom zadani napsany…

[guus]

Příloha:

• jak_na_matici.pdf (83051 bytů)

pockej kouknu na to ,jeste jsem sesmolil jak se dostane ta matice tak na to juknete v priloze

[guus]

Příloha:

• jak_na_matici.pdf (83051 bytů)

pockej kouknu na to jeste jsem sesmolil jak se dostane ta matice tak na to juknete v priloze

[guus]

promin tomu moc nerozumim ..predpokladam ze jsi napsal funkci ktera vypise 4 takze dle zadani ze bude bod uprostred. x y jsou indexy matice Aij resp. N^2-1=y a N^2=x? _nevim co je xsize ysize

[rADOn]

Pro matici velikosti xsize*ysize to vrati kolik ma bod x,y sousedu. Pomoci toho triku s implitnim pretypovanim boolu na int jak se popisuje v zadani: pokud je x na kraji tak x == 0, pokud je na opacnym kraji tak x == xsize etc. scitani vydrazdi pretypovani na cislo, int(True) je 1 cili to staci poscitat a mam pocet okraju.

[guus]

aha

no dobre...takze dejme tomu ze budu mit matici 8x9 ktera by mela odpovidat schematu pro N=3

zadam-li do tve funkce: def sousedi(x, y, xsize, ysize): return 4-(x == 0) + (x == xsize) + (y == 0) + (y == ysize)

print sousedi(0,0,8,9)

dostanu 4...pritom podle me bych mel dostat ze uzel o souradnicich 0,0 ma sousedy 2 dle schematu pro N=3

promin ja tomu proste nerozumim...

a co se tyka te konverze tak jsem uplne namydlenej

[rADOn]

To bude asi mit nejakou souvislost s tim ze tam byla chyba. Dve kdyz prictu spatnej koncovej index Nemuzes od toho cekam moc, pisu to z hlavy…

def sousedi(x, y, xsize, ysize): return 4 - int((x == 0) + (x == xsize-1) + (y == 0) + (y == ysize-1))

Vysvetleni te konverze je proste mily watsone, zkus si v konzoli zavolat int(True)

[guus]

jo aha...no je mi to jasny tak ze 70 procent...musim to poradne promyslet...co jsem teda pochopil tak ten mensi skriptik neudela nic jinyho nez ze pujde bod po bodu za predpokladu ze vhodne zacyklim x,y a vypise pocet sousedu...ok

ale jak tam ted zahrnu tu -1 tedy prechod z jednoho bodu do druheho a 0 "jinak"?

Priznam se ze jsem z toho trochu tumpachovej... Sedim nad tim uz tri tydny a na tuhle cast ne a ne prijit (asi jsem se uz taky zacyklil)...dekuji za cas mily Holmesi

[guus]

a mohl bych Te poprosit? nechtel by si mi prosim ten skript na generovani matice napsat? Ja vim je to drzost...