Diskretinio laiko sistemų testavimas

Matlabo terpėje įsisavinti motodiką diskretaus laiko sistemų tiesiškumui, invariantiškumui bei jų stabilumui nustatyti.
Darbo eiga : Pirmiausia laboratorinį darbą pradėsime nuo matlab terpės paruošimo, toliau vykdysime užduoties punktus ir bandysime išsiaiškinti ir patikrinti ar ši duotoji sistema yra tiesiška.
Sistemos tiesiškumo tikrinimas.
1 pav. Moduliatoriaus įėjimai : pjūklinė funkcija (a) ir impulsų seka (b)
Signalus, kai alfa=2,5 o beta= -1,5, kursime taip:
>>y1=x.*cos(pi*n/10);y2=y.*cos(pi*n/10);y3=2.5*y1+(-1.5*y2);
v=y1*2.5;w=y2*(-1.5);c1=2.5*x;q=c1.*cos(pi*n/10); c2=(-1.5)*y;g=c2.*cos(pi*n/10);
vaizduosime : >>subplot(2,2,1);stem(n,v);
%kai alfa =2.5 beta =-1.5 tikrinam ar homogenine
xlabel('vn');subplot(2,2,2);stem(n,w);xlabel('wn');subplot(2,2,3);stem(n,q);
xlabel('qn');subplot(2,2,4);stem(n,g);xlabel('gn');
pause
2 pav. Signalai : v(n), w(n), q(n) ir g(n) .
Įėjimai x1(n) ir x2(n) pateikti 1 pav. , o sekos 2 pav. Analizuodami 2 pav. Sekas darome išvadą, kad sekų v(n) ir w(n) bei q(n) ir g(n) atitinkamos reikšmės sutampa. Taigi nešiklio amplitudės moduliatorius – homogeninė sistema.
Dabar nustatykime, ar moduliatoriui galioja adityvumo sąlyga, tardami, kad alfa=beta=1. Signalus kursime taip :
y3=y1+y2; x1=x+y; z=x1.*cos(pi*n/10); subplot(2,1,1);stem(n,y3);xlabel('a');subplot(2,1,2);stem(n,z); xlabel('b');
pause
3 pav. Moduliatoriui taikomos adityvumo sąlygos tikrinimo rezultatai. Sekos: y(a), z(b).
3 pav. pavaizduotos sekos teigia mums, kad mūsų nagrinėjamas moduliatorius yra adityvioji sistema. Taigi moduliatorius keičiant signalo amplitudę – tiesinė sistema.
Sistemos invariantiškumo tikrinimas.
Nustatysime ar moduliatorius yra laike nekintanti sistema . Moduliatorių tikrinsime pagal LTI sistemos įėjimo ir išėjimo sekas , kai x(n) – 5 atskaitų pjūklinė funkcija.
>>n=[0:4];(n)>=0];x=n.*x;sublot(2,2,1);stem(n,x);xlabel(‘a’);
Patrauksime x(n) per vieną atskaitą į dešinę ir gausime x(n-1):
>>n1=n+1;x1=x; seka bus pavaizduota 4 pav. (b)
>>subplot(2,2,2);stem(n1,x1);xlabel(‚b‘);
Rasime: >>q=x.*cos(pi*n/10); Patrauksime q(n) per vieną atskaitą į...

Daugiau informacijos...

★ Klientai rekomenduoja

Šį rašto darbą rekomenduoja mūsų klientai. Ką tai reiškia?

Mūsų svetainėje pateikiama dešimtys tūkstančių skirtingų rašto darbų, kuriuos įkėlė daugybė moksleivių ir studentų su skirtingais gabumais. Būtent šis rašto darbas yra patikrintas specialistų ir rekomenduojamas kitų klientų, kurie po atsisiuntimo įvertino šį mokslo darbą teigiamai. Todėl galite būti tikri, kad šis pasirinkimas geriausias!