第一章等切面曲线和相似曲线 5 5*sin(t) Xs =-5*sin(t): Ys 5*cos(t) MATLAB提供两个M文件ade23和ode45,用于求积微分方程.在下面的主程序中,我们将调用 这些函数之一,并定义初始条件(注意,当t=0时,小孩在点(5,0)处和玩具在点(10,0)处) >>% main1 >>y0=[10o] >>[ty]=ode45(f’,[0100,yo) >> clf hold on; >>axis([-610-610]); >>axis >>plot(y(:,1),y(:,2)) 图1.3小孩在圆上行走 如果绘制两列y,则可得玩具的轨迹(参见图13),而且在同一图中,用语句加进小孩的曲线: >>t=0:0.05:6.3; >>[X, Xs, Y, Ys]= child(t) >>hold off 注意,在程序中,硬棒的长a没有明显出现;它由玩具的位置(初始条件),和t=0时小孩的位置(函 数chi1d)隐含地定义 我们用几个例子来结束本节.假设小孩沿sin函数的图形行走:X(t)=t和Y(t)=5sint,小 孩的曲线用虚线表示,初始条件c(0)=0和y(0)=10,我们得到图14 另一个例子,小孩在圆X(t)=5cost,Y(t)=5sint上行走,初始条件x(0)=0和y()=10, 我们得到一个象花的玩具轨迹(参见图1.5)