该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
这个是代码：
syms x y z
[deltafd]=xlsread('附件2.xls','FDOA','B2:B290');
[x1]=xlsread('附件3.xls','主星星历','B2:B290');
[x2]=xlsread('附件3.xls','邻星星历','B2:B290');
[y1]=xlsread('附件3.xls','主星星历','C2:C290');
[y2]=xlsread('附件3.xls','邻星星历','C2:C290');
[z1]=xlsread('附件3.xls','主星星历','D2:D290');
[z2]=xlsread('附件3.xls','邻星星历','D2:D290');
[v1x]=xlsread('附件3.xls','主星星历','E2:E290');
[v2x]=xlsread('附件3.xls','邻星星历','E2:E290');
[v1y]=xlsread('附件3.xls','主星星历','F2:F290');
[v2y]=xlsread('附件3.xls','邻星星历','F2:F290');
[v1z]=xlsread('附件3.xls','主星星历','G2:G290');
[v2z]=xlsread('附件3.xls','邻星星历','G2:G290');
f=6.372-2.225;
re=6378245;
c=299792458;
g=solve('deltafd=(1/lambda)*(((x-x2)*v2x+(y-y2)*v2y+(z-z2)*v2z)/r2)-(((x-x1)*v1x+(y-y1)*v1y+(z-z1)*v1z)/r1)','lambda=c/f','deltar=r2-r1','re=sqrt(x^2+y^2+z^2)','r2=sqrt((x-x2)^2+(y-y2)^2+(z-z2)^2)','r1=sqrt((x-x1)^2+(y-y1)^2+(z-z1)^2)','x1,y1,z1,x2,y2,z2,v1x,v2x,v1y,v2y,v1z,v2z,deltafd,r1,r2,deltar')
)
