| 5 | 1/1 | 返回列表 |
| 查看: 647 | 回復(fù): 3 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
promisewind金蟲(chóng) (著名寫(xiě)手)
|
[交流]
JJM算法的問(wèn)題 已有1人參與
|
||
|
自己編寫(xiě)了一個(gè)JJM算法,可是怎么出現(xiàn)的結(jié)果不能估計(jì)出弱信號(hào)的doa,得到的是與強(qiáng)信號(hào)對(duì)稱(chēng)的一個(gè)方向的信號(hào),不知道那里錯(cuò)了,還請(qǐng)大家賜教。 N=256; %設(shè)定采樣點(diǎn)數(shù) M =12; %設(shè)定陣元數(shù) fs=100; %設(shè)定采樣頻率 kd=1/2; %陣元間距與波長(zhǎng)的比值 f=[14 6]; %設(shè)兩個(gè)信號(hào)元的頻率,單位為MHz theta = [5, 34]*pi/180; %入射方向角向量 SNR=[-5, 40]; % 入射信號(hào)信噪比 AMP = 10.^(SNR./20); m=0:M-1; n=0:N-1; S = repmat(AMP',1,N).*exp(j*(2*pi*f.'/fs*n)); %產(chǎn)生兩個(gè)遠(yuǎn)場(chǎng)正弦信號(hào) A = exp(j*2*pi*kd*m'*sin(theta)); %方向矩陣 X=A*S; randn('seed',sum(100*clock)); noise = randn(M,N) + j*randn(M,N); X=X+noise; %接受數(shù)據(jù)矩陣 Rx=X*(X')/N; theta_stro = [34]*pi/180; T1=eye(M,M-1); for l=2:M T1(l,l-1)=-exp(-j*pi*sin(theta_stro(1))); end TT=T'; Rx1=TT*Rx*TT'; %計(jì)算協(xié)方差矩陣 %********************************************************************* [V1,D1]=eig(Rx1,2*TT*TT'); %自相關(guān)矩陣的特征分解,,得到特征值矩陣和與特征值一一對(duì)應(yīng)的特征向量矩陣 [Q1,index1]=sort(diag(D1)); %對(duì)特征值排序 for k=1:M-length(theta) %將與特征值對(duì)應(yīng)的特征向量也按照特征值的順序排列 En1(:,k)=V1(:,index1(k)); %取排了序的特征向量的前M-length(theta)個(gè)來(lái)構(gòu)造噪聲子空間En end a_theta=T'*exp(-j*2*pi*kd*[0:M-1].'*sin(theta_axis*pi/180)); for k=1:MD %進(jìn)行譜峰搜索 d1(k)=a_theta(:,k)'*a_theta(:,k)/(a_theta(:,k)'*En1*En1'*a_theta(:,k)); k=k+1; end plot(theta_axis,10*log10(abs(d1)),'k'); xlabel('到達(dá)角/\circ'); ylabel('歸一化空間譜/dB'); |

新蟲(chóng) (初入文壇)

金蟲(chóng) (著名寫(xiě)手)

金蟲(chóng) (著名寫(xiě)手)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 07化學(xué)280分求調(diào)劑 +4 | 722865 2026-03-23 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 081700 調(diào)劑 267分 +7 | 迷人的哈哈 2026-03-23 | 7/350 |
|
|
[考研] 一志愿國(guó)科過(guò)程所081700,274求調(diào)劑 +3 | 三水研0水立方 2026-03-23 | 3/150 |
|
|
[考研] 335分 | 材料與化工專(zhuān)碩 | GPA 4.07 | 有科研經(jīng)歷 +4 | cccchenso 2026-03-23 | 4/200 |
|
|
[考研] 291 求調(diào)劑 +4 | 化工2026屆畢業(yè)?/a> 2026-03-21 | 5/250 |
|
|
[考研] 0854電子信息求調(diào)劑 324 +3 | Promise-jyl 2026-03-23 | 3/150 |
|
|
[考研] 276求調(diào)劑。有半年電池和半年高分子實(shí)習(xí)經(jīng)歷 +8 | 材料學(xué)257求調(diào)劑 2026-03-23 | 9/450 |
|
|
[考研] 333求調(diào)劑 +6 | 87639 2026-03-21 | 10/500 |
|
|
[考研] 293求調(diào)劑 +3 | 濤濤Wjt 2026-03-22 | 5/250 |
|
|
[考研] 287求調(diào)劑 +8 | 晨昏線與星海 2026-03-19 | 9/450 |
|
|
[考研] 324求調(diào)劑 +6 | lucky呀呀呀鴨 2026-03-20 | 6/300 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 妮妮ninicgb 2026-03-21 | 4/200 |
|
|
[考研] 313求調(diào)劑 +4 | 肆叁貳壹22 2026-03-19 | 4/200 |
|
|
[考研] 296求調(diào)劑 +4 | www_q 2026-03-20 | 4/200 |
|
|
[考研] 299求調(diào)劑 +5 | shxchem 2026-03-20 | 7/350 |
|
|
[考研] 22 350 本科985求調(diào)劑,求老登收留 +3 | 李軼男003 2026-03-20 | 3/150 |
|
|
[考研] 296求調(diào)劑 +6 | www_q 2026-03-18 | 10/500 |
|
|
[考研] 22408 344分 求調(diào)劑 一志愿 華電計(jì)算機(jī)技術(shù) +4 | solanXXX 2026-03-20 | 4/200 |
|
|
[考研] 一志愿武漢理工材料工程專(zhuān)碩調(diào)劑 +9 | Doleres 2026-03-19 | 9/450 |
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|