请帮我看看这个M文件有问题吗？

美好

请各位学长帮我看看这个M文件有问题吗？谢谢！
x”+2*0.1*x’+f(x)=0.1+0.2*sinwt

         x-1      x>1
f(x)= 0       ︱x︱≤1
    x+1  x<-1

function y=xiti(t,x)

if x(1)>1

ff=x(1)-1;

elseif abs(x(1))<=1

ff=0;

else

ff=x(1)+1;

y=[x(2);0.1+0.2*sin(1.4*t)-2*0.1*x(2)-ff];

[ 本帖最后由 无水1324 于 2007-10-19 16:45 编辑 ]

无水1324

function y=xiti(t,x)
if x(1)>1
  ff=x(1)-1;
elseif abs(x(1))<=1
  ff=0;
else
  ff=x(1)+1;
end
y=[x(2);0.1+0.2*sin(1.4*t)-2*0.1*x(2)-ff];

无水1324

已经修改！

美好

哪里修改了呢？:@L

无水1324

:@L
就是后面差一个end
[t,x]=ode45('xiti',[0,100],[0 0])
plot(t,x)
运行这个就可以计算出来结果了，

VibrationMaster

这个题目是求刚度非线性二阶系统的强迫振动的响应. 就你目前的水平,没有三天时间肯定弄不出来.

美好

谢谢无水学长！