要对分钟数做一点小处理，比如12:40实际上是12.6667点，而9.5点实际上就是0.5*60=30，即9:30。
参考代码：
clc
clear all
close all
t = 9:0.01:15;
tr = abs(t-12) * 60; % 于12点相差的分钟数
m = 3.2097 ./ cos(tr/4);
figure
plot(t, m);
grid on
xlabel('t/h');
ylabel('m');
效果图
内容全部来源于网络收集，如有侵权，请联系网站删除：QQ：24596024