mihaispr
Administrator
 Inregistrat: acum 17 ani
Postari: 2142
|
|
In cele ce urmeaza vor fi prezentate anumite probleme care aplica script-urile FRECT,FTRIUNGHI,FSINUS,FCOSINUS in functie de natura semnalului(rectangular,triunghiular,sinusoidal,puls sinusoidal pt.script-ul FCOSINUS). Prin aceste script-uri se calculeaza integrala Fourier si vom avea si o reprezentare grafica.
Obs. Problemele in matlab se scriu in Matlab editor si se salveaza cu extensia .m
Realizati in matlab urmatoarele probleme:
Problema 1: Se da semnalul rectangular(aperiodic) de durata (k/100) ms(mili-secunde),simetric fata de originea timpului si de amplitudine wmax=100000/k.Utilizati scriptul aferent unui semnal rectangular si realizati si reprezentarea grafica in matlab. Fie k=2 pt.toate problemele de mai jos.
% Scriptul FRECT
% la noi a=k/100ms => a=k/100000s si wmax=100000/k, k=2
syms x w % declararea variabilele simbolice cu ajutorul cuvantului rezervat syms
wmax=50000; %initializare amplitudine maxima
int(1.0/a*exp(-j*w*x),-1/100000,1/100000); % se calculeaza integrala Fourier
ezplot(ans,[-wmax wmax]) % se reprezinta grafic
hold on %se retine graficul
u=-wmax:wmax:wmax; %u=val.initiala:pas:val finala;
y=0.0*u;
plot(u,y) % se traseaza orizontala y=0
hold off
Problema 2: Se da semnalul rectangular(aperiodic) de durata (k/100) ms(mili-secunde),simetric fata de originea timpului si de amplitudine wmax=100000/k.Utilizati scriptul aferent unui semnal rectangular si realizati si reprezentarea grafica in matlab. La acest semnal se aduna semnalul rectangular(aperiodic) de durata k/200ms,simtric fata de origine si de amplitudine 200000/k; k=2
% Scriptul FRECT
% la noi a1=k/100ms => a=k/100000s si wmax1=100000/k, k=2
% la noi a2=k/200ms => a=k/200000s si wmax2=200000/k, k=2
syms x w i1 i2 s % declararea variabilele simbolice cu ajutorul cuvantului rezervat syms
wmax1=50000; %initializare amplitudine maxima
i1=int(1.0/a*exp(-j*w*x),-1/100000,1/100000); % se calculeaza integrala Fourier
ezplot(ans,[-wmax1 wmax1]) % se reprezinta grafic
hold on %se retine graficul
u=-wmax1:wmax1:wmax1; %u=val.initiala:pas:val finala;
y=0.0*u;
plot(u,y) % se traseaza orizontala y=0
hold off
wmax2=100000;
i2=int(1.0/a*exp(-j*w*x),-1/200000,1/200000); % se calculeaza integrala Fourier
ezplot(ans,[-wmax2 wmax2]) % se reprezinta grafic
hold on %se retine graficul
u=-wmax2:wmax2:wmax2; %u=val.initiala:pas:val finala;
y=0.0*u;
plot(u,y) % se traseaza orizontala y=0
s=i1+i2;
hold off
|
|
