• Welcome to TUKE FÓRUM - Fórum pre študentov Technickej Univerzity v Košiciach.
 

Simulacne systemy

Started by 7ici7, 03.12.2008, 02:50:17

« predchdzajce - alie »

Sxx

podla mna je C spravne takto.....

nyquist(sin(omega*t)*SYS);
bode(sin(omega*t)*SYS);

:)

inak 3. priklad, ked mas schemu a robis maticu q=[.....]  a nemas rovnaky pocet vstupov pri vsetkych blokoch, tak chybajuce vstupy doplnis nulou?

a este jedna vec ma zaujima, jak sa kresli schema v simulinku zo zadanej dif. rovnice?


dabol

SUPER ;)

no, na prvu otazku ANO ked nema vstupy je nula. A na tu druhu netusim ;)
V spoločnosti chcem piť sám.

Sxx

a ten RLC obvod vies? respektive chapes? je zadane presne to iste co vyratane na webe len netreba ratat s cislami ale globalnymi premennymi


dabol

No tie globalne mas pouzite pri voziku, tj

function xdot=nabi(t,x);
global r,l,c,udc;
xdot=[...];

global r,l,c,udc;
r=input('zadaj hodnotu r');
.
.
udc=input('zadaj hodnotu udc');

a pokracujes ode45-kou
V spoločnosti chcem piť sám.

Sxx

nechce sa ti sem pastnut ten priklad jak to ma byt cele? :)


dabol

#30
nabi.m
-------------------
function xdot=nabi(t,x)
global R L C Udc;
xdot=[x(2)/c ; 1/L*(Udt-R*x(2)-x(1))];

------------------- pokracuje prikazovy riadok
X0=[0,0] % pociatocne podmineky mas zadane
global R L C Udc;
R=input('vloz R=');
L=input('vloz L=');
C=input('vloz C=');
Udt=input('vloz Udt=');
t0=input('vloz t0=');
tfin=input('vloz tfin=');
[t,x]=ode45('nabi',[t0 tfin],X0);
subplot(231), plot(t,x(:,1)), title('u v t');
subplot(232), plot(t,x(:,2)), title('i v t');
ic=x(:,2);
ur=R*ic;
subplot(233), plot(t,ur), title('ur v t');


Dajak taq ;) ale v matlabe som neskusal...
V spoločnosti chcem piť sám.

Sxx

mne v matlabe nejdu ani tie ukazkove :D furt dajaky error


dabol

V podstate ide o to aby si mal globalne premenne tie ktore pouzivas vo fkcii v tomto pripade nabi.m. Ostatne ako X0 uz nemusia byt globalne.
V spoločnosti chcem piť sám.

Sxx

a ten substitucny kanoicky tvar to je ktore


dabol

#34
No to je ta uprava ked pouzijes substituciu :

x1(t) = 1/C*x2(t)
x2(t) = 1/L*(Udt - R*x2(t) - x1(t))

toto potom prepisujes do xdot=[1/C*x(2), 1/L*(...)] prvy column pre x1(t) a druhy pre x2(t)
V spoločnosti chcem piť sám.

Sxx

no dik :), ja uz asi kladiem rano stanem a este mrknem potom


drobcek

Nemohol by niekto naskenovat prednaskyzo simulacnych a hodit to na web???dik... bp bp bp

dabol

Hmm... tak neviem, ale prednasky su snad tu nie ? ;)

http://matlab.fei.tuke.sk/ss/
V spoločnosti chcem piť sám.

drobcek

No tam su niketore prednasky len....

drobcek

Cize poprosim Vas ak sa tunajde daka dobra dus,ktora by naskenovala prednasky a hodila to tu na web.dakujem br

dabol

Otazka ku skuske...

Staci prist na termin, ktory je na stranke, alebo sa treba (da sa?) zapisat v IS? Lebo ja tam nemam nic len, ze mam udeleny zapocet. THNX
V spoločnosti chcem piť sám.

andy.gee

nazdar .. ked nemas v IS zapisany zapocet a mas ho spraveny tak sa musis ozvat jadlovskej alebo cviciacemu, ale skor asi jadlovskej,
lebo sa musis cez IS prihlasit na skusku nakolko tam su terminy vypisane.

dabol

No ja mam zapocet zapisany v IS, ze mam 28 bodov ale ked tuknem na terminy tak tam nic nemam ;(, totiz tento predmet opakujem takze to bude asi tym. Napisem teda Jadlovskej...
V spoločnosti chcem piť sám.

andy.gee

#43
no tak mame po prvom termine  :)  8.1.2009

najprv trosku pocity zo skusky .. mala pripravene 4 skupiny ale strasila, ze ich ma 8 takze sa nebudu opakovat tak neviem no ;) nam rozdala dve skupiny A a B, ja som mal A skupinu tak k tomu nieco napisem .. takze bola najprv pisomna cast (10 otazok), tu opravoval lonscak, ale aj ked opravoval dost mierne tak sa stalo, ze nieco ani neopravil. Aby vas pustila na ustnu cast tak ste museli mat tu NADPOLOVICNU vacsinu cize 31 bodov ale pustila aj tych co mali 29 a 30. Ostatnych poslala domov ale URCITE aj ked vas posle domov si popytajte a POZRITE PISOMKU ci mate dobre spocitane body a ci vam lonscak vsetko opravil.
Koho uz pusti na ustnu tak ten spravi urcite. Ustna cast vyzerala tak, ze si pozriete pisomku, a da vam este nejake tri otazky, ktore pisete na papier (konkretne napisem dolu) a potom s nou ich prejdete, tam mozte ziskat navyse nejake body (aj 6-9 bodov) a prejde s vami aj niektore zadanie (zvykla LDR, NDR (hlavne runge kuta ju zaujimala) a aj 4rte zadanie pozerala), tiez s vami prejde skusku tak zbezne iba.

Takze niektore otazky co si pamatam:
SKUPINA A1

1.Vypocitat urcity integral funkciou quad
integrovanu funkciu treba zadefinovaj aj pomocou m-fajlu aj pomocou niecoho ineho

2. nacitat zo suboru maticu hodnot
load merania.dat
zadefinovat x-ovu premennu
x=merania(:,1)
zadefinovat y-ovu premennu
y=meranie(:,2)

aproximovat vysledky x,y polynomom n-teho stupna (trebalo pouzit polyfit ale asi chcela este aj nieco ine sa mi vidi)
vykreslit do grafu povodne nacitane udaje
vykreslit aproximovany graf
porovnat chybu originalneho a aproximovaneho

3. diferencialnu rovnicu druheho radu riesit pomocou nejakeho symbolickeho toolboxu
(ja som pouzil dsolve)

4. pre zadanu funkciu f(x) vypisat tabulku 1000 hodnot funkcie v intervale <0,2pi>
vykresilt graf

5. Pre zadanu funkciu f(x,y) treba vypisat hodnoty v intervale <-pi,pi> so zadanym krokom
Funkciu treba interpolovat + vykreslit
nejaky prikaz meshgrid tu bol asi, ten prikaz generuje matice .. ale nevedel som co s tym dalej

6. Pomocou Simulinku riesit diferencialnu rovnicu
zadefinovat v MATLABE !! (prikaz setparameters) hodnoty premennych diferencialnej fcie
nakreslit simulinkovu schemu
tu nebola tazka diferencialka len sa netreba pomylit
uloha bola za 9 bodov, a ked ste mali iba cisto simulink tak 6 bodov istych

7. (za 6 bodov myslim) bol zadany nejaky dynamicky system  (pomocou zlomku)
   trebalo ho zadat do stavoveho priestou
   sys=tf(num,denum)

   problem bol ze num denum trebalo najprv roznasobit polynomi
   conv(a,b)

   vypocitat hodnoty a vykreslit graf so vstupnum jednotkovym skokom
   step(sys)

   vykreslit grafy s inym casovym vstupo lsim(sys,funkcia,cas)

   vykreslit nyquistovu charakteristiku nyquist(sys)

tu pripajam aj priklad podobny tomu co som spravil:

num=[5 50];
den=conv([7 1],[3 1]);
den2=conv(den,[2 1]);
sys=tf(num,den2);
[a b]=step(sys);
t=[0:0.1:10];
u=exp(-t);
[c d]=lsim(sys,u,t);
[e f]=bode(sys);
subplot(311);
step(sys);
subplot(312);
lsim(sys,u,t);
subplot(313);
nyquist(sys);

8. riesit diferencialnu rovnicu pomocou ode45 (subst. kan. tvar + m-fajly +zadefinovanie pociatocnych podmienok,cas .. proste vsetko ku tomu co sa len da) !!!
   uloha za NAJVIAC BODOV (15b)

9. Program co vypocita hodnoty vyrazu (bol zadany) pre 10 roznych hodnot a,b,c (Trebalo to robit cez for a este aj podmienku tam trebalo pouzit)

??. niekde bolo treba asi aj inline
    dalej bolo treba v niektorom priklade cykly spravit a vetvenia

10. este nejaka otazka bola len si nepamatam uz

HLAVNE SA ZAMERAJTE NA LDR, NDR, SIMULINK, blokove schemy, dynamicke systemy a kazdy urcite bude mat vetvenia a cykly.

poradie otazok bolo ine len som pisal co mi napadlo, a vela z toho co som napisal bolo aj minuly rok.
skupina B mala namiesto dynamickeho systemu zjednodusit blokovu schemu (so vsetkym cize maticu q, prikazy connect atd)
mozno ju tu niekto doplni ..

Hlavne sa nebojte, ona vas bude strasit vselijako ale potom vam pomoze na ustnej ;)

Este nejake otazky co boli na ustnej:  funkcie funkcii (quad, ode45, fzero atd .. chcela k nim aj priklady)
                                                   riadiace cykly ci ako sa to vola (vsetky for, while, switch, case, if, if else .. ku kazdemu priklad resp. ich vsetky zahrnut )
                                                  u nas chcela priklad na runge kutu, resp. aspon nejaku NDR, napisat k tomu priklad atd ..

Drzim palce  :bananawrite:

radix

Mohol by niekto hodit algoritmus na Runge-Kutta 4. radu? ??? ten co chce teta jadlovska ku skuske...

dabol

Ak sa najde dobra dusa, ktora ma vypracovane otazky zo skusky... PLZ nevahajte, Boh vam to vynahradi na detoch ;)
V spoločnosti chcem piť sám.

7ici7

#46
ludia z coho sa ucite, alebo ste sa ucili na skusku??? lebo neako na tej stranke matlab.fei...... je to dost rozhadzane a divne:) alebo sa to len mne zda:(ba

Treba si na skusku doniest aj popis k zadaniam???

radix

Quote from: radix on  18.01.2009, 01:54:32
Mohol by niekto hodit algoritmus na Runge-Kutta 4. radu? ??? ten co chcu ku skuske...

Sxx

ako je to s bodovanim? a vlastne tymi bodmi za ucast na prednaskach?


radix

10b ak si bol na vsetkych