untitled4.m

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

delta=1/128
    x0=pi/2
 i=1
for btcr=1:-delta:0
    btcrxq(i)=btcr*2^15     % Q15
    for n=1:16
    x1=x0-(2*pi-2*x0+sin(2*x0)-pi*btcr)/(2*cos(2*x0)-2)
    x0=x1
    end
alpha(i)=x1
alf1(i)=alpha(i)-pi/2       % 以电压峰值开始计时
T(i)=round((alf1(i)*10e-3)/(pi*8e-7))
alf1xq(i)=round(alf1(i)*2^13)          % Q13

i=i+1
end