[討論] 兩個矩陣運算

[JamesChen (業餘經濟學家)]

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 只是例子，實際資料大很多。

--

兩個我都 google 過，前者可能因為關鍵字不對找不太到。
後者搜尋 rebin 找到的方法也都還是 loop。
目前我暫時都是用 loop 在做的，想問問有沒有高手有解的。
感謝

--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 98.30.193.49
※ 文章網址: https://www.ptt.cc/bbs/MATLAB/M.1469840846.A.171.html
→ celestialgod: 1. arrayfun(@(a,b) b*ones(1,a), [200,300,100], [
→ celestialgod: 1,2,1],'uni',false)
→ celestialgod: 2看不懂
→ YoursEver: 樓上,他的2是巴斯卡三角形,b=[2,5]是指2項式展開的5次
→ YoursEver: [6, 25, 1]看起來是把某幾個冪次的係數加總,但規則不明
→ JamesChen: 抱歉  2 沒那麼複雜 b 只是 rebin 的 edge
→ JamesChen: a 是有 1個1  5個2  10個3 ... 的原始資料
→ JamesChen: 用 1 2 3 4 5 去 bin 過的 histogram
→ JamesChen: 我想要做更寬一點的 bin
→ JamesChen: 感謝一樓 arrayfunc 我太不熟了 慚愧
→ wsglu: 2. temp=cumsum(a);c=temp(b)  不曉得是不是原PO要的
→ JamesChen: 很接近了 再 diff 過後處理一下就是我要的了 感謝