ZB+K45 源代码

A = [1 2 3;4 5 6] % 生成2×3维矩阵 ZA = trapz(A) % 计算矩阵的1维积分 B = [4 5 6;7 8 9] % 生成2×3维矩阵 ZB = trapz(B,2) % 计算矩阵的2维积分 Y(…

clear n=5;%n=input('请输入节点数'); b=4;%b=input('请输入支路数'); %disp('请输入支路阻抗和节点功率矩阵，第一列存支路号，第二列存首节点号，第三列存尾节点号，第四列存支路自阻抗,第五列存尾节点给定功率'); …

%利用星敏感器输出矢量X、Z，确定惯性系下姿态矩阵A function A=AVV(X,Z,Xb,Zb) %X表示横轴矢量；Z轴表示光轴矢量；Xb表示横轴安装星体坐标；Zb表示光轴安装星体坐标 V2i=cross(X,Z); V2i=norm(V2i)…