%Ships pursuit
clc
clear
x0=zeros(1,1);y0=zeros(1,1);%违法船只初始位置
x1=zeros(1,1);y1=zeros(1,1);%值法船只初始位置
dt=0.01;%离散化时间,时间变成一个很小的间隔
v0=30;v1=40;%违法船只和执法船只的初始速度
x1(1,1)=20;%初始化执法船只的初始位置
i=1;whileabs(x0(1,i)-x1(1,i))>1e-4%当两者的距离小于极小值,代表追上了
d=sqrt((x1(1,i)-x0(1,i))^2+(y1(1,i)-y0(1,i))^2);%两船之间的距离
cosx=(x1(1,i)-x0(1,i))/d;%两船连线与x轴的余弦
sinx=(y0(1,i)-y1(1,i))/d;%两船连线与x轴的正弦
x0(1,1+i)=0;%更新一次数据,违法船只x不变
y0(1,1+i)=y0(1,i)+v0*dt;%违法船只y经过t后,y1=y0+vt
x1(1,1+i)=x1(1,i)-v1*dt*cosx;%执法船只朝着违法船只运动就是把xy分解,然后按照cos sin运动
y1(1,1+i)=y1(1,i)+v1*dt*sinx;
plot(x0(1,i),y0(1,i),'go',x1(1,i),y1(1,i),'ro')%go ro代表了绿色和红色
hold on %保持图像,让图像不被刷新
xlim([050])
ylim([-5579])%x y轴上限限制,也就是让最后显示图像出现的范围
pause(0.1)%延时函数 让图像延时0.1s再进行
i=i+1;%点数加1
end
d1 = sqrt((x1(1,i)-x1(1,1))^2+(y1(1,i)-y1(1,1))^2)
t = dt*i %单位小时
实现效果
最小二乘法直线拟合题目
代码及其详解
clc
clear
x =[1.3,2.0,3.6,4.2,5.3,6.6,7.8,8.3,9.1,10.2];
y =[1.2,3.2,5.9,6.5,10.3,11.2,13.0,14.6,16.9,18.3];plot(x,y,'s')
a1 =polyfit(x,y,1);% n=1 x一次拟合
%polyfit用来拟合散点图
x1 =[1.0:0.1:11.0];%x的取值范围
y1 =a1(1)*x1 +a1(2);%y的数值
%polyfit的1次拟合的第一个数据是斜率 第二个数据是常数
plot(x,y,'*')
hold on
plot(x1,y1,'ro')
k =a1(1)
b =a1(2)