海杂波分析的matlab程序

针对海杂波各种噪声分析，滤波，给除去杂后真实信号的matlab程序。

function [WB SumOfAll fd FirstOrder]=AllOrder（f0VpQQNQ0Reso）
%海浪散射系数仿真程序 个参数下面有说明 可以计算单基地雷达和双基地雷达的海浪散射系数

DopVal=zeros（1401）;
% f0=15; Vp=15; Q=0*pi/180;QN=90*pi/180.0;Q0=0*pi/180.0;Reso=2000;
g=9.8;
A=4/（3*pi）;
s=4;
K0=2*pi*f0/300;  %//无线电波长
WB=sqrt（2*g*K0*cos（Q0））; % //Bragg 基频

%/*     参数声明                                          
% Vp---输入风速19.5（m/s）	
% f0---雷达操作频率
% K0---雷达波数
% K1---一阶海洋波数
% K2---二阶海洋波数
% Wd---多普勒频率
% WB---Bragg多普勒频率
% Q ---风向（角度）
% QN---散射椭圆法线矢量方向
% Q0---半双基地角
% Q1---一阶海洋波方向
% Q2---二阶海洋波方向	
% s ---扩展因子 */    

% 	  //计算二阶散射截面积
% 	  //-3WBSingular1=0;
Singular2=0;
Singular3=0;
Singular4=0;
Singular5=0;
Singular6=0;
Singular7=0;
Singular8=0;
Singular9=0;
Singular10=0;
for N=1:160
    Wd=（-5.0+0.025*（N-1））*WB; 
    if（Wd+sqrt（2.0）*WB>0） %// Wd>=-sqrt（2）*WB部分
        if（N==145） 
            Wd=-sqrt（2）*WB; %// 奇异值
            Singular1=N;
        end
        if（abs（Wd+（2^0.75）*sqrt（（（1-sin（Q0））^0.5）/cos（Q0））*WB）<=0.0125*WB）
            Wd=-（2^0.75）*sqrt（（（1-sin（Q0））^0.5）/cos（Q0））*WB;
            Singular2=N;
            %//电场耦合系数产生的奇异值
        end
        DopVal（N）=dif（QN-piQN11WdWBf0VpQQNQ0ResoK0）+dif（QNQN+pi11WdWBf0VpQQNQ0ResoK0）;
    else  %//  Wd<-sqrt（2）WB
        if（abs（Wd+（2^0.75）*sqrt（（（1-sin（Q0））^0.5）/cos（Q0））*WB）<=0.0125*WB）
            Wd=-（2^0.75）*sqrt（（（1-sin（Q0））^0.5）/cos（Q0））*WB;
            Singular3=N;
            %//电场耦合系数产生的奇异值
        end
        if（abs（Wd+（2^0.75）*sqrt（（（1+sin（Q0））^0.5）/cos（Q0））*WB）<=0.0125*WB）
            Wd=-（2^0.75）*sqrt（（（1+sin（Q0））^0.5）/cos（Q0））*WB;
            Singular4=N;
        end
        DopVal（N）=dif（QN-piQN-acos（2*（WB/Wd）*（WB/Wd））11WdWBf0VpQQNQ0ResoK0）+dif（QN+acos（2*（WB/Wd）*（WB/Wd））QN+pi11WdWBf0VpQQNQ0ResoK0）;
    end
end
% //计算-WBfor N=162:200
    Wd=（-5+0.025*（N-1））*WB;
    if（abs（Wd+（2^0.75）*sqrt（（（1-sin（Q0））^0.5）/cos（Q0））*WB）<=0.0125*WB）
        Wd=-（2^0.75）*sqrt（（（1-sin（Q0））^0.5）/cos（Q0））*WB; 
        Singular5=N;
    end
    DopVal（N）=dif（QN-piQN-11WdWBf0VpQQNQ0ResoK0）+dif（QNQN+pi-11WdWBf0VpQQNQ0ResoK0）;
end
N=N+1;
if（N==201）  % //Wd=0
    Wd=（-5+0.025*（N-1））*WB;  
    DopVal（N）=dif（QN-0.5*piQN+0.5*pi-11WdWBf0VpQQNQ0ResoK0）; 
end
for N=202:240
    Wd=（-5+0.025*（N-1））*WB;
    if（abs（Wd-（2^0.75）*sqrt（（（1-sin（Q0））^0.5）/cos（Q0））*WB）<=0.0125*WB）
        Wd=（2^0.75）*sqrt（（（1-sin（Q0））^0.5）/cos（Q0））*WB;
        Singular6=N;
        %//电场耦合系数产生的奇异值
    end
    DopVal（N）=dif（QN-piQN1-1WdWBf0VpQQNQ0ResoK0）+dif（QNQN+pi1-1WdWBf0VpQQNQ0ResoK0）;
    %// printf（“%d%e\n“NDopVal[N]）;
end
%//计算WBfor N=242:401
    Wd=（-5+0.025*（N-1））*WB;
    if（Wd-sqrt（2.0）*WB<0） %// Wd>=sqrt（2）*WB部分
        if（N==257）
            Wd=sqrt（2）*WB; %// 奇异值
            Singular7=N;
        end
        if（abs（Wd-（2^0.75）*sqrt（（（1-sin（Q0））^0.5）/cos（Q0））*WB）<=0.0125*WB

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      12318  2008-11-09 09:56  AllOrder.m

----------- ---------  ---------- -----  ----

                12318                    1