24小时热门版块排行榜

返回列表

zhaomin80

铜虫 (小有名气)

应助: 4 (幼儿园)
金币: 175.6
红花: 2
帖子: 143
在线: 53.9小时
虫号: 1185832
注册: 2011-01-08
专业: 地理信息系统

[求助] 程序运行占用内存较高，造成机子死机如何简化程序？

syms b_s d_hh c1 d_vv k l
THETA=geotiffread('theta1111.tif');
HH=geotiffread('hh1111.tif');
VV=geotiffread('vv1111.tif');
p=nan(size(THETA));
i1=inversion1();
i2=d_hh-i1;
i3=[c1,d_hh,d_vv];
for n=1:size(THETA,1)
for o=1:size(THETA,2)
      theta=THETA(n,o);
      hh=HH(n,o);
      vv=VV(n,o);
      i4=[theta,hh,vv];
      i=subs(i2,i3,i4);
      j1=solve(i(1,1),b_s);
      j2=solve(i(2,1),b_s);
      j=[real(j1);real(j2)];

   if ~isnan(theta) && ~isnan(vv) && ~isnan(hh)
      if ~isempty(j1) && ~isempty(j2)
      l1=solve(j(1,1)==1+1.5995*k+1.861*k^3,k);
      l2=solve(j(2,1)==1+1.5995*k+1.861*k^3,k);
      l=[l1;l2];
      elseif ~isempty(j1) && isempty(j2)
      l=solve(j==1+1.5995*k+1.861*k^3,k);
      elseif isempty(j1) && ~isempty(j2)
      l=solve(j==1+1.5995*k+1.861*k^3,k);
      end
      for m=1:numel(l)
      if isreal(l(m,1)) && l(m,1)>=0.1 && l(m,1)<0.45
      p(n,o)=l(m,1);
      disp(num2str([n,o]));
      end
      end
   end
end
end

回复此楼

» 猜你喜欢

1楼 2014-09-22 23:38:23

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 zhaomin80 的主题更新

返回列表