testflvgznew.m

simulink和matlab的图形显示接口程序

clc;
clear;
t=0:20/128000:1.2-20/128000;
f=48;
xa=10+220*sqrt(2)*cos(2*pi*f*t+pi/6)+20*sqrt(2)*cos(6*pi*f*t+pi/6);
xb=10+221*sqrt(2)*cos(2*pi*f*t+pi/6-2*pi/3)+20*sqrt(2)*cos(6*pi*f*t+pi/6+2*pi/3);
xc=10+222*sqrt(2)*cos(2*pi*f*t+pi/6+2*pi/3)+20*sqrt(2)*cos(6*pi*f*t+pi/6-2*pi/3);
data=[xa;xb;xc];
count=0;
ah=0;

data(1,1250)=-1*data(1,1250);
data(1,1253)=-1*data(1,1253);
data(1,7650)=-1*data(1,7650);
data(1,7653)=-1*data(1,7653);

for j=1:1:7679
   if data(1,j)<0
        if data(1,j+1)>=0
               ah=[ah,j];
            
         end
    end
   
end

ah=ah(2:length(ah));

count=0;i=1;
while i<=length(ah)
    if (i+4<=length(ah))&(ah(i+4)-ah(i)<20)
        count=count+1;
        bh(count)=ah(i+3);
        i=i+5;
    elseif (i+3<=length(ah))&(ah(i+3)-ah(i)<20)
        count=count+1;
        bh(count)=ah(i+2);
        i=i+4;
    elseif (i+2<=length(ah))&(ah(i+2)-ah(i)<20)
        count=count+1;
        bh(count)=ah(i+1);
        i=i+3; 
    elseif (i+1<=length(ah))&(ah(i+1)-ah(i)<20)
        count=count+1;
        bh(count)=ah(i+1);
        i=i+2; 
    else
        count=count+1;
        bh(count)=ah(i);
        i=i+1;
    end
end

i=1;count=1;
while i<length(bh)
    a=bh(i+1)-bh(i);
    if a<100
        bh(i+1:length(bh)-1)=bh(i+2:length(bh));
    else
        i=i+1;
        count=count+1;
        if count==51
            bh=bh(1:51);
            break;
        end
    end
end
    
M=(max(bh)-bh(1));
fre=50/(M/6400);


for i=1:3
for j=1:1:128
    nm=fix((j-1)*50/fre);
   mm=(((j-1)/(128*fre))-nm/6400)*6400;
 dataf(i,j)=data(i,nm+1)+mm*(data(i,nm+2)-data(i,nm+1));
end
end

[fdat1,ftangle1]=fftfcn128(data(:,1:128),1,1,20);
[fdat,ftangle]=fftfcn128(dataf(:,1:128),1,1,20);
ftangle=180*ftangle/pi;
ftangle1=180*ftangle1/pi;