[問題] 請問如何將FOR迴圈改寫成矩陣運算已回收
小弟實力不夠深
厚請高手開示
程式碼如下
非常感謝
如果無法改寫 麻煩請推文一下
M1、U、N0皆為1x1的正整數，x為1x100的正整數，L1~L4為1xM1-1的正整數
最主要是計算出 D(a,b,c,d) 迴圈內的值
但是 D(a,b,c,d)矩陣牽連到 s4(d)、C(a,b,c)與 B(a,b)與 A(a)內的值
for i=1:100
for a=1:1:M1-1
s1(a)=1-exp(-N0.*x(i)./L1(a));
A(a)=binopdf(L1(a),M1-1,1/U)*s1(a);
for b=1:1:M1-1-L1(a)
s2(b)=1-exp(-N0.*x(i)./L2(b));
B(a,b)=binopdf(L2(b),M1-1-L1(a),1/(U-1))*s2(b)*A(a);
for c=1:1:M1-1-L1(a)-L2(b)
s3(c)=1-exp(-N0.*x(i)./L3(c));
C(a,b,c)=binopdf(L3(c),M1-1-L1(a)-L2(b),1/(U-2))*s3(c)*B(a,b);
for d=1:1:M1-1-L1(a)-L2(b)-L3(c)
sa(d)=(2*N0./L4(d)).^0.5;
sb(d)=(2*N0.*x(i)./L4(d)).^0.5;
s4(d)=1-marcumq(sa(d),sb(d),1);
D(a,b,c,d)=binopdf(L4(d),M1-1-L1(a)-L2(b)-L3(c),1/(U-3))*s4(d)*C(a,b,c);
end
end
end
end
R(i)=log2(U)/(U^3)*sum(sum(sum(sum(D))));
end
抱歉先前有些變數沒說明很清楚
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 122.116.245.77
→
06/23 22:13, , 1F
06/23 22:13, 1F
→
06/23 22:13, , 2F
06/23 22:13, 2F
→
06/23 22:15, , 3F
06/23 22:15, 3F
→
06/23 22:48, , 4F
06/23 22:48, 4F
→
06/23 23:11, , 5F
06/23 23:11, 5F
→
06/23 23:11, , 6F
06/23 23:11, 6F
→
06/23 23:27, , 7F
06/23 23:27, 7F
→
06/23 23:33, , 8F
06/23 23:33, 8F
→
06/23 23:40, , 9F
06/23 23:40, 9F
※ 編輯: hopebear 來自: 122.116.245.77 (06/24 00:21)
推
06/24 01:55, , 10F
06/24 01:55, 10F
→
06/24 01:55, , 11F
06/24 01:55, 11F
→
06/24 01:56, , 12F
06/24 01:56, 12F
→
06/24 01:57, , 13F
06/24 01:57, 13F
→
06/24 02:32, , 14F
06/24 02:32, 14F
→
06/24 02:34, , 15F
06/24 02:34, 15F
→
06/24 02:36, , 16F
06/24 02:36, 16F
推
06/24 03:07, , 17F
06/24 03:07, 17F
推
06/24 03:10, , 18F
06/24 03:10, 18F
→
06/24 03:11, , 19F
06/24 03:11, 19F
→
06/24 03:15, , 20F
06/24 03:15, 20F
→
06/24 03:18, , 21F
06/24 03:18, 21F
→
06/24 03:23, , 22F
06/24 03:23, 22F
→
06/24 03:25, , 23F
06/24 03:25, 23F
推
06/24 12:46, , 24F
06/24 12:46, 24F
討論串 (同標題文章)
完整討論串 (本文為第 1 之 3 篇)：
MATLAB 近期熱門文章
PTT數位生活區 即時熱門文章
