ISARMATLAB仿真程序

使用Matlab进行雷达信号处理系统仿真，能迅速建立起系统模型，设计理念可以在任何细节上得到体现，建模时间短，模型简单、清晰，计算精度高，同时在系统设计的任何阶段都能够很方便地修改模型、评估结果和验证系统行为。Matlab不仅是一个在科学研究及各类工程设计中便于使用的计算工具

function [wused] = BCS_fast（PHItsigma2eta）
[NM] = size（PHI）;
PHIt = PHI‘*t;
PHI2 = sum（PHI.^2）‘;
ratio = （PHIt.^2）./PHI2;
[maxrindex] = max（ratio）;
alpha = PHI2（index）/（maxr-sigma2）;
% compute initial mu Sig S Q
phi = PHI（:index）;
Hessian = alpha + phi‘*phi/sigma2;
Sig = 1/Hessian;
mu = Sig*PHIt（index）/sigma2;
left = PHI‘*phi/sigma2;
S = PHI2/sigma2-Sig*left.^2;
Q = PHIt/sigma2-Sig*PHIt（index）/sigma2*left;
%
for count = 1:1000

    s = S; q = Q;
    s（index） = alpha.*S（index）./（alpha-S（index））;
    q（index） = alpha.*Q（index）./（alpha-S（index））;
    theta = q.^2-s;

    % choice the next alpha that maximizes marginal likelihood
    ml = -inf*ones（1M）;
    ig0 = find（theta>0）;%寻找theta>0的索引
    % index for re-estimate
    [irefoowhich] = intersect（ig0index）;
    if ~isempty（ire）
        Alpha = s（ire）.^2./theta（ire）;
        delta = （alpha（which）-Alpha）./（Alpha.*alpha（which））;
        ml（ire） = Q（ire）.^2.*delta./（S（ire）.*delta+1）-log（1+S（ire）.*delta）;
    end
    % index for adding
    iad = setdiff（ig0ire）;%寻找theta>0但未被选中的索引
    if ~isempty（iad）
        ml（iad） = （Q（iad）.^2-S（iad））./S（iad）+log（S（iad）./（Q（iad）.^2））;
    end
    is0 = setdiff（[1:M]ig0）;%寻找theta<0的索引
    % index for deleting
    [idefoowhich] = intersect（is0index）;
    if ~isempty（ide）
        ml（ide） = Q（ide）.^2./（S（ide）-alpha（which））-log（1-S（ide）./alpha（which））;
    end

    [ML（count）idx] = max（ml）;%最大化的位置
    % check if terminates?
    if count > 2 & abs（ML（count）-ML（count-1）） < abs（ML（count）-ML（1））*eta
        break;
    end

    % update alphas
    which = find（index==idx）;
    if theta（idx） > 0
        if ~isempty（which） % re-estimate 不是空集则重估
            Alpha = s（idx）^2/theta（idx）;
            Sigii = Sig（whichwhich）; mui = mu（which）; Sigi = Sig（:which）;
            delta = Alpha-alpha（which）;
            ki = delta/（1+Sigii*delta）;
            mu = mu-ki*mui*Sigi;
            Sig = Sig-ki*Sigi*Sigi‘;
            comm = PHI‘*（phi*Sigi）/sigma2;
            S = S + ki*comm.^2;
            Q = Q + ki*mui*comm;
            %
            alpha（which） = Alpha;
        else % adding
            Alpha = s（idx）^2/theta（idx）;
            phii = PHI（:idx）; Sigii = 1/（Alpha+S（idx））; mui = Sigii*Q（idx）;
            comm1 = Sig*（phi‘*phii）/sigma2;
            ei = phii-phi*comm1;
            off = -Sigii*comm1;
            Sig = [Sig+Sigii*comm1*comm1‘ off; off‘ Sigii];
            mu = [mu-mui*comm1; mui];
            comm2 = PHI‘*ei/sigma2;
            S = S - Sigii*comm2.^2;
            Q = Q - mui*comm2;
            %
            index = [index;idx];
            alpha = [alpha;Alpha];
            phi = [phiphii];
        end
    else
        if ~isempty（which） % deleting
            Sigii = Sig（whichwhich）; mui = mu（which）; Sigi = Sig（:which）;
            Sig = Sig-Sigi*Sigi‘/Sigii; Sig（:which） = []; Sig（which:） = [];
            mu  = mu-mui/Sigii*Sigi; mu（which） = [];
  

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

     文件       5479  2009-12-12 13:22  米波一维两球高频\Algorithms\BCS_fast_rvm.m

     文件       1133  2009-12-12 13:22  米波一维两球高频\Algorithms\CoSaMP.m

     文件      10040  2017-03-11 16:50  米波一维两球高频\Algorithms\ExCoV.m

     文件       8090  2009-12-12 13:22  米波一维两球高频\Algorithms\ExCoVapp.m

     文件       7322  2009-12-29 09:36  米波一维两球高频\Algorithms\ExCoVls.m

     文件      21647  2009-12-12 13:22  米波一维两球高频\Algorithms\GPSR_BB.m

     文件        666  2018-05-08 16:56  米波一维两球高频\Algorithms\gui_fig51.m

     文件      14446  2009-12-12 13:54  米波一维两球高频\Algorithms\hard_l0_Mterm_mod.m

     文件       2763  2009-12-12 13:22  米波一维两球高频\Algorithms\l1qc_logbarrier.m

     文件       4179  2009-12-12 13:22  米波一维两球高频\Algorithms\l1qc_newton.m

     文件       1597  2009-12-13 10:55  米波一维两球高频\Algorithms\SBL.m

     文件          0  2018-05-08 16:59  米波一维两球高频\Algorithms\Untitled6.m

     文件       3392  2017-06-30 11:13  米波一维两球高频\BCS_fast.m

     文件        676  2017-11-20 11:16  米波一维两球高频\cx1.m

     文件        664  2018-01-20 19:47  米波一维两球高频\fft_1d.m

     文件        538  2018-01-20 10:33  米波一维两球高频\fft_isar.m

     文件      36057  2018-01-20 09:51  米波一维两球高频\mibo1d_qiu_gaopin_imag.txt

     文件      36057  2018-01-20 09:51  米波一维两球高频\mibo1d_qiu_gaopin_real.txt

     文件       2646  2018-01-20 16:42  米波一维两球高频\OneD.m

     文件        490  2017-12-27 09:56  米波一维两球高频\Solve_OMP.m

     文件       4436  2018-05-04 08:57  米波一维两球高频\TwoD.m

     目录          0  2018-05-08 16:57  米波一维两球高频\Algorithms

     目录          0  2018-01-20 19:47  米波一维两球高频

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

               162318                    23