遍历.m

多基时差雷达定位算法。以雷达测量误差分布规律入手

for m=1:5   %因为从图形可看出5年后，药性减弱,不会是我们的优化结果
    [t,x]=ode45('shu4',[0,m],[7,1.7,0.001]);   %加药
    laoshu1(m)=x(length(t),2);              %m年后老鼠量
    cao1(m)=x(length(t),1);           %m年后草的量
    cz1=x(length(t),1)-7;           %m年后草的增加量
    q2=cz1*1*1000;        %m年后该方法带来的草的收益
    q1=m*1000;            %每年买药花的钱
    zongxiaohao1(m)=q1-q2;    %为了治理草原m年总共的花费，包含草达不到7吨时损耗的钱
    for n=1:15                 %m年后开始引入天敌
        [t,x]=ode45('shu2',[0,n],[7,laoshu1(m),0.023]);    %引入天敌
        laoshu2(m,n)=x(length(t),2);     %加药m年，又引入天敌n年后的老鼠量
        cao2(n)=x(length(t),1);       %加药m年，又引入天敌n年后的草量
        cz2=x(length(t),1)-7;                %草的增加量
        p2=cz2*1*1000;             %该方法带来的草的收益
        p1=2.3*1000;               %引入天敌的花费
        zongxiaohao2(m,n)=p1-p2;  %为了治理草原引入天敌总共的花费
    end
end
zongxiaohao1
zongxiaohao2
laoshu2