untitled2.m

MATLAB6P5下的GUI设计－实时参数二阶时域曲线的自动绘制.rar

DR=0,DI=0
delt=0.02/32
%u[]=0
for n=0:31
y=1.414*sin(2*pi*50*delt*n+4*pi/3)+0.2*cos(2*pi*150*delt*n)+1.414
DR=DR+y*sin(2*pi*n/32)
DI=DI+y*cos(2*pi*n/32)
temp(n+1)=DR
temp1(n+1)=DI
temp2(n+1)=y*sin(2*pi*n/32)
temp3(n+1)=y*cos(2*pi*n/32)
%波形里加一个常数不影响DR和DI的值
end
Dr=DR/32
Di=DI/32
D=sqrt(2)*sqrt(Dr*Dr+Di*Di)
ph=atan(Di/Dr)
ph1=ph+3.1415926
for n=0:31
    x(n+1)=1.414*sin(2*pi*50*delt*n+4*pi/3)+0.5*cos(2*pi*150*delt*n)+0.5*cos(2*pi*250*delt*n)+1.414
    if(x(n+1)<=0)
        x(n+1)=0
    end
    if(x(n+1)>=3.3)
        x(n+1)=3.3
    end
    
    temp4(n+1)=x(n+1)*sin(2*pi*n/32)  %这个算例的结果表明x(n+1)*sin或cos的值不会超过2*A
    temp5(n+1)=x(n+1)*cos(2*pi*n/32)  % A是基波幅值。x(n+1) 的幅值会超过2*A，但不会大于4*A
end

plot(x)
%for n=1:32
 %   u(n)=1.414*sin(2*pi*50*delt*(n-1)+pi/3)%+0.5*sin(2*pi*150*(n-1)*delt+pi/4)
 %end

%for i=3:32
 %  coswt=(u(i-2)+u(i))/(2*u(i-1))
  % sinwt=sqrt(1-coswt*coswt)
  % ph(i-2)=atan((u(i-2)*sinwt)/(u(i-1)-u(i-2)*coswt))
  % E(i-2)=sqrt(u(i-2)*u(i-2)+u(i-1)*u(i-1)-2*u(i-2)*u(i-1)*coswt)/(1.414*sinwt)
  %end