美好 发表于 2007-10-19 16:09

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

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

         x-1      x>1
f(x)= 0       ︱x︱≤1
    x+1x<-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=;


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

无水1324 发表于 2007-10-19 16:44

回复 #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;
end
y=;

无水1324 发表于 2007-10-19 16:45

回复 #1 美好 的帖子

已经修改!

美好 发表于 2007-10-19 16:54

回复 #3 无水1324 的帖子

哪里修改了呢?:@L

无水1324 发表于 2007-10-19 16:58

回复 #4 美好 的帖子

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

VibrationMaster 发表于 2007-10-19 16:58

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

美好 发表于 2007-10-20 09:47

谢谢无水学长!
页: [1]
查看完整版本: 请帮我看看这个M文件有问题吗?