function park0try(fname,srate,nc)
%park vector and its extended method
%fname:filename
%srate:sample rate(kHz)
%nc:date select,0-only 3 columns;1-first 3 column of 4 columns;2-back 3
%column of 4 column
fid=fopen(fname);
if abs(nc)<=0.001
[a1,count]=fscanf(fid,'%g,%g,%g',[3,inf]);
ct=3;a=a1';
elseif abs(nc-1)<=0.001
[a1,count]=fscanf(fid,'%g,%g,%g,%g',[4,inf]);
a1=a1';a=a1(:,1:3);ct=4;
elseif abs(nc-2)<=0.001
[a1.count]=fscanf(fid,'%g,%g,%g,%g',[4,inf]);
a1=a1';a=a1(:,2:4);ct=4;
end
fclose(fid);
count
if count/ct<=5000
npoint=12;
end
if (count/ct<=10000)+(count/ct>=5000)
npoint=13;
end
if (count/ct>=10000)+(count/ct<=20000)
npoint=14;
end
aa=a(:,1);bb=a(:,2);cc=a(:,3);
npoint
np=2^npoint
srate
nc
t=(0:1:count/ct-1)/srate;
t=t';
id=2/3*(aa.*cos(100*pi*t)+bb.*cos(100*pi*t-2/3*pi)+cc.*cos(100*pi*t+2/3*pi));
iq=-2/3*(aa.*sin(100*pi*t)+bb.*sin(100*pi*t-2/3*pi)+cc.*sin(100*pi*t+2/3*pi));
% id=sqrt(2)/sqrt(3)*aa.*cos(100*pi*t)-1/sqrt(6)*bb.*cos(100*pi*t-2/3*pi)+1/sqrt(6)*cc.*cos(100*pi*t+2/3*pi);
% iq=1/sqrt(2)*bb.*cos(100*pi*t-2/3*pi)-1/sqrt(2)*cc.*cos(100*pi*t+2/3*pi);
figure;
plot(id,iq);title('PARK图像');
评论2