要对分钟数做一点小处理,比如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');效果图
你根据什么说它错了.
你的参数这么多,(除了自变量以外,还有s,r,theta0)
你要人家怎么给你简化.
你试一下
simple(l),这里给出它的化简式,没有一个是短的.
你可以先用一些值代替s,r,theta0
看一下.
建议你自己把dy和dx求出来,然后用数值积分计算,命令是quad