How to plot angle (in radians) versus phase shift?

Question

Rahul 2023-6-20

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1985524-how-to-plot-angle-in-radians-versus-phase-shift

评论： David Goodmanson 2023-6-20

采纳的回答： Naman

在 MATLAB Online 中打开

Hi,

I want a profile with phase angle plotted against angle in radians.

Presently I have a code with phase angle plotted against a linear distance which is as below.

Need your valuable advice how to do the necessary changes in this code.

Also I need to determine its slope.

Sincerely,

rc

clc;
clear;
phi_1 = [0.1724*pi  0.1472*pi  2.18*pi  2.25*pi  4*pi  3.588*pi  1.012*pi  4.112*pi  4.31*pi  4.25*pi  0.832*pi];
x1 = linspace(0,1,11);
coefficients1=polyfit(x1,phi_1,1);
slope1=coefficients1(1);
intercept1=coefficients1(2);
figure(1)
hold on
plot(x1,phi_1,'o')
plot([0 1],intercept1 + slope1*[0  1])

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

David Goodmanson 2023-6-20

在 MATLAB Online 中打开

Hi Rahul,

you can shorten up the code by eliminating the slope and intercept calculations and plotting the line with

plot([0 1],polyval(coefficients1,[0 1])).

It's not possible to say anything about plotting using angle without knowing the lower and upper limits of the angle. Anyway, any kind of straight line fit to this data looks pretty dubious.

请先登录，再进行评论。

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

Answer 1

Naman 2023-6-20

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1985524-how-to-plot-angle-in-radians-versus-phase-shift#answer_1259069

在 MATLAB Online 中打开

Hi Rahul,

To plot phase angle against angle in radians instead of linear distance, you will need to replace the x-axis values with angle values. Assuming you want to plot angle values from 0 to 2*pi, you can generate them using the linspace function and modify the plot accordingly. Here's the modified code:

clc;

clear;

phi_1 = [0.1724*pi 0.1472*pi 2.18*pi 2.25*pi 4*pi 3.588*pi 1.012*pi 4.112*pi 4.31*pi 4.25*pi 0.832*pi];

theta = linspace(0, 2*pi, numel(phi_1)); % generate angle values

coefficients1 = polyfit(theta, phi_1, 1);

slope1 = coefficients1(1);

intercept1 = coefficients1(2);

figure(1)

hold on

plot(theta, phi_1, 'o') % plot phase angle against angle in radians

plot([0 2*pi], intercept1 + slope1 * [0 2*pi]) % plot linear fit

xlabel('Angle (radians)') % update x-axis label