追击走私船题目
代码及其详解
%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;
while abs(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([0 50])
ylim([-55 79])%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)
实现效果
