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

%當權向量的方向角選擇得和發送信號都不符合時，方向圖零點的個數等于M-1；
%當權向量的方向角選擇得和發送信號之一符合時，另一個信號的強烈干擾置為0
%信號中不再含有強干擾
%保證輸出功率最小的同時，約束條件為W‘*a（theta0）=1）

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

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