Elektroda.pl
Elektroda.pl
X
IGE-XAO
Proszę, dodaj wyjątek dla www.elektroda.pl do Adblock.
Dzięki temu, że oglądasz reklamy, wspierasz portal i użytkowników.

opór elektryczny - zadanie

14 Mar 2005 19:37 4388 28
  • IGE-XAO
  • Pomocny post
    Poziom 29  
    P2P3 = 1/n oraz P1P2 = 1/m
    Rab² = (1/n)² + (1/m)² (wzór kwadratów trójkąta)

    jeżeli n=m to

    Rab² = 2/n² => Rab = √2/n² => Rab = √2/n

    czyli:
    a) 0,14
    b) 0,014
    c) 0,0028
    d) 0

    Pozdro.
  • Poziom 10  
    Dzięki za pomoc. Wykorzystanie po prostu twierdzenia Pitagorasa do rozwiazania tego zadanka wydawalo mi sie za proste... Zajrzałem do zeszytu i tam zobaczyłem, ze niepotrzebnie dorysowałem odcinek łączący punkty A i B. To znacznie uprościło zadanko. Mea culpa. Link u góry nie będzie działał juz. Błąd poprawiłem i własciwe zadanko teraz umieściłem pod:

    http://ciapek.uci.agh.edu.pl/~rafalz/zadanie_z_TO_ok.doc

    Teraz wszystko jest ok. Nie wydaje mi się, żeby istniało jakieś proste rozwiązanie tego problemu, skoro prowadzący ćwiczenia z teorii obwodów dał na niego rozwiązanie aż 3 miesiące.

    Jeżeli ktoś zna rozwiązanie tego problemu lub choć ma pomysł jak to by można ugryźć [np. jak i w czym napisac program], proszę o kontakt [namiary w pliku].
    Pozdrawiam. I mimo mojej pomyłki thx, john_t!!
  • IGE-XAO
  • Pomocny post
    Poziom 12  
    Witam
    Program nazywa sie MatLab, a nie MathLab.
    Zadanie masz juz rozwiazane, a jak tylko zdobedziesz MatLaba to napisanie skryptu zajmie Ci 1 minute.
    , ja tu widze z 10 linii kodu.
    Pozdrawiam.
  • Pomocny post
    Poziom 29  
    Jeśli dobrze rozumiem to tak jakby ustawić 4 rezystory połączone w rogach kwadratu i zmierzyć rezystancję po przekątnej.
    Jeśli tak to rozwiązanie jest jeszcze prostsze.

    1/Rab = 1/(P1P2+P2P3) + 1/(P3P4+P4P1)
    Jeżeli P1P2=P2P3=P3P4=P4P1=R to
    1/Rab = 1/(2R) + 1/(2R)
    1/Rab = 1/R
    Rab = R
    W odniesieniu do n => Rab(n) = 1/n
    i tak:

    a) 0,1
    b) 0,01
    c) 0,002
    d) Dla n (dążącego do nieskończoności) Rab(n) = 1/n = 0
  • Poziom 10  
    dobap napisał:
    Witam
    Program nazywa się MatLab, a nie MathLab.
    Zadanie masz juz rozwiazane, a jak tylko zdobedziesz MatLaba to napisanie skryptu zajmie Ci 1 minute.
    , ja tu widze z 10 linii kodu.
    Pozdrawiam.


    Hej! Właśnie idąc za twoja radą zacząłem ściagac triala MatLaba. Plik jest spory także licze, że 2 dni mi wystarczy na jego zassanie. Dzieki za nadzieje, ze rozwiazanie tego zadania jest osiagalne. :o

    Chcialbym sie jeszcze dowiedziec, co rozumiesz poprzez "Zadanie masz juz rozwiazane". Chodzi o to, że w przykładach, demo lub helpie w MatLabie masz juz to zadanie rozwiązane?

    Nie znam wogóle progsa ani jego jezyka, takze watpie zeby mi to zajelo minutke. Najpierw pewnie bede musial przebrnac przez tutoriale itp zeby sie do tego zabrac wogole. Jesli mowisz ze napisanie skryptu zajmuje ci jedna minutke, to bylbys tak uprzejmy napisać chociaż jego zarys? Bardzo to skróciłoby mi czas na opracowanie problemu. A dla ciebie to tylko chwilka!

    Pozdrawiam!

    Dodano po 40 [minuty]:

    john_t napisał:
    Jeśli dobrze rozumiem to tak jakby ustawić 4 rezystory połączone w rogach kwadratu i zmierzyć rezystancję po przekątnej.
    Jeśli tak to rozwiązanie jest jeszcze prostsze.

    1/Rab = 1/(P1P2+P2P3) + 1/(P3P4+P4P1)
    Jeżeli P1P2=P2P3=P3P4=P4P1=R to
    1/Rab = 1/(2R) + 1/(2R)
    1/Rab = 1/R
    Rab = R
    W odniesieniu do n => Rab(n) = 1/n
    i tak:

    a) 0,1
    b) 0,01
    c) 0,002
    d) Dla n (dążącego do nieskończoności) Rab(n) = 1/n = 0


    Hmm, jutro mam ćwiczenia z Teorii Obwodów, to pokaże prowadzacemu Twoje rozwiazanie. Tyle ze obawiam sie ze skoro rozwiązanie jest faktycznie tak proste, to jest ono złe. Prowadzacy wspominal, ze dla n=m=10 da sie to jeszcze na papierze policzyc, ale dla 100 i 500 trzeba ponoc juz napisac program - on nie zna szybkiej metody "ręcznej".
    Dam znać czy twoja metoda okaże się dobra.
    Dzieki i pozdrawiam! :spoko:
  • Pomocny post
    Poziom 12  
    Witam

    Wybacz, ze wyrazilem sie nie dosc jasno, mowiac, ze masz rozwiazane, mialem na mysli rozwiazania kolegi John_t.
    Z tego co widze to myslisz, ze jest nie poprawne, ja nie mam zdania.
    Na szczescie pozeganlem sie z Teoria Obwodow juz ponad rok temu i tez nie obylo sie bez krwi:).
    Przyjze sie temu zadaniu, ale jesli wchodza w gre te wszystkie magiczne, genialne, maciezowe algorytmyto boje sie, ze moge nie dojsc do porozumienia z tym zadaniem-sprobuje.
    Co do MatLaba sluze pomoca.

    P.S.
    Jesli rozwiazanie jest iteracyjne....to masz racje dla 5 czy 10 krokow zrobic to recznie jest niezlym wysilkem....a co tu mowic o 100 krokach.....

    Jak juz bedziesz mial MatLaba to mozesz zajrzec na www.mathworks.com/
    Znajdziesz tam kupe gotowych skryptow.

    Pozdrawiam.
  • Poziom 10  
    dobap napisał:
    Witam

    Wybacz, ze wyrazilem się nie dosc jasno, mowiac, ze masz rozwiazane, mialem na mysli rozwiazania kolegi John_t.
    Z tego co widze to myslisz, ze jest nie poprawne, ja nie mam zdania.
    Na szczescie pozeganlem się z Teoria Obwodow juz ponad rok temu i tez nie obylo się bez krwi:).
    Przyjze się temu zadaniu, ale jesli wchodza w gre te wszystkie magiczne, genialne, maciezowe algorytmyto boje się, ze moge nie dojsc do porozumienia z tym zadaniem-sprobuje.
    Co do MatLaba sluze pomoca.

    P.S.
    Jesli rozwiazanie jest iteracyjne....to masz racje dla 5 czy 10 krokow zrobic to recznie jest niezlym wysilkem....a co tu mowic o 100 krokach.....

    Pozdrawiam.


    Dzięki. MatLaba tylko dwa razy używałem bo do tej pory mialem tylko dwie laborki z Metod Numerycznych. Także tylko znam najogolniejsze podstawy. Jakby Ci sie udało zrobić to zadanie, skakałbym z radości!Zależy mi na tym, bo po paru zajeciac z Teorii Obwodow tez juz czuje zapach krwi w powietrzu :] Takze jesli to zadanie gwarantuje polowe minimum zaliczenia to genialnie i sprawa jest warta swieczki. W Tobie moja nadzieja! :cunning:
    Pozdrawiam! :spoko:
  • Pomocny post
    Poziom 29  
    _druss_ chyba w końcu wiem o co tu chodzi.
    To co przedstawiłem powyżej to stosunek rezystancji pojedynczej komórki do układu liniowego. A tu chodzi o to że jest to sieć, kwadratowa macierz rezystorów o wymiarze (n x m) i należy wyznaczyc rezystancię całej sieci w punktach AB. I tu rzeczywiście rozwiązanie nie jest takie proste !!!
    Ale pomyślę nad tym.
  • Poziom 10  
    john_t napisał:
    _druss_ chyba w końcu wiem o co tu chodzi.
    To co przedstawiłem powyżej to stosunek rezystancji pojedynczej komórki do układu liniowego. A tu chodzi o to że jest to sieć, kwadratowa macierz rezystorów o wymiarze (n x m) i należy wyznaczyc rezystancię całej sieci w punktach AB. I tu rzeczywiście rozwiązanie nie jest takie proste !!!
    Ale pomyślę nad tym.


    No właśnie, skoro prowadzący dał na to 2 miesiące, to to nie moze byc proste. Zreszta biorac to na logike, kto byl dal polowe minimum zaliczenia za cos prostego.
    Dziękuje! Mam nadzieję że uda sie cos wykombinowac. Ja juz sciagam MatLaba, liczac ze ten program pomoze rozwiazac problem.
  • Pomocny post
    Pomocny dla użytkowników
    Wrzuciłem układ na symulator (dzięki czemu mam dość wprowadzania rezystorów na 3 miesiące i 22 dni) i wyszło:

    Dla ułatwienia wszystkie rezystory = 1 kΩ

    Układ 4x4 węzły -> Rab=714.3 Ω
    Układ 6x6 węzłów -> Rab=668Ω
    Układ 8x8 węzłów -> Rab = 654.2Ω
    Układ 10x10 węzłów -> Rab=648.6Ω
  • Poziom 10  
    Paweł Sujko napisał:
    Wrzuciłem układ na symulator (dzięki czemu mam dość wprowadzania rezystorów na 3 miesiące i 22 dni) i wyszło:

    Dla ułatwienia wszystkie rezystory = 1 k?

    Układ 4x4 węzły -> Rab=714.3 ?
    Układ 6x6 węzłów -> Rab=668?
    Układ 8x8 węzłów -> Rab = 654.2?
    Układ 10x10 węzłów -> Rab=648.6?


    Dzięki! Czyli już mam wynik do podpunktu a) mojego zadania, a mianowicie Rab=648.6?. Jednak opisywana przez Ciebie metodą nie da się zrobić zadania dla n=m=100 a tym bardziej 500! Chyba że ktos ma cierpliwosc wprowadzać 499? oporników do kompa jesli dobrze licze. Poza tym prowadzacy żąda wydruku procedury, ktora zostala uzyta do obliczeń. Ale postęp już jest! Może wspólnymi siłami sie coś nam wszyskiem uda wykombinować.

    A tak na marginesie, jakiego symulatora użyłes do tych obliczeń?
    Czy to zadanie da się zrobić w Workbenchu?

    Pozdrawiam.
  • Pomocny post
    Pomocny dla użytkowników
    W Circuit Makerze

    W układzie 4x4 udało mi się wyprowadzić, że Rab=(5/7)*Rgałęzi

    Układ ma symetrię po przekątnej przechodzącej przez punkty A i B

    Na drugiej głównej przekątnej wartości napięć wynoszą 0 dla wszystkich przebadanych układów od 4x4 do 10x10 więc pewnie dla wyższych będzie to samo.

    Układ 4x4 dał się rozciąć na dwie równoległe elektrycznie części idąc po przekątnej AB.

    Przez narożniki na przekątnej AB nie płynie prąd - można je wyciąć.

    Niektóre punkty dają się połaczyć, bo mają ten sam potencjał (niestety Circ Maker) nie skopiował obliczonych napięć węzłowych.

    Przy badaniu trzeba zasilać układ symetrycznie względem masy (wyprowadzonej ze środka baterii) wtedy w napięciach węzłowych widać symetrię układu.

    A swoją drogą to należałoby psami poszczuć wymyślaczy takich zadań, nie mogą sobie sami wziąć omomierza, lutownicy i pastwić się nad sobą ? ;)
  • Pomocny post
    Pomocny dla użytkowników
    W Circuit Makerze

    W układzie 4x4 udało mi się wyprowadzić, że Rab=(5/7)*Rgałęzi

    Układ ma symetrię po przekątnej przechodzącej przez punkty A i B

    Na drugiej głównej przekątnej wartości napięć wynoszą 0 dla wszystkich przebadanych układów od 4x4 do 10x10 więc pewnie dla wyższych będzie to samo.

    Układ 4x4 dał się rozciąć na dwie równoległe elektrycznie części idąc po przekątnej AB.

    Przez narożniki na przekątnej AB nie płynie prąd - można je wyciąć.

    Niektóre punkty dają się połaczyć, bo mają ten sam potencjał (niestety Circ Maker) nie skopiował obliczonych napięć węzłowych.

    Przy badaniu trzeba zasilać układ symetrycznie względem masy (wyprowadzonej ze środka baterii) wtedy w napięciach węzłowych widać symetrię układu.
  • Poziom 10  
    <wow> dzięki Paweł Sujko! Jako że dopiero sie wkrecam w teorie obwodow, zajmie mi troche czasu zanim to rozkminie co napisales. Ale wszystko jest przejrzyscie bardzo zrobione przez Ciebie, także na pewno to ulatwi sprawe.

    Opierając sie na tych wszystkich informacjach przez Ciebie zaprezentowanych, postaram się - mam nadzieje z pomoca uzytkownikow forum - rozgryzc algorytm rozwiazywania tego typu zagadnien.

    Próbowałem dzis z pomoca OriginLaba znalezc zaleznosc przedstawionych przez Ciebie danych eksperymentalnych od ilości wezłów n [n=m], ale z marnym skutkiem. MatLab mi sie konczy sciagac, moze z jego pomoca cos wykmini sie.

    Dzieki jeszcze raz! Pozdro!
  • Poziom 12  
    witam

    Pamietaj, ze MatLab to jest jezyk programowania, z przyblizeniem powiedzmy jak C++, on nie znajdzie zaleznosci, on wykona skrypt, czyli rownania, algorytmy, ktore kazesz TY mu wykonac.
    Z tulbox-ami rozwiazuje algorytmy, matematyczne zalezosci....ale on nic nie znajdzie, nie zrobi nic poza tym co mu rozkazales.

    Pozdrawiam i powodzenia.
  • Poziom 10  
    Hej dobap!
    Wiadomo, że musze mu to zadac w jego jezyku. Mówiłes, ze ty umiesz poslugiwac sie MatLabem jako tako. Mniej wiecej masz pomysl jak ten skrypt powinien wygladac/jakiego toolboxa wykorzystać? Bylbym wdzieczny gdybys przyblizyl mi jak mam to ugryzc. Bez pomocy z zewnatrz napisanie tego skryptu zajmie mi pewnie kupe czasu takze prosze o pomoc, wskazówki itp.

    MathLaba juz skonczylem sciagac, teraz koncze instalacje, jutro sie zabieram do rozkminy.

    Pozdro!
  • Pomocny post
    Pomocny dla użytkowników
    Dla układu 6x6 rezystancja zastępcza wynosi:

    Rab=(331/495)*Rgałęzi

    nawet się ładnie ten układ redukuje, tylko to dodawanie ułamków zwykłych ... ;)

    Do obliczeń komputerowych trzeba użyć metody potencjałów węzłowych.
    Buduje się N równań prądowych z N niewiadomymi potencjałami węzłowymi a potem rozwiązuje się tak powstały układ równań.

    Jeden węzeł musi być węzłem masy - proponuję wybrać węzeł A. Jeżeli węzeł B zasilisz źródłem prądowym o wydajności jeden to po rozwiązaniu n równań z n niewiadomymi potencjałami węzłowymi otrzymasz rezystancję zastępczą sieci wyrażoną potencjałem w punkcie B

    Polecam książkę:

    Białko Michał:
    "Analiza układów elektronicznych wspomagana mikrokomputerem"
  • Poziom 10  
    Dzięki. Teraz wlasciwie to po tym co mi Paweł napisałes to teraz tylko zostaje zaprogramować to w matlabie. A sadze juz ze wprowadzenie petli liczacej 500 rownan juz bedzie do wykminienia.
    Sprobuje zdobyc ta ksiazke ktora polecasz w uczelnianej bibliotece.
    Dzieki jesczez raz za zaangazowanie!
  • Pomocny post
    Pomocny dla użytkowników
    Proponuję zautomatyzować budowanie macierzy admitancyjnej, bo ręczne wprowadzanie może być "letko" upierdliwe. Można też ograniczyć ilość węzłów korzystając z symetrii układu.
  • Poziom 10  
    Spoko. Wyczaje tylko co to jest macierz admitacyjna :] Dopiero wchodze w powazniejsza teorie obwodow, takze nie mam w tym wszystkim bieglosci. Mam nadzieję że z pomocą tej ksiazki co mi poleciles uda sie cos wykombinowac. Sprawdzalem - sa w bibliotece 3 sztuki takze nie bedzie problemow z wypozyczeniem. Pozdro!
  • Poziom 10  
    Hej, Paweł Sujko! Sorki ze cie tak męczę ale znalazlbys chwile zeby napisac jak doszedles do tego wyniku 5/7 dla 4x4 i 331/445 dla 6x6? Dalbys rade walnac skany ze swoimi obliczeniami recznymi czy cos w tym stylu? Pozyczylem ta ksiazke ktora polecales i faktycznie w 4. rozdziale mniej wiecej jest napisane jak tworzyc macierze do zadan z rezystorami. Tyle ze nigdzie nie moge znalezc opisu tej metody potencjalow wezlowych o ktorej pisales. Takze bylbym bardzi wdzieczny jakbys mi pokazal na tych przykladach 4x4 i 6x6 jak przebiegaja obliczenia, to wtedy szybko sczailbym algorytm jaki trzeba zastosowac.

    Sorry, ze moze zadaje infantylne pytania ale jesli chodzi o te sprawy to jestem żółtodziobem.

    A propos tych wcześniejszych obliczeń oporow wykonanych przez Ciebie czyli 714,3; 668; 654,2; i 648,6, to wyglada na to ze sa ok bo dla n dazacego do nieskonczonosci Rab=2/pi - kumpel znalazl analityczne rozwiazanie tego problemu, tyle ze go nie rozumie, bo zastosowana jest tam jakas hardcorowa matma wyższa. A dodam ze facet jest najwiekszym mozgiem w grupie. Takze mysle ze istnieje jakas wykladnicza badz ekspotencjalna zaleznosc miedzy Rab a n [lub n^2]. Takze majac 4-5 punktow sprobuje dopasowac mniej wiecej jakas funkcje dla tej zaleznosci. Moze sie ta moja hipoteza potwierdzi, zwlaszcza, ze funkcja powinna miec asymptote pozioma o rownaniu y=2/pi a zbior wartosci to (1;2/pi). Ale generalnie nawet jesli znajde taka funkcje aproksymującą, to bedzie ona przydatna tylko do sprawdzenia poprawnosci wynikow... :/

    Pozdrawiam i jesczez raz prosze o wrzucenie na forum sposobow obliczeniowych jakie ty wykorzystales do ukladow 4x4 i 6x6.
    Pozdrawiam!
  • Pomocny post
    Pomocny dla użytkowników
    Nie wiem czy uda się uogólnić metodę przeze mnie użytą na wyższe n.

    Istota obliczeń jest zawarta w zamieszczonych już schematach.

    1. Ponieważ zasilasz układ po przekątnej A_B to widać, że cały układ jest symetryczny względem tej osi - składa się z dwóch trójkątów rezystorowych górnego i dolnego połączonych równolegle na osi A-B.
    Tak więc wycinamy elementy powyżej osi A-B pozostawiając oczywiście połączenia rezystorów z dolnego trójkąta leżące na osi.

    Tak uzyskany układ ma rezystancję dwa razy większą niż cały układ.

    2. Na osi prostopadłej do A-B (druga główna przekątna) wszystkie punkty mają ten sam potencjał więc możemy je zewrzeć.


    Kolejne przekształcenia układu 4x4 pokazane są na kolejnych rysunkach

    Uc = Ud = Ue = Uf -> wynika z symetrii układu i jego zasilania.

    Ug=Uh -> pozwala wyciąć rezystory narożne
    Ui=Uj -> jak wyżej

    Po wycięciu góry mamy Rab'=2*Rab (rys. krata4x4_3)

    Na ostatnim po zwarciu punktów o tym samym potencjale mamy dwa identyczne układy szeregowe.

    Każdy podukład składa się z rezystora R i dołączonego do niego rezystora (5/2)*R

    $$Rab=\frac{\frac{5}{2}*R*R}{\frac{5}{2}*R+\frac{2}{2}R}=\frac{5}{7}*R$$


    Dla ułatwienia rezystancja podukładu wynosi Rab całego układu (ot taka ciekawostka wskazująca na marnację rezystorów w reszcie układu :)


    W układach wyższych rzędów postępuje się tak samo ale jest to upierdliwsze, bo trzeba stosować przekształcenie gwiazda-trójkąt by zwijać kolejne oczka.

    Poniższe 2 układy są równoważne elektrycznie.

    Z gwiazdy na trójkąt:

    R5=R1+R2+(R1*R2)/R3

    R6=R2+R3+(R2*R3)/R1

    R4=R1+R3+(R1*R3)/R2
  • Pomocny post
    Poziom 12  
    Czolem

    Sprawa nie jest tak hop-siup. Zaimplementowac wezlowke nie powinno byc tak najgorzej - stosunkowa prosta metoda - ale jak juz zapewne wiesz w MatLabie chodzi o to zebys nie musial pisac - jak C++ - jak sie liczy macierze itp itd, jak chcesz szybka transformate Furiera to wpisujesz "fft()" i Ci programik wszystko sam zrobi, a jak go poprosisz to jeszcze narysuje.
    Nie wiem jakiego toolboxa bedziesz potrzebowal, ale z tego co pamietam o wezlowce to raczej nie powinno byc zadnego klopotu z tym, prawdopodobnie wystarczy podstawowy MatLab(zawsze siec p2p jest otwarta:)).
    Ja na Twoim miejscu bym jeszcze zahaczyl o PSpica - tak profilaktycznie sprawdzic sobie wyniki:)

    Pozdrawiam.
  • Poziom 10  
    Dziekówa. MatLaba juz mam i sporo info od Was. Postaram się coś wykombinować. Mam troche nauki biezacej na weekend takze jak nie zdaze to w prerwie swiatecznej cos poprobuje!
    Pozdrawiam!
  • Pomocny post
    Poziom 17  
    wyniki dla n=6,8,10 sprawdziłem z programem symulacyjnym i jest ok
    tylko dla n=4 liczy błędnie ale to jest na żywca poprawione na sam koniec.
    korzystam z metody prądów oczkowych .w załączniku jest rysunek
    musze przyznac ze to dosc trudne zadanie . dla n=100 jeszcze mój program policzy , ale dość długo
    z pewnością istnieje prostsze rozwiązanie
  • Poziom 10  
    Dzięki za rozwiązanie. Wczytałem do Matlaba twój program i chyba jest ok tylko ja jestem początkujący jeśli chodzi o ten program [jeśli chodzi o teorię obwodow tez :)].

    Mam pare pytanek co do niego.
    Jak określa sie ilość punktów w boku kwadratu? Po wczytaniu tego programu wpisuje np. "n=100", nastepnie "R_Z" a program wskazuje cały czas tę sama liczbe: 0,6478 [sprawdzałem tez wpisując 'format long' - to nie błąd przybliżenia]. Pewnie moje wątpliwości wynikają z braku pracy z programem więc byłbym wdzięczny, gdybyś mi to pokrótce wyjaśnił.

    Domyślam się jedynie, że R_Z to opór zastępczy całej kwadratowej kratownicy, n - dlugość boku kwadratu [ilość punktów], ans - wynik programu; nie wiem natomiast co oznaczyłes przez wypiszg i wypiszd oraz jakie napiecie i natęzenie poprzez U oraz I.

    Byłbym wdzięczny, gdybyś spróbował mniej więcej wyjaśnić te moje wątpliwości. Byłoby super gdybys streścił na czym polega stosowana przez Ciebie metoda prądów oczkowych. Jestem tutaj nowy także prosze o wyrozumiałość.

    Pozdrawiam!!
  • Pomocny post
    Poziom 17  
    powinieneś n zmieniać w programie (w pliku) . zasejwować i włączyć .nie dziwię się ,że Ci cały czas ten sam wynik wywala.
    te wypiszd i wypiszg służyły mi tylko jako indeksy pomocnicze przy wypełnianiu tabeli (macierzy) oporności R.

    otwórz sobie rysunek - tu jest 8 na 8 punktów , zapuść program wpiusąc n=8
    w oknie matlaba wpisz wypiszd i wypiszg - poatrz na rysunek - teraz już wiesz co oznaczają te wektory

    w oknie matlaba wpisz R - zobaczysz jak wygląda macierz oporności

    wektor I - to wektor prądów oczkowych (na rysunku zaznaczone numerami od 1..13) . cała zabawa polega na tym żeby znaleźć I(r) r=13 na moim rysunku.
    układ jest symetryczny to wiadomo ,że drugi taki prąd będzie pobierany z naszego wyimaginowanego źródła napięciowego zasilającego całość - stąd linijka nr 54 wygląda tak : R_Z=1/I(r)/2;

    układ tak na dobrą sprawę jest powójnie symetryczny - osie symetrii - w kształcie litery X . popatrz teraz na rysunek - pod 12, a obok 10 też płynie jakiś prąd oczkowy . ze względu na symetrię ma on taką samą wartość jak 9 !!! wpisz teraz w oknie matlaba R - w niektórych miejscach masz dwójki - to dlatego ,że tam "siedzą" prądy - I9 oraz jego bliźniak .

    tak na dobrą sprawę powinieneś sam spróbować na kartce wypełnić jakąś macierz R np dla przypadku n=6 i sprawdzić z progamem

    wektor U - wektor napięć oczkowych
  • Poziom 10  
    Dzieki. Ten algorytm faktycznie liczy tylko mniej wiecej dla n do 150, pozniej juz komp wywala 'out of memory'. Korzystajac z fachowej pomocy zoptyamalizowalem twoj skrypt z uzyciem macierzy rzadkiej i teraz dla n=500 liczy mniej wiecej w 10 min.
    Jeszcze raz dzieki, bez Ciebie bym nie ruszyl z miejsca.