wuchachuli.m

这是一些自己几个月的劳动成果

clear
R=[3 6 9 10 12 15 18 20 30 40 50 60];;
W1=[4.36 6.91 7.32 7.75 9.21 10.34 11.6 12.28 15.47 17.35 19.47 21.85];
r=3:0.01:60;
t1=interp1(R,W1,r,'spline');
figure(1);
plot(R,W1,':',R,W1,'+',r,t1)
xlabel('距离(km)','FontWeight','bold');
ylabel('误差(mils)','FontWeight','bold');
title('等效热噪声引起角误差');
R1=3:3:18;
r1=3:0.1:18;
W3=[0.2 0.2 0.2 0.2 0.2 0.2];
t3=interp1(R1,W3,r1,'spline');
figure(2);
subplot(211);plot(R1,W3,':',R1,W3,'+',r1,t3)
xlabel('(a)         距离(km)','FontWeight','bold');
ylabel('误差(mils)','FontWeight','bold');
title('伺服噪声引起角误差')
W4=[0.2 0.2 0.2 0.2 0.2 0.2];
t4=interp1(R1,W4,r1,'spline');
subplot(212);plot(R1,W4,':',R1,W4,'+',r1,t4)
xlabel('(b)        (距离(km)','FontWeight','bold');
ylabel('误差(mils)','FontWeight','bold');
W5=[0.2 0.2 0.2 0.2 0.2 0.2];
t5=interp1(R1,W3,r1,'spline');
figure(3);
subplot(211);plot(R1,W5,':',R1,W5,'+',r1,t5)
xlabel('(a)        距离(km)','FontWeight','bold');
ylabel('误差(mils)','FontWeight','bold');
title('伺服机械误差引起角误差')
W6=[0.2 0.2 0.2 0.2 0.2 0.2];
t6=interp1(R1,W6,r1,'spline');
subplot(212);plot(R1,W6,':',R1,W6,'+',r1,t6)
xlabel('(b)        距离(km)','FontWeight','bold');
ylabel('误差(mils)','FontWeight','bold');
W7=[0.72 0.44 0.38 0.38 0.42 0.49];
t7=interp1(R1,W7,r1,'spline');
figure(4);
subplot(211);plot(R1,W7,':',R1,W7,'+',r1,t7)
xlabel('(a)        距离(km)','FontWeight','bold');
ylabel('误差(mils)','FontWeight','bold');
title('随机噪声引起角误差')
W8=[0.35 0.3 0.3 0.34 0.42 0.5];
t8=interp1(R1,W8,r1,'spline');
subplot(212);plot(R1,W8,':',R1,W8,'+',r1,t8)
xlabel('(b)        距离(km)','FontWeight','bold');
ylabel('误差(mils)','FontWeight','bold');
W9=[4.44 6.97 7.34 9.17 10.36 11.62];
t9=interp1(R1,W9,r1,'spline');
figure(5);
subplot(211);plot(R1,W9,':',R1,W9,'+',r1,t9)
xlabel('(a)        距离(km)','FontWeight','bold');
ylabel('总误差（均方根误差）(mils)','FontWeight','bold');
title('综合角误差')
W10=[4.41 7.07 7.33 9.22 10.36 11.62];
t10=interp1(R1,W10,r1,'spline');
subplot(212);plot(R1,W10,':',R1,W10,'+',r1,t10)
xlabel('(b)        距离(km)','FontWeight','bold');
ylabel('总误差（均方根误差）(mils)','FontWeight','bold');