draw.m

来自「画圆」· M 代码 · 共 30 行

function kalman()
alpha=0:pi/250:2*pi;%角度[0,2*pi] 
R=20;               %半径 
x=R*cos(alpha); 
y=R*sin(alpha);
subplot(1,2,1);
plot(x,y,'-') 
axis equal 
ynoi=randn(1,501); 
ynoi=ynoi/std(ynoi); 
ynoi=ynoi-mean(ynoi); 
ay=0; 
by=sqrt(0.06); 
ynoi=ay+by*ynoi; 
xnoi=randn(1,501); 
xnoi=xnoi/std(xnoi); 
xnoi=xnoi-mean(xnoi); 
ax=0; 
bx=sqrt(0.05); 
xnoi=ax+bx*xnoi; 
x1=x+xnoi;
y1=y+ynoi;
subplot(1,2,2);
plot(x1,y1,'-') 
axis equal