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)
; Z5 P1 o- o  d) G" H  cdx1=diff(x)
" Q- V/ q! C- ^0 C$ }+ `dy1=diff(y)
dx11=diff(x,2)
" E* v3 {5 n- h9 M3 a  |# y& Edy11=diff(y,2)
/ B& c" G7 d1 ?8 O; F7 ?3 ~1 P
p=abs(sqrt((dx1.^2+dy1.^2).^3)/(dx1.*dy11-dy1.*dx11))
hold on
, [: c& U# I3 Uplot(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)
1 E/ I! L7 R; d7 A" j6 oy=(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)

3 {% j# A# j7 q6 t1 q  K) k% x
dx1=dx1(1:length(dx11))
dy1=dx1(1:length(dy11))
# g1 x5 h' ]$ k7 Z% A" p4 v* Z2 q
) c- P& l% t6 a5 P- g
p=abs(sqrt((dx1.^2+dy1.^2).^3)/(dx1.*dy11-dy1.*dx11))
8 Q5 y7 m$ H& H' p3 B3 V3 R- h1 Fhold on
plot(t,p)
这个时候代码无报错，但绘图只有坐标轴没有图像
! i% q; I7 v# ~# p) P" h* o
6 m! D, O9 {# ~/ N5 O1 _希望大神们可以给小弟指一指错误，谢谢大家
/ T6 _; ]) }* B7 b6 R
( \1 r  `- o! c) B

regngfpcb

不太会啊

mutougeda

按照你的代码p是定值，你少写了一个点，正确代码如下
9 L; v( T, H; }, A5 np=abs(sqrt((dx1.^2+dy1.^2).^3)./(dx1.*dy11-dy1.*dx11))

thinkfunny

不是点运算的问题，是离散向量求差分的时候，长度会比原向量少阶数个元素（一阶差分少一个，二阶少两个），你先设置的t向量不可能和后面用差分求出来的p有相同维度的。既然原始的x和y都是t的简单函数，那就直接用解析式求精确微分（而不是用离散向量求差分）。
) V' N0 X! z/ L改成
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);
- a* Y. f* V, }' T. Q0 }* m" Rdx1 = diff( x );
% Y/ |" L6 i+ w: b, [# udy1 = diff( y );
- Z  h6 F9 v6 a, ?dx11 = diff( x, 2 );
% {* {! i) [* [dy11 = diff( y, 2 );
p = abs(sqrt((dx1.^2+dy1.^2).^3)./(dx1.*dy11-dy1.*dx11));
  `7 c" ?& l; C2 L* d( @8 b* EPFunction = matlabFunction( p );
t = 0 : 0.001 : (13*pi/18);
p = PFunction( t );
. d1 F- v$ V- N: Gplot( t, p )

greensmile

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

按照你的回答我运行了一下，结果得出的图像感觉有点奇怪。。不太像是曲率半径的图像。
然后我又按照解析式求导，把得出来的公式输入进去，代码可以正常运算了，但是还是只有坐标轴，绘制不出来图像是怎么回事呀。。
这是改正后的代码
9 ^, E% n6 V2 x8 c6 qt=0:0.001: (13*pi/18)
8 p$ O; U" C- D; \8 I7 k7 xx=-(50+32.5*(1-cos(t*18/13))).*sin(t)-15*cos(t)
" K" x; H# [* M* `$ @y=(50+32.5*(1-cos(t*18/13))).*cos(t)-15*sin(t)
! W% W+ x- O' }# Y$ N& p* }dx1=-(50+32.5*(1-cos(t*18/13))).*cos(t)-(45*sin(t*18/13)-15).*sin(t)
  d2 U8 F$ q% \+ Y. d4 Ody1=-(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)
4 @  w% T! u. udy11=(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)
3 x+ F2 a, o/ b! @1 z/ [1 n3 J" \p=abs(((dx1.^2+dy1.^2).^(3/2))/(dx1.*dy11-dy1.*dx11))
3 Q" b( h, ^1 w! n2 x. P) S6 S- xhold on
plot(t,p)

thinkfunny

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

除法没用点运算，导致p的维度不对，必须是与t维度相同的向量才能拿来跟t一起绘图。