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

greensmile

t=0:0.001: (13*pi/18)
& e  r- j7 M: w0 Ex=-(50+32.5*(1-cos(t*18*pi/(13*pi)))).*sin(t)-15*cos(t)
, k2 L1 [. c. v' h% \. oy=(50+32.5*(1-cos(t*18*pi/(13*pi)))).*cos(t)-15*sin(t)
dx1=diff(x)
! M! V! f; X* K1 s, ndy1=diff(y)
5 b5 `  W, P# W  L# zdx11=diff(x,2)
dy11=diff(y,2)
: n0 }6 d' L& F8 C% u$ k
p=abs(sqrt((dx1.^2+dy1.^2).^3)/(dx1.*dy11-dy1.*dx11))
hold on
% G0 r& `" B$ u3 M0 i! s: jplot(t,p)
9 L3 B  c) \4 v/ {5 C* X7 l运算时报错矩阵维度必须一致，然后我尝试改了一下。。
' z) D* u5 h" n6 |% yt=0:0.001: (13*pi/18)
8 ~/ c9 O2 o' ~: m4 A6 ^9 @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)
( f. x! W/ Q/ M+ u) vdx1=diff(x)
1 ^8 ^% G- E3 M. ~7 Zdy1=diff(y)
dx11=diff(x,2)
dy11=diff(y,2)

. W: {% v) v8 W
' Y; i9 E. _7 @8 B, z2 Vdx1=dx1(1:length(dx11))
dy1=dx1(1:length(dy11))

% o) K! N) D% w# a
p=abs(sqrt((dx1.^2+dy1.^2).^3)/(dx1.*dy11-dy1.*dx11))
hold on
plot(t,p)
这个时候代码无报错，但绘图只有坐标轴没有图像

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

9 n; V; H$ S+ V- L

regngfpcb

不太会啊

mutougeda

4 V7 H4 X) D! N0 H) j4 W按照你的代码p是定值，你少写了一个点，正确代码如下
p=abs(sqrt((dx1.^2+dy1.^2).^3)./(dx1.*dy11-dy1.*dx11))

thinkfunny

不是点运算的问题，是离散向量求差分的时候，长度会比原向量少阶数个元素（一阶差分少一个，二阶少两个），你先设置的t向量不可能和后面用差分求出来的p有相同维度的。既然原始的x和y都是t的简单函数，那就直接用解析式求精确微分（而不是用离散向量求差分）。
+ ~8 F, y, P* E& Y+ Q改成
, t& R! x- p2 s5 m1 osyms 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);
: }4 A4 u2 n! H# xdx1 = diff( x );
dy1 = diff( y );
dx11 = diff( x, 2 );
dy11 = diff( y, 2 );
p = abs(sqrt((dx1.^2+dy1.^2).^3)./(dx1.*dy11-dy1.*dx11));
$ c3 g% a; U' v7 x# s6 I: P( s" r4 I- {) kPFunction = matlabFunction( p );
t = 0 : 0.001 : (13*pi/18);
/ z& o$ c; |" I  Pp = PFunction( t );
plot( t, p )

greensmile

thinkfunny 发表于 2020-7-27 18:32
* {1 x9 J& W7 S* g  t+ ?9 O不是点运算的问题，是离散向量求差分的时候，长度会比原向量少阶数个元素（一阶差分少一个，二阶少两个）， ...

% G; B% N$ L5 s: e- p按照你的回答我运行了一下，结果得出的图像感觉有点奇怪。。不太像是曲率半径的图像。
然后我又按照解析式求导，把得出来的公式输入进去，代码可以正常运算了，但是还是只有坐标轴，绘制不出来图像是怎么回事呀。。
这是改正后的代码
; `, v9 f3 ?& k* _  A- {9 a9 jt=0:0.001: (13*pi/18)
x=-(50+32.5*(1-cos(t*18/13))).*sin(t)-15*cos(t)
5 y. e% K1 }$ Z! o, P+ Zy=(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)
' M# l& D- {0 E$ b6 i7 A3 Tdx11=-(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 ~. x3 `2 f1 fdy11=(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)
( A9 B. g0 o- C0 B5 y, R4 ^- f1 ^p=abs(((dx1.^2+dy1.^2).^(3/2))/(dx1.*dy11-dy1.*dx11))
/ Q0 t! d* u: B2 uhold on
7 X  `, c) M2 K# J- Q2 z. m# [plot(t,p)

thinkfunny

greensmile 发表于 2020-7-27 18:34
' t5 Q2 P% C5 `+ T  k# K2 E" n按照你的回答我运行了一下，结果得出的图像感觉有点奇怪。。不太像是曲率半径的图像。
/ \6 A1 p4 v. g然后我又按照解析 ...

除法没用点运算，导致p的维度不对，必须是与t维度相同的向量才能拿来跟t一起绘图。
4 Q! ]( Q. a; c/ ^/ x( u+ n; r