syms gamma_s alpha_s psi_a psi_t tau_1 mu nMN nMP
gamma_s = sym(45)*(pi/sym(180))%太阳方位角(单位:度)
alpha_s = sym(30)*(pi/sym(180))%太阳高度角(单位:度)
psi_a = sym(0)%方向轴倾斜方位角(单位:度)
psi_t = sym(0)%方向轴倾角(单位:度)
tau_1 = sym(0)%俯仰轴与方向轴垂直度(单位:度)
mu = sym(0)%电池板与俯仰轴平行度(单位:度)
nMN = [sym(1) sym(0) sym(0)]
nMP = 型盯返[sym(0) sym(0) sym(1)]
iS = [cos(alpha_s).*sin(gamma_s) ,cos(alpha_s).*cos(gamma_s), sin(alpha_s)]
R_1 = 则猛[cos(gamma) , -sin(gamma) , sym(0) sin(gamma) , cos(gamma) , sym(0) sym(0) , sym(0) , sym(1)]
R_2 = [cos(psi_a) , -sin(psi_a) , sym(0) sin(psi_a) , cos(psi_a) , sym(0) sym(0) , sym(0) , sym(1)]
R_3 = [cos(psi_t) , sym(0) , sin(psi_t) sym(0) , sym(1) , sym(0) -sin(psi_t) , sym(0) , cos(psi_t)]
R_4 = [cos(tau_1) , sym(0) , sin(tau_1) sym(0) , sym(1) , sym(0) -sin(tau_1) , sym(0) , cos(tau_1)]
R_5 = [sym(1) , sym(0) , sym(0) sym(0) , cos(alpha) , -sin(alpha) sym(0) , sin(alpha) , cos(alpha)]
R_6 = [cos(mu) , -sin(mu) , sym(0) 卜饥 sin(mu) , cos(mu) , sym(0) sym(0) , sym(0) , sym(1)]
R = R_1 * R_2 * R_3 * R_4 * R_5 * R_6
iMN = R * nMN
iMP = R * nMP
s = solve([iMN'*iS'==0, iMP'*iS'==0],[gamma,alpha])
错误的最大原因是, 最后一行, 你的向量相乘, 你应该是想要一个数!!! 不是一个矩阵
我把所有的东西弄成了sym, 有没有必要, 我不知道, 你自己可以去掉试试
最后别问我为什么有两个一模一样的答案, 我不知道, 你自己弄吧
????在地球上任何地点,只要有只表,有根竹竿,一根卷尺,就可知道当地经纬度。但表必须与陪绝该国标准时校对。????方法如下: 1、先迟乱前算两分日????比如在中国某地,杆影最短时是中午13点20分,且杆长与影长之比为1,则可知该地是北纬45°(tgα=1),东经100°(从120°里1小时减15°,4分钟减1°)杆长与影长之比需查表求α,码清这里用了特殊角。????2、再算两至日经度的算法不变 纬度在北半球冬至α+23.5°,夏至α-23.5°在任意一天加减修正值即可。 ????3、修正值算法:就是距两分或两至日的天欢迎分享,转载请注明来源:内存溢出
评论列表(0条)