%capon波束形成器，也稱MVDR波束形成器。
%它試圖使噪聲以及來自非θ方向的任何干擾所貢獻(xiàn)的功率為最小但又能保持在觀測(cè)方向θ上的信號(hào)功率不變

%當(dāng)權(quán)向量的方向角選擇得和發(fā)送信號(hào)都不符合時(shí)，方向圖零點(diǎn)的個(gè)數(shù)等于M-1；
%當(dāng)權(quán)向量的方向角選擇得和發(fā)送信號(hào)之一符合時(shí)，另一個(gè)信號(hào)的強(qiáng)烈干擾置為0
%信號(hào)中不再含有強(qiáng)干擾
%保證輸出功率最小的同時(shí)，約束條件為W‘*a（theta0）=1）

close?all;????%關(guān)閉所有正在運(yùn)行的窗口
clear?all;????%清空緩存
clc;??????????%清屏?命令窗口
M=8;??????????%假設(shè)空間?陣列的?陣元數(shù)目?M=8
%X=[];?????????%陣元接收信號(hào)矢量為x（t）。某一固定時(shí)刻，接收信號(hào)是一個(gè)列向量X=[];
theta=[-20?60?80];???????%設(shè)置仿真所用的信號(hào)?入射角度，單位是度。仿真中，構(gòu)造信號(hào)時(shí)作為已知量，估計(jì)的時(shí)候是不知道的
theta=theta.*pi/180;?????%對(duì)上述單位為?度?的交度進(jìn)行轉(zhuǎn)換，化為?弧度?形式
d_=0.5;???????????????????%因?yàn)榫嚯xd一般設(shè)置為波長(zhǎng)λ的一半，所以此處直接令：d/λ=0.5
W=2*pi*d_.*sin（theta）;????%?導(dǎo)向矢量中，三個(gè)信號(hào)源的?空間相位

%%%%??????構(gòu)造發(fā)送信號(hào)??
N=1000;????????????????????%假設(shè)快拍數(shù)??N=1000
n=1:N;?????????????????????%?n從1開始，以步長(zhǎng)1?增長(zhǎng)到?1000
s1=（10）*cos（2*pi*0.01*n）;??%不考慮載波，此處設(shè)置三個(gè)信號(hào)s1s2s3?為幅度，頻率都不相同的信號(hào)，
???????????????????????????%因?yàn)樾盘?hào)頻率不同，因此一定是非相干信號(hào)，其中的?n?代表該采樣信號(hào)的采樣時(shí)刻
???????????????????????????%?n?從1?到1000?，每個(gè)信號(hào)Si?分別采樣1000次，即每個(gè)信號(hào)有1000個(gè)樣本點(diǎn)
s2=（17）*cos（2*pi*0.15*n）;
s3=（28）*cos（2*pi*0.35*n）;
noise=wgn（MN0）;??????????%產(chǎn)生一個(gè)M*N?的高斯白噪聲矩陣，噪聲能量為?0dB。因?yàn)?，假設(shè)已知噪聲功率為1所以：10*log1=0
???????????????????????????%?注意：每個(gè)采樣點(diǎn)在每個(gè)天線陣元上都要加入高斯白噪聲
%%%%%%%????接收信號(hào)X?
?????????????????????????????????%?注意，此處求接收信號(hào)矩陣X，是先求該矩陣中第m個(gè)陣元，q時(shí)刻的接收信號(hào)向量X（mq）
?for?m=1:M???????????????????????%?大循環(huán)為?天線陣元M
?????for?q=1:N???????????????????%?小循環(huán)?為每個(gè)天線陣元上接收到的1000個(gè)采樣信號(hào)點(diǎn)的值，q從1?到1000，默認(rèn)步長(zhǎng)為1
?????????Y=[s1（q）?s2（q）?s3（q）];??%?構(gòu)造