| 查看: 524 | 回复: 1 | ||
[求助]
MATLAB元包机编程算法出错求指导
|
|
clc;clear all; %初始化参数 N=100; % 矩阵是一个 N×N的 T=1000; % 参数值 S=500; % 元胞的可能取值为1-500 Q=300000; % 参数值 R=831; % 参数值 p0=exp(-Q/(R*T)); % 计算概率P0值 kk=1000; % 循环次数即总的时间步数 A=round(rand(N)); %生成随机整数矩阵499=500-1 imh = image(cat(3,A,A,A)); ti=title('CAS=0', 'Fontsize',14,'Fontname','Times New Roman' );%显示标题 % A1=round(rand(n)); % 循环部分 for k=100:110; k Am=zeros(N+2); %生成一个全为0的矩阵Am Am(2:end-1,2:end-1)=A; %将A赋值给Am的中间部分 Am(1,2:end-1)=A(end,  ; %将A最下面一行补到矩阵Am的上边Am(end,2:end-1)=A(1, ; %将A第一行补到矩阵Am的下面Am(2:end-1,1)=A(:,end); %将A最右边一列补到矩阵Am的最左边 Am(2:end-1,end)=A(:,1); %将A的第一列补到矩阵Am的最右边 Am(1,1)=A(end,end); %将A的右下角补到Am左上角 Am(end,end)=A(1,1); %将A的左上角补到Am右下角 Am(end,1)=A(1,end); %将A的右上角补到Am左下角 Am(1,end)=A(end,1); %将A的左下角补到Am右上角 rand(1); %生成一个0-1之间的随机数，改 A1= Am(1:end-2,1:end-2); %A5的左上近邻 A2= Am(2:end-1,1:end-2); %A5的上近邻 A3 = Am(3:end,1:end-2); %A5的右上近邻 A4 = Am(1:end-2,2:end-1); %A5的左近邻 A6 = Am(3:end,2:end-1); %A5的右近邻 A7 = Am(1:end-2,3:end); %A5的上近邻 A8 = Am(2:end-1,3:end); %A5的下近邻 A9= Am(3:end,3:end); %A5的右下近邻 A5=Am(end-1,end-1); %Am的中间部分A5 % index=round(1+7* rand(1)); %生成1-8的参数值便于提取 index=round(1+7* rand(1)); B={A1;A2;A3;A4;A6;A7;A8;A9}; if A1==A2==A3==A4==A6==A7==A8==A9==A5& (rand(1)<P0); A5=A1; if A2==A4==A6|A4==A6==A8| A2==A6==A8& (rand(1)<P0); A5=A2; elseif A1==A3==A7|A1==A3==A9|A1==A7==A9 & (rand(1)<P0); A5=A1; else A5=B(index)& (rand<P0); %改 end A=A5; % set(Ii, 'CData',[]); %更新灰度值的显示 % set(ti,'string',['CAS=',num2str(k)]); %显示更新的CAS时间步数 set(imh, 'cdata', cat(3,A,A,A) ); pause(1); %暂停1秒，显示动画效果 end end @beefly |
回复此楼» 猜你喜欢
【复旦大学】二维材料方向招收2026年博士研究生1名
已经有0人回复
-
北京纳米能源与系统研究所 王中林院士/曹南颖研究员课题组2026级硕/博/博后招生
已经有10人回复
-
物理学I论文润色/翻译怎么收费?
已经有97人回复
-
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
已经有23人回复
-
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
2楼2017-01-16 09:11:07