testxiangxu.m

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

clc;
clear;

t=0:20/128000:0.02-20/128000;
xa=10+220*sqrt(2)*cos(2*pi*50*t+pi/6)+20*sqrt(2)*cos(2*pi*50*t+pi/6);
xb=10+220*sqrt(2)*cos(2*pi*50*t+pi/6-2*pi/3)+20*sqrt(2)*cos(2*pi*50*t+pi/6+2*pi/3);
xc=10+220*sqrt(2)*cos(2*pi*50*t+pi/6+2*pi/3)+20*sqrt(2)*cos(2*pi*50*t+pi/6-2*pi/3);
data=[xa;xb;xc];
[fdat,ftangle]=fftfcn128(data,1,1,20);

jk=1;
for i=1:1
    u01(i,jk)=(fdat(((i-1)*3+1),2)*cos(ftangle(((i-1)*3+1),2))+fdat(((i-1)*3+2),2)*cos(ftangle(((i-1)*3+2),2))+fdat(((i-1)*3+3),2)*cos(ftangle(((i-1)*3+3),2)))/3;
    u02(i,jk)=(fdat(((i-1)*3+1),2)*sin(ftangle(((i-1)*3+1),2))+fdat(((i-1)*3+2),2)*sin(ftangle(((i-1)*3+2),2))+fdat(((i-1)*3+3),2)*sin(ftangle(((i-1)*3+3),2)))/3;
   u0(i,jk)=sqrt(u01(i,jk)*u01(i,jk)+u02(i,jk)*u02(i,jk))
   if u01(i,jk)>0
       aa0(i,jk)=180*atan( u02(i,jk)/u01(i,jk))/pi;
   else
       aa0(i,jk)=180+180*atan( u02(i,jk)/u01(i,jk))/pi;
   end
   u11(i,jk)=(fdat(((i-1)*3+1),2)*cos(ftangle(((i-1)*3+1),2))+fdat(((i-1)*3+2),2)*cos(ftangle(((i-1)*3+2),2)+2*pi/3)+fdat(((i-1)*3+3),2)*cos(ftangle(((i-1)*3+3),2)+4*pi/3))/3;
   u12(i,jk)=(fdat(((i-1)*3+1),2)*sin(ftangle(((i-1)*3+1),2))+fdat(((i-1)*3+2),2)*sin(ftangle(((i-1)*3+2),2)+2*pi/3)+fdat(((i-1)*3+3),2)*sin(ftangle(((i-1)*3+3),2)+4*pi/3))/3;
   u1(i,jk)=sqrt(u11(i,jk)*u11(i,jk)+u12(i,jk)*u12(i,jk))
   if u11(i,jk)>0
       aa1(i,jk)=180*atan( u12(i,jk)/u11(i,jk))/pi;
   else
       aa1(i,jk)=180+180*atan( u12(i,jk)/u11(i,jk))/pi;
   end
   u21(i,jk)=(fdat(((i-1)*3+1),2)*cos(ftangle(((i-1)*3+1),2))+fdat(((i-1)*3+2),2)*cos(ftangle(((i-1)*3+2),2)+4*pi/3)+fdat(((i-1)*3+3),2)*cos(ftangle(((i-1)*3+3),2)+2*pi/3))/3;
   u22(i,jk)=(fdat(((i-1)*3+1),2)*sin(ftangle(((i-1)*3+1),2))+fdat(((i-1)*3+2),2)*sin(ftangle(((i-1)*3+2),2)+4*pi/3)+fdat(((i-1)*3+3),2)*sin(ftangle(((i-1)*3+3),2)+2*pi/3))/3;
   u2(i,jk)=sqrt(u21(i,jk)*u21(i,jk)+u22(i,jk)*u22(i,jk))
   if u21(i,jk)>0
       aa2(i,jk)=180*atan( u22(i,jk)/u21(i,jk))/pi;
   else
       aa2(i,jk)=180+180*atan( u22(i,jk)/u21(i,jk))/pi;
   end
end