马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
dt=0.0001; T=15; t=0:dt:T;
g=9.8; L=1; m=10; k=5; theta0=3.1; v0=0;
v=zeros(size(t)); theta=zeros(size(t));
v(1)=v0; theta(1)=theta0;
for n=1:length(t)
v(n+1)=v(n)+(g*sin(theta(n))-k./m.*v(n)).*dt;
theta(n+1)=theta(n)-1./L.*v(n).*dt;
end
[AX,H1,H2]=plotyy(t,v(1:length(t)),t,theta(1:length(t)),'plot');
set(H1,'linestyle','-'); set(H2,'linestyle','-.');
set(AX(2),'ylim',[-6 6]);
set(get(AX(1),'ylabel'),'string','线速度v(t) m/s');
set(get(AX(2),'ylabel'),'string','角位移 \theta(t)/rad'); xlabel('时间 t/s');
legend(H1,'线速度 v(t)',2); legend(H2,'角位移 \theta(t)',1);
(我的电脑有问题,传不了图片,抱歉)
但我想要的是右侧的y轴显示的是-6:2:6 (就像左侧y轴显示的刻度一样)
[ 本帖最后由 ChaChing 于 2009-5-30 20:21 编辑 ] |