| 查看: 825 | 回复: 0 | ||
[求助]
程序运行占用内存较高，造成机子死机如何简化程序？
|
|
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 |
回复此楼» 猜你喜欢
假如你的研究生提出不合理要求
已经有12人回复
-
实验室接单子
已经有7人回复
-
全日制（定向）博士
已经有5人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
对氯苯硼酸纯化
已经有3人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
-
北核录用
已经有3人回复