改变SNR的music算法

clear all;
close all;
clc;
j=sqrt(-1);
seita=[3,2,6]’;
N=size(seita,1); %求出信号个数
M=11; %阵子个数
kp=1000;
lmda=0.24; %信号波长
d=0.12; %阵子的间距
jthresh=1000; %给出一个比较的较大值
jdbeg=1; %采样的第一个点
jdend=180; %采样的最后一个点
jdnum=180; % 采样的数目
a9=40; %信噪比
a9=10^(a9/20); %a9=100
for m=1:M
for n=1:N
A(m,n)=exp(j2pi*(m-1)*(d/lmda)*sin(seita(n)*pi/180));
end
end

% a1=pi/180;
% for n=1:N
% a2=jpisin(seita(n)a1); %导向矩阵赋值
% for m=1:M
% A(m,n)=exp(a2(m-1));
% end
% end
a_jd=zeros(1,jdnum);
drt=(jdend-jdbeg)/(jdnum-1); %步进
t1=[0:0.005:25]; %最大快拍5000个
out1=exp(j2pi10t1);
out2=exp(j2pi20t1); %信号模型
out3=exp(j2pi40t1);
cd1=10;