| 查看: 864 | 回复: 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 |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
-
浙江师范大学是怎么坑我的
已经有7人回复
-
求教各位大神：2026教育部人文社科青年基金项目何时公示呀？
已经有6人回复
-
面上有专家说实验设备不是我们单位的
已经有8人回复
-
科研求助
已经有4人回复
-
博士去二本高校当辅导员，如何调整心态？
已经有11人回复
-
江苏省自然基金 什么时候出结果
已经有6人回复
-
闲聊
已经有3人回复