- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 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
您可能关注的文档
最近下载
- 新人教版高中数学必修第二册统计全套课件.pptx VIP
- 台球厅消防安全应急预案.docx VIP
- 海外代理协议合同协议.docx VIP
- 初中教科研课题:《初中语文预习方法研究》课题研究工作报告.doc VIP
- 2025至2030年中国新疆维吾尔自治区建筑市场运行态势及行业发展前景预测报告.docx
- 简述10KV 高压配电柜安装.doc VIP
- GB50148-2010 电气装置安装工程电力变压器油浸电抗器、互感器施工及验收规范.pdf VIP
- 2025航天恒星科技有限公司招聘80+人笔试历年参考题库附带答案详解.pdf
- RB∕T 174-2021 司法鉴定法庭科学机构能力专业要求.pdf
- CP-717安装指南.doc VIP
文档评论(0)