newton_method_btcr.m

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

delta=1/128
    x0=pi/2
 i=1
for btcr=1:-delta:0
    btcrxq(i)=btcr*2^15     % Q15
    btcrx(i)=btcr
    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
t1(i)=10*(alpha(i)-pi/2)/pi   % 转换为时间，单位毫秒，从电压峰值开始计时
alphaxq(i)=alpha(i)*2^13        % Q13
b(i)=(2*pi-2*x1+sin(2*x1))/pi   % 验算用
i=i+1
%for n=1:8
%    x1=x0-(x0*exp(x0)-1)/(exp(x0)+x0*exp(x0))
%    x0=x1
%    end
end
for n=2:129
    temp(n)=b(n)-b(n-1);     % 结果是正确的
    separation(n)=(alpha(n)-alpha(n-1))*180/pi;
end