2012年全国数学建模竞赛优秀选

附录1计算运动轨迹程 ,- .631471,- plot3(Y(:,1),Y(:,3),Y(:,5),'r*'),grid%t0=50;%%初始时间为50秒k0=(t0-0)/bc+1;%50秒标号tend=170;%%初始时间为50秒提取50秒到170秒初始坐标状态，间隔0.2秒，依次为X,Y,Zdy(6)=-附录二步长取值分析 -5379.544693-407.095342 .715--407.095342%%Y0(1:(k+1):hh,5)).^2);附录间配准程%%把06号观测时间统一到50:0.2:170clear;clc;closeallfprintf('拟合求解06号A的一元二次多项式系数') 预测%%对比预测结果准确plot(zs6(:,1),YCa6,'ro'),holdplot(zs6(:,1),zs6(:,2),'g*'),gridxlabel('时间plot(zs6(:,1),YCb6,'r*',zs6(:,1),zs6(:,3),'g*'),holdontitle('06号b=ZS/XS预测值和真实值对比图 求解50:0.2:170a和bYCA6=Aa6(1)*t.^2+Aa6(2)*t+Aa6(3);%a预测 预测YC6T1AB=[t',ones((170-50)/bc+1,1),YCA6',YCB6'];%%依次为时间、1,、a、YC6_1AB=[ones((170-50)/bc+1,1),YCA6',YCB6'];%%依次为时间、1,、a、%%对比预测结果准确%%%%title('06号50:bc:170秒a预测值和真实值对比图%%%title('06号50:bc:170秒b预测值和真实值对比图%%%%为把09号观测时间统一到%%zs9为09号数据真实fprintf('拟合求解09号A的一元二次多项式系数') 预测%%对比预测结果准确plot(zs9(:,1),zs9(:,2),'g*'),gridxlabel('时间plot(zs9(:,1),zs9(:,3),'g*'),gridxlabel('时间title('09号b=ZS/XS预测值和真实值对比图 求解50:0.2:170a和bYCA9=Aa9(1)*t.^2+Aa9(2)*t+Aa9(3);%a预测 预测YC9_1AB=[ones((170-50)/bc+1,1),YCA9',YCB9'];%%依次为1,、a、b对比预测结果准确%%%%title('06号50:bc:170秒a预测值和真实值对比图%%%title('06号50:bc:170秒b预测值和真实值对比图%%%%06号三维位 ,- .771852,- %t0=50;%%初始时间为50秒k0=(t0-0)/bc+1;%50秒标号tend=170;%%初始时间为50秒title('06号50:0.2:170秒三维位置图提取50秒到170秒初始坐标状态，间隔0.2秒，依次为T,X,Y,Z%%09号三维位 ,- .631471,- %t0=50;%%初始时间为50秒k0=(t0-0)/bc+1;%50秒标号tend=170;%%初始时间为50秒提取50秒到170秒初始坐标状态，间隔0.2秒，依次为T,X,Y,Z附录四逐点交汇定轨%%%%%%坐标转DATA96=DATA9(:,2:4)-YC9_1AB=[ones((170-50)/bc+1,1),YCA9',YCB9';];%%依次 1,、a、K=(tend-%%当输入YC9_1AB(k,:)XS后，C9表示观测飞行器转换后的坐标，A9表示飞行器转换后的矩阵Afork=1:K%%当输入YC6_1AB(k,:)XS后，C6表示观测飞行器转换后的坐标，A6表示飞行器转换后的矩阵Afork=1:K%%XS6b的第一行，XS9bfork=1:K%%求解6号和9号的观测值校正函%C69=CC6-CC9;%%坐标差%%%YCAX=AX(1)*tt.^4+AX(2)*tt.^3+AX(3)*tt.^2+AX(4)*tt.^1+AX(5);%X预测%%YCAY=AY(1)*tt.^4+AY(2)*tt.^3+AY(3)*tt.^2+AY(4)*tt.^1+AY(5);%Y%%YCAZ=AZ(1)*tt.^4+AZ(2)*tt.^3+AZ(3)*tt.^2+AZ(4)*tt.^1+AZ(5);%Z预测%aaa=[YCAX,YCAY,YCAZ]-附录五坐标转换调用function%%要求不能取YZ%%a为坐标，b为观测点相对坐标%%输出c为观测点相对地球坐标function%%c为观测点相对地球坐标,a为坐%%输出b为观测点相对坐标附录六基于多项式模型计算飞行轨迹的%%CC0为空间飞行物的坐fori=1:length(CC0(:,1))-1fori=1:length(dd(:,1))-1fori=1:length(CC0(:,1))-1fori=1:length(cos)-m(i)=F.*(-g*cos(i)+sqrt(g^2*cos(i)^2+sqrt(sum(a(i,:).^2))-g^2))./(sqrt(sum(a(i,:).^2))-%%CC0的初始数%%50秒对应的数%[T,Y]=ode45('weifen2',[50.4:bc:170],[- %%50.4秒对应的数 fprintf('基于多项式的累积误差为wucha2=sum((Y(:,1)-CC0(3:601,1)).^2)+sum((Y(:,3)-CC0(3:601,2)).^2)+sum((Y(:,5)-%%反求a，b,K%%BB9,BB6分别为9号和6号的观测到的飞行物的坐标分别相对于9号和号的坐%%CC9,CC6分别为9号和6号的观测到的飞行物的坐标分别相对于9号和fork=1:K 附录七基于指数模型计算飞行fori=1:length(CC0(:,1))-1fori=1:length(dd(:,1))-1fori=1:length(CC0(:,1))-1fori=1:length(cos)-m(i)=F.*(-g*cos(i)+sqrt(g^2*cos(i)^2+sqrt(sum(a(i,:).^2))-g^2))./(sqrt(sum(a(i,:).^2))-plot(shijian,m,'b-',shijian,y,'r-[T,Y]=ode45('weifen9',[50.4:bc:170],[- plot3(Y(:,3),Y(:,1),Y(:,5),'r.'),holdwucha1=sum((Y(:,1)-CC0(3:601,1)).^2)+sum((Y(:,3)-CC0(3:601,2)).^2)+sum((Y(:,5)-附录八基于指数模型系统误差%%6号数据X=[ones(599,1),-X1=[ones(599,1),-附录星测轨程%%A99为9号观测值A，B，weixing为9号三维坐%%校正%%坐标转DATA96=DATA9-YC9_1AB=[ones(K-1,1),A99(1:K-1,1),A99(1:K-1,1)];%%依次 1,、a、YC6_1AB=[ones(K-1,1),A99(2:K,2)