MATLAB实例:构造网络连接图(Network Connection)及计算图的代数连通度(Algebraic Connectivity)
MATLAB实例:构造网络连接图(Network Connection)及计算图的代数连通度(Algebraic Connectivity)
作者:凯鲁嘎吉 - 博客园 http://www.cnblogs.com/kailugaji/
1. 图的代数连通度(Algebraic Connectivity)
图的代数连通度:Laplace图谱的次小特征值。
2. 网络连接图(Network Connection)的构造
随机生成一个具有50个节点的传感器网络。节点随机放置在3.5 x 3.5方形区域内,通信距离为0.8。如下图所示,共有159条边,其代数连通度为:0.3007。
3. MATLAB程序
demo_Create_Network_Connection.m
%创建无向图 网络连接图 Network Connection.
clc;
close all;
clear; Conf.Square = 3.5; %方形区域的边长
Conf.NodeNumber = 50; %节点个数
Conf.CommDist = 0.8; %最大通信距离 is_create_network = 1;
if is_create_network == 1
[ Network, Dists ] = CreateNetworksFunc(Conf);
save Network_1.mat Network
else
load Network_1.mat
end nodenum = size(Network.Nodes.loc,1); %节点个数
lap_matrix = zeros(nodenum); %节点数*节点数 图的Laplace矩阵:diag(d1,d2,...dn)-邻接矩阵,di为节点i的度
for i=1:nodenum
idx = Network.Nodes.neighbors{i}; %邻接节点的id
lap_matrix(i,idx) = -1; %负的邻接矩阵
lap_matrix(i,i) = length(idx); %对角线元素为节点的度
end
eig_val = eig(lap_matrix); %lap_matrix的特征值
eig_val = sort(eig_val,'ascend'); %从小到大排序,最小特征值为0
algeb_conn = eig_val(2) % algebraic connectivity 代数连通度:lap_matrix的第二小特征值>0,连通图
avg_deg = sum(diag(lap_matrix))/nodenum % average values 节点度的均值 DrawNetworks(Network);
% DrawNetworks(Network, Dists); %把所有的边的长度(通信距离)都标出来了
print(gcf,'-dpng','Network_1.png'); %保存图片
CreateNetworksFunc.m
function [ Network, Dists ] = CreateNetworksFunc(Conf)
% 创建无向图 网络连接图 Network Connection.
num = Conf.NodeNumber; %节点个数
square = Conf.Square; %方形区域的边长
maxDist = Conf.CommDist; %最大通信距离 loc = square*rand(num,2) - square/2; %num*2的随机数 节点坐标
Dists = Euclid_Dist(loc(:,1),loc(:,2)); %节点数*节点数,对角线元素为0 % without self-loop 不存在节点自己到自己的路径,对角线上的元素为无穷大
Dists = Dists + 10*maxDist*eye(num); Neighbors = cell(num,1);
maxDegree = 0; %节点的最大度,与节点相邻的最大边数
edges = 0; %图的总边的个数,无向图的度/2
for i=1:num
Neighbors{i} = find(Dists(i,:)<=maxDist); %找邻接节点的id
if length(Neighbors{i}) > maxDegree
maxDegree = length(Neighbors{i}); %节点的最大度
end
edges = edges + length(Neighbors{i});
end Nodes.loc = loc;
Nodes.neighbors = Neighbors; Network.maxDegree = maxDegree;
Network.edges = edges/2; %% undirected graph
Network.Conf = Conf;
Network.Nodes = Nodes;
end function dist = Euclid_Dist(X,Y)
% 求两两节点之间的距离,输出[节点*节点]的矩阵,距离矩阵
len = length(X);
xx = repmat(X,1,len); %节点数*节点数
yy = repmat(Y,1,len);
dist = sqrt((xx-xx').^2+(yy-yy').^2); %节点数*节点数
end
DrawNetworks.m
function fig = DrawNetworks( Network )
%画无向图 网络连接图 Network Connection.
% function fig = DrawNetworks( Network, Dists ) %把所有的边的长度(通信距离)都标出来了 num = Network.Conf.NodeNumber; %节点个数
loc = Network.Nodes.loc; %节点坐标
square = Network.Conf.Square; %方形区域的边长
Neighbors = Network.Nodes.neighbors; %邻接节点的id fig = figure;
plot(loc(:,1),loc(:,2),'ro','MarkerSize',8,'LineWidth',2); %节点是红色圆圈
side=ceil(square/2);
axis([-side,side,-side,side]);
for i=1:num
for k = 1:length(Neighbors{i})
j = Neighbors{i}(k);
% c = num2str(Dists(i,j),'%.2f');
% text((loc(i,1) + loc(j,1))/2,(loc(i,2) + loc(j,2))/2,c,'Fontsize',10); %把所有的边的长度(通信距离)都标出来了
% hold on;
line([loc(i,1),loc(j,1)],[loc(i,2),loc(j,2)],'LineWidth',0.8,'Color','b'); %线是蓝色
end
end
set(gcf, 'Color', 'w'); %白色
end
4. 参考文献
[1] Hua J, Li C. Distributed variational Bayesian algorithms over sensor networks[J]. IEEE Transactions on Signal Processing, 2015, 64(3): 783-798.
[2] 肖恩利, 束金龙, 闻人凯. 图的代数连通度及其点连通度[J]. 华东师范大学学报(自然科学版), 2003, 2003(4):1-4.
[3] Junhao Hua. Distributed Variational Bayesian Algorithms. Github, 2017.
MATLAB实例:构造网络连接图(Network Connection)及计算图的代数连通度(Algebraic Connectivity)的更多相关文章
- MATLAB实例:聚类网络连接图
MATLAB实例:聚类网络连接图 作者:凯鲁嘎吉 - 博客园 http://www.cnblogs.com/kailugaji/ 本文给出一个简单实例,先生成2维高斯数据,得到数据之后,用模糊C均值( ...
- 使用putty连接虚拟机上的centos提示Network:connection refused
转自:https://yeyuan.iteye.com/blog/1266484 今天早上开机之后,像往常一样使用putty连接linux的时候,突然提示Network:connection refu ...
- MATLAB实例:求相关系数、绘制热图并找到强相关对
MATLAB实例:求相关系数.绘制热图并找到强相关对 作者:凯鲁嘎吉 - 博客园 http://www.cnblogs.com/kailugaji/ 用MATLAB编程,求给定数据不同维度之间的相关系 ...
- MATLAB实例:散点密度图
MATLAB实例:散点密度图 作者:凯鲁嘎吉 - 博客园http://www.cnblogs.com/kailugaji/ MATLAB绘制用颜色表示数据密度的散点图 数据来源:MATLAB中“fit ...
- MATLAB实例:绘制折线图
MATLAB实例:绘制折线图 作者:凯鲁嘎吉 - 博客园 http://www.cnblogs.com/kailugaji/ 条形图的绘制见:MATLAB实例:绘制条形图 用MATLAB将几组不同的数 ...
- 初次使用git就遭遇不测,提示没有这个服务连接和需要配置git的一个http参数 NO network connection,SSl host could not be verified ...
第一次使用git 拉取服务上的项目到本地,结果,在拿到访问的url地址后,输入用户名密码,失败了. --eclispe 4.5.3 继承了git客户端插件的版本 ------下一步后,报错 NO n ...
- Delaunay三角剖分及MATLAB实例
https://blog.csdn.net/piaoxuezhong/article/details/68065170 一.原理部分 点集的三角剖分(Triangulation),对数值分析(如有限元 ...
- 处于同一个域中的两台Sql server 实例无法连接
处于同一个域中的两台Sql server 实例无法连接,报的错误信息如下: A network-related or instance-specific error occurred while es ...
- openstack(liberty): 简单网络连接图
openstack起初的网络部分是和计算核心nova合在一起的,后来被拆分出来,独立成为一个模块, 现在名为Neutron. 本博文是学习记录,记录的是基于GRE tunnel技术的neutron和计 ...
随机推荐
- gitlab如何从Github导入项目
本文简单演示如何Github上导入项目到私人搭建的Gitlab中,搭建过程参考:CentOS7 搭建gitlab服务器. Gitlab版本是gitlab-ce-12.0.2,界面可能稍有差异,但应该影 ...
- Git多个远程仓库不同步时的补救办法
git本地仓库是可以与多个远程仓库关联的,如果想知道怎么配置,请参考Git如何使用多个托管平台管理代码 . 当git remote关联了多个远程仓库时,总会遇到一些问题.今天就遇到了两个远程仓库不一致 ...
- LeetCode 将一个按照升序排列的有序数组,转换为一棵高度平衡二叉搜索树
第108题 将一个按照升序排列的有序数组,转换为一棵高度平衡二叉搜索树. 本题中,一个高度平衡二叉树是指一个二叉树每个节点 的左右两个子树的高度差的绝对值不超过 1. 示例: 给定有序数组: [-10 ...
- C#基础之多线程与异步
1.基本概念 多线程与异步是两个不同概念,之所以把这两个放在一起学习,是因为这两者虽然有区别,但也有一定联系. 多线程是一个技术概念,相对于单线程而言,多线程是多个单线程同时处理逻辑.例如,假如说一个 ...
- tesseract-OCR + pytesseract安装
1. tesseract-OCR下载安装 地址:https://digi.bib.uni-mannheim.de/tesseract/ 选择一个版本下载,下载完成点击**.exe进行安装,若无其他需求 ...
- Delphi CreateProcess 创建一个新的进程和它的主线程
Delphi CreateProcess WIN32API函数CreateProcess用来创建一个新的进程和它的主线程,这个新进程运行指定的可执行文件 CreateProcess百科名片 WIN32 ...
- MySQL 字符串分隔成行&子串统计
利用help_topic表把字符串转换成行(分隔符号',') ),) as `id` FROM mysql.`help_topic`; 统计字符串每个子串出现次数(分隔符号',') ),) AS `s ...
- 小白到大神,你需要了解的 sqlite 最佳实践
本文微信公众号「AndroidTraveler」首发. 背景 本文是对一篇英文文档的翻译,原文请见文末链接. 并发数据库访问 假设你实现了自己的 SQLiteOpenHelper. public cl ...
- springboot配置文件(一)
一.YAML语法 1.基本语法 k 空格 v 表示一对键值对(必须有空格),以空格的缩进来控制层级关系,只要是左对齐的一列数据,都表示同一个层级.属性和值大小写敏感 server: port: 808 ...
- 26.异常检测---孤立森林 | one-class SVM
novelty detection:当训练数据中没有离群点,我们的目标是用训练好的模型去检测另外发现的新样本 outlier dection:当训练数据中包含离群点,模型训练时要匹配训练数据的中心样 ...