{[ promptMessage ]}

Bookmark it

{[ promptMessage ]}

# ass ans - Direction Finding d = 0.4 ff = 400 ww = 2*pi*ff c...

This preview shows pages 1–6. Sign up to view the full content.

Sheet1 Page 1 3.1 d z1 = 10 * exp(-1i*2*pi/3) z2 = -5+5j zvect(z1,'r') hold zvect(z2, 'b') zvect(z1*z2, 'g') zprint(z1*z2) Z = X + jY Magnitude Phase Ph/pi Ph(deg) 68.3 18.3 70.71 0.262 0.083 15.00 3.1.h z1=10*exp(-j*2*pi/3) z2=-5+j*5 subplot(221) subplot(222) subplot(223) subplot(224) N = 200 r = sqrt((1:N)/50.*(1:N)/50+2.25) sig = exp (j*2*pi*r) plot((1:N)/50, r, '-mo') hold plot((1:N)/50, sig, '-mo') function [x,t] = one_cos(A,w,phi,dur) t = (-dur/2): 1/(20*(w/(2*pi))): (dur/2) x = A*cos(w*t+phi) plot(t,x),grid function [xx,tt] = syn_sin(fk, Xk, fs, dur, tstart) if nargin<5, tstart=0 end if length(fk)~=length(Xk), error('The length of fk and Xk are unequal') end xx = (1:(dur-tstart)*fs+1)*0 for k=1:length(Xk) tt=tstart:1/fs:dur xa(k,:) = Xk(k)*exp(1i*2*pi*fk(k)*tt) xx=xx+xa(k,:)

This preview has intentionally blurred sections. Sign up to view the full version.

View Full Document
Sheet1 Page 2 end plot(tt,xx) tt0 = o:.01:1 [xx0,tt0] = syn_sin([0,100,250],[10,14*exp(-pi/3),8*j],10000,0.1,0) plot(xx0,tt0) syn_sin([1,1,1],[2,2*exp(-j*pi*1.25),1-j],10000,10,-0.5)

This preview has intentionally blurred sections. Sign up to view the full version.

View Full Document

This preview has intentionally blurred sections. Sign up to view the full version.

View Full Document
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: Direction Finding d = 0.4 ff = 400 ww = 2*pi*ff c = 1000/3 lambda = c/ff yr = 100 xx = -50:1:300 [X1,X2,theta_true] = DF_gen(xx) Xout = X1.*conj(X2) phaseout = angle(Xout) theta = asin(c/d*(phaseout/ww)) plot(xx,theta*180/pi, 'r-') hold on, plot( xx, theta_true, 'g-'), hold off The process of decomposing any periodic function into its constituent sine or cosine waves is called Fourier analysis. Once the harmonic content of a sustained wave from Fourier analysis are known, the capability of synthesizing that sound fro m Sheet1 Page 3 hold zvect(z1,'r') zvect(z2,'b') hold zvect(conj(z1),'r'),zvect(conj(z2),'b'),ucplot hold zvect(conj(1/z1),'r') zvect(conj(1/z2),'b'),ucplot hold zvect(z1*z2) ucplot %-- gives 20 samples per period Sheet1 Page 4 Sheet1 Page 5 ucplot Sheet1 Page 6...
View Full Document

{[ snackBarMessage ]}