Re: Matlab 6-problemy z uzyskaniem charakterystyk czasowych danej
Masz problem? Zapytaj na forum elektroda.pl
From: Waldemar Krzok <waldemar.krzok_at_nospam_charite.de>
Subject: Re: Matlab 6-problemy z uzyskaniem charakterystyk czasowych danej
Date: Thu, 20 Jan 2005 13:21:28 +0100
Piotr wrote:
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
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