chuli.m

数模中的钻井问题的matlab源码

function z=chuli(x,y,xita)
X=x+y.*i;
X=X.*(cos(xita)+sin(xita)*i);
a=real(X);
b=imag(X);
a=mod(a,1);
b=mod(b,1);
X=a+b.*i;
for m=1:12
    if a(m)<.1&b(m)<.1
        X=[X a(m)+1+b(m)*i a(m)+(b(m)+1)*i a(m)+1+(b(m)+1)*i];
    elseif a(m)<.1
        X=[X a(m)+1+b(m)*i];
    elseif b(m)<.1
        X=[X a(m)+(b(m)+1)*i];
    end
end
z=X;