作者ljuyentintho (小劉)
看板MATLAB
標題[問題]一個bug求解
時間Wed Jun 1 00:27:16 2016
如題
程式碼如下:
clc
clear all
R=2;
A=2;
J=3;
Fn=kron([1:R],ones(1,J*A))';
Sn=repmat(kron([1:A],ones(1,J))',R,1);
Tn=repmat([1:J]',R*A,1);
x=[repmat('X',R*A*J,1),num2str(Fn),num2str(Sn),num2str(Tn)];
X=mat2cell(x,ones(R*A*J,1),4);
randnum=8.*rand(R*A,1);
for i=1:R*A
temp=randperm(J);
X{(i-1)*J+temp(1),2}=randnum(i);
end
X(cellfun(@isempty,X))={0};
idx=find(cell2mat(X(:,2)));
Y1=cumsum(randnum(1:2:end));
Y2=cumsum(randnum(2:2:end));
y=[repmat('Y',R*A*J,1),num2str(Fn),num2str(Sn),num2str(Tn)];
Y=mat2cell(y,ones(R*A*J,1),4);
Y(idx(1:2:end),2)=mat2cell(Y1,ones(1,length(idx(1:2:end))));
Y(idx(2:2:end),2)=mat2cell(Y2,ones(1,length(idx(2:2:end))));
Y(cellfun(@isempty,Y))={0};
finalx=cell2mat(X(:,2));
finaly=cell2mat(Y(:,2));
if (finaly==0)
probabilityY=0;
else
probabilityY=0.125*exp(-(0.125*finaly));
end
finaly
probabilityY
這是其中一組輸出:
finaly =
7.8749
0
0
0
0
3.7263
9.0695
0
0
6.7079
0
0
probabilityY =
0.0467
0.1250
0.1250
0.1250
0.1250
0.0785
0.0402
0.1250
0.1250
0.0540
0.1250
0.1250
我希望當finalY這個矩陣的值為0時
probabilityY對應到的值也是0
不為0時
probabilityY=0.125*exp(-(0.125*finaly))
也寫了if跟else判斷式
但是為何finaly原本為0的地方還是被代入probabilityY=0.125*exp(-(0.125*finaly)) ?
謝謝好心的大大
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 140.114.54.170
※ 文章網址: https://webptt.com/m.aspx?n=bbs/MATLAB/M.1464712038.A.F59.html
1F：→ ejialan: if 後面接向量要全部都是1才會是true 所以逐一判斷才對 06/01 09:20
2F：→ ejialan: 或是把if end那五行改成 06/01 09:21
3F：→ ejialan: probabilityY=0.125*exp(-(0.125*finaly)); 06/01 09:21
4F：→ ejialan: probabilityY(~finaly)=0; 06/01 09:22
5F：→ ljuyentintho: 已解決 感謝 06/01 19:08