clear all;
clc;
t=0:0.01:4*pi;
w=1;
n=4;
fs=zeros(1,n);
a=zeros(1,n);
b=zeros(1,n);
c=[1 0 0 1];
d=[1 1 0 0];
for k=1:n
fs(k)=panduanab11(c(k),d(k));
a(k)=cos(fs(k));
b(k)=sin(fs(k));
m=a(k)*cos(w*t)-b(k)*sin(w*t);
if k==1
h1=m;
elseif k==2
h2=m;
elseif k==3
h3=m;
else h4=m;
end
end
plot(t,h1,t,h2,t,h3,t,h4);
E=1;
cc=zeros(1,n);
dd=zeros(1,n);
aa=zeros(1,n);
bb=zeros(1,n);
for k=1:n
aa(k)=E*cos(fs(k));
bb(k)=E*sin(fs(k));
if aa(k)>0
cc(k)=1;
else cc(k)=0;
end
if bb(k)>0
dd(k)=0;
else dd(k)=1;
end
end
c
d
fs
cc
dd