| 查看: 553 | 回复: 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 |
回复此楼» 猜你喜欢
桂林理工大学物理学专业招收调剂
已经有18人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有67人回复
-
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
-
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）
已经有6人回复
-
【修正版】单环用磁约束低速电子实验方案（简化版）
已经有0人回复
-
桂林理工大学物理学专业招收调剂，还有三个名额！！！
已经有20人回复
-
考博自荐
已经有1人回复
-
山东大学第二批博士研究生招生
已经有0人回复
2楼2017-01-16 09:11:07
5