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

greensmile

t=0:0.001: (13*pi/18)
# I8 h8 c. k1 |( cx=-(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)
7 k. ?8 X* Y' p* w( p: ^3 I+ sdy1=diff(y)
dx11=diff(x,2)
dy11=diff(y,2)
* @: y/ C+ l% c3 O4 M& a/ M
5 @: S% ~7 F% j+ u& _; lp=abs(sqrt((dx1.^2+dy1.^2).^3)/(dx1.*dy11-dy1.*dx11))
hold on
# R. o- j5 k5 U) F" m" bplot(t,p)
运算时报错矩阵维度必须一致，然后我尝试改了一下。。
t=0:0.001: (13*pi/18)
( ~; W. [, }" A/ 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)
dx1=diff(x)
dy1=diff(y)
( u# B- S% L3 m5 U& Odx11=diff(x,2)
dy11=diff(y,2)


) i" @! f! v$ q( Q9 ldx1=dx1(1:length(dx11))
+ `% v, ~9 ~9 ^/ o5 ^dy1=dx1(1:length(dy11))
, U3 D: ]6 D+ C3 t* R6 i3 |

. F. X' B: s! w% C5 o( h; Lp=abs(sqrt((dx1.^2+dy1.^2).^3)/(dx1.*dy11-dy1.*dx11))
8 s0 }! u/ u$ G9 chold on
& L  i$ ?3 t0 B* N- h2 H3 l% q! ?plot(t,p)
' }9 T4 ^2 h% g. R' _这个时候代码无报错，但绘图只有坐标轴没有图像

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

% I+ p8 b/ r9 g% I7 k8 `  r+ f, u+ X

regngfpcb

不太会啊

mutougeda

4 l* }9 g, r0 E$ l- Z按照你的代码p是定值，你少写了一个点，正确代码如下
5 k9 |6 E' W' M: lp=abs(sqrt((dx1.^2+dy1.^2).^3)./(dx1.*dy11-dy1.*dx11))

thinkfunny

不是点运算的问题，是离散向量求差分的时候，长度会比原向量少阶数个元素（一阶差分少一个，二阶少两个），你先设置的t向量不可能和后面用差分求出来的p有相同维度的。既然原始的x和y都是t的简单函数，那就直接用解析式求精确微分（而不是用离散向量求差分）。
- [/ ~; w+ E( i9 C- `/ n0 C改成
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 );
dy1 = diff( y );
. r' w5 x& R3 K$ Q1 udx11 = diff( x, 2 );
dy11 = diff( y, 2 );
( A3 N7 L6 ~( Y( q2 i# f" [p = abs(sqrt((dx1.^2+dy1.^2).^3)./(dx1.*dy11-dy1.*dx11));
PFunction = matlabFunction( p );
t = 0 : 0.001 : (13*pi/18);
! [. Q8 }3 r5 ^# m8 @8 I8 _* v4 sp = PFunction( t );
plot( t, p )

greensmile

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

按照你的回答我运行了一下，结果得出的图像感觉有点奇怪。。不太像是曲率半径的图像。
然后我又按照解析式求导，把得出来的公式输入进去，代码可以正常运算了，但是还是只有坐标轴，绘制不出来图像是怎么回事呀。。
9 O+ Y* e4 t. k* c1 L4 ~! ~这是改正后的代码
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)
0 H; x7 B" {, K5 L  Fdx1=-(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)
" `" A, z  m5 \) J8 edx11=-(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)
. Z/ p3 ~( c1 L$ f* C& Ydy11=(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)
! T' I7 {+ Y( ?. g7 c3 ep=abs(((dx1.^2+dy1.^2).^(3/2))/(dx1.*dy11-dy1.*dx11))
hold on
plot(t,p)
; j5 I* b' I  S- R) b% ?' h

thinkfunny

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

除法没用点运算，导致p的维度不对，必须是与t维度相同的向量才能拿来跟t一起绘图。
! _1 r/ j+ W/ z( }3 t) N- Z