📄 cacmack2.m
字号:
m=8;
nb=8;
zz=(m*(nb-1)+1)^2;
t=m*nb^2;
w=zeros(1,t);
f=zeros(1,t);
cycles=26;
%x=-3.14:6.28/64:3.14;
%y=x ;
%[sx,sy]=size(x);
sx=57;
sy=57;
for p=1:cycles
p1(p)=0;
for i=1:sy
for j=1:sy
x(i)=-3.14+(i-1)*3.14/28;
y(j)=-3.14+(j-1)*3.14/28;
%z(i,j)=((x(i)-3)^2+(y(j)-3)^2)^0.5;
z(i,j)=(sin(x(i)))*(cos(y(j)));
sk(1)=1;
for j1=1:m
i1=m-j1+1;
ax=fix((i+i1)/m);
ay=fix((j+i1)/m);
sk(j1)=ax+ay*nb+(j1-1)*nb^2+1;
end
for j1=1:m
w1(j1)=w(sk(j1));
f(sk(j1))=f(sk(j1))+1;
d(i,j)=sum(w1);
end
for j2=1:m
f1(j2)=1/((f(sk(j2)))^1.0);
fh=sum(f1);
end
for j1=1:m
w2(sk(j1))=(1/((f(sk(j1))^1.0)))*(z(i,j)-d(i,j))/fh;
w(sk(j1))=w(sk(j1))+w2(sk(j1));
end
%p1(p)=p1(p)+abs(d(i,j)-z(i,j));
end
end
i
j
%p1(p)=p1(p)+abs(d-z);
ter(p)=sum(sum(abs(d(i,j)-z(i,j))));
%rt(p)=sum(sum(abs(d(i,j)-z(i,j))))/zz;
end
p=1:cycles;
plot(p,ter,'g');
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -