close?all
clear?all
clc
%%?輸入參數
syms?y（t）?x（t）??????????????????????????%?創建符號變量/函數
a?=?input（‘請輸入緯度=‘）;
x0?=?input（‘請輸入x方向初始坐標x0=‘）;
y0?=?input（‘請輸入y方向初始坐標y0=‘）;
u0?=?input（‘請輸入x方向初始速度u0=‘）;
v0?=?input（‘請輸入y方向初始速度v0=‘）;
T?=?input（‘請輸入總時間T=‘）;

%%?常數設置
tPause?=?0;?????????%?繪圖間隔時間
c?=?a*pi/180;???????%?弧度制轉換
w?=?pi/50;??????????%?地球自轉角速度，實際約為2pi/86400，為效果明顯，暫設為2pi/100
l?=?67;?????????????%?傅科擺擺長
g?=?9.8;????????????%?重力加速度

%%?求解微分方程組
eqns?=?[diff（xt2）-2*w*diff（yt）*sin（c）+x*g/l==0?diff（yt2）+2*w*diff（xt）*sin（c）+y*g/l==0];
Dx?=?diff（xt）;
Dy?=?diff（yt）;
cond?=?[x（0）?==?x0y（0）?==?y0?Dx（0）?==?u0?Dy（0）?==?v0];
[x?y]?=?dsolve（eqnscond）;

t?=?0:0.1:T;????????%?時間

x?=?eval（x）;????????%?運行符號解，得到解向量
y?=?eval（y）;????????%?同上

%%?可視化
set（gcf‘color‘[1?1?1]）;
axis（[-0.6?0.6?-1?0.2]）;
axis?off
axis?equal
z?=?-（l^2-x.^2-y.^2）.^0.5;

%?豎軸
verticalaxis?=?line（[00][00][0min（z）]‘color‘‘k‘‘linestyle‘‘-‘‘linewidth‘1.5）;
%?擺線初始化
pendulumline?=?line（[0x（1）][0y（1）][0z（1）]‘color‘‘y‘‘linestyle‘‘-‘‘linewidth‘2）;??
%?擺錘初始化
pendulumbob?=?line（x（1）y（1）z（1）‘color‘‘r‘‘marker‘‘.‘‘markersize‘?25）;???????????????
%?擺錘起點
pendulumtra?=?line（[x（1）x（1）][y（1）y（1）][z（1）z（1）]）;????????????????????????????????????

axis?square
axis（[min（x）?max（x）?min（y）?max（y）?min（z）?0]）
hold?on
grid?on

for?i?=?1:size（t2）
????%?更新擺線位置
????set（pendulumline‘xdata‘[0x（i）]‘ydata‘[0y（i）]‘zdata‘[0z（i）]）;?
????%?更新擺錘位置
????set（pendulumbob‘xdata‘x（i）‘ydata‘y（i）‘zdata‘z（i））;?
????%?更新軌跡
????set（pendulumtra‘xdata‘x（1:i）‘ydata‘y（1:i）‘zdata‘z（1:i））?
????%?繪制
????drawnow;
????%?暫停
????pause（tPause）;
end?

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-02-24?15:00??2_3_foucaultpendulum\
?????文件????????1965??2018-02-24?14:56??2_3_foucaultpendulum\dsolve_pendulum.m
?????目錄???????????0??2018-02-24?21:06??__MACOSX\
?????目錄???????????0??2018-02-24?21:06??__MACOSX\2_3_foucaultpendulum\
?????文件?????????210??2018-02-24?14:56??__MACOSX\2_3_foucaultpendulum\._dsolve_pendulum.m
?????文件?????????162??2018-02-24?14:29??2_3_foucaultpendulum\myODE.m
?????文件?????????210??2018-02-24?14:29??__MACOSX\2_3_foucaultpendulum\._myODE.m
?????文件????????1857??2018-02-24?14:59??2_3_foucaultpendulum\ode45_pendulum.m
?????文件?????????120??2018-02-24?14:59??__MACOSX\2_3_foucaultpendulum\._ode45_pendulum.m