追击走私船题目

在这里插入图片描述 2.将曲线画为动态曲线

代码及其详解

%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+vtx1(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 %单位小时

动态曲线代码

%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;
writerObj = VideoWriter('Pursuit.avi')      %写入文件
open(writerObj)         %打开文件
while abs(x0(1,i)-x1(1,i)) > 1e-4   %当两者的距离小于极小值，代表追上了xlabel('东')ylabel('北')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+vtx1(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代表了绿色和红色axis([-10,50,-55,79])            %设置座标的xmin xmax ymin ymaxframe = getframe                 %获得一帧图像writeVideo(writerObj,frame)      %把这一帧显示到图像里面hold off    %保持图像，让图像不被刷新%如果把这里改为hold on 那么轨迹路线也会显示出来pause(0.1)%延时函数 让图像延时0.1s再进行i=i+1;    %点数加1
end
xlim([-10 50])
ylim([-55 79])%x y轴上限限制，也就是让最后显示图像出现的范围
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)