pwangeng 发表于 2007-5-24 17:56

程序错误：不知道如何定义

function d=fun(t,y,w)

N=14;
y=zeros(N,1);
w=2100;
m1=4;
m2=32.1;
m3=50.0;
g=9.81;
e=0.00005;
c=0.00011;
k=2.5e7;
delta2=0.6e-3;
c1=1050;
c2=2100;
k1=7.5e7;
k2=2.5e9;
cb1=350;
cb2=500;
ox1=y(1);
oy1=y(2);
odx1=y(8);
ody1=y(9);
ox2=y(3);
oy2=y(4);
odx2=y(10);
ody2=y(11);
ox3=y(5);
oy3=y(6);
odx3=y(12);
ody3=y(13);
oy4=y(7);
ody4=y(14);
if oy4<0
    cb=cb2;
    kb=k2;
elseif(oy4>=0)&(oy4<=delta2)
      cb=0;
      kb=0;
    else
         cb=cb1;
         kb=k1;
end
M=[m1   0   0    0    0   0    0;
   0    m1    0    0    0   0    0;
   0    0   m2   0    0   0    0;
   0    0   0    m2   0   0    0;
   0    0   0    0    m1    0    0;
   0    0   0    0    0   m1   0;
   0    0   0    0    0   0    m3;];
C=[c1   0   0    0    0   0    0;
   0    c1    0    0    0   0    0;
   0    0   c2   0    0   0    0;
   0    0   0    c2   0   0    0;
   0    0   0    0    c1    0    0;
   0    0   0    0    0   c1   0;
   0    0   0    0    0   0    cb;];

K=[k   0   -k    0    0   0    0;
   0    k    0    -k    0   0    0;
   -k   0    2*k    0    -k    0    0;
   0    -k   0    2*k    0    -k    0;
   0    0    -k    0   k    0    0;
   0    0    0    -k   0    k    0;
   0    0    0   0   0    0    kb;];

fx=oilx( ox1, oy1, odx1, ody1, w);
fy=oily( ox1, oy1, odx1, ody1, w);
fx1=oilx( ox3,oy3-oy4,odx3,ody3-ody4,w);
fy1=oily( ox3,oy3-oy4,odx3,ody3-ody4,w);

F=[   fx;
      fy-m1*g;
      m2*e*w^2*cos(w*t);
      m2*e*w^2*sin(w*t)-m2*g;
      fx1;
      fy1-m1*g;
      -fy1-m3*g ];
   
C=C/w;
K=K/w^2;
F=F/c/w^2;

for i=1:1:N/2
    y1(i,1)=y(i);
    y2(i,1)=y(i+N/2);
end

yy2=inv(M)*(F-C*y2-K*y1);
d=zeros(N,1);
for i=1:1:N/2
    d(i)=y2(i,1);
    d(i+N/2)=yy2(i,1);
end



错误提示：

??? Input argument "t" is undefined.
Error in ==> fun at 85
F=[   fx;


我弄不清楚t应该在哪里该怎么样定义。

zhlong 发表于 2007-5-24 18:09

楼主还是那个松动故障的程序吗,t不是函数的自变量吗?

咕噜噜 发表于 2007-5-24 19:43

回复 #1 pwangeng 的帖子

fx=oilx( ox1, oy1, odx1, ody1, w);
fy=oily( ox1, oy1, odx1, ody1, w);
fx1=oilx( ox3,oy3-oy4,odx3,ody3-ody4,w);
fy1=oily( ox3,oy3-oy4,odx3,ody3-ody4,w);
这个里面的oilx在那里，没定义吧

咕噜噜 发表于 2007-5-24 19:45

回复 #1 pwangeng 的帖子

t在F前面定义就可以

eight 发表于 2007-5-24 20:01

原帖由 pwangeng 于 2007-5-24 17:56 发表 http://www.chinavib.com/forum/images/common/back.gif
function d=fun(t,y,w)

N=14;
y=zeros(N,1);
w=2100;
m1=4;
m2=32.1;
m3=50.0;
g=9.81;
e=0.00005;
c=0.00011;
k=2.5e7;
delta2=0.6e-3;
c1=1050;
c2=2100;
k1=7.5e7;
k2=2.5e9;
cb1=350;
...

请先看看置顶帖：聚宝盆， 然后看看里面提到的精华贴

pwangeng 发表于 2007-5-25 12:20

不好意思，其实这个程序并不存在错误问题，只是昨天我拿函数M文件直接运行了，毕业设计做头晕了，呵呵。
页: [1]
查看完整版本: 程序错误：不知道如何定义