I want to plot graph for my intersatellite range rate. How can i do that

Question

Toshi 2024-1-12

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2069206-i-want-to-plot-graph-for-my-intersatellite-range-rate-how-can-i-do-that

评论： Sam Chak 2024-2-25

采纳的回答： Hassaan

here is the given code :

recPos = [-7500 12500 15000];

recVel = [112.5 100.66 -50];

transPos = [-8337.5 7330.13 -12000];

transVel = [112.5 100.66 -40];

startTime = datetime(2021,4,25);

stopTime = datetime(2021,4,26);

sampleTime = 60;

sc = satelliteScenario(startTime,stopTime,sampleTime);

gs = groundStation(sc);

[p,pdot] = pseudoranges(recPos,transPos,recVel,transVel);

disp('Range between Satellites:');

disp(pseudoranges(recPos,transPos));

plot(pdot);

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Sam Chak 2024-2-25

Hey @Toshi, I wanted to inquire if the issue regarding plotting the intersatellite range rate has been resolved. It's an intriguing problem, and I'm curious to know if there have been any recent developments or breakthroughs.

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Hassaan 2024-1-12

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2069206-i-want-to-plot-graph-for-my-intersatellite-range-rate-how-can-i-do-that#answer_1388496

在 MATLAB Online 中打开

% Define the initial positions and velocities

recPos = [-7500, 12500, 15000];

recVel = [112.5, 100.66, -50] / 3600; % Convert to km/s

transPos = [-8337.5, 7330.13, -12000];

transVel = [112.5, 100.66, -40] / 3600; % Convert to km/s

% Define the time parameters

startTime = datetime(2021, 4, 25);

stopTime = datetime(2021, 4, 26);

sampleTime = 60; % Sampling time in seconds

% Calculate the number of samples

numSamples = hours(stopTime - startTime) * (3600 / sampleTime);

% Initialize arrays to store the range rate and time vector

rangeRates = zeros(numSamples, 1);

timeVec = startTime:seconds(sampleTime):stopTime - seconds(sampleTime);

% Calculate the range rate at each sample

for i = 1:numSamples

% Update the positions based on the velocities

recPos = recPos + recVel * sampleTime;

transPos = transPos + transVel * sampleTime;

% Calculate the range as the Euclidean distance between recPos and transPos

range = norm(recPos - transPos);

% Calculate the range rate as the dot product of the relative velocity

% and the unit vector in the direction from recPos to transPos

relativeVelocity = recVel - transVel;

rangeRates(i) = dot(relativeVelocity, (transPos - recPos) / range);

end

% Plot the range rate

figure;

plot(timeVec, rangeRates);

title('Intersatellite Range Rate Over Time');

xlabel('Time');

ylabel('Range Rate (km/s)');

grid on;

This script calculates the range rate at each time step and then plots it. It assumes a linear motion model for the satellites, which might not be accurate in a real-world scenario, but it's sufficient for demonstrating the concept with dummy data.

Make sure to replace the pseudoranges function call with actual code to compute the range and range rate, as the provided snippet is just a placeholder based on a simplified model. If you have a more accurate way to calculate the range and range rate, use that instead.

-----------------------------------------------------------------------------------------------------------------------------------------------------------------

If you find the solution helpful and it resolves your issue, it would be greatly appreciated if you could accept the answer. Also, leaving an upvote and a comment are also wonderful ways to provide feedback.

Professional Interests