小世界网络简介及及MATLAB建模

小世界网络MATLAB建模

1.简介

小世界网络存在于数学、物理学和社会学中，是一种数学图的模型。在这种图中大部份的结点不与彼此邻接，但大部份结点可以通过任一其它节点经少数几步就可以产生联系。若将一个小世界网络中的点代表一个人，而联机代表人与人之间是相互认识的，则这小世界网络可以反映陌生人通过彼此共同认识的人而起来产生联系关系的小世界现象。

在日常生活中，有时你会发现，某些你觉得与你隔得很“遥远”的人，其实与你“很近”。小世界网络就是对这种现象的数学描述。用数学中图论的语言来说，小世界网络就是一个由大量顶点构成的图，其中任意两点之间的平均路径长度比顶点数量小得多。除了社会人际网络以外，小世界网络的例子在生物学、物理学、计算机科学等领域也有出现。许多经验中的图可以用小世界网络来作为模型。因特网、公路交通网、神经网络都呈现小世界网络的特征。小世界网络最早是由邓肯·瓦茨（Duncan Watts）和斯蒂文·斯特罗加茨（Steven Strogatz）在1998年引进的，将高聚合系数和低平均路径长度作为特征，提出了一种新的网络模型，一般就称作瓦茨-斯特罗加茨模型（WS模型），这也是最典型的小世界网络的模型。

由于WS小世界模型构造算法中的随机化过程有可能破坏网络的连通性，纽曼(Newman)和瓦茨(Watts)提出了NW小世界网络模型，该模型是通过用“随机化加边”模式来取代WS小世界网络模型构造中的“随机化重连”。

在考虑网络特征的时候，使用两个特征来衡量网络：特征路径长度和聚合系数。

特征路径长度（characteristic path length）：在网络中，任选两个节点，连同这两个节点的最少边数，定义为这两个节点的路径长度，网络中所有节点对的路径长度的平均值，定义为网络的特征路径长度。这是网络的全局特征。

聚合系数(clustering coefficient)：假设某个节点有k个边，则这k条边连接的节点之间最多可能存在的边的个数为k(k-1)/2，用实际存在的边数除以最多可能存在的边数得到的分数值，定义为这个节点的聚合系数。所有节点的聚合系数的均值定义为网络的聚合系数。聚合系数是网络的局部特征，反映了相邻两个人之间朋友圈子的重合度，即该节点的朋友之间也是朋友的程度。

我们可以发现规则网络具有很高的聚合系数，大世界（large world，意思是特征路径长度很大），其特征路径长度随着n(网络中节点的数量)线性增长，而随机网络聚合系数很小，小世界（small world，意思是特征路径长度小），其特征路径长度随着log(n)增长中说明，在从规则网络向随机网络转换的过程中，实际上特征路径长度和聚合系数都会下降，到变成随机网络的时候，减少到最少。但这并不是说大的聚合系数一定伴随着大的路径长度，而小的路径长度伴随着小的聚合系数，小世界网络就具有大的聚合系数，而特征路径长度很小。试验表明，少量的short cut的建立能够迅速减少特征路径长度，而聚合系数变化却不大，因为某一个short cut的建立，不仅影响到所连接的节点的特征路径长度，而且影响到他们邻居的路径长度，而对整个网络的聚合系数影响不大。这样，少量的short cut的建立就能使整个网络不知不觉地变成小世界网络。

实际的社会、生态、等网络都是小世界网络，在这样的系统里，信息传递速度快，并且少量改变几个连接，就可以剧烈地改变网络的性能，如对已存在的网络进行调整，如蜂窝电话网，改动很少几条线路，就可以显著提高性能。

2.小世界网络构成原则

WS小世界网络的构成原则为：从一个环状的规则网络开始，网络含有N个结点，每个结点向与它最近邻的K个结点连出K条边，并满足N>>K>>In(N)>>1。随后进行随机化重连，以概率p随机地重新连接网络中的每个边，即将边的一个端点保持不变，而另一个端点取为网络中随机选择的一个节点。其中规定，任意两个不同的节点之间至多只能有一条边，并且每一

个节点都不能有边与自身相连。这样就会产生pNK/2条长程的边把一个结点和远处的结点联系起来。改变p值可以实现从规则网络(p=0)向随机网络(p=1)转变。

NW小世界网络的构成原则为：从一个环状的规则网络开始，网络含有N个结点，每个结点向与它最近邻的K个结点连出K条边，并满足N>>K>>In(N)>>1。随后进行随机化加边，以概率p在随机选取的一对节点之间加上一条边。其中，任意两个不同的节点之间至多只能有一条边，并且每一个节点都不能有边与自身相连。改变p值可以实现从最近邻耦合网络(p=0)向全局耦合网络(p=1)转变。在p足够小和N足够大时，NW小世界模型本质上等同于WS 小世界模型。

3.MATLAB建模

建立一个初始节点数为20的NW网络。MATLAB程序如下：

function matrix = SW()

%By 201121250314

tic

N=20;m=4;% 初始化网络数据

p=0.1;% 以概率p=0.1在随机选取的一对结点之间加上一条边

matrix=sparse([],[],[],20,20,0);% 创建一个20*20的全0稀疏矩阵

%建立初始的环状的规则网络

%结点网络有N个节点

%每个结点向与它最近邻的m个结点连出边

%求出邻接矩阵

for i=m+1:N-m

for j=i-m:i+m

matrix(i,j)=1;

end

end

for i=1:m

for j=1:i+m

matrix(i,j)=1;

end

end

for i=N-m+1:N

for j=i-m:N

matrix(i,j)=1;

end

end

for i=1:m

for j=N-m+i:N

matrix(i,j)=1;

matrix(j,i)=1;

end

end

%逆时针的边重连,从节点到N-m-1

for i=1:N-m-1

for j=i+1:i+m

r=rand(1);% 随机选取一个数

if r<=p

unconect=find(matrix(i,:)==0);% 取出邻接矩阵中的非0元素位置M=length(unconect);% 求出非0元素个数

r1=ceil(M*rand(1));% 正向取整

matrix(i,unconect(r1))=1;

matrix(unconect(r1),i)=1;% 连接这一对点

%matrix(i,j)=0; matrix(j,i)=0;% 加上这个是SW小世界网络

end

end

end

%逆时针的边重新连接,从节点N-m到N-1

for i=N-m+1:N-1

for j=[i+1:N 1:i- N+m]

r=rand(1);

if r<=p

unconect=find(matrix(i,:)==0);

r1=ceil(length(unconect)*rand(1));

matrix(i,unconect(r1))=1;

matrix(unconect(r1),i)=1;

%matrix(i,j)=0;matrix(j,i)=0;

end

end

end

%逆时针的边重新连接,节点N

for i=N

for j=1:m

r=rand(1);

if r<=p

unconect=find(matrix(i,:)==0);

r1=ceil(length(unconect)*rand(1));

matrix(i,unconect(r1))=1;

matrix(unconect(r1),i)=1;

matrix(i,j)=0;matrix(j,i)=0;

end

end

end

%恢复小世界网络的邻接矩阵

for m=1:N

matrix(m,m)=0;% 去掉自身节点形成的环

end

%存储邻接矩阵

%save data matrix;

toc %计算程序耗时

end

上述程序建立了一个NW小世界网络，求出了其邻接矩阵，用tu_plot()函数画出邻接矩阵的图，就得出了该小世界网络的图形。

function tu_plot(rel,control)

%由邻接矩阵画连接图，输入为邻接矩阵rel，必须为方阵；

%control为控制量，0表示画出的图为无向图，1表示有向图。默认值为0

r_size=size(rel);

%a=size(x)返回的是一个行向量，该行向量第一个元素是

%x的行数，第2个元素是x的列数

if nargin<2 %nargin是用来判断输入变量个数的函数

control=0; %输入变量小于2，即只有一个，就默认control为0

end

if r_size(1)~=r_size(2)%行数和列数不相等，不是方阵，不予处理

disp('Wrong Input! The input must be a square matrix!');

return;

end

len=r_size(1);

rho=10;%限制图尺寸的大小

r=2/1.05^len;%点的半径

theta=0:(2*pi/len):2*pi*(1-1/len);

[pointx,pointy]=pol2cart(theta',rho);

theta=0:pi/36:2*pi;

[tempx,tempy]=pol2cart(theta',r);

point=[pointx,pointy];

hold on

for i=1:len

temp=[tempx,tempy]+[point(i,1)*ones(length(tempx),1),point(i,2)*ones(length(tempx),1)];

plot(temp(:,1),temp(:,2),'r');

text(point(i,1)-0.3,point(i,2),num2str(i));%画点

end

for i=1:len

for j=1:len

if rel(i,j)

link_plot(point(i,:),point(j,:),r,control); %连接有关系的点

end

end

end

set(gca,'XLim',[-rho-r,rho+r],'YLim',[-rho-r,rho+r]);

axis off

function link_plot(point1,point2,r,control)%连接两点

temp=point2-point1;

if (~temp(1))&&(~temp(2))

return;%不画子回路

end

theta=cart2pol(temp(1),temp(2));

[point1_x,point1_y]=pol2cart(theta,r);

point_1=[point1_x,point1_y]+point1;

[point2_x,point2_y]=pol2cart(theta+(2*(theta

point_2=[point2_x,point2_y]+point2;

if control

arrow(point_1,point_2);

else

plot([point_1(1),point_2(1)],[point_1(2),point_2(2)]);

end

得出NW小世界网络的图像如下：

4.分析由于采用的是随机加边的模式，故，每次得到的图形细节都有所不同。程序中，如果加入去除随机加边的原边，即随机重连，就可以得到SW小世界网络。

(完整版)BP神经网络matlab实例(简单而经典).doc

p=p1';t=t1'; [pn,minp,maxp,tn,mint,maxt]=premnmx(p,t); % 原始数据归一化 net=newff(minmax(pn),[5,1],{'tansig','purelin'},'traingdx'); %设置网络,建立相应的BP 网络net.trainParam.show=2000; % 训练网络 net.trainParam.lr=0.01; net.trainParam.epochs=100000; net.trainParam.goal=1e-5; [net,tr]=train(net ,pn,tn); %调用TRAINGDM 算法训练BP 网络 pnew=pnew1'; pnewn=tramnmx(pnew,minp,maxp); anewn=sim(net,pnewn); anew=postmnmx(anewn,mint,maxt); %对 BP 网络进行仿真%还原数据 y=anew'; 1、 BP 网络构建 （1）生成 BP 网络 net newff ( PR,[ S1 S2...SNl],{ TF1 TF 2...TFNl }, BTF , BLF , PF ) PR ：由R 维的输入样本最小最大值构成的R 2 维矩阵。 [ S1 S2...SNl] ：各层的神经元个数。 {TF 1 TF 2...TFNl } ：各层的神经元传递函数。 BTF ：训练用函数的名称。 （2）网络训练 [ net,tr ,Y, E, Pf , Af ] train (net, P, T , Pi , Ai ,VV , TV ) （3）网络仿真 [Y, Pf , Af , E, perf ] sim(net, P, Pi , Ai ,T ) {'tansig','purelin'},'trainrp' BP 网络的训练函数 训练方法 梯度下降法 有动量的梯度下降法 自适应 lr 梯度下降法 自适应 lr 动量梯度下降法弹性梯度下降法训练函数traingd traingdm traingda traingdx trainrp Fletcher-Reeves 共轭梯度法traincgf Ploak-Ribiere 共轭梯度法traincgp

神经网络建模及Matlab中重要的BP网络函数

神经网络建模及Matlab中重要的BP网络函数一、神经组织的基本特征 1．细胞体是一个基本的初等信号处理器，轴突是信号的输出通路，树突是信号的输入通路。信号从一个神经细胞经过突触传递到另一个细胞。 2．不同的神经元之间有不同的作用强度，称为联接强度。当某细胞收到信号时，它的电位发生变化，如果电位超过某一阈值时，该细胞处于激发态，否则处于抑制状态。 3．两神经元之间的联接强度随其激发与抑制行为相关性的时间平均值正比变化，也就是说神经元之间的联接强度不是一成不变的。这就是生物学上的Hebb律。

∑t j ij t S w )(二、人工神经元的M-P 模型（McCulloch 、Pitts,1943） 1．构造一个模拟生物神经组织的人工神经网络的三要素： （1）．对单个神经元给出定义； （2）．定义网络结构：决定神经元数量及连接方式； （3）．给出一种方法，决定神经元之间的联接强度。 2．M-P 模型 其中，t 表示时间 S i (t)表示第i 个神经元在t 时刻的状态，S i (t)=1表示处于激发态，S i (t)=0表示处于抑制态 w ij 表示第j 个神经元到第i 个神经元的联接强度，称之为权，可正可负 表示第i 个神经元在t 时刻所接收到的所有信号的线性迭加。 μi 表示神经元i 的阈值, 可以在模型中增加一个S k (t)=1神经元k ，并且w ik =-μi ,则阈值可归并到和号中去。 注： 1．M-P 神经元虽然简单，但可以完成任何计算。 2．神经元的状态可以取[0，1]中的连续值，如用以下函数代替θ(x): ???<≥=-=+∑0 0011x x x t S w t S i j j ij i )() )(()(θμθ

小世界网络

4.2 小世界网络 4.2.1 小世界网络简介 1998年, Watts和Strogatz 提出了小世界网络这一概念，并建立了WS模型。实证结果表明，大多数的真实网络都具有小世界特性（较小的最短路径）和聚类特性（较大的聚类系数）。传统的规则最近邻耦合网络具有高聚类的特性，但并不具有小世界特性；而随机网络具有小世界特性但却没有高聚类特性。因此这两种传统的网络模型都不能很好的来表示实际的真实网络。Watts和Strogatz建立的小世界网络模型就介于这两种网络之间，同时具有小世界特性和聚类特性，可以很好的来表示真实网络。 4.2.2 小世界模型构造算法 1、从规则图开始：考虑一个含有N个点的最近邻耦合网络，它们围成一个环，其中每个节点都与它左右相邻的各K/2节点相连，K是偶数。 2、随机化重连：以概率p随机地从新连接网络中的每个边，即将边的一个端点保持不变，而另一个端点取为网络中随机选择的一个节点。其中规定，任意两个不同的节点之间至多只能有一条边，并且每一个节点都不能有边与自身相连。 在上述模型中，p=0对应于完全规则网络，p=1则对应于完全随机网络，通过调节p 的值就可以控制从完全规则网络到完全随机网络的过渡。 相应程序代码（使用Matlab实现） ws_net.m （位于“代码”文件夹内） function ws_net() disp('小世界网络模型') N=input('请输入网络节点数'); K=input('请输入与节点左右相邻的K/2的节点数'); p=input('请输入随机重连的概率'); angle=0:2*pi/N:2*pi-2*pi/N; x=100*cos(angle); y=100*sin(angle); plot(x,y,'r.','Markersize',30); hold on; %生成最近邻耦合网络； A=zeros(N); disp(A); for i=1:N if i+K<=N for j=i+1:i+K A(i,j)=1; end else for j=i+1:N A(i,j)=1; end for j=1:((i+K)-N) A(i,j)=1; end

无标度网络模型构造

课题：无标度网络模型构造 姓名赵训 学号201026811130 班级实验班1001

一、源起 无标度网络（或称无尺度网络）的概念是随着对复杂网络的研究而出现的。 “网络”其实就是数学中图论研究的图，由一群顶点以及它们之间所连的边构成。在网络理论中则换一套说法，用“节点”代替“顶点”，用“连结”代替“边”。复杂网络的概念，是用来描述由大量节点以及这些节点之间错综复杂的联系所构成的网络。这样的网络会出现在简单网络中没有的特殊拓扑特性。 自二十世纪60年代开始，对复杂网络的研究主要集中在随机网络上。随机网络，又称随机图，是指通过随机过程制造出的复杂网络。最典型的随机网络是保罗·埃尔德什和阿尔弗雷德·雷尼提出的ER模型。ER模型是基于一种“自然”的构造方法：假设有个节点，并假设每对节点之间相连的 可能性都是常数。这样构造出的网络就是ER模型网络。科学家们最初使用这种模型来解释现实生活中的网络。 ER模型随机网络有一个重要特性，就是虽然节点之间的连接是随机形成的，但最后产生的网络的度分布是高度平等的。度分布是指节点的度的分布情况。在网络中，每个节点都与另外某些节点相连，这种连接的数目叫做这个节点的度。在网络中随机抽取一个节点，它的度是多少呢？这个概率分布就称为节点的度分布。 在一般的随机网络（如ER模型）中，大部分的节点的度都集中在某个特殊值附近，成钟形的泊松分布规律（见下图）。偏离这个特定值的概率呈指数性下降，远大于或远小于这个值的可能都是微乎其微的，就如一座城市中成年居民的身高大致的分布一样。然而在1998年，Albert-László Barab ási、Réka Albert等人合作进行一项描绘万维网的研究时，发现通过超链接与网页、文件所构成的万维网网络并不是如一般的随机网络一样，有着均匀的度分布。他们发现，万维网是由少数高连接性的页面串联起来的。 绝大多数（超过80%）的网页只有不超过4个超链接，但极少数页面（不到总页面数的万分之一）却拥有极多的链接，超过1000个，有一份文件甚至与超过200万个其他页面相连。与居民身高的例子作类比的话，就是说大多数的节点都是“矮个子”，而却又有极少数的身高百丈的“巨人”。Barab ási等人将其称为“无标度”网络。

WS小世界网络模型的程序代码(matlab)

程序仿真实例 例一、请输入最近邻耦合网络中节点的总数N:30 请输入最近邻耦合网络中每个节点的邻居数K:4 请输入随机化重连的概率p:0.9 例二、请输入最近邻耦合网络中节点的总数N:40 请输入最近邻耦合网络中每个节点的邻居数K:2 请输入随机化重连的概率p:0.7

Matlab的m文件代码如下： N=input('请输入最近邻耦合网络中节点的总数N:'); K=input('请输入最近邻耦合网络中每个节点的邻居数K:'); if K>floor(N-1)|mod(K,2)~=0; disp('参数输入错误：K值必须是小于网络节点总数且为偶数的整数'); return ; end angle=0:2*pi./N:2*pi-2*pi/N; angle=0:2*pi/N:2*pi-2*pi/N; x=100*sin(angle); y=100*cos(angle); plot(x,y,'ro','MarkerEdgeColor','g','MarkerFaceColor','r','MarkerSize',8); hold on; A=zeros(N); for i=1:N for j=i+1:i+K/2 jj=j; if j>N jj=mod(j,N); end A(i,jj)=1; A(jj,i)=1; end end %WS小世界网络的代码 p=input('请输入随机化重连的概率p:'); for i=1:N for j=i+1:i+K/2

jj=j; if j>N jj=mod(j,N); end p1=rand(1,1); if p1

无标度网络matlab建模

复杂系统无标度网络研究与建模 XXX 南京信息工程大学XXXX系，南京 210044 摘要：21世纪是复杂性的世界，基于还原论的世界观与方法论已经无法满足当前人们对作为一个整体系统的自然界和人类社会的认识和研究，利用系统科学的方法对科学重新审视已近变为迫切的需要。现实生活中众多复杂网络都具有无标度性，这种无标度网络的增长性和择优连接性很好的解释了富者越富的“马太效应”。对无标度网络的深入研究，让人们深刻的认识到其在Internet、地震网、病毒传播和社会财富分布网中的理论与现实意义。本文通过对复杂网络中的无标度网络的分析与研究，介绍了无标度网络区别于一般随机网络的特性与现实意义，并利用了Matlab生成了一个无标度网络。 关键词：无标度网络，幂律特性，模型建立 1 引言 任何一种网络都可以看作是由一些节点按某种方式连接在一起而构成的一个系统，曾经关于网络结构的研究常常着眼于包含几十个到几百个节点的网络，而近几年关于复杂网络的研究中则常常可以见上万个节点的网络，网络规模尺度上的改变也促使网络分析方法做相应的改变，而复杂网络是近年来随着网络规模、理论和计算机技术的飞速发展而出现的一个新的研究方向。它的出现不仅顺应了现代科技的发展趋势，而且反映了在以信息科学为支柱的新世纪中，各学科理论及应用交叉、渗透和融合的发展趋势[1]。复杂系统主要研究其个体之间相互作用所产生的系统的整体性质与行为“复杂系统的复杂性体现在系统的整体性质与行为往往不是系统各个个体的状态的简单综合”因此，复杂系统的研究不能采用还原论的方法，而要从整体上进行研究。 在对复杂系统的研究中，美国物理学家Barabasi和Albert通过对万维网的研究，发现万维网中网页连接的度分布服从幂律分布，而万维网中少数网页(Hub点)具有非常大的连接，大多数网页的连接数甚小Barabasi等把度分布为幂律分布(Power law)的复杂网络称为无标度网络（scale-free net）[2]。 经过众多的科研工作者的努力，已经证实了现实世界中无论是自然界还是人类社会都广泛的存在着具有度分布符合幂律分布的无标度网络，如生物网络、Internet网、WWW网、演员合作网、科学研究合作网、财富分布网、地震网、电站供电网、科技引文网和病毒传播网等。Newman将这些复杂网络粗略地分成四类:社会网络、信息网络、技术网络和生物网络[3]。

复杂网络模型的matlab实现

function [DeD,aver_DeD]=Degree_Distribution(A) %% 求网络图中各节点的度及度的分布曲线 %% 求解算法：求解每个节点的度，再按发生频率即为概率，求P(k) %A————————网络图的邻接矩阵 %DeD————————网络图各节点的度分布 %aver_DeD———————网络图的平均度 N=size(A,2); DeD=zeros(1,N); for i=1:N % DeD(i)=length(find((A(i,:)==1))); DeD(i)=sum(A(i,:)); end aver_DeD=mean(DeD); if sum(DeD)==0 disp('该网络图只是由一些孤立点组成'); return; else figure; bar([1:N],DeD); xlabel('节点编号n'); ylabel('各节点的度数K'); title('网络图中各节点的度的大小分布图'); end figure; M=max(DeD); for i=1:M+1; %网络图中节点的度数最大为M,但要同时考虑到度为0的节点的存在性 N_DeD(i)=length(find(DeD==i-1)); % DeD=[2 2 2 2 2 2] end P_DeD=zeros(1,M+1); P_DeD(:)=N_DeD(:)./sum(N_DeD); bar([0:M],P_DeD,'r'); xlabel('节点的度 K'); ylabel('节点度为K的概率 P(K)'); title('网络图中节点度的概率分布图'); 平均路径长度 function [D,aver_D]=Aver_Path_Length(A) %% 求复杂网络中两节点的距离以及平均路径长度 %% 求解算法：首先利用Floyd算法求解出任意两节点的距离，再求距离的平均值得平均路

BP神经网络地设计实例(MATLAB编程)

神经网络的设计实例(MATLAB编程) 例1 采用动量梯度下降算法训练BP 网络。训练样本定义如下： 输入矢量为 p =[-1 -2 3 1 -1 1 5 -3] 目标矢量为t = [-1 -1 1 1] 解：本例的MATLAB 程序如下： close all clear echo on clc % NEWFF——生成一个新的前向神经网络% TRAIN——对BP 神经网络进行训练 % SIM——对BP 神经网络进行仿真pause % 敲任意键开始 clc % 定义训练样本 P=[-1, -2, 3, 1; -1, 1, 5, -3]; % P 为输入矢量T=[-1, -1, 1, 1]; % T 为目标矢量

clc % 创建一个新的前向神经网络 net=newff(minmax(P),[3,1],{'tansig','purelin'},'traingdm') % 当前输入层权值和阈值 inputWeights=net.IW{1,1} inputbias=net.b{1} % 当前网络层权值和阈值 layerWeights=net.LW{2,1} layerbias=net.b{2} pause clc % 设置训练参数 net.trainParam.show = 50; net.trainParam.lr = 0.05; net.trainParam.mc = 0.9; net.trainParam.epochs = 1000; net.trainParam.goal = 1e-3; pause clc % 调用TRAINGDM 算法训练BP 网络 [net,tr]=train(net,P,T);

基于Matlab的无标度网络仿真

无标度网络 无标度网络（或称无尺度网络）的概念是随着对复杂网络的研究而出现的。“网络”其实就是数学中图论研究的图，由一群顶点以及它们之间所连的边构成。在网络理论中则换一套说法，用“节点”代替“顶点”，用“连结”代替“边”。复杂网络的概念，是用来描述由大量节点以及这些节点之间错综复杂的联系所构成的网络。ER模型随机网络有一个重要特性，就是虽然节点之间的连接是随机形成的，但最后产生的网络的度分布是高度平等的。度分布是指节点的度的分布情况。在网络中，每个节点都与另外某些节点相连，这种连接的数目叫做这个节点的度。在网络中随机抽取一个节点，它的度是多少呢？这个概率分布就称为节点的度分布。 自二十世纪60年代开始，对复杂网络的研究主要集中在随机网络上。随机网络，又称随机图，是指通过随机过程制造出的复杂网络。最典型的随机网络是保罗·埃尔德什和阿尔弗雷德·雷尼提出的ER模型。ER模型是基于一种“自然”的构造方法：假设有n个节点，并假设每对节点之间相连的可能性都是常数。这样构造出的网络就是ER模型网络。 Matlab程序如下： SFNG： function SFNet = SFNG(Nodes, mlinks, seed) seed = full(seed); pos = length(seed); rand('state',sum(100*clock)); Net = zeros(Nodes, Nodes, 'single'); Net(1:pos,1:pos) = seed; sumlinks = sum(sum(Net)); while pos < Nodes pos = pos + 1; linkage = 0; while linkage ~= mlinks rnode = ceil(rand * pos); deg = sum(Net(:,rnode)) * 2; rlink = rand * 1; if rlink < deg / sumlinks && Net(pos,rnode) ~= 1 && Net(rnode,pos) ~= 1 Net(pos,rnode) = 1; Net(rnode,pos) = 1; linkage = linkage + 1; sumlinks = sumlinks + 2; end end end clear Nodes deg linkage pos rlink rnode sumlinks mlinks

课题：WS小世界网络模型构造

课题：WS小世界网络模型构造 姓名赵训 学号 2 班级计算机实验班

一、WS 小世界网络简介 1998年, Watts和Strogatz 提出了小世界网络这一概念,并建立了WS模型。实证结果表明,大多数的真实网络都具有小世界特性(较小的最短路径) 和聚类特性(较大的聚类系数) 。传统的规则最近邻耦合网络具有高聚类的特性,但并不具有小世界特性;而ER 随机网络具有小世界特性但却没有高聚类特性。因此这两种传统的网络模型都不能很好的来表示实际的真实网络。 Watts 和Strogatz建立的WS小世界网络模型就介于这两种网络之间,同时具有小世界特性和聚类特性,可以很好的来表示真实网络。 二、WS小世界模型构造算法 1、从规则图开始：考虑一个含有N个点的最近邻耦合网络，它们围成一个环，其中每个节点都与它左右相邻的各K/2节点相连，K是偶数。 2、随机化重连：以概率p随机地从新连接网络中的每个边，即将边的一个端点保持不变，而另一个端点取为网络中随机选择的一个节点。其中规定，任意两个不同的节点之间至多只能有一条边，并且每一个节点都不能有边与自身相连。 在上述模型中，p=0对应于完全规则网络，p=1则对应于完全随机网络，通过调节p的值就可以控制从完全规则网络到完全随机网络的过渡，如图a所示。 图a 相应程序代码（使用Matlab实现） ws_net.m （位于“代码”文件夹内） function ws_net() disp('WS小世界网络模型') N=input('请输入网络节点数'); K=input('请输入与节点左右相邻的K/2的节点数'); p=input('请输入随机重连的概率'); angle=0:2*pi/N:2*pi-2*pi/N; x=100*cos(angle); y=100*sin(angle); plot(x,y,'r.','Markersize',30); hold on; %生成最近邻耦合网络； A=zeros(N);

Matlab训练好的BP神经网络如何保存和读取方法(附实例说明)

Matlab训练好的BP神经网络如何保存和读取方法（附实例说明） 看到论坛里很多朋友都在提问如何存储和调用已经训练好的神经网络。 本人前几天也遇到了这样的问题，在论坛中看了大家的回复，虽然都提到了关键的两个函数“save”和“load”，但或多或少都简洁了些，让人摸不着头脑（呵呵，当然也可能是本人太菜）。通过不断调试，大致弄明白这两个函数对神经网络的存储。下面附上实例给大家做个说明，希望对跟我有一样问题的朋友有所帮助。 如果只是需要在工作目录下保到当前训练好的网络，可以在命令窗口 输入：save net %net为已训练好的网络 然后在命令窗口 输入：load net %net为已保存的网络 加载net。 但一般我们都会在加载完后对网络进行进一步的操作，建议都放在M文件中进行保存网络和调用网络的操作 如下所示： %% 以函数的形式训练神经网络 functionshenjingwangluo() P=[-1,-2,3,1; -1,1,5,-3]; %P为输入矢量 T=[-1,-1,1,1,]; %T为目标矢量 net=newff(minmax(P),[3,1],{'tansig','purelin'},'traingdm') %创建一个新的前向神经网络 inputWeights=net.IW{1,1} inputbias=net.b{1} %当前输入层权值和阀值 layerWeights=net.LW{2,1} layerbias=net.b{2} net.trainParam.show=50; net.trainParam.lr=0.05; net.trainParam.mc=0.9;

小世界网络综述

关于小世界网络的文献综述 一、小世界网络概念方面的研究 Watts和Strogatz开创性的提出了小世界网络并给出了WS小世界网络模型。小世界网络的主要特征就是具有比较小的平均路径长度和比较大的聚类系数。所谓网络的平均路径长度，是指网络中两个节点之间最短路径的平均值。聚类系数被用来描述网络的局部特征，它表示网络中两个节点通过各自相邻节点连接在一起的可能性，以及衡量网络中是否存在相对稳定的子系统。规则网络具有大的特征路径长度和高聚类系数，随机网络则有短的特征路径长度和比较小的聚类系数[1]。 Guare于1967年在《今日心理学》杂志上提出了“六度分离”(Six Degrees of Separation) 理论，即“小世界现象”。该理论认为，在社交网络中存在短路径，即人们只要知道自己认识的人，就能很快地把信息传递到任何远方目标[2]。 .Stanleymilgram的邮件试验，后来的“培根试验”，以及1998年《纽约时代周刊》的关于莱温斯基的讽刺性游戏，都表现出：似乎在庞大的网络中各要素之间的间隔实际很“近”，科学家们把这种现象称为小世界效应[3]。研究发现，世界上任意两个人可以平均通过6个人联系在一起，人们称此现象为“六度分离” [2]。 二、小世界网络模型方面的研究 W-S模型定义了两个特征值：a.特征路径的平均长度L。它是指能使网络中各个结点相连的最少边长度的平均数，也就是上面说的小世界网络平均距离。b.集团化系数C。网络结点倾向于结成各种小的集团,它描述网络局部聚类特征。 稍后，Newman和WattS对上述的WS模型作了少许改动，提出了另一个相近但较好的(NW)小世界网络模型[5]，其做法是不去断开原来环形初始网络的任何一条边、而只是在随机选取的节点对之间增加一条边(这时，新连接的边很可能是长程边)。这一模烈比WS模型容易分析，因为它在形成过程中不会出现孤立的竹点簇。 其次，还有Monasson小世界网络模型[6]以及一些其它的变形模型包括BW 小世界网络模型等等[7]。 三、小世界网络应用方面的研究 ①、在生物学领域的应用 Wdt怡和StrogatZ证明疾病全球传播所需的时间和特征路径长度非常相似，只要在传播网络中加人一些捷径就可以使传播速度明显加快。运用病毒在小世界网络中的传播性质可推出信息在一个平均分离度为6的网络中传播要比在平均分离度为一百或一百万的网络中快得多[8]。 许多知名的生物网络表现出了小世界网络节点间的关连性。一般的小世界网络模型,也利用了网络的无向和无标度特性来展示网络中各节点之间的联系。这种网络模型能模拟一些神经网络的重要性质，例如，染色体结合的方向和标度。 [9][10]。 有学者研究了基于神经网络的有小世界结构的联想记忆模型。这一网络检索某一存储的模型的有效性展示了混乱的有限价值的阶段转换。更加常规化的网络很难恢复这个模型，而对混合的不对称的状态更有效。[11]。 ②、在博弈论方面的应用

无标度网络及MATLAB建模

无标度网络 1.简介 传统的随机网络（如ER模型），尽管连接是随机设置的，但大部分节点的连接数目会大致相同，即节点的分布方式遵循钟形的泊松分布，有一个特征性的“平均数”。连接数目比平均数高许多或低许多的节点都极少，随着连接数的增大，其概率呈指数式迅速递减。故随机网络亦称指数网络。 现实世界的网络大部分都不是随机网络，少数的节点往往拥有大量的连接，而大部分节点却很少，一般而言他们符合zipf定律，(也就是80/20马太定律)。人们给具有这种性质的网络起了一个特别的名字——无标度网络。这里的无标度是指网络缺乏一个特征度值（或平均度值），即节点度值的波动范围相当大。 现实中的交通网，电话网和Internet都是无标度网络，在这种网络中，存在拥有大量连接的集散节点。分布满足幂律的无标度网络还具有一个奇特的性质—“小世界”特性。虽然万维网中的页面数已超过80亿，但平均来说，在万维网上只需点击19次超链接，就可从一个网页到达任一其它页面。 无标度网络具有严重的异质性，其各节点之间的连接状况（度数）具有严重的不均匀分布性：网络中少数称之为Hub点的节点拥有极其多的连接，而大多数节点只有很少量的连接。少数Hub点对无标度网络的运行起着主导的作用。从广义上说，无标度网络的无标度性是描述大量复杂系统整体上严重不均匀分布的一种内在性质。 1999 年, Albert、Jeong和Barabs发现万维网网页的度分布不是通常认为的Poisson 分布,而是重尾特征的幂律分布,而且万维网基本上是由少数具有大量超链接的网页串连起来的, 绝大部分网页的链接很少,他们把网络的这个特性称为无标度性(Scale-free nature, SF)。1999 年Barabs和Albert考察了实际网络的生成机制, 发现增长和择优连接是实际网络演化过程的两个基本要素, 他们创造性地构建了能够产生无标度特性的第一个网络模型——BA 模型。 BA 网络主要具有以下特性: 具有幂律度分布, 是一个无标度网络; 具有小世界特征。幂律度分布的重尾特征导致无标度网络中有少数具有大量连接边的中枢点, 择优连接必然产生“富者愈富”的现象。BA 网络同时具有鲁棒性和脆弱性，面对结点的随机失效, 网络具有鲁棒性；但面对蓄意攻击时, 由于中枢点的存在, 网络变得十分脆弱, 很容易陷于瘫痪。 特别地, 网络传染性疾病在无标度网络中不存在传播阈值, 疾病一旦产生就在网络上迅速传播并达到稳定状态。如果没有人为干预, 疾病将在网络中永远存在, 不会自动灭绝。这对制定无标度网络上的网络疾病防控策略提出了重大挑战。 2.BA无标度网络构成原则 ( 1) 增长: 网络开始于少数几个结点(初始设定为m0个) , 每个相等时间间隔增加一个新点, 新点与m个(m小于等于m0)不同的已经存在于网络中的旧点相连产生m条新边。 （2）择优连接：新点与旧点i相连的概率P取决于结点i的度数ki。

小世界网络及其性质

小世界网络及其性质 复杂网络是多主体系统的一个子集，对它的研究是计算经济学的一个重要研究领域。而贸易网络又是复杂网络在社会科学领域的一个子集，因此要讨论贸易网络，首先要关注关于复杂网络的一般性研究。在本章中，我们先介绍复杂网络研究的一个重要成果，即小世界模型，这为后面讨论贸易网络做准备。因为小世界模型只是一个纯粹的数学问题，并不包含某一学科的特殊含义，因此要首先讨论贸易网络的经济学意义，这是从一个简单的分工协调问题开始的。之后，我们进一步的追问，贸易网络是否也具有小世界这样一个普遍存在的性质。 1、多主体系统中的复杂网络 我们生活在各种各样的网络之中。在与同学、朋友、老师交往的时候，我们处于一个人际关系网络中；在我们使用各种电器的时候，我们处于一个电力网络中，很多发电厂、变电所、输电线构成了这个网络；我们去银行取钱的时候，便处于一个银行网络之中，而银行又可以与各种投资者，贷款人联系，这又是一个更大的网络。复杂网络是最近几年新兴起来的一个研究方向，如图 3.1 所示，复杂网络是多主体系统的一个子集，它里面还包括了社会关系网络（社会学）、神经网络（生物学）、计算机网络（计算机科学）、贸易网络（经济学）等等诸多的网络类别。复杂网络是多主体系统研究的一个重要分支，对社会网络、贸易网络的分析是计算经济学（ACE）的一个主要领域。 ACE 在社会科学角度对复杂网络的研究主要关注以下几个方面：（1）市场或 人际关系网络的拓扑结构究竟是什么样的。（2）这些网络结构的微观基础是什么，即如何从个体的行为出发，通过自下而上的建模涌现出这样的网络。（3）社会科学领域的网络与其它领域的网络能否找到一致的共性特征。社会科学可

小世界网络简介及MATLAB建模

小世界网络简介及MATLAB建模 1.简介 小世界网络存在于数学、物理学和社会学中，是一种数学图的模型。在这种图中大部份的结点不与彼此邻接，但大部份结点可以通过任一其它节点经少数几步就可以产生联系。若将一个小世界网络中的点代表一个人，而联机代表人与人之间是相互认识的，则这小世界网络可以反映陌生人通过彼此共同认识的人而起来产生联系关系的小世界现象。 在日常生活中，有时你会发现，某些你觉得与你隔得很“遥远”的人，其实与你“很近”。小世界网络就是对这种现象的数学描述。用数学中图论的语言来说，小世界网络就是一个由大量顶点构成的图，其中任意两点之间的平均路径长度比顶点数量小得多。除了社会人际网络以外，小世界网络的例子在生物学、物理学、计算机科学等领域也有出现。许多经验中的图可以用小世界网络来作为模型。因特网、公路交通网、神经网络都呈现小世界网络的特征。 小世界网络最早是由邓肯·瓦茨（Duncan Watts）和斯蒂文·斯特罗加茨（Steven Strogatz）在1998年引进的，将高聚合系数和低平均路径长度作为特征，提出了一种新的网络模型，一般就称作瓦茨-斯特罗加茨模型（WS模型），这也是最典型的小世界网络的模型。 由于WS小世界模型构造算法中的随机化过程有可能破坏网络的连通性，纽曼(Newman)和瓦茨(Watts)提出了NW小世界网络模型，该模型是通过用“随机化加边”模式来取代WS小世界网络模型构造中的“随机化重连”。 在考虑网络特征的时候，使用两个特征来衡量网络：特征路径长度和聚合系数。 特征路径长度（characteristic path length）：在网络中，任选两个节点，连同这两个节点的最少边数，定义为这两个节点的路径长度，网络中所有节点对的路径长度的平均值，定义为网络的特征路径长度。这是网络的全局特征。 聚合系数(clustering coefficient)：假设某个节点有k个边，则这k条边连接的节点之间最多可能存在的边的个数为k(k-1)/2，用实际存在的边数除以最多可能存在的边数得到的分数值，定义为这个节点的聚合系数。所有节点的聚合系数的均值定义为网络的聚合系数。聚合系数是网络的局部特征，反映了相邻两个人之间朋友圈子的重合度，即该节点的朋友之间也是朋友的程度。 我们可以发现规则网络具有很高的聚合系数，大世界（large world，意思是特征路径长度很大），其特征路径长度随着n(网络中节点的数量)线性增长，而随机网络聚合系数很小，小世界（small world，意思是特征路径长度小），其特征路径长度随着log(n)增长中说明，在从规则网络向随机网络转换的过程中，实际上特征路径长度和聚合系数都会下降，到变成随机网络的时候，减少到最少。但这并不是说大的聚合系数一定伴随着大的路径长度，而小的路径长度伴随着小的聚合系数，小世界网络就具有大的聚合系数，而特征路径长度很小。试验表明，少量的short cut的建立能够迅速减少特征路径长度，而聚合系数变化却不大，因为某一个short cut的建立，不仅影响到所连接的节点的特征路径长度，而且影响到他们邻居的路径长度，而对整个网络的聚合系数影响不大。这样，少量的short cut的建立就能使整个网络不知不觉地变成小世界网络。 实际的社会、生态、等网络都是小世界网络，在这样的系统里，信息传递速度快，并且少量改变几个连接，就可以剧烈地改变网络的性能，如对已存在的网络进行调整，

复杂网络理论及其研究现状

复杂网络理论及其研究现状 复杂网络理论及其研究现状 【摘要】简单介绍了蓬勃发展的复杂网络研究新领域，特别是其中最具代表性的是随机网络、小世界网络和无尺度网络模型；从复杂网络的统计特性、复杂网络的演化模型及复杂网络在社会关系研究中的应用三个方面对其研究现状进行了阐述。 【关键词】复杂网络无标度小世界统计特性演化模型 一、引言 20世纪末，以互联网为代表的信息技术的迅速发展使人类社会步入了网络时代。从大型的电力网络到全球交通网络，从Internet 到WWW，从人类大脑神经到各种新陈代谢网络，从科研合作网络到国际贸易网络等，可以说，人类生活在一个充满着各种各样的复杂网络世界中。 在现实社会中，许多真实的系统都可以用网络的来表示。如万维网（WWW网路）可以看作是网页之间通过超级链接构成的网络；网络可以看成由不同的PC通过光缆或双绞线连接构成的网络；基因调控网络可以看作是不同的基因通过调控与被调控关系构成的网络；科学家合作网络可以看成是由不同科学家的合作关系构成的网络。复杂网络研究正渗透到数理科学、生物科学和工程科学等不同的领域，对复杂网络的定性与定量特征的科学理解，已成为网络时代研究中一个极其重要的挑战性课题，甚至被称为“网络的新科学”。 二、复杂网络的研究现状 复杂网络是近年来国内外学者研究的一个热点问题。传统的对网络的研究最早可以追溯到18世纪伟大数学家欧拉提出的著名的“Konigsberg七桥问题”。随后两百多年中，各国的数学家们一直致力于对简单的规则网络和随机网络进行抽象的数学研究。规则网络过于理想化而无法表示现实中网络的复杂性，在20世纪60年代由Erdos和Renyi（1960）提出了随机网络。进入20世纪90年代，人们发现现实世界中绝大多数的网络既不是完全规则，也不是完全随机

BP神经网络matlab实例

神经网络Matlab p=p1';t=t1'; [pn,minp,maxp,tn,mint,maxt]=premnmx(p,t); %原始数据归一化net=newff(minmax(pn),[5,1],{'tansig','purelin'},'traingdx');%设置网络,建立相应的BP网络 net.trainParam.show=2000; % 训练网络 net.trainParam.lr=0.01; net.trainParam.epochs=100000; net.trainParam.goal=1e-5; [net,tr]=train(net ,pn,tn); %调用TRAINGDM算法训练BP网络 pnew=pnew1'; pnewn=tramnmx(pnew,minp,maxp); anewn=sim(net,pnewn); %对BP网络进行仿真 anew=postmnmx(anewn,mint,maxt); %还原数据 y=anew'; 1、BP网络构建 （1）生成BP网络 = net newff PR S S SNl TF TF TFNl BTF BLF PF (,[1 2...],{ 1 2...},,,) R?维矩阵。 PR：由R维的输入样本最小最大值构成的2

S S SNl：各层的神经元个数。 [1 2...] TF TF TFNl：各层的神经元传递函数。 { 1 2...} BTF：训练用函数的名称。 （2）网络训练 = [,,,,,] (,,,,,,) net tr Y E Pf Af train net P T Pi Ai VV TV （3）网络仿真 = [,,,,] (,,,,) Y Pf Af E perf sim net P Pi Ai T {'tansig','purelin'},'trainrp' BP网络的训练函数 训练方法训练函数 梯度下降法traingd 有动量的梯度下降法traingdm 自适应lr梯度下降法traingda 自适应lr动量梯度下降法traingdx 弹性梯度下降法trainrp Fletcher-Reeves共轭梯度法traincgf Ploak-Ribiere共轭梯度法traincgp Powell-Beale共轭梯度法traincgb 量化共轭梯度法trainscg 拟牛顿算法trainbfg 一步正割算法trainoss Levenberg-Marquardt trainlm

MATLAB的建模和仿真

课程设计说明书 题目：基于Matlab的IIR滤波器设计与仿真班级：2012 级电气五班 姓名：王璐 学号：201295014178 指导教师：张小娟 日期：2015年 1 月12日

课程设计任务书

基于MATLAB的IIR滤波器设计与仿真 前言 数字信号处理（digital signal processing，DSP）是从20世纪60年代以来，随着信息学科和计算机学科的高速发展而迅速发展起来的一门新兴学科。数字信号处理是把信号用数字或符号表示的序列，通过计算机或通用（专用）信号处理设备，用数字的数值计算方法处理（例如滤波、变换、压缩、增强、估计、识别等），以达到提取有用信息便于应用处理的目的。数字信号处理系统有精度高、灵活性高、可靠性高、容易大规模集成、时分复用、可获得高性能指标、二维与多维处理等特点。正是由于这些突出的特点，使得它在通信、语音、雷达、地震测报、声呐、遥感、生物医学、电视、仪器中得到愈来愈广泛的应用。在数字信号处理中起着重要的作用并已获得广泛应用的是数字滤波器（DF，Digital Filter），根据其单位冲激响应函数的时域特性可分为两类：无限冲激响应IIR（Infinite Impulse Response）滤波器和有限冲激响应FIR（Finite Impulse Response）滤波器。MATLAB的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来结算问题要比用C，FORTRAN等语言完成相同的事情简捷得多，并且MATLAB也吸收了像Maple等软件的有点，使MATLAB成为一个强大的数学软件，在新的版本中也加入了对C，FORTRAN，C++，JA V A的支持。可以直接调用，用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用。 1 数字滤波器概述 数字滤波器是对数字信号实现滤波的线性时不变系统。数字滤波实质上是一种运算过程，实现对信号的运算处理。输入数字信号（数字序列）通过特定的运算转变为输出的数字序列，因此，数字滤波器本质上是一个完成特定运算的数字计算过程，也可以理解为一台计算机。描述离散系统输出与输入关系的卷积和差分方程只是给数字信号滤波器提供运算规则，使其按照这个规则完成对输入数据的处理。时域离散系统的频域特性：Y(eωj)=X(eωj)H(eωj) 其中Y(eωj)、X(eωj)分别是数字滤波器的输出序列和输入序列的频域特性（或称为

无标度网络matlab建模

无标度网络m a t l a b建模Last revision on 21 December 2020

复杂系统无标度网络研究与建模 XXX 南京信息工程大学XXXX系，南京 210044 摘要：21世纪是复杂性的世界，基于还原论的世界观与方法论已经无法满足当前人们对作为一个整体系统的自然界和人类社会的认识和研究，利用系统科学的方法对科学重新审视已近变为迫切的需要。现实生活中众多复杂网络都具有无标度性，这种无标度网络的增长性和择优连接性很好的解释了富者越富的“马太效应”。对无标度网络的深入研究，让人们深刻的认识到其在Internet、地震网、病毒传播和社会财富分布网中的理论与现实意义。本文通过对复杂网络中的无标度网络的分析与研究，介绍了无标度网络区别于一般随机网络的特性与现实意义，并利用了Matlab 生成了一个无标度网络。 关键词：无标度网络，幂律特性，模型建立 1 引言 任何一种网络都可以看作是由一些节点按某种方式连接在一起而构成的一个系统，曾经关于网络结构的研究常常着眼于包含几十个到几百个节点的网络，而近几年关于复杂网络的研究中则常常可以见上万个节点的网络，网络规模尺度上的改变也促使网络分析方法做相应的改变，而复杂网络是近年来随着网络规模、理论和计算机技术的飞速发展而出现的一个新的研究方向。它的出现不仅顺应了现代科技的发展趋势，而且反映了在以信息科学为支柱的新世纪中，各学科理论及应用交叉、渗透和融合的发展趋势[1]。复杂系统主要研究其个体之间相互作用所产生的系统的整体性质与行为“复杂系统的复杂性体现在系统的整体性质与行为往往不是系统各个个体的状态的简单综合”因此，复杂系统的研究不能采用还原论的方法，而要从整体上进行研究。 在对复杂系统的研究中，美国物理学家Barabasi和Albert通过对万维网的研究，发现万维网中网页连接的度分布服从幂律分布，而万维网中少数网页(Hub点)具有非常大的连接，大多数网页的连接数甚小Barabasi等把度分布为幂律分布(Power law)的复杂网络称为无标度网络（scale-free net）[2]。