CRB_克拉美罗界_MUSIC算法CRB_阵列信号处理中克拉美罗界函数CRB_克拉美罗_CRB

%%%%%%双平行线阵 CRB 仿真 model = gen_model; truth= gen_truth(model); meas= gen_meas(model,truth); syms DOA1 phi1 DOA2 phi2 DOA3 phi3 signal = [DOA1 DOA2 DOA3; phi1 phi2 phi3]; Ax1 = exp(j*model.twpi*model.d.'*(sin(signal(1,1)*model.derad).*cos(signal(2,1)*model.derad)));%%%% direction matrix Ax Ax2 = exp(j*model.twpi*model.d.'*(sin(signal(1,2)*model.derad).*cos(signal(2,2)*model.derad)));%%%% direction matrix Ax Ax3 = exp(j*model.twpi*model.d.'*(sin(signal(1,3)*model.derad).*cos(signal(2,3)*model.derad)));%%%% direction matrix Ax Ax = [Ax1 Ax2 Ax3]; Ay1 = exp(j*model.twpi*model.D.'*(sin(signal(1,1)*model.derad).*sin(signal(2,1)*model.derad))); Ay2 = exp(j*model.twpi*model.D.'*(sin(signal(1,2)*model.derad).*sin(signal(2,2)*model.derad))); Ay3 = exp(j*model.twpi*model.D.'*(sin(signal(1,3)*model.derad).*sin(signal(2,3)*model.derad))); Ay = [Ay1 Ay2 Ay3]; clum = size(Ax,2); for i=1:clum A(:,i) = kron(Ay(:,i),Ax(:,i)); end a1 = A(:,1); a2 = A(:,2); a3 = A(:,3); D = [diff(a1,DOA1),diff(a2,DOA2),diff(a3,DOA3),diff(a1,phi1),diff(a2,phi2),diff(a3,phi3)]; CRB = countCRB(A,D,model,meas,truth);