如何選擇在MatLab中碰到矩形的隨機線

我正在創建隨機的起點和終點。我想繪製與放置在原點的矩形相交/相交的那些人。我發現我的代碼錯過了一些行，如圖所示。在那之後，我想要數一下軌道撞上矩形。對於expample軌道從正面來了，從右側等如何選擇在MatLab中碰到矩形的隨機線

退出

我的代碼是

function [hot, cold, h] = MuonTracks(tracks)%#eml 

    % NOTE: no variables larger than 1x1 are initialized 

    width = 1e-4; 
    height = 2e-4; 

    % constant used for Laplacian noise distribution 
    bL = 15/sqrt(2); 

    % Loop through all tracks 
    X = []; 
    bhit= []; 
    hot = 0; 
    ii = 0; 
    TopBottom= 0; 
    LeftRight= 0; 
    TopRight= 0; 
    TopLeft= 0; 
    LeftBottom= 0; 
    RightBottom= 0; 
    ihit= 0; 
    while ii <= tracks 

     ii = ii + 1; 

     % Note that I've inlined (== copy-pasted) the original laprnd() 
     % function call. This was necessary to work around limitations 
     % in loops in Matlab, and prevent the nececessity of those HUGE 
     % variables. 
     % 
     % Of course, you can still easily generalize all of this: 

     % the new data 
     u = rand-0.5; 

     Ystart = 15; 
     Xstart = 80*rand-40; 
     Xend = Xstart - bL*sign(u)*log(1-2*abs(u)); 
     %Xend=laprnd(tracks,1,Xstart,15); 

     b = (Ystart*Xend)/(Xend-Xstart); 


     % the test 
     if ((b < height && b > 0)) ||... 
      (Xend < width/2 && Xend > -width/2) 

      hot = hot+1; 

      % growing an array is perfectly fine when the chances of it 
      % happening are so slim 
      X = [X [Xstart; Xend]]; %#ok 
      bhit=b; 

     end 
    end 

    % This is trivial to do here, and prevents an 'else' 
    cold = tracks - hot; 

    % Now plot the chosen ones 
    h = figure; 
    hold all  
    %Y = bsxfun(@times, 15, ones(size(X))); 
    if (size(X)==0) 
     %Disp('No hits were found'); 
     msgbox('No tracks were found','Result','warn'); 
    elseif (size(X)>1) 
     Y = bsxfun(@times, [15; 0], ones(size(X))); 
     plot(X, Y, 'r'); 
     msgbox([num2str(hot) ' tracks were found'],'Result','help',num2str(hot)); 
    else 
     Y = bsxfun(@times, [15; 0], ones(size(X))); 
     plot(X, Y, 'r'); 
     msgbox([num2str(hot) ' track was found'],'Result','help',num2str(hot)); 
    end 
    %X; 
    %Y; 
    %size(X,2) 
    while ihit<size(X,2) 

     ihit=ihit+1 
     %X(2,ihit) 

     if ((X(2,ihit)==-width && (bhit<=0 && bhit<=height))&&(bhit==0 && (X(2,ihit)>=-width && X(2,ihit)>=width))) 
      LeftBottom=LeftBottom+1; 
     elseif ((bhit==height && (X(2,ihit)>=-width && X(2,ihit)>=width)) && (X(2,ihit)==width && (bhit<=0 && bhit<=height))) 
      TopRight=TopRight+1; 
     elseif ((bhit==height && (X(2,ihit)>=-width && X(2,ihit)>=width)) && (bhit==0 && (X(2,ihit)>=-width && X(2,ihit)>=width))) 
      TopBottom=TopBottom+1; 
     elseif ((X(2,ihit)==-width && (bhit<=0 && bhit<=height)) && (X(2,ihit)==width && (bhit<=0 && bhit<=height))) 
      LeftRight=LeftRight+1; 
     elseif ((X(2,ihit)==-width && (bhit<=0 && bhit<=height)) && (bhit==height && (X(2,ihit)>=-width && X(2,ihit)>=width))) 
      TopLeft=TopLeft+1; 
     elseif ((X(2,ihit)==width && (bhit<=0 && bhit<=height)) && (bhit==0 && (X(2,ihit)>=-width && X(2,ihit)>=width))) 
      RightBottom=RightBottom+1; 
     else 
      display('sth is wrong'); 
     end 
     X(2,ihit) 
    end 
    %X(1,:); 
    %Y(1,:); 
    LeftBottom 
    TopRight 
    TopBottom 
    LeftRight 
    TopLeft 
    RightBottom 
    %display('sdfghjk'); 
end

任何想法會更受歡迎！

來源

2013-04-14 Thanos

那麼你想隨機生成一條與你的矩形相交的線，或者只是從給定的一組線中選擇一條呢？ –

@EitanT：非常感謝您的評論！我已經產生了起點和終點。所以我想要的是繪製與rectange相交的線，並計算每個交點在哪一側出現的次數（即：總共10個交點，從上到下2個，從左到右3個，從上到下4個左邊，右邊1，右邊0，左邊0，右上邊0）。 – Thanos

Here你有一個功能，能夠相交整個組的段並返回交點。我希望它有幫助。

來源

2013-04-14 15:45:34

如何選擇在MatLab中碰到矩形的隨機線

回答

相關問題