Changeset 781

Show
Ignore:
Timestamp:
01/12/10 23:09:15 (14 years ago)
Author:
smidl
Message:

example finished

Files:
1 modified

Legend:

Unmodified
Added
Removed
  • applications/bdmtoolbox/tutorial/mpdm/dist_ctrl_example.m

    r774 r781  
    2828C1.class = 'LQG_ARX'; 
    2929C1.ARX = A1; 
    30 C1.Qu = 0.1*eye(1); 
     30C1.Qu = 0.01*eye(1); 
    3131C1.Qy = 1*eye(1); 
    3232C1.yreq = 1; 
     
    6666 
    6767    loss_non_coop(i) = M.DS_y'*C1.Qy*M.DS_y + M.DS_u1'*C1.Qu*M.DS_u1 + M.DS_u2'*C1.Qu*M.DS_u2; 
     68    if loss_non_coop(i)>100 
     69        %keyboard 
     70    end 
    6871end 
    6972mean(loss_non_coop) 
     73 
     74loss_coop = zeros(1,Ntrials); 
     75for i=1:Ntrials 
     76    P1.neighbours = {'P2'}; 
     77    P2.neighbours = {'P1'}; 
     78    M= arena(DS,{P1,P2},exper); 
     79 
     80    loss_coop(i) = M.DS_y'*C1.Qy*M.DS_y + M.DS_u1'*C1.Qu*M.DS_u1 + M.DS_u2'*C1.Qu*M.DS_u2; 
     81    if loss_coop(i)>100 
     82        %keyboard 
     83    end 
     84end 
     85mean(loss_coop) 
     86