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)
2 d3 [( M7 J8 X/ Z; T; f: N  wy=(50+32.5*(1-cos(t*18*pi/(13*pi)))).*cos(t)-15*sin(t)
" H5 _& b( |- X' [4 ^- tdx1=diff(x)
dy1=diff(y)
dx11=diff(x,2)
! x5 c5 \' L' P2 n+ cdy11=diff(y,2)
1 e4 X  \/ D' a% @2 C* X+ [+ o
$ r4 }5 B1 M+ F9 z! i7 Z& Pp=abs(sqrt((dx1.^2+dy1.^2).^3)/(dx1.*dy11-dy1.*dx11))
/ L( ]. T$ x+ jhold on
plot(t,p)
运算时报错矩阵维度必须一致，然后我尝试改了一下。。
  r) i' \8 b0 J; M: ~; At=0:0.001: (13*pi/18)
0 [* s# B# R1 gx=-(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)
0 m5 [* l' H+ u0 l& e4 Vdy1=diff(y)
dx11=diff(x,2)
dy11=diff(y,2)

5 k, d( e$ S4 q% P0 V: a" R
dx1=dx1(1:length(dx11))
dy1=dx1(1:length(dy11))
4 C/ ~  M4 A2 @4 V4 W
9 T$ g3 g$ I+ a. J
* }* ^9 N2 d0 E$ ~3 gp=abs(sqrt((dx1.^2+dy1.^2).^3)/(dx1.*dy11-dy1.*dx11))
hold on
plot(t,p)
这个时候代码无报错，但绘图只有坐标轴没有图像

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

regngfpcb

不太会啊

mutougeda

( q) |/ ]* ~* T9 z+ M  A按照你的代码p是定值，你少写了一个点，正确代码如下
p=abs(sqrt((dx1.^2+dy1.^2).^3)./(dx1.*dy11-dy1.*dx11))

thinkfunny

不是点运算的问题，是离散向量求差分的时候，长度会比原向量少阶数个元素（一阶差分少一个，二阶少两个），你先设置的t向量不可能和后面用差分求出来的p有相同维度的。既然原始的x和y都是t的简单函数，那就直接用解析式求精确微分（而不是用离散向量求差分）。
改成
) M- a, ]  G& ]3 t: A/ _syms t
x = -(50+32.5*(1-cos(t*18*pi/(13*pi)))).*sin(t)-15*cos(t);
! s/ u) ]& {3 y! `$ @* sy = (50+32.5*(1-cos(t*18*pi/(13*pi)))).*cos(t)-15*sin(t);
dx1 = diff( x );
3 ?* E# f' ]- u, o& Rdy1 = diff( y );
4 U5 O6 j3 H3 H1 q1 B% W  adx11 = diff( x, 2 );
5 v6 y( ?6 _) {2 I0 l+ M6 s" Jdy11 = diff( y, 2 );
p = abs(sqrt((dx1.^2+dy1.^2).^3)./(dx1.*dy11-dy1.*dx11));
PFunction = matlabFunction( p );
% F8 e+ h7 C- r: }: Nt = 0 : 0.001 : (13*pi/18);
p = PFunction( t );
* G! N8 |# Q4 Z3 T" gplot( t, p )

greensmile

thinkfunny 发表于 2020-7-27 18:32
. d/ O4 X" Q- |# a, [7 M4 X1 K6 c6 ]不是点运算的问题，是离散向量求差分的时候，长度会比原向量少阶数个元素（一阶差分少一个，二阶少两个）， ...

* `1 q  F1 A, y* A3 V% i按照你的回答我运行了一下，结果得出的图像感觉有点奇怪。。不太像是曲率半径的图像。
然后我又按照解析式求导，把得出来的公式输入进去，代码可以正常运算了，但是还是只有坐标轴，绘制不出来图像是怎么回事呀。。
% V+ Q9 y* m. C. ~: ^. {这是改正后的代码
t=0:0.001: (13*pi/18)
x=-(50+32.5*(1-cos(t*18/13))).*sin(t)-15*cos(t)
/ j9 u9 s' |' [1 `9 p1 R5 Hy=(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)
$ {9 C  Z3 H, h7 k3 D" A! ndy1=-(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)
; V7 z  w; X% \- b* q* c) i$ V* q% Ndy11=(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))
3 @; s3 j& m: `hold on
+ @0 T+ c) e. }: D/ [7 `8 }+ cplot(t,p)

thinkfunny

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

3 i- \8 G8 T1 k. \; s- f1 [# t9 [除法没用点运算，导致p的维度不对，必须是与t维度相同的向量才能拿来跟t一起绘图。