Code showing standing wave due to terminator absorbing wave power
I've discussed how reflection results in a standing wave. It can be useful to visualise how complete and partial reflection look.At some point I'll make this into something more fun and user-friendly. For now, if you are really curious to visualise how a standing wave looks when 96% of the energy has been absorbed, then paste the following into Octave (Matlab will do), and change Ar to 0.04:
clear;
close;
Ai = 1; %incident amplitude
Ar = 0.9; %amplitude of combined radiated and reflected waves - try varying between 0 (full absorption) and 1(no absorption).
NumWaves = 2;
wtmax = NumWaves*2*pi; %maximum time to plot to give NumWaves # of cycles
kx= 1:100;
wt = 1:wtmax/(20*NumWaves):wtmax %discretisation: 20 points per wave
for ii = 1:length(wt)
wi(ii,:) = Ai*sin(wt(ii)+kx);
wr(ii,:) = Ar*sin(wt(ii)-kx);
end
ws = wi+wr;
plot(ws)
No comments:
Post a Comment
Comments
Note: only a member of this blog may post a comment.