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

greensmile

t=0:0.001: (13*pi/18)
/ W, W. M: j# u- p* E, u: Sx=-(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)
. [+ n+ d6 U* ~% p5 H& c: Ydx1=diff(x)
dy1=diff(y)
9 Z+ g/ w$ |* O! [$ ^6 }! Rdx11=diff(x,2)
dy11=diff(y,2)
( I! S' T6 ]. V8 h4 s
p=abs(sqrt((dx1.^2+dy1.^2).^3)/(dx1.*dy11-dy1.*dx11))
6 M  l8 Z$ y$ w) ahold on
plot(t,p)
运算时报错矩阵维度必须一致，然后我尝试改了一下。。
t=0:0.001: (13*pi/18)
x=-(50+32.5*(1-cos(t*18*pi/(13*pi)))).*sin(t)-15*cos(t)
, j3 {9 R0 \; Z1 ~+ Jy=(50+32.5*(1-cos(t*18*pi/(13*pi)))).*cos(t)-15*sin(t)
& r/ ?- F+ L2 Y( D: P3 Zdx1=diff(x)
dy1=diff(y)
( S8 N! c) G! p6 G$ r$ ]dx11=diff(x,2)
dy11=diff(y,2)
- p- b9 I4 G+ O  `# d
- R! a; z' y2 k
dx1=dx1(1:length(dx11))
dy1=dx1(1:length(dy11))


p=abs(sqrt((dx1.^2+dy1.^2).^3)/(dx1.*dy11-dy1.*dx11))
8 o1 [# Y* R* Ehold on
! G! a) D2 T; u2 o7 a+ oplot(t,p)
  d- V$ @4 [5 r8 h& `  L6 X这个时候代码无报错，但绘图只有坐标轴没有图像
( n& G! Y- ~& [# E# s
$ M1 N/ _  M7 A! \# F6 @! ]9 F: e6 I希望大神们可以给小弟指一指错误，谢谢大家

regngfpcb

不太会啊

mutougeda

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

thinkfunny

不是点运算的问题，是离散向量求差分的时候，长度会比原向量少阶数个元素（一阶差分少一个，二阶少两个），你先设置的t向量不可能和后面用差分求出来的p有相同维度的。既然原始的x和y都是t的简单函数，那就直接用解析式求精确微分（而不是用离散向量求差分）。
改成
syms t
" I$ N6 Z  r; c  J- d" ~# M5 m$ ?, wx = -(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 );
dx11 = diff( x, 2 );
dy11 = diff( y, 2 );
: R. f& Z- o0 h6 b- _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 );
plot( t, p )

greensmile

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

按照你的回答我运行了一下，结果得出的图像感觉有点奇怪。。不太像是曲率半径的图像。
然后我又按照解析式求导，把得出来的公式输入进去，代码可以正常运算了，但是还是只有坐标轴，绘制不出来图像是怎么回事呀。。
这是改正后的代码
t=0:0.001: (13*pi/18)
x=-(50+32.5*(1-cos(t*18/13))).*sin(t)-15*cos(t)
. x, `4 g1 a; K$ Ky=(50+32.5*(1-cos(t*18/13))).*cos(t)-15*sin(t)
& ~/ m3 y  T, [+ \& ]: d  c  ^dx1=-(50+32.5*(1-cos(t*18/13))).*cos(t)-(45*sin(t*18/13)-15).*sin(t)
4 G  u* o# [3 }- [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)
; B( Z0 t, k% S1 {. K# Cdy11=(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)
, F2 H' h' P) w9 O7 W' x1 lp=abs(((dx1.^2+dy1.^2).^(3/2))/(dx1.*dy11-dy1.*dx11))
8 ?' D  Y9 F9 D. uhold on
plot(t,p)
$ u! g, W% _( ]  O$ B

thinkfunny

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

3 Y4 ]4 X7 m, u: y/ @7 o) h除法没用点运算，导致p的维度不对，必须是与t维度相同的向量才能拿来跟t一起绘图。
9 C8 t' R% l* ^. z% e