已知可视卫星位置 （X:--14687383.4,Y:494666.7,Z:21815130.7) (X:-21196834.3,Y:10285023.5,Z:12008060.9) (X:4786048.2,Y:25852523.5,Z:3431362.8) (X:-10594640.3,Y:21523420.2,Z:11244612.1) (X:-1499577.0,Y:15272138.8,Z:22099342.4) (X:8198094.4,Y:16682058.1,Z:18863012.1) (X:-17457205.9,Y:-8604158.0,Z:18989327.4) (X:17055061.1,Y:3556726.5,Z:19795781.2) (X:-1487664.0,Y:24096277.1,Z:-10650122.6) 伪距（m） (21933314.42,21243122.33,22462827.50,20239446.11,20998877.44,21813415.63,24858805.55,24851219.61,24432291.70) ，设接收机钟差为0，不考虑接收机钟差时的接收机三维坐标MATLAB程序

好的，根据提供的可视卫星位置和伪距数据，可以使用 MATLAB 求解接收机的三维坐标。以下是一个简单的程序：

% 可视卫星位置
satellite_pos = [-14687383.4, 494666.7, 21815130.7;
                 -21196834.3, 10285023.5, 12008060.9;
                  4786048.2, 25852523.5, 3431362.8;
                 -10594640.3, 21523420.2, 11244612.1;
                  -1499577.0, 15272138.8, 22099342.4;
                   8198094.4, 16682058.1, 18863012.1;
                 -17457205.9, -8604158.0, 18989327.4;
                  17055061.1, 3556726.5, 19795781.2;
                  -1487664.0, 24096277.1, -10650122.6];

% 伪距
pseudoranges = [21933314.42, 21243122.33, 22462827.50, ...
                20239446.11, 20998877.44, 21813415.63, ...
                24858805.55, 24851219.61, 24432291.70];

% 求解接收机位置
A = [-2 * satellite_pos(1, 1) + 2 * satellite_pos(2:9, 1), ...
     -2 * satellite_pos(1, 2) + 2 * satellite_pos(2:9, 2), ...
     -2 * satellite_pos(1, 3) + 2 * satellite_pos(2:9, 3)];
b = pseudoranges(2:9).^2 - pseudoranges(1)^2 - ...
    satellite_pos(2:9, 1).^2 + satellite_pos(1, 1).^2 - ...
    satellite_pos(2:9, 2).^2 + satellite_pos(1, 2).^2 - ...
    satellite_pos(2:9, 3).^2 + satellite_pos(1, 3).^2;
x = A \ b;
receiver_pos = [x(1), x(2), x(3)];

% 输出结果
fprintf('接收机位置：\n');
fprintf('X：%f\n', receiver_pos(1));
fprintf('Y：%f\n', receiver_pos(2));
fprintf('Z：%f\n', receiver_pos(3));

需要注意的是，该程序仅仅是一个简单的示例，实际应用中可能需要考虑更多的因素，如测量误差的处理、卫星钟差的影响等。