matlab求图形的聚类系数,求助,为什么画不出来聚类系数的图?一直为0啊
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
format long;
N=1000;
k=5;
L=zeros(14,20);
C=zeros(14,20);
for i=1:14
p(15-i,1)=1/2^(i-1);
end
A=zeros(N);
for i=1:N
for j=i+1:i+k
jj=j;
if j>N
jj=mod(j,N);
end
A(i,jj)=1;A(jj,i)=1;
end
end
D1=A;
D1(find(D1==0))=inf;
for i=1:N
D1(i,j)=0;
end
m=1;
while m<=N
for i=1:N
for j=1:N
if D1(i,j)>D1(i,m)+D1(m,j)
D1(i,j)=D1(i,m)+D1(m,j);
end
end
end
m=m+1;
end
L0=sum(sum(D1))/(N*(N-1)); %平均路径长度
Ci0=zeros(N,1);
for i=1:N
aa1=find(D1(i,:)==1); %寻找子图的邻居节点
if isempty(aa1)
Ci0(i)=0;
else
m1=length(aa1)
if m1==1
Ci0(i)=0;
else
B1=D1(aa1,aa1); % 抽取子图的邻接矩阵
Ci0(i)=length(find(B1==1))/(m1*(m1-1));
end
end
end
C0=mean(Ci0);
for z=1:14
for g=1:20
B=zeros(N)
for i=1:N
for j=i+1:i+k
jj=j;
if j>N
jj=mod(j,N);
end
B(i,jj)=1; B(jj,i)=1; %%生成最近邻耦合网络的邻接矩阵
end
end
for i=1:N
for j=1:k
p_rand=rand(1,1);
if p_rand
bb=randint(1,1,[1,N]);
if B(i,bb)==0&&B(bb,i)==0&&bb~=i
j2=j+i;
if j2>N
j2=mod(j2,N);
end
B(i,j2)=0;
B(j2,i)=0;
B(i,bb)=1;
B(bb,i)=1;
end
end
end
end
N1=size(A,2)
D=B;
D(find(D==0))=inf;
for i=1:N
D(i,i)=0;
end
m2=1;
while m2<=N
for i=1:N
for j=1:N
if D(i,j)>D(i,m2)+D(m2,j)
D(i,j)=D(i,m2)+D(m2,j); %
end
end
end
m2=m2+1;
end
Ci=zeros(N,1);
for i=1:N
aa=find(D(i,:)==1); %寻找子图的邻居节点
if isempty(aa)
Ci(i)=0;
else
m3=length(aa);
if m3==1
Ci(i)=0;
else
BB=D(aa,aa); % 抽取子图的邻接矩阵
Ci(i)=length(find(BB==1))/(m3*(m3-1));
end
end
end
C(z,g)=mean(Ci);
end
end
figure
LWS=mean(L,2);
CWS=mean(C,2);
semilogx(p,LWS/L0,'ro');
hold on;
semilogx(p,CWS/C0,'b*');
matlab求图形的聚类系数,求助,为什么画不出来聚类系数的图?一直为0啊相关推荐
- matlab求函数的无穷积分,[求助]贝塞尔函数的无穷积分怎么求解?
回复: [求助]贝塞尔函数的无穷积分怎么求解? 对于这种积分格式,上限是正无穷,下限是zero,我编了个函数文件,回头你自己把函数的表达式代入就可以了. function I = IntGaussLa ...
- MATLAB求图片两圆圆心,求助:如何求此图中两圆的圆心距?
对于你给出的图片,我根据之前提到的思路编写了以下程序.感觉还蛮准的.这个程序只能用于两圆相交的情况.当小圆在大圆内部时,需要在其中加以判定.仅供参考. clear;clc c = imread('1. ...
- Matlab 求纳什效率系数(NSE)
Matlab求NSE系数 function nse = NSE(sim, obs) %输入参数分别为模拟值.实测值(序列) ave_obs = sum(obs(:))/numel(obs); %实测数 ...
- 利用MATLAB求系统响应
例子 给定一个低通滤波器的差分方程为: 输入信号: a) 分别求出系统对x1和x2的响应序列,并画出其波形. b) 求出系统的单位冲响应,画出其波形. (3)给定系统的单位脉冲响应为: 用线性卷积法分 ...
- matlab求xk符号解,Matlab学习手抄(6)——非线性方程的数值解法
1. 二分法 定理: 函数f(x)在区间[a,b]上连续,如果f(a)*f(b)<0,则f(x)=0在[a,b]内至少有一个实根x'. [a,b]对半平分为两个小区间,判断f改变正负号的子区间, ...
- 【 MATLAB 】使用 MATLAB 求某输入的稳态响应
稳态响应是在本科,更确切地说是在考研时候经常接触的,例如电路基础.信号与系统中就常有稳态响应这个词. 今天在补习MATLAB数字信号处理时,突然碰到稳态响应,一脸懵逼,感觉并不理解什么是稳态响应了,也 ...
- matlab求二元函数极值算法_最优化计算与matlab实现(3)——进退法
参考资料 <精通MATLAB最优化计算(第二版)> 数值实现 Matlab 2019a 目录 石中居士:最优化计算与Matlab实现--目录zhuanlan.zhihu.com 进退法 ...
- matlab求两向量夹角_高考数学解题技巧:专题21 平面向量中最值、范围问题高中数学黄金解题模板...
↑ 点击上方"红色笔记"关注我们 高中各科解题技巧下载 回复数字:800009 高考数学解题技巧:专题24 数列求和方法高中数学黄金解题模板 高考数学解题技巧:专题26 含参不等 ...
- matlab求解不等式方程组解集,matlab求不等式方程组的可行解
求高手用matlab解一个三元的方程组 我算的你这个解都是无穷大.具体过程如下:第一步:建立M文件函数myf.mfunctionf=myf(x)f(1)=600+x(1)*cos(10)+x(2)*s ...
最新文章
- 第六篇:Feedforward Networks 前向网络
- 函数图像在图形计算机的应用,浅谈图形计算器在高中函数教学中的应用
- Spring之DAO二
- hbase实践之写流程拾遗
- The test form is only available for requests from the local machine解决方法
- jAvA中deprecate,在Java中使用Deprecated方法或类是错误的吗?
- netstat查询端口被占用
- AngularJS与服务器交互(4)
- PDF加页码怎么设置?这里有你想知道的答案
- LVGL使用华为鸿蒙字体
- 51单片机之共阳极静态数码管
- STM32F103C8T6+LD3320语音识别模块智能灯控
- GraphX入门简介
- ubuntu安装和使用labellmg(含pyqt4的安装教程)
- 茅台抢购特产平台(葫芦娃)
- 小学数学解题思维窍门
- pycharm参数设置;解决 from imutils import contours
- 【Algorithm】Prim
- JDK 8的新特性-Lambda表达式 精品文章总结
- 焚烧系统里的木马 教你彻底清除木马病毒
热门文章
- Acwing第 22 场周赛【未完结】
- core identity mysql_Asp.Net Core Identity 4 改成 MySql/MariaDB
- DS1819 对应版本的FFMPEG_Qt音视频开发33-ffmpeg安卓版
- 今天悄悄的给你说几个HashCode的破事。
- 牛X!Eclipse 开始支持 Java 14~
- org.activiti.bpmn.exceptions.XMLException: src-resolve: 无法将名称 'extension' 解析为 'element declaration'
- Java11安装及环境变量设置
- java网络编程(七)
- 算法-Valid Anagram
- Android以当前Activity为基准进行截屏