function?newEdges=addrandedges（numbAddedgesN）
%Function?newEdges=addrandedges（numbAddedgesN）?adds?numbAdd?random?
%edges?to?the?current?edge?set.??Generally?used?to?“small-worldify“?a?
%graph?by?adding?some?long?range?connections?（i.e.?drastically?decrease?
%the?graph?diameter）
%
%Inputs:????numbAdd?-?Number?of?additional?edges?to?add
%???????????edges?-?An?Mx2?vector?indicating?the?edges?of?the?graph
%???????????N?-?Optional?number?of?nodes?in?the?graph?if?different?than
%???????????????N?=?max（max（edges））?（i.e.?there?are?isolated?nodes）
%
%Outputs:???newEdges?-?An?（M?+?numbAdd）x2?vector?of?edges?for?the?new?
%???????????????graph
%
%
%Note:?If?duplicate?edges?are?generated?randomly?the?duplicates?are
%???ignored?-?Therefore?the?actual?number?of?edges?added?may?be?less?than
%???numbAdd?（e.g.?if?the?input?graph?were?fully?connected）
%
%
%5/22/03?-?Leo?Grady

%?Copyright?（C）?2002?2003?Leo?Grady?
%???Computer?Vision?and?Computational?Neuroscience?Lab
%???Department?of?Cognitive?and?Neural?Systems
%???Boston?University
%???Boston?MA??02215
%
%?This?program?is?free?software;?you?can?redistribute?it?and/or
%?modify?it?under?the?terms?of?the?GNU?General?Public?License
%?as?published?by?the?Free?Software?Foundation;?either?version?2
%?of?the?License?or?（at?your?option）?any?later?version.
%
%?This?program?is?distributed?in?the?hope?that?it?will?be?useful
%?but?WITHOUT?ANY?WARRANTY;?without?even?the?implied?warranty?of
%?MERCHANTABILITY?or?FITNESS?FOR?A?PARTICULAR?PURPOSE.??See?the
%?GNU?General?Public?License?for?more?details.
%
%?You?should?have?received?a?copy?of?the?GNU?General?Public?License
%?along?with?this?program;?if?not?write?to?the?Free?Software
%?Foundation?Inc.?59?Temple?Place?-?Suite?330?Boston?MA??02111-1307?USA.
%
%?Date?-?$Id:?addrandedges.mv?1.2?2003/08/21?17:29:29?lgrady?Exp?$
%========================================================================%

%Determine?number?of?points
if?nargin?????N=max（max（edges））;
end

%Randomly?generate?new?edges
tmpEdges=round（（N-1）*rand（numbAdd2）+1）;

%Remove?duplicates
W=adjacency（[edges;tmpEdges]ones（length（edges）+numbAdd1）N）;
newEdges=adjtoedges（W）;