Comparison of MUSIC and ESPRIT using MATLAB. Besides the algorithm development a

Question

Comparison of MUSIC and ESPRIT using MATLAB.

Besides the algorithm development and the MATLAB code, your output should also contain the
following five parts:

(a) In MATLAB, retrieve your signal set using command: ,-, num2str(naam)]; filename = [ 'signal set, num2 str (nstu) load (filename) 'nstu' is your designated student number, 'ndam' is a flag indicating whether you are using the signal set that contains a damped (ndam=1) signal or an undamped (ndam=0) one. In the signal set, vector 'v' is the complex-valued sampled voltage sequence, and 'fs' is the sampling frequency. Then, analyze the frequency mode of the voltage sequence 'v' using MUSIC and ESPRIT Please name your program hw4.m and start your MATLAB code with function hw4 (nstu, ndam) [HINT: You should work with both the damped case and the undamped one. Each data set also contains a vector 'fr,' you should NOT need use that for this assignment. Please also use the complex magnitude estimation method instead of the absolute value one. You can also use K -7 directly without having to estimate the number of modes in both methods.] (b) You can directly obtain the signal magnitudes |Ak], initial phases frequencies wk, and the damping factor Bk (ESPRIT only) for each of the K frequency modes. From the estimated | Akl, k, wk, and 3k values, describe how you can obtain the corresponding states a [n] Then implement this in MATLAB HINT: For MUSIC, you can assume that 3k = 0,Yk.

Explanation / Answer

N = 10;
ula = phased.ULA('NumElements',N,'ElementSpacing',0.5)
fc = 300e6;
fs = 8192;
lambda = physconst('LightSpeed')/fc;
pos = getElementPosition(ula)/lambda;
ang1 = [90;73]; ang2 = [90;68];
angs = [ang1 ang2];
Nsamp = 1024;
noisePwr = 0.01;
rs = rng(2012);
signal=sensorsig(pos,Nsamp,angs,noisePwr);
ang_true=az2broadside(angs(1,:),angs(2,:))
rootmusicangle=phased.RootMUSICEstimator('SensorArray',ula,...'OperatingFrequency',fc,...'NumSignalsSource','Property','NumSignals',2)
rootmusicangle.ForwardBackwardAveraging = true;
ang = rootmusicangle(signal)

espritangle=phased.ESPRITEstimator('SensorArray',ula,...'OperatingFrequency',fc,'ForwardBackwardAveraging',true,..'NumSignalsSource','Property','NumSignals',2)

ang = espritangle(signal)
musicangle = phased.MUSICEstimator('SensorArray',ula,...'OperatingFrequency',fc,'ForwardBackwardAveraging',true,...'NumSignalsSource','Property','NumSignals',2,...'DOAOutputPort',true)

[~,ang] = musicangle(signal)
plotSpectrum(musicangle)

Related Questions

Navigate

• Browse (All)
• Browse C
• Subjects

---

---

Integrity-first tutoring: explanations and feedback only — we do not complete graded work. Learn more.