Plot +1 +1 -1 -1 rectangular pulse

17 visualizaciones (últimos 30 días)
nasz13
nasz13 el 12 de En. de 2018
Comentada: nasz13 el 9 de Feb. de 2018
Hi,
I would like to ask for advice to help me to plot the rectangular pulse having +1 +1 -1 -1 in amplitude. Before, I plot with +1 -1 +1 -1, and this is fine because I use remainder, so it goes by odd and even. However, I don't know how to plot it with +1 +1 -1 -1 code. Can anyone please advice me on this? Below is my code;
for m = 0:Npulses
if rem(m,2) == 0 %i is even
Tx1 = exp(1j*pi*gamma*(t-m*PRI).^2);
W=ones(size(t));
W(t<m*PRI)=0;
W(t>(m*PRI+tau))=0;
s1 = Tx1.*W+s1;
else
Tx2 = -exp(1j*pi*gamma*(t-m*PRI).^2);%
W=ones(size(t));
W(t<m*PRI)=0;
W(t>(m*PRI+tau))=0;
s1 = Tx2.*W+s1;
end
end
where Npulses is number of pulses and PRI is Pulse Repetitive Interval.
  1 comentario
nasz13
nasz13 el 9 de Feb. de 2018
Hi,
At this moment, I need to change the code to +1 -1 -1 +1 and also -1 +1 +1 -1. Can anyone please help me with this, based on my code above?

Iniciar sesión para comentar.

Respuesta aceptada

Surya Pullakhandam
Surya Pullakhandam el 12 de En. de 2018
Editada: Surya Pullakhandam el 12 de En. de 2018
Previously the pattern was repeating after every 2 occasions (+1 -1)(+1 -1). So, you used rem function with respect to 2. Now the pattern has to repeat after every 4 occasions (+1 +1 -1 -1)(+1 +1 -1 -1). So, let us use rem function with respect to 4. I think the changes as done below will suffice your requirement.
for m = 0:Npulses
if rem(m,4) == 0 | rem(m,4) == 1 %1st and 2nd term
Tx1 = exp(1j*pi*gamma*(t-m*PRI).^2);
W=ones(size(t));
W(t<m*PRI)=0;
W(t>(m*PRI+tau))=0;
s1 = Tx1.*W+s1;
else
Tx2 = -exp(1j*pi*gamma*(t-m*PRI).^2);
W=ones(size(t));
W(t<m*PRI)=0;
W(t>(m*PRI+tau))=0;
s1 = Tx2.*W+s1;
end
end
  1 comentario
nasz13
nasz13 el 13 de En. de 2018
Thank you so much!! I really appreciate it, Surya.

Iniciar sesión para comentar.

Más respuestas (0)

Categorías

Más información sobre Matrix Indexing en Help Center y File Exchange.

Etiquetas

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by