matlab输入两个点的坐标共四个参数 （x1,y1,x2,y2 ）计算 pI-P2 的方位角，要求结果显示 degree min second

可以使用以下代码实现：

% 输入两个点的坐标
x1 = input('请输入第一个点的 x 坐标：');
y1 = input('请输入第一个点的 y 坐标：');
x2 = input('请输入第二个点的 x 坐标：');
y2 = input('请输入第二个点的 y 坐标：');

% 计算方位角
dx = x2 - x1;
dy = y2 - y1;
angle = atan2(dy, dx);

% 将弧度转换为度分秒
deg = rad2deg(angle);
deg_int = floor(deg);
min_float = (deg - deg_int) * 60;
min_int = floor(min_float);
sec_float = (min_float - min_int) * 60;

% 输出结果
fprintf('点(%g, %g)到点(%g, %g)的方位角为：%d° %d'' %.2f" \n', x1, y1, x2, y2, deg_int, min_int, sec_float);

使用示例：

请输入第一个点的 x 坐标：0
请输入第一个点的 y 坐标：0
请输入第二个点的 x 坐标：1
请输入第二个点的 y 坐标：1
点(0, 0)到点(1, 1)的方位角为：45° 0' 0.00"

相关问题

用matlab写一段代码，输入两个坐标系的四个点坐标，输出bursha七参数

由于Bursha七参数的计算涉及到矩阵运算，因此需要使用matlab中的矩阵操作函数。下面是一段简单的matlab代码，可以计算出两个坐标系之间的Bursha七参数。

% 输入两个坐标系的四个点坐标
A1 = [x1_1, y1_1, z1_1; x1_2, y1_2, z1_2; x1_3, y1_3, z1_3; x1_4, y1_4, z1_4];
A2 = [x2_1, y2_1, z2_1; x2_2, y2_2, z2_2; x2_3, y2_3, z2_3; x2_4, y2_4, z2_4];

% 构造增广矩阵
n = size(A1, 1);
M = [A1, ones(n, 1), zeros(n, 3); zeros(n, 3), A1, ones(n, 1)];
N = [A2; A2];

% 求解Bursha七参数
p = pinv(M) * N;
B = p(1:3, :);
t = p(4:6, :);
s = p(7, :);

% 输出结果
disp('Bursha七参数为：');
disp('旋转矩阵：');
disp(B);
disp('平移向量：');
disp(t);
disp('尺度因子：');
disp(s);

在上述代码中，输入的是两个坐标系的四个点坐标，其中A1表示第一个坐标系的四个点坐标，A2表示第二个坐标系的四个点坐标。这里使用的是增广矩阵法求解Bursha七参数，将两个坐标系的点坐标组成增广矩阵，然后通过矩阵运算求解出旋转矩阵B、平移向量t和尺度因子s，最后将结果输出。

用matlab编写函数式文件计算下式：y1=2x1+x2,y2=5x1-x2

可以使用以下的MATLAB代码实现：

% 读取输入文件
data = load('input.txt');
% 计算 y1 和 y2
y1 = 2 * data(:,1) + data(:,2);
y2 = 5 * data(:,1) - data(:,2);
% 输出结果
fprintf('%d %d\n', [y1, y2]');

假设输入文件为`input.txt`，格式为：

x1 x2
1 2
3 4
5 6

运行代码，输出结果如下：

4   8  12
-3  -7 -11

其中第一行为 y1 的计算结果，第二行为 y2 的计算结果。