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; % Operating frequency
fs = 8192; % Sampling frequency
lambda = physconst('LightSpeed')/fc; % Wavelength
pos = getElementPosition(ula)/lambda; % Element position in wavelengths
ang1 = [90;73]; ang2 = [90;68]; % Direction of the signals
angs = [ang1 ang2];
Nsamp = 1024; % Number of snapshots
noisePwr = 0.01; % Noise power
rs = rng(2012); % Set random number generator
signal = sensorsig(pos,Nsamp,angs,noisePwr);
