dt=1/20; % 设定dt为步长时间
A=zeros(300,3); % 先分配300个点
for k=1:3; % 计算10个周期
for m=1:20: (150+k); % 计算每一个周期内每点的增量情况
[dx,dy]=Doubleheightdeposit(A(m,2),A(m,3)-A(150+k-1,3),dt); % 调用函数
A(m,2)=A(m,2)+dx;
A(m,3)=A(m,3)+dy;
end;
A(150+k+1,3)=A(150+k,3); %增加一个新点，其高度为前一点的高度(为什么这一句会出错？)
m=1:20: (150+k+1);
plot(A(m,2),A(m,3),'-k*');
end
***====*****=====******
function [dx,dy]=Doublem(A(m,2),A(m,3),dt) % 要调用的函数
dx=1;
dy=2;
|