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

greensmile

t=0:0.001: (13*pi/18)
. Z6 R+ r& I8 O+ L' m8 `# nx=-(50+32.5*(1-cos(t*18*pi/(13*pi)))).*sin(t)-15*cos(t)
( y4 U& ?  ]8 K4 ^y=(50+32.5*(1-cos(t*18*pi/(13*pi)))).*cos(t)-15*sin(t)
dx1=diff(x)
9 L2 H# j9 E( {: y0 I6 Ndy1=diff(y)
0 u' ^0 D) K& d3 `; t. ~$ g$ a. ldx11=diff(x,2)
dy11=diff(y,2)

, \  E+ W: w* M& Z- np=abs(sqrt((dx1.^2+dy1.^2).^3)/(dx1.*dy11-dy1.*dx11))
hold on
- z# z1 F$ H4 J: Fplot(t,p)
1 L9 ^1 ~5 [# O  u3 H, u运算时报错矩阵维度必须一致，然后我尝试改了一下。。
4 z' H: j! z3 e( p4 f! Y, e8 c$ Vt=0:0.001: (13*pi/18)
x=-(50+32.5*(1-cos(t*18*pi/(13*pi)))).*sin(t)-15*cos(t)
1 g  [  t: p, Q$ ~1 @1 Sy=(50+32.5*(1-cos(t*18*pi/(13*pi)))).*cos(t)-15*sin(t)
dx1=diff(x)
dy1=diff(y)
5 y( R' n0 Z- S% Y# A( e& Xdx11=diff(x,2)
dy11=diff(y,2)
( O6 ]8 X" t; b2 J+ Z
$ x4 @, v* a6 l  `8 t! P
dx1=dx1(1:length(dx11))
( F  P5 p  z2 H9 o; }dy1=dx1(1:length(dy11))
4 x& S% x1 d8 x; `3 Z
& z5 M. e# [* Y7 T  N) p, V2 Q/ ^$ d/ f
) G5 l9 F7 @* U1 {9 J# [& Yp=abs(sqrt((dx1.^2+dy1.^2).^3)/(dx1.*dy11-dy1.*dx11))
( F8 {. B- e$ C7 n' [0 Dhold on
& x( w1 K- H) B" m8 A2 splot(t,p)
4 T' K& Z) T, ]- _5 g这个时候代码无报错，但绘图只有坐标轴没有图像
2 e# M$ X0 F- `: p
4 |1 b5 x7 j3 ~7 t希望大神们可以给小弟指一指错误，谢谢大家
5 X$ F( }0 p. K2 h6 e
+ o) `$ Y& X, W; }

regngfpcb

不太会啊

mutougeda

8 r8 {% b3 E8 k- L( {( r按照你的代码p是定值，你少写了一个点，正确代码如下
0 v# W+ o! `! @( qp=abs(sqrt((dx1.^2+dy1.^2).^3)./(dx1.*dy11-dy1.*dx11))

thinkfunny

不是点运算的问题，是离散向量求差分的时候，长度会比原向量少阶数个元素（一阶差分少一个，二阶少两个），你先设置的t向量不可能和后面用差分求出来的p有相同维度的。既然原始的x和y都是t的简单函数，那就直接用解析式求精确微分（而不是用离散向量求差分）。
改成
syms t
! V" `! d8 G3 _, ]x = -(50+32.5*(1-cos(t*18*pi/(13*pi)))).*sin(t)-15*cos(t);
  z* p9 U9 o' X/ u3 I/ C$ [y = (50+32.5*(1-cos(t*18*pi/(13*pi)))).*cos(t)-15*sin(t);
dx1 = diff( x );
dy1 = diff( y );
) o2 z; E  ^+ p3 D) N; jdx11 = diff( x, 2 );
dy11 = diff( y, 2 );
9 G' |% x% F& z9 y/ U4 W; Dp = abs(sqrt((dx1.^2+dy1.^2).^3)./(dx1.*dy11-dy1.*dx11));
, `5 E( }. y4 C* {2 l5 ?. A9 dPFunction = matlabFunction( p );
t = 0 : 0.001 : (13*pi/18);
1 L( ~* _9 R& O1 p. k0 |  R- e. Gp = PFunction( t );
plot( t, p )

greensmile

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

2 J$ z5 z; a# N! T; A按照你的回答我运行了一下，结果得出的图像感觉有点奇怪。。不太像是曲率半径的图像。
% N2 \+ O& L4 C& S% }然后我又按照解析式求导，把得出来的公式输入进去，代码可以正常运算了，但是还是只有坐标轴，绘制不出来图像是怎么回事呀。。
这是改正后的代码
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)
& E0 o% l% \4 Z% J) Pdx1=-(50+32.5*(1-cos(t*18/13))).*cos(t)-(45*sin(t*18/13)-15).*sin(t)
, C) ?: T# `1 g, S# @; vdy1=-(50+32.5*(1-cos(t*18/13))).*sin(t)+(45*sin(t*18/13)-15).*cos(t)
. T2 \% c5 \2 c) T8 M7 kdx11=-(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))
" \- ?' g9 |. V9 Z- B/ [9 Khold on
plot(t,p)
+ A4 i5 K+ B  B' x( ?+ L/ I

thinkfunny

greensmile 发表于 2020-7-27 18:34
按照你的回答我运行了一下，结果得出的图像感觉有点奇怪。。不太像是曲率半径的图像。
  K+ p! P+ H+ \$ X! u8 F9 W  k然后我又按照解析 ...

4 X- [" k. @5 t# `! M% C* e8 G除法没用点运算，导致p的维度不对，必须是与t维度相同的向量才能拿来跟t一起绘图。
' j2 D) D( q/ x! S