Matlab绘图，代码无报错，但绘图只有坐标轴没有图像，请大神指导

greensmile

t=0:0.001: (13*pi/18)
  r  o2 d3 o* \9 n, sx=-(50+32.5*(1-cos(t*18*pi/(13*pi)))).*sin(t)-15*cos(t)
- F+ P) x! P; w) b* ty=(50+32.5*(1-cos(t*18*pi/(13*pi)))).*cos(t)-15*sin(t)
dx1=diff(x)
2 S& Y4 {9 H0 Q; @( u4 Mdy1=diff(y)
dx11=diff(x,2)
8 Z6 H/ Z4 \1 I" [% Pdy11=diff(y,2)
% S. M9 T8 q8 T! Y; ~2 B! Z
; {* U+ R' G! \$ E2 Z2 op=abs(sqrt((dx1.^2+dy1.^2).^3)/(dx1.*dy11-dy1.*dx11))
/ Y) z  f  y. c7 |3 s4 s4 vhold on
: @  E/ v2 n% s7 z3 p1 c# @  u2 U. O9 ?plot(t,p)
( N3 n7 w+ o  \4 E+ b2 q运算时报错矩阵维度必须一致，然后我尝试改了一下。。
/ _8 L. v! X4 \" B' Kt=0:0.001: (13*pi/18)
x=-(50+32.5*(1-cos(t*18*pi/(13*pi)))).*sin(t)-15*cos(t)
y=(50+32.5*(1-cos(t*18*pi/(13*pi)))).*cos(t)-15*sin(t)
dx1=diff(x)
1 R' q  j- J; A; X6 pdy1=diff(y)
& W. V: V; ]+ _dx11=diff(x,2)
. n* p5 c7 d6 j8 ]dy11=diff(y,2)
" ~# j" C  A- e& }2 \

dx1=dx1(1:length(dx11))
/ f% i: ]( T2 G6 G7 Fdy1=dx1(1:length(dy11))


* f1 G9 o+ x- d# b6 D) Np=abs(sqrt((dx1.^2+dy1.^2).^3)/(dx1.*dy11-dy1.*dx11))
hold on
! L8 O1 e! {& \. o8 Q& s( Iplot(t,p)
3 y7 Y' O" g7 t+ J3 a3 @. Y这个时候代码无报错，但绘图只有坐标轴没有图像
; Y8 B- m) W+ _- P( w. G
希望大神们可以给小弟指一指错误，谢谢大家

thinkfunny

greensmile 发表于 2020-7-27 18:34
% {# h1 d& O. m& B" s按照你的回答我运行了一下，结果得出的图像感觉有点奇怪。。不太像是曲率半径的图像。
  T2 J( t% f4 Z5 F' y8 U+ j然后我又按照解析 ...

除法没用点运算，导致p的维度不对，必须是与t维度相同的向量才能拿来跟t一起绘图。
4 [6 T3 o7 T+ R0 u. R0 e1 b" w+ \

greensmile

thinkfunny 发表于 2020-7-27 18:32
不是点运算的问题，是离散向量求差分的时候，长度会比原向量少阶数个元素（一阶差分少一个，二阶少两个）， ...

* Q0 h+ I0 h4 M# g7 R) U按照你的回答我运行了一下，结果得出的图像感觉有点奇怪。。不太像是曲率半径的图像。
9 I. C2 ~: k2 a然后我又按照解析式求导，把得出来的公式输入进去，代码可以正常运算了，但是还是只有坐标轴，绘制不出来图像是怎么回事呀。。
4 a, y8 E- K( |1 A这是改正后的代码
t=0:0.001: (13*pi/18)
0 K/ ?* a- w1 Yx=-(50+32.5*(1-cos(t*18/13))).*sin(t)-15*cos(t)
y=(50+32.5*(1-cos(t*18/13))).*cos(t)-15*sin(t)
dx1=-(50+32.5*(1-cos(t*18/13))).*cos(t)-(45*sin(t*18/13)-15).*sin(t)
dy1=-(50+32.5*(1-cos(t*18/13))).*sin(t)+(45*sin(t*18/13)-15).*cos(t)
dx11=-(2*45*sin(t*18/13)-15).*cos(t)-(810*cos(t*18/13)/13-50-32.5*(1-cos(t*18/13))).*sin(t)
/ V7 F3 @1 }  w) I7 m3 z. ody11=(810*cos(t*18/13)/13-50-32.5*(1-cos(t*18/13))).*cos(t)-(2*45*sin(t*18/13)-15).*sin(t)
2 t+ G6 a2 i/ F6 j9 _1 kp=abs(((dx1.^2+dy1.^2).^(3/2))/(dx1.*dy11-dy1.*dx11))
  [# ~- F( v1 Y8 K: \, P: Nhold on
1 `' _$ z3 n  I: h: x" L# ^plot(t,p)

thinkfunny

不是点运算的问题，是离散向量求差分的时候，长度会比原向量少阶数个元素（一阶差分少一个，二阶少两个），你先设置的t向量不可能和后面用差分求出来的p有相同维度的。既然原始的x和y都是t的简单函数，那就直接用解析式求精确微分（而不是用离散向量求差分）。
& B* q: F2 X& S0 F改成
) `  t, X$ u3 F5 n' H! Dsyms t
x = -(50+32.5*(1-cos(t*18*pi/(13*pi)))).*sin(t)-15*cos(t);
' J% J( e# V; ]y = (50+32.5*(1-cos(t*18*pi/(13*pi)))).*cos(t)-15*sin(t);
0 m9 N+ t8 X0 e5 n+ F2 Tdx1 = diff( x );
dy1 = diff( y );
# O1 N, [2 Q1 \/ s* B9 f( r4 }- [$ Ndx11 = diff( x, 2 );
, |! h  [% c) ^4 V1 ]' z9 sdy11 = diff( y, 2 );
: a- L2 a+ d* K8 ~+ j0 D0 Dp = abs(sqrt((dx1.^2+dy1.^2).^3)./(dx1.*dy11-dy1.*dx11));
8 x6 `/ s/ E6 J1 PPFunction = matlabFunction( p );
t = 0 : 0.001 : (13*pi/18);
/ C& m# D7 E9 F5 u9 u" y4 S  S( fp = PFunction( t );
% j& m( v2 ]" N9 ~  Bplot( t, p )

mutougeda

按照你的代码p是定值，你少写了一个点，正确代码如下
p=abs(sqrt((dx1.^2+dy1.^2).^3)./(dx1.*dy11-dy1.*dx11))

regngfpcb

不太会啊