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

greensmile

t=0:0.001: (13*pi/18)
( S  |: e; g/ J$ }  v3 kx=-(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)
7 a% x2 I3 z: B( b3 W; y0 R* sdx1=diff(x)
dy1=diff(y)
dx11=diff(x,2)
+ x" V& T# D/ H" Fdy11=diff(y,2)
! R: s4 Q& l& y
p=abs(sqrt((dx1.^2+dy1.^2).^3)/(dx1.*dy11-dy1.*dx11))
- w. x5 S1 o1 P7 \5 `: g* W* fhold on
plot(t,p)
3 }2 i+ A( t, ^运算时报错矩阵维度必须一致，然后我尝试改了一下。。
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)
) }% N& x( ~' M2 P. _4 p+ n9 ndx1=diff(x)
dy1=diff(y)
9 s! L+ T" m$ i! {" |5 ^dx11=diff(x,2)
dy11=diff(y,2)
' o  _- K) {) M7 U
. p0 B2 {0 Z- `) A% I) \% D
, W4 k. |* G8 K. M, Hdx1=dx1(1:length(dx11))
dy1=dx1(1:length(dy11))
/ ^% a" p0 c* [6 v& g
" L4 I( x" d- Z. Q, k
& J; T: b. x# Q- {! Q" g4 }- ^p=abs(sqrt((dx1.^2+dy1.^2).^3)/(dx1.*dy11-dy1.*dx11))
' V% W& t$ Z7 e. L, ~5 j7 r7 mhold on
plot(t,p)
5 {, K: S6 l# F! h' \$ s8 Z这个时候代码无报错，但绘图只有坐标轴没有图像

希望大神们可以给小弟指一指错误，谢谢大家
1 q9 r! j' P; [9 R0 g6 n. v

regngfpcb

不太会啊

mutougeda

( R% y3 w# h2 k5 S* j" X+ k: h按照你的代码p是定值，你少写了一个点，正确代码如下
% P: q, ^% L( T% a# Bp=abs(sqrt((dx1.^2+dy1.^2).^3)./(dx1.*dy11-dy1.*dx11))

thinkfunny

不是点运算的问题，是离散向量求差分的时候，长度会比原向量少阶数个元素（一阶差分少一个，二阶少两个），你先设置的t向量不可能和后面用差分求出来的p有相同维度的。既然原始的x和y都是t的简单函数，那就直接用解析式求精确微分（而不是用离散向量求差分）。
; F8 X* W2 p. R% Y* N0 M; _2 C改成
syms t
, k0 T# t& U: Y5 x) Px = -(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 );
/ \" n! J) }- ldx11 = diff( x, 2 );
& ^9 I0 }8 m: Z% v3 P$ Bdy11 = diff( y, 2 );
, E( B2 y" H% b' J! s$ gp = abs(sqrt((dx1.^2+dy1.^2).^3)./(dx1.*dy11-dy1.*dx11));
PFunction = matlabFunction( p );
4 A7 D# ~/ }3 \1 T7 st = 0 : 0.001 : (13*pi/18);
p = PFunction( t );
plot( t, p )

greensmile

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

按照你的回答我运行了一下，结果得出的图像感觉有点奇怪。。不太像是曲率半径的图像。
然后我又按照解析式求导，把得出来的公式输入进去，代码可以正常运算了，但是还是只有坐标轴，绘制不出来图像是怎么回事呀。。
$ c1 p7 c; h" [+ ^" `/ {( @  a- w这是改正后的代码
t=0:0.001: (13*pi/18)
x=-(50+32.5*(1-cos(t*18/13))).*sin(t)-15*cos(t)
y=(50+32.5*(1-cos(t*18/13))).*cos(t)-15*sin(t)
6 _2 x2 B) L' t; i) A& Y( cdx1=-(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)
dy11=(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))
- r7 j& N( A( P/ P+ j; Y! Bhold on
7 E7 l- x) l  N. F, u6 Xplot(t,p)

thinkfunny

greensmile 发表于 2020-7-27 18:34
6 W5 |; G/ [# L+ f6 y按照你的回答我运行了一下，结果得出的图像感觉有点奇怪。。不太像是曲率半径的图像。
7 K2 m% }/ `% ^  l  V# Y然后我又按照解析 ...

  b8 A$ D$ m" M0 a& R, N8 n8 E+ U除法没用点运算，导致p的维度不对，必须是与t维度相同的向量才能拿来跟t一起绘图。