空间前方交会程序.docVIP

  1. 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
  2. 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  3. 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  4. 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  5. 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  6. 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  7. 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
空间前方交会程序

#includestdio.h #includemath.h const int n=3; void inverse (double c[n][n]); templatetypename T1,typename T2void transpose (T1*mat1,T2*mat2,int a,int b); templatetypename T1,typename T2void multi(T1*mat1,T2 * mat2,T2 * result,int a,int b,int c); templatetypename Tvoid input (T*mat,int a,int b); templatetypename Tvoid output(T*mat,char*s,int a,int b); void main() { double X=0,Y=0,Z=0; double AT[3][4],ATA[3][3],XG[3][1],ATl[3][1]; double a1,a2,a3,b1,b2,b3,c1,c2,c3,l1,l2,l3,l4,l5,l6,l7,l8,l9,l10,l11,l12,lx1,ly1,lx2,ly2; double a10,a20,a30,b10,b20,b30,c10,c20,c30; float r,w,k; printf(输入r:\n); scanf(%f,r); printf(输入w:\n); scanf(%f,w); printf(输入k:\n); scanf(%f,k); a1=cos(r)*cos(k)-sin(r)*sin(w)*sin(k); a2=-cos(r)*sin(k)-sin(r)*sin(w)*cos(k); a3=-sin(r)*cos(w); b1=cos(w)*sin(k); b2=cos(w)*cos(k); b3=-sin(w); c1=sin(r)*cos(k)+cos(r)*sin(w)*sin(k); c2=-sin(r)*sin(k)+cos(r)*sin(w)*cos(k); c3=cos(r)*cos(w); float r1,w1,k1; printf(输入r1:\n); scanf(%f,r1); printf(输入w1:\n); scanf(%f,w1); printf(输入k1:\n); scanf(%f,k1); a10=cos(r)*cos(k)-sin(r)*sin(w)*sin(k); a20=-cos(r)*sin(k)-sin(r)*sin(w)*cos(k); a30=-sin(r)*cos(w); b10=cos(w)*sin(k); b20=cos(w)*cos(k); b30=-sin(w); c10=sin(r)*cos(k)+cos(r)*sin(w)*sin(k); c20=-sin(r)*sin(k)+cos(r)*sin(w)*cos(k); c30=cos(r)*cos(w); double x0=0.0,y0=0.0,f=0.15324; double x1=-86.15,y1=-68.99,x2=0.0,y2=0.0; double Xs=393795.45,Ys=27476.46,Zs=7572.69; l1=f*b1+(x1-x0)*a3; l2=f*a1+(x1-x0)*b3; l3=f*c1+(x1-x0)*c3; lx1=f*a1*Xs+f*b1*Ys+f*c1*Zs+(x1-x0)*a3*Xs+(x1-x0)*b3*Ys+(x1-x0)*c3*Zs; l4=f*a2+(y1-y0)*a3; l5=f*b2+(y1-y0)*b3; l6=f*c2+(y1-y0)*c3; ly1=f*a2*Xs+f*b2*Ys+f*c2*Zs+(y1-y0)*a3*Xs+(y1-y0)*b3*Ys+(y1-y0)*c3*Zs; l7=f*b10+(x2-x0)*a30; l8=f*a10+(x2-x0)*b30; l9=f*c10+(x2-x0)*c30; lx2=f*a10*Xs+f*b10*Ys+f*c10*Zs+(x2-x0)*a30*Xs+(x2-x0)*b30*Ys+(x2-x0)*c30*Zs; l10=f*a20+(y2-y0)*a30; l11=f*b20+(y2-y0)*b30; l12=f*c20+(y2-y0)*c30; ly2=f*a20*Xs+f*b20*Ys+f*c20*Zs+(y2-y0)*a30*Xs+(y2-y0)*b30*Ys+(y2-y0)*c30*Zs; double A[4][3]={l1,l

文档评论(0)

xcs88858 + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

版权声明书
用户编号:8130065136000003

1亿VIP精品文档

相关文档