[讨论] 两个矩阵运算

楼主: JamesChen (James)   2016-07-30 09:07:23
1.
ex: a = [200 300 100];
b = [1 2 1];
要产生 c = [ones(1,200) 2*ones(1,300) ones(1,100)];
请问要怎么用非 loop 的方式快速产生?
a b 只是举例 实际上数字大很多,且要做很多次
2.
ex: a = [1 5 10 10 5 1];
b = [2 5];
要产生 c = [6 , 25 , 1];
就是 a 已经是 binned 过的资料了 b 是换一种方式去 bin a。
请问有没有不用 loop 的方式。
一样 a b 只是例子,实际资料大很多。
作者: celestialgod (天)   2016-07-30 09:12:00
1. arrayfun(@(a,b) b*ones(1,a), [200,300,100], [1,2,1],'uni',false)2看不懂
作者: YoursEver (猪是妳)   2016-07-30 09:22:00
楼上,他的2是巴斯卡三角形,b=[2,5]是指2项式展开的5次[6, 25, 1]看起来是把某几个幂次的系数加总,但规则不明
楼主: JamesChen (James)   2016-07-30 09:35:00
抱歉 2 没那么复杂 b 只是 rebin 的 edgea 是有 1个1 5个2 10个3 ... 的原始资料用 1 2 3 4 5 去 bin 过的 histogram我想要做更宽一点的 bin感谢一楼 arrayfunc 我太不熟了 惭愧
作者: wsglu   2016-07-30 10:09:00
2. temp=cumsum(a);c=temp(b) 不晓得是不是原PO要的
楼主: JamesChen (James)   2016-07-30 10:24:00
很接近了 再 diff 过后处理一下就是我要的了 感谢

Links booklink

Contact Us: admin [ a t ] ucptt.com