如何用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]