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

pwangeng

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

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

咕噜噜

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在那里，没定义吧

咕噜噜

t在F前面定义就可以

eight

原帖由 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;
...

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

pwangeng

不好意思，其实这个程序并不存在错误问题，只是昨天我拿函数M文件直接运行了，毕业设计做头晕了，呵呵。