应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 关于Matlab中for循环的问题，求达人指点，谢谢了
51/1返回列表
查看: 1288  |  回复: 7
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

tao317594856

金虫 (小有名气)

[求助] 关于Matlab中for循环的问题，求达人指点，谢谢了已有1人参与

c1=[];
c2=[];
N=100000;
a1=-pi/2+2*(pi/2)*rand(N,1);
a2=-pi/3+2*(pi/3)*rand(N,1);
a4=-pi/2+2*(pi/2)*rand(N,1);
a5=-pi/3+2*(pi/3)*rand(N,1);
b=cos(a1+a4)./10 + cos(a1)./10 + (cos(a1 + a4).*cos(a5))./5 + (cos(a1).*cos(a2))./5;
c=sin(a1+a2)./10+sin(a1+a4)./10+sin(a1)./10+sin(a1-a2)./10+(sin(a1+ a4).*cos(a5))./5;
d=sin(a2)./5 + sin(a5)./5;
c1=[c1 d]; %将Z轴数值赋值到C1中
c2=[c2 b c]; %将x、y轴数值赋值到C2中
c11=max(c1);
c12=min(c1) %求出Z轴的最大值和最小值
% 定义Z轴分隔间隔的大小
syms derta1
% 定义Z轴分隔的数目
k1=3;
xl=abs(c11-c12);    %求解Z坐标最大值与最小值的差
derta1=xl/k1;
c3=[c1 c2];
c4=sortrows(c3,1);
for j=1:k1-1
    c5{j}=[];
for i=1:N
if c4(i,1)<=(c12+(xl/k1)*(j+1))&&c4(i,1)>(c12+(xl/k1)*(j))
c5{j}=[c5{j};c12+derta1*(j-1) c4(i,2) c4(i,3)]; %%%%%%%%%%% 不同点，改其横坐标为分段数值
end
end
end

以上是我写的程序，问题是从第一个for循环开始的。我想要在判断出来if 条件后，将满足条件的C4中的行矩阵赋值为[] 空集，这样就可以大大的减小循环时间，可是我没有写出来 ，写成如下的程序后
k3=0
for j=1:k1-1
    c5{j}=[];
    k2=1
   for i=k2:N-k3
if c4(i,1)<=(c12+(xl/k1)*(j+1))&&c4(i,1)>(c12+(xl/k1)*(j))
c5{j}=[c5{j};c12+derta1*(j+0.5) c4(i,2) c4(i,3)] %%%%%%%%%%% 不同点，改其横坐标为分段数值
c4(i,=[]
k2=i
k3=k3+1
else
    k2=k2+1
end
end
end
发现，for 循环后面跟的循环体在循环过程中是不能变化的，因此没办法再每次循环到符合条件的C4矩阵中的行抽离出来，导致将数据分多少层就需要进行N(数据点数)*K(分层数)次计算。
不知道有什么别的方法可以实现每循环一次，遇到符合条件的变量就抽离出来，在进行下一次循环的时候就可以扫描N-1个数据了，然后还能继续运行呢？求大神指教啊，
求大神指教！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2014-10-31 23:36:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tao317594856

金虫 (小有名气)
引用回帖:
3楼: Originally posted by zhswee at 2014-11-01 00:51:31
不清楚你要实现什么。

就是将三维散点数据按Z轴分段啊
赞一下 回复此楼
7楼2014-11-01 14:07:21
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

zhswee

铁杆木虫 (著名写手)
不清楚你要实现什么。

[ 发自小木虫客户端 ]
赞一下 回复此楼
3楼2014-11-01 00:51:31
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

johnliu1983

至尊木虫 (著名写手)

【答案】应助回帖

感谢参与，应助指数 +1
你可以把内层循环改成这样的试试,把for i=1:N,改成 for i=1:length(c4)
然后再删数据。这样每次循环到内循环的时候就会改变循环次数了。有什么不明白的再问我。
赞一下 回复此楼
4楼2014-11-01 06:14:24
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tao317594856

金虫 (小有名气)
引用回帖:
4楼: Originally posted by johnliu1983 at 2014-11-01 06:14:24
你可以把内层循环改成这样的试试,把for i=1:N,改成 for i=1:length(c4)
然后再删数据。这样每次循环到内循环的时候就会改变循环次数了。有什么不明白的再问我。

首先，我是想内循环的时候，每循环判定符合条件1次，c4就减少符合条件的相应的一行，而您所说的是在最外部每次取K值的时候，减少响应符合上一次k值条件的所有的行。
其次，即使把for i=1:N,改成 for i=1:length(c4)，将程序写成如下：
for j=1:k1-1
    c5{j}=[];
for i=1:length(c4)
if c4(i,1)<=(c12+(xl/k1)*(j+1))&&c4(i,1)>(c12+(xl/k1)*(j))
c5{j}=[c5{j};c12+derta1*(j+0.5) c4(i,2) c4(i,3)]; %%%%%%%%%%% 不同点，改其横坐标为分段数值
c4(i,=[];
end
end
end
也是无法执行的，因为for 后面的 i=1:length(c4)，已经确定i是从1执行到c4的初始长度了，在该循环下c4的长度是不能改变的

最后，还请您再帮忙思考思考，怎么能在内循环层里，每次遇到满足条件的就剔除，然后还能在循环体 里直接执行新C4的行数次数的循环呢，是不是有别的循环循环函数之类的可以满足要求呢，因为这样可以大量减小计算量，有很大优势的。
赞一下 回复此楼
5楼2014-11-01 11:19:30
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭