ball.m

来自「matlab的两个小球碰撞的模拟程序」· M 代码 · 共 30 行

function ball()
x=64;
y=64;
a=6;
b=7;
r=10;
for i=1:200
    B=ones(128)*256;
    for j=1:128
        for k=1:128
            if sqrt((x-j)^2+(y-k)^2)<r
                B(j,k)=0;
            end
        end
    end
    imshow(B);
    %%pause;close;
    x=x+a;
    y=y+b;
    if (x<=r)
        a=abs(a);
    elseif (x>=128-r)
        a=abs(a)*(-1);
    elseif (y<=r)
        b=abs(b);
    elseif (y>=128-r)
        b=abs(b)*(-1);
    end
end