642712189 发表于 2014-4-28 13:54

如何用MATLAB模拟傅科摆运动轨迹

其中地球重力加速度g=9.8。绳长67。w为地球自转角速度
程序如下

a=input(‘输入纬度=‘);
q=input(‘按次格式依次输入=')
c=a*pi/;
=ode45(‘fkbfun’,,q,[].c);
xlabel(‘x’)
ylabel(‘y’);
comet(x(:,1),x(:,3))
函数文件是一个独立的文件,文件名为 fkbfun.m.
functiontt=fkb(t,x,flag,c)]
a=(2*pi*sin(c))/100;(放大了地球的角速度)
b=9.8/67;
tt=[x(2);
2*a*x(4)-b*x(1);
x(4);
-2*a*x(2)-b*x(3)];

为什么运行错误呢   求高手讲解!
页: [1]
查看完整版本: 如何用MATLAB模拟傅科摆运动轨迹