function?[powerpowerdBdoasignal_out]=Beamforming_CircleArray（stylemethodscandthitasigRNf0fsW）
%--------------2012-02-17------------------------------
%?圓弧陣波束形成，360度全方位掃描
%?輸入:????????style?：信號形式?????1---窄帶信號（單頻信號）????2---寬帶信號（LFM）
%?????????????method?:?波束形成算法類型??
%?????????????????????????????????“1”---?時域的時延波束形成
%?????????????????????????????????“2”---?相移波束形成
%?????????????????????????????????“3”---?頻域波束形成
%?????????????????????????????????“4”---?經過三次FFT的頻域波束形成，延遲濾波器和輸入信號在頻域相乘
%??????????????scan?:?開角范圍（單位：度）
%???????????dthita?:?掃描步長?（單位：度）
%????????????????sig?:?圓陣輸入信號
%??????????????????R?:?圓陣半徑
%??????????????????N?:?陣元個數
%?????????????????f0?:?窄帶的中心頻率或者寬帶的最小頻率
%?????????????????fs?：采樣率
%??????????????????W?:?加權值
%?輸出：???????power?：能量
%????????????powerdB?:?能量歸一化分貝輸出
%??????????????doa?:?信號波達方向（?單位：度）
%?????????signal_out?:?波束在預成方向上疊加后的輸出信號，每一行代表一個方向上的輸出
%---------------------------------------------------------
c=1500;??%%默認聲速為1500m/s
[sig_Nsig_len]=size（sig）;
if（?sig_N?~=?N?）
????error（‘信號和圓陣個數不匹配‘）;
end
if?（?style?~=1?&?style?~=2??）
????error（‘請正確選擇輸入信號類型‘）;
end
if?（?method?~=1?&?method?~=2?&?method?~=3?&?method?~=4?）
????error（‘請正確選擇波束形成算法‘）;
end

psi?=?2?*?pi?*（0:N-1）?/?N;?%%圓心角
thita=[0:dthita:359]*pi/180;
power=zeros（1length（thita））;
scan_begin?=?thita?-?scan?*pi/180/?2;
scan_begin（?scan_begin?scan_end?=?thita?+?scan?*?pi/180/?2;
scan_end（scan_end?>=?2*pi）?=?scan_end（scan_end?>=?2*pi）-2*pi;
flag=zeros（length（thita）N）;?%%?標志位，對應位為1時，表示該號陣元參與計算
%%----------??確定預成方向上參與計算的陣元號存入flag中??----
for?a=thita
????????????index=round（a/（dthita*pi/180）+1）;
????????????if（?scan_end?（index）?/?（scan_begin?（index）+eps）?<=?（1.0+0.0001））
????????????????flag（indexfind（scan_begin?（index）./?（psi+eps）?<=（1.0+eps）））=1;
????????????????flag（indexfind（（psi-eps）?./（scan_end?（index）+eps）?<=?（1.0+eps）））=1;?
????????????else
????????????????flag（indexfind（scan_begin?（index）./（psi+eps）?<=?（1.0+eps）?&?psi?./（scan_end?（index）+eps）<=?（1.0+eps）））=1;
????????????end?
end
%%-------------------------------------------------------------
if（?method?==?1）????????????????%%----??時延波束形成對入射信號不限制，窄帶、寬帶均可??（?method?==?1）?----?
????beishu=1;
????%%----?升采樣（可選）?----
%?????????beishu=10;
%?????????sig=resample（sig.‘beishu1）;
%?????????sig=sig.‘;
%?????????fs=fs*beishu;
%?????????sig_len=sig_len*beishu;
????%%------------------------
????MaxTao=round（2*R/c*fs）;
????for?a=thita
%?????????a*180/pi
????????sig_buchang=zeros（Nsig_len+MaxTao）;
????????index=round（a/（dthita*pi/180）+1）;
????????Number_Array_InUse=find（flag（index:）==1）;
????????Ntao0=round（R*（1-cos（a-psi（Number_Array_InUse）））/c*fs）;
????????for?kk=1:length（Number_Array_InUse）
????????????sig_buchang（Number_Array_InUse（kk）MaxTao-Ntao0（kk）+1:MaxTao-Ntao0（kk）+sig_len）=sig（Number_Array_InUse（kk）:）;
????????end
????????sig_buchang=resample（sig_buchang.‘1beishu）;
????????sig_buchang=sig_buc

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2014-09-11?15:42??CircleBeamforming\
?????文件????????7420??2014-09-11?15:34??CircleBeamforming\Beamforming_CircleArray.m
?????文件????????3830??2014-02-28?16:41??CircleBeamforming\CreateSignal_CircleArray.m
?????文件????????2040??2014-09-11?15:44??CircleBeamforming\Simulation1.m