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

greensmile

t=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)
- e/ X: `' J# cdy1=diff(y)
& H5 B8 ^# s; E8 Bdx11=diff(x,2)
dy11=diff(y,2)

# x1 G; r9 Q8 n4 _0 Lp=abs(sqrt((dx1.^2+dy1.^2).^3)/(dx1.*dy11-dy1.*dx11))
hold on
$ x, c- t, i- {# S6 a6 Dplot(t,p)
' y* c0 v* r) c, U6 c  ]6 o运算时报错矩阵维度必须一致，然后我尝试改了一下。。
+ B4 Q% k6 Y" t0 Vt=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)
! g- y) z5 H. u4 Udx1=diff(x)
7 X' ^- p& \7 a" i4 ~: f: Mdy1=diff(y)
dx11=diff(x,2)
dy11=diff(y,2)
/ D5 m$ ^8 X6 T1 w( h( A

9 J' N5 V. Q$ P# f: Y  Ddx1=dx1(1:length(dx11))
dy1=dx1(1:length(dy11))
& R  O% r9 o4 W" W
& M* m: x. |# d) e4 l# T0 Q
0 z7 P/ @+ U( L' X; ip=abs(sqrt((dx1.^2+dy1.^2).^3)/(dx1.*dy11-dy1.*dx11))
/ `2 N! y9 H6 w4 w% w, T' Q6 [* Zhold on
plot(t,p)
" U1 D  ~( |2 I这个时候代码无报错，但绘图只有坐标轴没有图像

希望大神们可以给小弟指一指错误，谢谢大家

" P% l8 E* s0 N1 t% n

regngfpcb

不太会啊

mutougeda

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

thinkfunny

不是点运算的问题，是离散向量求差分的时候，长度会比原向量少阶数个元素（一阶差分少一个，二阶少两个），你先设置的t向量不可能和后面用差分求出来的p有相同维度的。既然原始的x和y都是t的简单函数，那就直接用解析式求精确微分（而不是用离散向量求差分）。
改成
syms t
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 );
; k' ]) F' `8 r$ [# Hdy1 = diff( y );
2 ?' b9 f# W; [: F0 v2 v& pdx11 = diff( x, 2 );
dy11 = diff( y, 2 );
. z0 u* I7 |' k) wp = abs(sqrt((dx1.^2+dy1.^2).^3)./(dx1.*dy11-dy1.*dx11));
6 T. ?% E1 d3 x. D9 K- SPFunction = matlabFunction( p );
8 W' \3 I0 l* D* ~5 [! t  t9 k) `t = 0 : 0.001 : (13*pi/18);
8 w% k# T. d" L" R# j8 {9 jp = PFunction( t );
plot( t, p )

greensmile

thinkfunny 发表于 2020-7-27 18:32
3 [2 M+ Z. x6 M8 g# f* @不是点运算的问题，是离散向量求差分的时候，长度会比原向量少阶数个元素（一阶差分少一个，二阶少两个）， ...

2 p8 w6 a, r$ k! B  K8 k5 S按照你的回答我运行了一下，结果得出的图像感觉有点奇怪。。不太像是曲率半径的图像。
% `5 K! ?/ r" N# |" p  x$ M( y然后我又按照解析式求导，把得出来的公式输入进去，代码可以正常运算了，但是还是只有坐标轴，绘制不出来图像是怎么回事呀。。
这是改正后的代码
t=0:0.001: (13*pi/18)
1 r0 p3 E; s6 W) y  J  ~) Q# nx=-(50+32.5*(1-cos(t*18/13))).*sin(t)-15*cos(t)
! x3 p7 Y( m( x: Y" dy=(50+32.5*(1-cos(t*18/13))).*cos(t)-15*sin(t)
$ `% R7 p$ s3 d3 K2 [$ |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)
9 z- ]2 G2 ], A: @" ndy11=(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
plot(t,p)

thinkfunny

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

8 C& f! Y# S8 w7 B  {7 Z" b除法没用点运算，导致p的维度不对，必须是与t维度相同的向量才能拿来跟t一起绘图。