阵列信号处理对角加载算法matlab程序
《阵列信号处理对角加载算法matlab程序》由会员分享,可在线阅读,更多相关《阵列信号处理对角加载算法matlab程序(2页珍藏版)》请在装配图网上搜索。
1、%- 对角加载( LSMI 和 SMI)方向图 -%总结:这种算法主要给出了一种对角加载值的计算方法,对误差具有一定的稳健性,研究发现%当数据协方差矩阵中含有信号分量会影响算法的性能。clear all;clear all;clc;ratio_d_and_w=0.5;N_array=20;%阵列数N_signal=60;% 样本数ang1=0*pi/180;% 所需信号的方向SNR=5;%信噪比ASd=sqrt(10.(SNR/10);ang2=40*pi/180;% 干扰信号的方向INR=45;%干噪比ASi=sqrt(10.(INR/10);Sd=ASd*(randn(1,N_signal
2、)+i*randn(1,N_signal);%Sd 为所需信号 Si=ASi*(randn(1,N_signal)+i*randn(1,N_signal);%Si 为干扰信号Ni=randn(N_array,N_signal)+i*randn(N_array,N_signal);%Ni 内噪声Desired_Array=zeros(N_array,N_signal);Interferential_Array=zeros(N_array,N_signal);for LL=1:N_signalInterferential_Array(:,LL)=Si(LL)*test(ang2,N_array,r
3、atio_d_and_w).;Desired_Array(:,LL)=Sd(LL)*test(ang1,N_array,ratio_d_and_w).;endX=zeros(N_array,N_signal);X= Interferential_Array +Ni;Rx=X*X/N_signal;mm=std(diag(Rx);%对角加载值的确定下限%mm=trace(Rx)/N_array;% 对角加载值的确定上限R1=Rx+mm*eye(size(Rx);R=inv(R1);A_est=test(ang1,N_array,ratio_d_and_w);C= A_est;w_SMI=R*C/
4、(C*R*C);% 对角加载w_LSMI=inv(Rx)*C/(C*inv(Rx)*C);% 普通的 Capon 算法x_axis=-90:0.5:90;ang=(x_axis)*pi/180;for j=1:length(ang)steer=test(ang(j),N_array,ratio_d_and_w);f1(j)=w_SMI*steer;f2(j)=w_LSMI*steer;endf1=10*log10(abs(f1)/(max(max(abs(f1);f2=10*log10(abs(f2)/(max(max(abs(f2);subplot(2,1,1);plot(x_axis,f1
5、),grid on;hold on;plot(x_axis,f2, -),grid on;xlabel(theta/ );ylabel(F/dB);title(a);%*X=Desired_Array +Interferential_Array +Ni;Rx=X*X/N_signal;mm=std(diag(Rx);%对角加载值的确定下限%mm=trace(Rx)/N_array;% 对角加载值的确定上限R1=Rx+mm*eye(size(Rx);R=inv(R1);A_est=test(ang1,N_array,ratio_d_and_w);C= A_est;w_SMI=R*C/(C*R*C
6、);% 对角加载w_LSMI=inv(Rx)*C/(C*inv(Rx)*C);% 普通的 Capon 算法x_axis=-90:0.5:90;ang=(x_axis)*pi/180;for j=1:length(ang)steer=test(ang(j),N_array,ratio_d_and_w);f1(j)=w_SMI*steer;f2(j)=w_LSMI*steer;endf1=10*log10(abs(f1)/(max(max(abs(f1);f2=10*log10(abs(f2)/(max(max(abs(f2);subplot(2,1,2);plot(x_axis,f1),grid on;hold on;plot(x_axis,f2, -),grid on;xlabel(theta/ );ylabel(F/dB);title(b);
- 温馨提示:
1: 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
2: 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
3.本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
5. 装配图网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。