how can i plot this code ?

Question

Tomer Segev 2020-9-30

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/602530-how-can-i-plot-this-code

评论： Rik 2020-10-1

hey, I want to plot this code the ri veriable is a scalar and not a vector, also how can I define Z1 and Z2 without the operator @ ?

s= [0,5*10^-2,1*10^-1,1.5*10^-1,2*10^-1,2.5*10^-1,3*10^-1,3.5*10^-1,4*10^-1,4.5*10^-1,5*10^-1,5.5*10^-1,6*10^-1,6.5*10^-1,7*10^-1,7.5*10^-1,8*10^-1,8.5*10^-1,9*10^-1,9.5*10^-1,9.9*10^-1,9.99*10^-1,1,1,1,1,1];
omega= [3.74*10^-1, 3.8*10^-1,3.87*10^-1,3.94*10^-1,4.02*10^-1,4.11*10^-1,4.2*10^-1,4.29*10^-1,4.4*10^-1,4.51*10^-1,4.64*10^-1,4.78*10^-1,4.94*10^-1,5.12*10^-1,5.33*10^-1,5.57*10^-1,5.86*10^-1,6.23*10^-1,6.72*10^-1,7.46*10^-1,8.71*10^-1,9.56*10^-1,9.68*10^-1,9.72*10^-1,9.75*10^-1,9.8*10^-1,9.86*10^-1];
plot(s,omega,'.');
hold on
Z1=@(s) (1+ ((1-s.^2)^(1/3))*((1+s).^1/3+(1-s).^1/3));
Z2=@(s,Z1) ((3*s.^2+Z1.^2).^1/2);
p=2;
w=3;
c= 3*10^8;
G=6.67384*10^-11;
M=1;
ri= (G*M./c^2).*(3.+Z2+sqrt(9+6.*Z2-(Z1).^2-2*(Z2.*Z1)));
plot(ri,omega,'r');

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Alan Stevens 2020-9-30

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/602530-how-can-i-plot-this-code#answer_502966

在 MATLAB Online 中打开

You need to plot them separately as they have completely different 'x' values (indeed the ri values seem unbelievably small!):

s= [0,5*10^-2,1*10^-1,1.5*10^-1,2*10^-1,2.5*10^-1,3*10^-1,3.5*10^-1,4*10^-1,4.5*10^-1,5*10^-1,5.5*10^-1,6*10^-1,6.5*10^-1,7*10^-1,7.5*10^-1,8*10^-1,8.5*10^-1,9*10^-1,9.5*10^-1,9.9*10^-1,9.99*10^-1,1,1,1,1,1];
omega= [3.74*10^-1, 3.8*10^-1,3.87*10^-1,3.94*10^-1,4.02*10^-1,4.11*10^-1,4.2*10^-1,4.29*10^-1,4.4*10^-1,4.51*10^-1,4.64*10^-1,4.78*10^-1,4.94*10^-1,5.12*10^-1,5.33*10^-1,5.57*10^-1,5.86*10^-1,6.23*10^-1,6.72*10^-1,7.46*10^-1,8.71*10^-1,9.56*10^-1,9.68*10^-1,9.72*10^-1,9.75*10^-1,9.8*10^-1,9.86*10^-1];
Z1= (1+ ((1-s.^2).^(1/3)).*((1+s).^1/3+(1-s).^1/3));
Z2= ((3*s.^2+Z1.^2).^1/2);
p=2;
w=3;
c= 3*10^8;
G=6.67384*10^-11;
M=1;
ri= (G*M./c^2).*(3.+Z2+sqrt(9+6.*Z2-(Z1).^2-2*(Z2.*Z1)));
subplot(2,1,1);
plot(s,omega,'.');
xlabel('s'),ylabel('omega')
subplot(2,1,2)
plot(ri,omega,'r');
xlabel('ri'),ylabel('omega')