| 查看: 854 | 回复: 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 |
回复此楼» 猜你喜欢
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有11人回复
-
重要来源：本周末出结果
已经有5人回复
-
FileCode能看出啥？
已经有23人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有36人回复
-
关于Filecode分析方法
已经有9人回复
-
不应该看fileCode
已经有6人回复
-
好奇怪的filecode
已经有6人回复
-
长年满屏的广告，版主太不责任了。基金也等的急
已经有4人回复
-
filecode
已经有15人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有4人回复
