lagelangri.rar_lagrange multiplier_saddlefx4_乘子法约束_拉格朗日_拉格朗日乘数

%函数名netplot %使用方法输入请help netplot %无返回值 %函数只能处理有向图 %作者：tiandsp %最后修改：2012.12.26 function netplot(A,flag) %调用方法输入netplot(A,flag)，无返回值 %A为邻接矩阵或关联矩阵 %flag=1时处理邻接矩阵 %flag=2时处理关联矩阵 %函数只能处理有向图 if flag==1 %邻接矩阵表示有向图 D_netplot(A); return; end if flag==2 %关联矩阵表示有向图 [m n]=size(A); %关联矩阵变邻接矩阵 W=zeros(m,m); for i=1:n a=find(A(:,i)~=0); if A(a(1),i)==1 W(a(1),a(2))=1; else W(a(2),a(1))=1; end end D_netplot(W); return; end function D_netplot(A) [n n]=size(A); w=floor(sqrt(n)); h=floor(n/w); x=[]; y=[]; for i=1:h %使产生的随机点有其范围，使显示分布的更广 for j=1:w x=[x 10*rand(1)+(j-1)*10]; y=[y 10*rand(1)+(i-1)*10]; end end ed=n-h*w; for i=1:ed x=[x 10*rand(1)+(i-1)*10]; y=[y 10*rand(1)+h*10]; end plot(x,y,'r*'); title('网络拓扑图'); for i=1:n for j=1:n if A(i,j)~=0 c=num2str(A(i,j)); %将A中的权值转化为字符型 text((x(i)+x(j))/2,(y(i)+y(j))/2,c,'Fontsize',10); %显示边的权值 arrow([x(i) y(i)],[x(j) y(j)]); %带箭头的连线 end text(x(i),y(i),num2str(i),'Fontsize',14,'color','r'); %显示点的序号 hold on; end end end end