trans.m

camlaw的ｍａｔｌａｂ程序 [e-mail]

function [xer,yer]=trans(xe,ye,lab,fi1,fi2)

N=size(xe);
N(1)=[];
for n=1:N
    xe(n)=xe(n)-lab*cos(fi1(n));
    ye(n)=ye(n)-lab*sin(fi1(n));
    
    xer(n)=xe(n)*cos(fi2(n))+ye(n)*sin(fi2(n));
    yer(n)=xe(n)*sin(-fi2(n))+ye(n)*cos(fi2(n));
end