dsp_simulation.m

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

% 定点DSP模拟
DR=0,DI=0
delt=0.02/32
for n=0:31
    x(n+1)=1.414*sin(2*pi*50*delt*n+4*pi/3)+0.2*sin(2*pi*100*delt*n)+0.5*cos(2*pi*150*delt*n)+0.3*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
    xq(n+1)=x(n+1)*1023/3.3         % Q0,AD的3.3伏对应1023
    sintab(n+1)=sin(2*pi*n/32)*32768    %Q15
    costab(n+1)=cos(2*pi*n/32)*32768    %Q15
    DR=DR+xq(n+1)*sintab(n+1)/2^15       % Q0,15+0-0=15
    DI=DI+xq(n+1)*costab(n+1)/2^15
    temp(n+1)=DR
    temp1(n+1)=DI
    
    %波形里加一个常数不影响DR和DI的值
    
end
Dr=DR/32
Di=DI/32
D=sqrt(2)*sqrt(Dr*Dr+Di*Di)           % Q0
ph=atan(DI/DR)
ph1=ph+pi
plot(x)