声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 857|回复: 4

[综合讨论] 为什么运行结果为0?

[复制链接]
发表于 2008-9-9 21:30 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
被调用的函数
function ss=Rugkut3(t,y,flag,A,B,u )
ss=A*y+B*u;

运行的程序代码
clf
clear
clc;
kaxi=0.01; mu=0.02;
kaxit=0.002;  
h=0.1;
NT=50;
t=0:h:NT;
N=length(t);
ft=sqrt(1+0.5*mu)/(1+mu);
dalt=0.6;
fe=1.0;  
t0=0;
e=0.6;
zz0(:,1)=[0;0;0;0];
for I=1:N
   
     Xg(I)=cos(fe*t(I));
   if abs(zz0(2,1)-zz0(1,1))-dalt>=1*0.0000001
   A=[0 0 1 0;0 0 0 1;0 mu*ft.^2 -2*kaxi 2*mu*kaxit*ft*e;0 -(1+mu)*ft.^2 2*kaxi
   -(1+mu)*2*kaxit*ft*e];
  B=[0;0;1;0];
  u=Xg(I);
   k1=Rugkut3(t0,zz0,A,B,Xg(I));
   k2=Rugkut3(t0+h/2,zz0+(h/2)*k1,A,B,Xg(I));
   k3=Rugkut3(t0+h/2,zz0+(h/2)*k2,A,B,Xg(I));
   k4=Rugkut3(t0+h,zz0+h*k3,A,B,Xg(I));
    zz1=zz0+(1/6)*h*(k1+2*k2+2*k3+k4);
  
   else
       zz1(1,1)=zz0(1,1);
       zz1(2,1)=zz0(2,1);
       zz1(3,1)=(e+mu)/(e*(1+mu))*zz0(3,1)+mu*(e+1)/(e*(1+mu))*zz0(4,1);
       zz1(4,1)=(1+e)/(e*(1+mu))*zz0(3,1)-(mu*e-1)/(e*(1+mu))*zz0(4,1);
   
    end
    Z12(:,I)=zz1(:,1);
    t0=t0+h;
    zz0=zz1;  
end
figure(1)
         plot(t,(Z12(1,:)));
         title({'系统位移-时间关系图'})
         xlabel('t')
         ylabel(' x_1')
      hold on;

figure(3),
plot(Z12(1,:),Z12(3,:))
       title({'系统位移-速度关系图'})
       xlabel(' x_1')
       ylabel('dx_1/dt')
      
        hold on;
figure(2),
  plot(t,Z12(1,:)+dalt,'g:',t,(Z12(2,:)+Z12(1,:)),'r-',t,Z12(1,:)-dalt,'g: ');
       title({'碰撞效果图'})
       xlabel('t')
       ylabel('x_1-\delta              x_2               x_1+\delta')
       hold on
我刚学matlab,有很多问题不懂,望各位同学能指点指点,我急需要解决这个问题,谢谢!!!!!!!!!!
回复
分享到:

使用道具 举报

发表于 2008-9-10 08:31 | 显示全部楼层

回复 楼主 dearchina 的帖子

你说的是哪个结果
 楼主| 发表于 2008-9-10 16:10 | 显示全部楼层
每个图像中的每一时刻z12的值均为零或相等
发表于 2008-9-10 16:43 | 显示全部楼层

回复 板凳 dearchina 的帖子

你可以设置断点,一步步运行

评分

1

查看全部评分

 楼主| 发表于 2008-9-10 19:43 | 显示全部楼层
能否指点一下怎么设置断点
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-12 03:40 , Processed in 0.057860 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表