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

greensmile

t=0:0.001: (13*pi/18)
& r  h7 m3 `1 E  sx=-(50+32.5*(1-cos(t*18*pi/(13*pi)))).*sin(t)-15*cos(t)
) v! z7 |$ v, U2 ^2 Xy=(50+32.5*(1-cos(t*18*pi/(13*pi)))).*cos(t)-15*sin(t)
  m! X& w7 b, m' w$ jdx1=diff(x)
dy1=diff(y)
! o& o1 S1 d' o+ @8 ?2 bdx11=diff(x,2)
dy11=diff(y,2)

p=abs(sqrt((dx1.^2+dy1.^2).^3)/(dx1.*dy11-dy1.*dx11))
hold on
plot(t,p)
! O8 z, ^! o+ ?" t! w5 ?! K运算时报错矩阵维度必须一致，然后我尝试改了一下。。
t=0:0.001: (13*pi/18)
$ m( F/ [. \; g/ A8 v0 M' ax=-(50+32.5*(1-cos(t*18*pi/(13*pi)))).*sin(t)-15*cos(t)
5 M/ l! r- j' G" h# g7 Uy=(50+32.5*(1-cos(t*18*pi/(13*pi)))).*cos(t)-15*sin(t)
: D* `( v  s* Wdx1=diff(x)
& l8 P' z2 r/ Q5 @6 l% {' xdy1=diff(y)
5 z. Q7 v# Q3 Ndx11=diff(x,2)
) S' F' c, D/ M2 Q0 e/ ]% {$ Qdy11=diff(y,2)


/ q3 N9 k5 ]2 H" ]* P& j8 jdx1=dx1(1:length(dx11))
dy1=dx1(1:length(dy11))
' I" s3 a! M- R9 B: B9 D5 W
9 `1 `# B+ `" g2 G4 k) y' w5 h- f
3 R' g$ Y' j1 H4 U  T& rp=abs(sqrt((dx1.^2+dy1.^2).^3)/(dx1.*dy11-dy1.*dx11))
hold on
4 X# T: X+ C, v- B: Zplot(t,p)
这个时候代码无报错，但绘图只有坐标轴没有图像

/ y: r! T# a4 ~% ]0 t6 u希望大神们可以给小弟指一指错误，谢谢大家

regngfpcb

不太会啊

mutougeda

& Y, l1 M5 g$ o7 {; P& z按照你的代码p是定值，你少写了一个点，正确代码如下
$ {' b+ ]1 M" b' Lp=abs(sqrt((dx1.^2+dy1.^2).^3)./(dx1.*dy11-dy1.*dx11))

thinkfunny

不是点运算的问题，是离散向量求差分的时候，长度会比原向量少阶数个元素（一阶差分少一个，二阶少两个），你先设置的t向量不可能和后面用差分求出来的p有相同维度的。既然原始的x和y都是t的简单函数，那就直接用解析式求精确微分（而不是用离散向量求差分）。
+ ?9 \, m) ]9 S9 K+ E- V$ j改成
  C$ Q# L8 k. j2 X/ u, Q3 Csyms t
" ]: j1 c/ q, ]3 j% @5 W' f/ fx = -(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 );
dy1 = diff( y );
: U5 x/ n- v$ |7 R1 Fdx11 = diff( x, 2 );
% ]' Z/ |9 M* N3 P/ n. Pdy11 = diff( y, 2 );
p = abs(sqrt((dx1.^2+dy1.^2).^3)./(dx1.*dy11-dy1.*dx11));
PFunction = matlabFunction( p );
t = 0 : 0.001 : (13*pi/18);
p = PFunction( t );
, s& C: O& {9 ^  y2 y* N4 Wplot( t, p )

greensmile

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

按照你的回答我运行了一下，结果得出的图像感觉有点奇怪。。不太像是曲率半径的图像。
然后我又按照解析式求导，把得出来的公式输入进去，代码可以正常运算了，但是还是只有坐标轴，绘制不出来图像是怎么回事呀。。
; Y( i4 b  T7 U. c2 `  ?这是改正后的代码
( n, o7 Q3 O* _& ^  St=0:0.001: (13*pi/18)
, x- `" h  o# ]1 ex=-(50+32.5*(1-cos(t*18/13))).*sin(t)-15*cos(t)
+ n+ ?9 a+ k' G$ L' |% |1 ty=(50+32.5*(1-cos(t*18/13))).*cos(t)-15*sin(t)
+ E: k# J8 g, v3 `/ O* _* t. ]dx1=-(50+32.5*(1-cos(t*18/13))).*cos(t)-(45*sin(t*18/13)-15).*sin(t)
7 V# @  p9 L; Q/ X8 Kdy1=-(50+32.5*(1-cos(t*18/13))).*sin(t)+(45*sin(t*18/13)-15).*cos(t)
4 N, j) W$ S5 f9 E7 ?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)
% R. h' g! G& i7 Rdy11=(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)
p=abs(((dx1.^2+dy1.^2).^(3/2))/(dx1.*dy11-dy1.*dx11))
hold on
. G# K! a5 V( I- u' Lplot(t,p)
  z! @% O7 F  K8 W5 V- E

thinkfunny

greensmile 发表于 2020-7-27 18:34
按照你的回答我运行了一下，结果得出的图像感觉有点奇怪。。不太像是曲率半径的图像。
然后我又按照解析 ...

$ ?0 U, z7 U1 `8 k除法没用点运算，导致p的维度不对，必须是与t维度相同的向量才能拿来跟t一起绘图。