 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 611 | 回復(fù): 1 | ||
[求助]
程序運(yùn)行結(jié)果不對，自己查不出來錯(cuò)，請各位大俠幫忙
|
|
下面是寫的一段UCA ROOT-MUSIC程序，但是運(yùn)行結(jié)果不對，不是設(shè)置的[20,50]。同樣是本程序的噪聲子空間，用譜峰搜索的MUSIC就是對的，但是用ROOT-MUSIC就不對。應(yīng)該就是root-music那段程序有錯(cuò)，但是自己查也查不出來錯(cuò)誤，請各位大俠幫幫忙，非常感謝！ sig_num = 1000;%%信號(hào)的快拍數(shù) snr=10; rou=0; x1=randint(sig_num,1,4);%產(chǎn)生sig_sum*1的矩陣，數(shù)值在0-3之間 x2=randint(sig_num,1,4); x3=randint(sig_num,1,4); x222=rou*x1+sqrt(1-(abs(rou))^2)*x3;%產(chǎn)生相關(guān)的兩路信號(hào) x2=round(x222); for i=1:sig_num if x2(i)>=4 x2(i)=3; elseif x2(i)<=0 x2(i)=0; end end y1=dmod(x1,1,0.1,4,'psk',4); y1=y1(1:sig_num); y2=dmod(x2,1,0.1,4,'psk',4); y2=y2(1:sig_num); x=[y1';y2']; s=x; q=2; p=12;%智能天線陣元數(shù) fc=10.7*10^6; lemda=3*10^8/fc; banjing=20; jh=pi./180; k0=2*pi/lemda; h=ceil(k0*banjing); P=2*h+1; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %均勻圓陣轉(zhuǎn)化為均勻線陣 alpha=[90,90]; theta=[20,50]; kesei=2*pi*banjing./lemda;%% j=sqrt(-1); n=[0: (p-1)]'; gama=2*pi*n/p; for i=1:q for n=0: (p-1) a(n+1,i)=exp(j*kesei*sin(theta(i)*jh)*cos(alpha(i)*jh-gama(n+1))); end end x11=a*s; vn=var(x11.')./(10^(snr/10)); n11=randn(p,sig_num)+j*randn(p,sig_num); a1=(vn./var(n11.')); b1=((n11.'-ones(size(n11,2),1)*mean(n11.'))).'; for i=1:p c1(i, : )=a1(i)*b1(i, : ); end xn=x11+c1; for i=-h:h cv1(i+h+1)=j^(-abs(i)); end cv=diag(cv1); for m=-h:h for i=0:p-1 w(i+1,m+h+1)=exp(-j*m*gama(i+1)); end end V=1./sqrt(p)*w; Fe=cv*V'; for i=-h:h alph(i+h+1)=2*pi*i/P; end for i=-h:h for m=-h:h v(m+h+1,i+h+1)=exp(j*2*pi*i*m./P); end end W=(1/sqrt(P))*v; Fr=W'*Fe; rx=xn*xn'./sig_num; ry=Fr*rx*Fr'; r=real(ry); % r=Fe*rx*Fe'; [u,sig,v]=svd(r);%%svd奇異值分解，特征值降序排列 signalspace=u(:,q); noisespace=u(:,q+1:7); % vv=conj(noisespace)*noisespace'; UU=noisespace*noisespace.'; %%%%%%root-music%%%%%%%%%%%%%%%%% J1=besselj(1:h,k0*banjing*sin(90*jh)); J11=diag([fliplr(J1),besselj(0,k0*banjing*sin(90*jh)),J1]); Q=J11*W*UU*W'*J11; num1=1; num2=1; A1(1,1: P-1)=0; A2(1,1: P-1)=0; for l=P-1:-1:1 j=P; jj=j-l; for n=1:jj i=j-l; A1(1,num1)=A1(1,num1)+Q(i,j); A2(1,num2)=A2(1,num2)+Q(j,i); j=j-1; j; end num1=num1+1; num2=num2+1; end A1; A2; A=[A1, trace(Q),fliplr(A2)]; ra=roots([A]); rb=ra(abs(ra)<1); [dumm,I]=sort(abs(abs(rb)-1)); w=angle(rb(I(1:q))); fi=-w/jh [ Last edited by 影月水瀾 on 2013-7-2 at 21:53 ] |
回復(fù)此樓» 猜你喜歡
085601 329分調(diào)劑
已經(jīng)有3人回復(fù)
-
本2一志愿C9-333分，材料科學(xué)與工程，求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿大連理工大學(xué)材料求調(diào)劑
已經(jīng)有5人回復(fù)
-
材料與化工調(diào)劑一志愿大連海事085600，349
已經(jīng)有6人回復(fù)
-
一志愿中海洋320化學(xué)工程與技術(shù)學(xué)碩求調(diào)劑
已經(jīng)有4人回復(fù)
-
材料科學(xué)與工程調(diào)劑
已經(jīng)有11人回復(fù)
-
311求調(diào)劑一志愿合肥工業(yè)大學(xué)
已經(jīng)有4人回復(fù)
-
一志愿食品科學(xué)與工程083200求調(diào)劑
已經(jīng)有4人回復(fù)
-
材料科學(xué)與工程求調(diào)劑
已經(jīng)有10人回復(fù)
-
一志愿北化085600材料專碩275|有文章專利｜求調(diào)劑
已經(jīng)有11人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- 求Matlab高手解決線性方程組的迭代求解問題
已經(jīng)有17人回復(fù)
- q2r.x運(yùn)行出錯(cuò)
已經(jīng)有12人回復(fù)
- packmol如何構(gòu)建初始結(jié)構(gòu)供dl_poly使用？
已經(jīng)有9人回復(fù)
- 急求fortran運(yùn)行錯(cuò)誤原因，在線等
已經(jīng)有7人回復(fù)
- 高斯安裝成功，運(yùn)行提示出錯(cuò)。
已經(jīng)有10人回復(fù)
- 從excel讀取數(shù)據(jù)多項(xiàng)式擬合出錯(cuò)
已經(jīng)有3人回復(fù)
- C++編程，關(guān)于循環(huán)結(jié)構(gòu)的，大家看看我這程序哪里錯(cuò)了？
已經(jīng)有21人回復(fù)
- 一個(gè)區(qū)間搜索的c語言程序 出問題了
已經(jīng)有8人回復(fù)
- 回答一些關(guān)于Multiwfn的疑問以及未來Multiwfn的發(fā)展打算
已經(jīng)有37人回復(fù)
- 【求助】此程序沒有錯(cuò)誤，怎么運(yùn)行不出來？
已經(jīng)有6人回復(fù)
- 【求助】fortran程序未出錯(cuò)，但運(yùn)行不出結(jié)果，中途中斷！
已經(jīng)有6人回復(fù)
- 【求助】運(yùn)行出錯(cuò) 虛擬內(nèi)存
已經(jīng)有10人回復(fù)
- 【求助】MS4.3在運(yùn)行，卻一直沒有CASTEP這個(gè)程序在工作，計(jì)算結(jié)果也遲遲不出現(xiàn)。。。?/a>
已經(jīng)有15人回復(fù)
- 【求助】合成氨過程模擬 結(jié)果錯(cuò)誤，怎么辦？并問下如何修改收斂次數(shù)！【已解決】
已經(jīng)有11人回復(fù)
- 【求助】修改MATLAB支持向量機(jī)核函數(shù)程序遇到的疑問？求助
已經(jīng)有5人回復(fù)
2樓2018-06-24 12:41:00
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 影月水瀾 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 285求調(diào)劑 +5 | AZMK 2026-03-30 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 理學(xué)07化學(xué) 303求調(diào)劑 +4 | 睿08 2026-03-27 | 4/200 |
|
|
|
[有機(jī)交流]
 甲基亞磺磺酸鈉和甲基磺酸酯反應(yīng)機(jī)理
 10+3
|
kaobao456 2026-03-29 | 4/200 |
|
|
|
[考研] 085701環(huán)境工程求調(diào)劑 +11 | 多久上課 2026-03-27 | 12/600 |
|
|
|
[考研] 材料工程專碩求調(diào)劑 +8 | hyl3153942 2026-03-29 | 8/400 |
|
|
|
[考研] 346求調(diào)劑 一志愿070303有機(jī)化學(xué) +7 | 蘿卜燉青菜 2026-03-28 | 8/400 |
|
|
|
[考研] 材料學(xué)碩333求調(diào)劑 +14 | 北道巷 2026-03-24 | 14/700 |
|
|
|
[考研] 334分 一志愿武理 材料求調(diào)劑 +16 | 李李不服輸 2026-03-26 | 16/800 |
|
|
|
[考研] 298求調(diào)劑 +3 | 什么是胖頭魚 2026-03-30 | 5/250 |
|
|
|
[考研] 348求調(diào)劑 +6 | 小懶蟲不懶了 2026-03-28 | 6/300 |
|
|
|
[考研] 337求調(diào)劑 +6 | 《樹》 2026-03-29 | 6/300 |
|
|
|
[考研] 279求調(diào)劑 +4 | 蝶舞輕繞 2026-03-29 | 4/200 |
|
|
|
[考研] 356求調(diào)劑 +3 | gysy?s?a 2026-03-28 | 3/150 |
|
|
|
[考研] 0703化學(xué)調(diào)劑，求導(dǎo)師收 +9 | 天天好運(yùn)來上岸?/a> 2026-03-24 | 10/500 |
|
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-27 | 6/300 |
|
|
|
[考研] 材料求調(diào)劑一志愿哈工大324 +7 | 閆旭東 2026-03-28 | 9/450 |
|
|
|
[考研] 331環(huán)境科學(xué)與工程求調(diào)劑 +3 | 熠然好運(yùn)氣 2026-03-27 | 3/150 |
|
|
|
[考研] 一志愿南師大0703化學(xué) 275求調(diào)劑 +4 | Ripcord上岸 2026-03-27 | 4/200 |
|
|
|
[考研] 調(diào)劑推薦 +5 | 清酒714 2026-03-26 | 6/300 |
|
|
|
[考研] 求b區(qū)院校調(diào)劑 +4 | 周56 2026-03-24 | 5/250 |
|