请教个避免用循环的编程技巧
程序：a=;
b=;
c=;
aa=10:10:70;
for i=1:length(aa)-1
j=find(a<=aa(i) & b>aa(i));
cc(i)=c(j);%由于实际应用中的特殊设置，此处的j只可能为一个数
end
cc(end+1)=c(end);
怎么不用循环实现这个过程，谢谢！
[ 本帖最后由 花如月 于 2007-12-13 22:35 编辑 ] maybe
cc= c(1)*(aa>=a(1)&aa<=b(1))...
+ c(2)*(aa>=a(2)&aa<=b(2))...
+ c(3)*(aa>=a(3)&aa<=b(3));
原帖由 hyl2323 于 2007-12-13 09:39 发表 http://www.chinavib.com/forum/images/common/back.gif
a=;
b=;
c=;
aa=a(1):10:b(end);
cc=如果aa(i)处于a(i)与b(i)之间，则cc(i)取c(i)的值。
怎么不用循环实现这个过程，谢谢！
aa有10个数据,a,b有3个,那aa后面那7个该怎么比啊?? 经过测试应该没有问题：“如果aa(i)处于a(i)与b(i)之间，则cc(i)取c(i)的值”边界没有取等于的情况，因此cc的第一个数总是为空，如果需要自己加上clear,clc
a=;
b=;
c=;
aa=a(1):10:b(end);
index=1:length(c);
cc=c((a(index)<aa(index))&(aa(index)<b(index)));没刷新已经有人回复了，我和2楼的思路一致。不过判断的时候可以完全矢量化，一个句子就够了。多余的7个就自动截断了，因为结果由c的长度决定。看这个例子：
>> a=;c=1:3;
>> c(a)
??? Subscript indices must either be real positive integers or logicals.
>> c(logical(a))
ans =
2 3
以前一直注意过那个错误提示。这个例子说明逻辑0可以做下标，得到的是空数相当于删除操作
[ 本帖最后由 eight 于 2007-12-13 11:02 编辑 ]
回复 #2 sugarcane 的帖子
a,b,c我只简单列了三个数，其实都是长度为n的行向量。所以你说的cc= c(1)*(aa>=a(1)&aa<=b(1))...
+ c(2)*(aa>=a(2)&aa<=b(2))...
+ c(3)*(aa>=a(3)&aa<=b(3));
其中的1 2 3，还是得用循环阿？ 结果是由aa的长度决定，对不起，我的意思没有表达清楚,问题重新修改编辑过了。
[ 本帖最后由 hyl2323 于 2007-12-13 13:38 编辑 ] 问题还没有解决，楼上所说的我看明白了，谢谢。问题关键是aa与a b c不同长度。
回复 #7 hyl2323 的帖子
又多了个j，问题似乎描述的有些问题，看不太懂。[ 本帖最后由 花如月 于 2007-12-13 13:53 编辑 ] 现在应该看得懂了，我琢磨了好久更正了描述了，我把实际应用说说吧，我是用来处理声音的3rd倍频程的，我想把每个3rd倍频程再细分，例如250Hz中心频率的倍频程范围为220-280Hz，声压级为L，我要将这个频带细分为220：10：280，那么各细分频带的声压级均为L-10lg((280-220)/10)。有30几个3rd倍频程。现在应该表述清楚了吧。可能不了解这个专业的难以明白物理意思。那就看主题帖吧。
[ 本帖最后由 hyl2323 于 2007-12-13 14:07 编辑 ] 我用带for循环可以实现如下：帮助理解。
for i=1:length(aa)-1
j=find(a<=aa(i) & b>aa(i));
cc(i)=c(j);%由于实际应用中的特殊设置，此处的j只可能为一个数
end
cc(end+1)=c(end);
如何用向量计算代替这个for？值得思考！---最后编辑已经改对了！
[ 本帖最后由 hyl2323 于 2007-12-13 17:56 编辑 ] 当然，我此处的for只要循环大概2000次，不会很耗资源的，只是学究一下，毕竟matlab最好少用循环，如何化循环为矩阵计算，确为不俗之技巧。高手们快支招，我也在冥思苦想。
回复 #10 hyl2323 的帖子
这个程序根本无法运行，赋值的矩阵维数都不一样???In an assignmentA(I) = B, the number of elements in B and
I must be the same. j是个向量，不可能放到cc(i)里
因此的表述还是有问题，你的a，b是用来确定区间起始，c应该是代表的分的区间段数。
cc的结果是aa里每个元素所在的区间索引吧？这样理解对不？
[ 本帖最后由 花如月 于 2007-12-13 15:47 编辑 ] 我身边没有matlab，没有调试过，可能是处理最后一个aa(end)返回的j为空矩阵[]的原因？我下班后再来弄。
[ 本帖最后由 hyl2323 于 2007-12-13 16:00 编辑 ]
回复 #12 花如月 的帖子
a，b是区间，c是这个区间的一个函数总值，不是段数，cc相当于是一个插值后的c。只是插值是按照那个不等式来操作的。回复 #14 hyl2323 的帖子
最好先给个常规的解法
页:
[1]
2