btcr.m

来自「MATLAB6P5下的GUI设计－实时参数二阶时域曲线的自动绘制.rar」· M 代码 · 共 16 行

% 计算BTCR，忽略w和L。Btcr=(2*pi-2*a+sin2a)/pi

delwt=pi/256
alpha=pi/2
% wt=90~180,共128点
for n=1:128
    btcr1(n)=(2*pi-2*alpha+sin(2*alpha))/pi
     btcrq(n)=round(btcr1(n)*32768)     % Q15
    alpha=alpha+delwt
    if(n>1)
        delta(n-1)=btcrq(n)-btcrq(n-1)
    end
end
%plot(btcrq)
plot(btcr1,'k')