Matlab 6-problemy z uzyskaniem charakterystyk czasowych danej transmitancji

Witam. Wykonuje projekt z Podstaw Sterowania Cyfrowego, mam za zadanie przy zadanej transmitancji wykonac charakterystyki, wyprobowac roznego rodzaje regulatorow etc I teraz mam problem. a mianowicie napisany przeze mnie kod programu w matlabie daje mi bledne wyniki przy charakterystyce impulsowej i skokowej. Wyglada to dosc dziwnie tak jakby nie dzialala suma....zamieszczam kod

CHOICE = MENU('projekt', 'log amp', 'log fazy', 'log amp-fazowa','amp-fazowa', 'rzeczywista', 'urojona','impulsowa','skokowa','koniec') k=1; u=0.4; t2=0.1; t1=0.2; t0=0.5; w=0.01:0.01:1000; t=0.01:0.01:10; h=0.0; (...) elseif CHOICE==7

syms s; g=(k.*(1+u.*s))./((1+t1.*s).*(1+t0.*t0.*s.*s).*(1+t2.*s)); ilaplace(g); t=0.01:0.01:10; h=(-40/29*exp(-5*t)+15/13*exp(-10*t)+85/377*cos(2*t)+875/377*sin(2*t)); plot(t,h); TITLE('charakterystyka impulsowa'); XLABEL('t'); YLABEL('Amplituda'); " i teraz przy takim 'h' rysuje mi jedynie sinusoide!!! taka sama gdy h jest rowne tylko"85/377*cos(2*t)+875/377*sin(2*t));" temu czlonowi...jesli zostawie tylko "h=(-40/29*exp(-5*t)+15/13*exp(-10*t)" to rysuje ladne kombinacje e potencjalnych funkcji...Jako , ze z matlabem stykam się po raz pierwszy prosze was o pomoc Z gory dziekuje za pomoc

Reply to
Piotr
Loading thread data ...

to co prawda kompletnie ntg, ale mogę spróbować. Po pierwsze, po co liczysz transformatę odwrotną laplace'a, jak nie używasz wyniku do niczego. Po drugie, napisz, jak mają wyglądać równania prozą, to sprawdzę, czy dobrze zakodowałeś. Człony z exp są dość szybko w okolicach zera, więc to, że widzisz sinus (trochę kosinusa też tam się pałęta(, nie dziwota.

Waldek

Reply to
Waldemar Krzok

ElectronDepot website is not affiliated with any of the manufacturers or service providers discussed here. All logos and trade names are the property of their respective owners.