Determinarea sensibilitatii

DOWNLOAD REFERAT

Preview referat: Determinarea sensibilitatii

%program pentru reprezentarea grafica a functiilor de sensibilitate ale lui y(t) in raport cu parametrii m, KR si KV.

t=0:.0001:0.1;
a=cos(100*t); %intrare armonicã.
%a=sin(100*t); %intrare sinusoidala.
%a=1+0*t; %intrare treapta.

%matricile sistemului inainte de aplicarea functiilor de sensibilitate
A=[0 0;-10^3 -10]; %matricea A a sistemului
B=[0;1]; %matricea B a sistemului
C=[1 0;0 1]; %matricea C a sistemului
D=[0;0]; %matricea D a sistemului

%determinarea matricilor noului sistem dupa aplicarea functiilor de sensibilitate.
AM=[0 0;10^5 10^3]; %derivata lui A in raport cu parametrul m
AKR=[0 0;-10 0]; %derivata lui A in raport cu parametrul KR
AKV=[0 0;0 -10]; %derivata lui A in raport cu parametrul KV
B1=[0;0]; %derivatele lui B in raport cu parametrii m, KR, KV
X=[0 0;0 0]; % notatii folosite pentru calcularea matricilor
Y=[1 1;1 1]; %notatii folosite pentru calcularea matricilor

Atot=[A X X AM; X A X AKR; X X A AKV; X X X A]; %matricea A a noului sistem
Btot=[B1; B1;B1;B]; %matricea B a noului sistem
Ctot=[Y X X X; X Y X X; X X Y X; X X X Y]; %matricea C a noului sistem
Dtot=[B1;B1;B1;B1]; %matricea D a noului sistem

%reprezentarea grafica a functiilor de sensibilitate si a iesirii noului sistem.
y=lsim(Atot,Btot,Ctot,Dtot,a,t);

subplot(211);
plot(t,y(:,1),'k');grid;
xlabel('Timp[secunde]');
ylabel('Amplitudine');
title('Functia de sensibilitate in raport cu m');

subplot(212);
plot(t,y(:,2),'k');grid;
xlabel('Timp[secunde]');
ylabel('Amplitudine');
title('Functia de sensibilitate in raport cu KR');

pause;
clf;
subplot(211);
plot(t,y(:,3),'k');grid;
xlabel('Timp[secunde]');
ylabel('Amplitudine');
title('Functia de sensibilitate in raport cu KV');

subplot(212);
plot(t,y(:,4),'k',t,y(:,4)+y(:,1).*100,'r',t,y(:,4)+y(:,2).*8.0e+3,'g',t,y(:,4)+y(:,3).*8.0e-7,'b');grid;
xlabel('Timp[secunde]');
ylabel('Amplitudine');
title('Iesirea circuitului'); « mai multe referate din Informatica