版塊導航: 正在加載中...

登錄注冊

應《網絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯網跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

北京石油化工學院2026年研究生招生接收調劑公告

返回列表

影月水瀾

新蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 358.5
紅花: 1
帖子: 53
在線: 61.1小時
蟲號: 1826817
注冊: 2012-05-20
專業(yè): 敏感電子學與傳感器

[求助] 程序運行結果不對，自己查不出來錯，請各位大俠幫忙

下面是寫的一段UCA ROOT-MUSIC程序，但是運行結果不對，不是設置的[20,50]。同樣是本程序的噪聲子空間，用譜峰搜索的MUSIC就是對的，但是用ROOT-MUSIC就不對。應該就是root-music那段程序有錯，但是自己查也查不出來錯誤，請各位大俠幫幫忙，非常感謝！

sig_num = 1000;%%信號的快拍數
snr=10;
rou=0;
x1=randint(sig_num,1,4);%產生sig_sum*1的矩陣，數值在0-3之間
x2=randint(sig_num,1,4);
x3=randint(sig_num,1,4);
x222=rou*x1+sqrt(1-(abs(rou))^2)*x3;%產生相關的兩路信號
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;%智能天線陣元數
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;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%均勻圓陣轉化為均勻線陣
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 ]

回復此樓

» 猜你喜歡

學碩274求調劑已經有11人回復
一志愿鄭大材料工程290求調劑已經有12人回復
環(huán)境工程 085701，267求調劑已經有8人回復
083000學碩274求調劑已經有10人回復
一志愿085600中科院寧波所276分求調劑已經有4人回復
材料專業(yè)申博已經有3人回復
307求調劑已經有20人回復
285求調劑已經有8人回復
材料與化工（0856）304求B區(qū)調劑已經有8人回復
370求調劑已經有3人回復

» 本主題相關價值貼推薦，對您同樣有幫助:

求Matlab高手解決線性方程組的迭代求解問題已經有17人回復
q2r.x運行出錯已經有12人回復
packmol如何構建初始結構供dl_poly使用？已經有9人回復
急求fortran運行錯誤原因，在線等已經有7人回復
高斯安裝成功，運行提示出錯。已經有10人回復
從excel讀取數據多項式擬合出錯已經有3人回復
C++編程，關于循環(huán)結構的，大家看看我這程序哪里錯了？已經有21人回復
一個區(qū)間搜索的c語言程序出問題了已經有8人回復
回答一些關于Multiwfn的疑問以及未來Multiwfn的發(fā)展打算已經有37人回復
【求助】此程序沒有錯誤，怎么運行不出來？已經有6人回復
【求助】fortran程序未出錯，但運行不出結果，中途中斷！已經有6人回復
【求助】運行出錯虛擬內存已經有10人回復
【求助】MS4.3在運行，卻一直沒有CASTEP這個程序在工作，計算結果也遲遲不出現。。。?/a> 已經有15人回復

【求助】合成氨過程模擬結果錯誤，怎么辦？并問下如何修改收斂次數！【已解決】已經有11人回復
【求助】修改MATLAB支持向量機核函數程序遇到的疑問？求助已經有5人回復

1樓 2013-07-01 20:04:08

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

suisuimumu

新蟲 (初入文壇)

應助: 0 (幼兒園)
帖子: 1
在線:
蟲號: 9182627
注冊: 2018-06-24

樓主，能不能給看一下使用譜峰搜索的代碼呢？非常感謝！

贊一下

回復此樓

2樓2018-06-24 12:41:00

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關版塊跳轉我要訂閱樓主影月水瀾的主題更新