matlab大作业Word文档格式.docx
- 文档编号:22276726
- 上传时间:2023-02-03
- 格式:DOCX
- 页数:9
- 大小:74.46KB
matlab大作业Word文档格式.docx
《matlab大作业Word文档格式.docx》由会员分享,可在线阅读,更多相关《matlab大作业Word文档格式.docx(9页珍藏版)》请在冰豆网上搜索。
wz0=20;
%偏航轴初始旋转角速度
tf=2*pi/((Ix/It-1)*wx0);
ts=[0,tf];
w0=[wx0;
wy0;
wz0];
%初始条件
[t,w]=ode45('
dwdt'
ts,w0);
[m,n]=size(w);
plot3(w(:
2),w(:
3),w(:
1));
axis([-30,30,-30,30,0,60]);
holdon
gg(w);
theta=atan((It/Ix)*(sqrt(wy0^2+wz0^2)/wx0));
%beta=atan(w(:
3)./w(:
2));
beta=0.35;
gama=atan(sqrt(wy0^2+wz0^2)/wx0);
k1=tan(theta-gama)/tan(gama);
k2=cos(theta-gama)/cos(gama);
H(:
1)=k2*w(:
1);
2)=k2*k1*w(:
2);
3)=k2*k1*w(:
3);
gg1(H);
xlabel('
z'
),ylabel('
y'
),zlabel('
x'
)
gridon,holdon
x=0:
80;
y=0*x;
z=0*x;
x1=0:
y1=0*x1;
%z=(sqrt(wy0^2+wz0^2)/wx0)*x;
zz=(It/Ix)*(sqrt(wy0^2+wz0^2)/wx0)*x1;
%plot3(-z,y,x);
%plot3(y1,y1,x1,'
--'
);
%plot3(y,z,x);
%plot3(y,-z,x);
plot3(y,z,x,'
m-'
'
Linewidth'
plot3(y1,zz,x1,'
holdon%动量矩H的方向
text(0,10,70),text(30,0,70),text(-10,10,60),holdoff
其中包含如下函数体:
(1)functionomegad=dwdt(t,w)
globalIxItwx0Omega
Omega=(Ix-It)*wx0/It;
wx0=w
(1);
wy=w
(2);
wz=w(3);
omegad=[0;
-Omega*wz;
Omega*wy];
(2)functiongg(a)
gama=atan(sqrt(a(1,2)^2+a(1,3)^2)/a(1,1));
fori=0:
0.1:
1
b=i*a;
plot3(b(:
2),b(:
3),b(:
end
[m,n]=size(a);
forj=1:
30
k=fix(1+(j-1)*0.5*m/30);
alpha=atan(a(k,2)/a(k,3));
x=0:
a(1,1);
y=sin(alpha)*tan(gama)*x;
z=cos(alpha)*tan(gama)*x;
plot3(y,z,x);
plot3(-y,-z,x);
(3)functiongg1(a)
A=[cos(-theta),0,-sin(-theta);
0,1,0;
sin(-theta),0,cos(-theta)];
%转换矩阵
b=i*a*A;
1),'
g'
XX(:
1)=x'
;
XX(:
2)=y'
3)=z'
YY(:
YY(:
2)=-y'
3)=-z'
XX=XX*A;
YY=YY*A;
plot3(XX(:
2),XX(:
3),XX(:
c'
plot3(YY(:
2),YY(:
3),YY(:
逆进动
2、正进动
Ix=60;
It=80;
正进动
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- matlab 作业