Question

如何用Matlab求参数方程的导数？

已知参数方程x = a * cos(t) ^ 3，y = b * sin(t) ^ 3，用Matlab求dy / dx。如果直接使用：——————————————————————syms a b t;x = a * cos(t) ^ 3;y = b * sin(t) ^ 3;diff(y, 'x')——————————————————————进行计算的话，会得到：——————————————————————ans = 0——————————————————————

Answer 1

1. clear;clc;

2. syms a b t;

3. x = a * cos(t) ^ 3;

4. y = b * sin(t) ^ 3;

5. dydx=diff(y,t)/diff(x,t)

7. %%以下为图形验证

8. a=2;b=3;

9. x=str2func(['@(t)',vectorize(subs(x))]);

10. y=str2func(['@(t)',vectorize(subs(y))]);

11. k=str2func(['@(t)',vectorize(subs(dydx))]);

12. t=linspace(-pi,pi);

13. X=[-a,a];

14. n=0;

15. f=figure;

16. while isvalid(f)

17. n=mod(n,99)+1;

18. T=t(n);

19. hold off

20. plot(x(t),y(t));

21. hold on

22. Y=k(T)*(X-x(T))+y(T);

23. plot(X,Y);

24. plot(x(T),y(T),'ko','markerfacecolor','w');

25. xlim([-a,a])

26. ylim([-b,b])

27. pause(0.02)

28. end

追问

大佬，Matlab是你开发的吧？这么6的吗？
不过我的是Matlab R2017b，用实时脚本（*.mlx）执行到第15行之后就卡在while那了，那个动图一直画不出来……

追答

你好，根据我的试验，实时脚本中的图形本质上是一个在程序执行完毕后的截图，而真正的图像被默认设置为不可见，这样在循环执行完毕前，无法出现这个图形。
解决方案是，将f=figure;修改为f=figure('visible','on');

Answer 2

syms a b t;
x = a * cos(t)^3;
y = b * sin(t)^3;

dx1dt = diff(x,t);

% 一阶导数
dy1dx = diff(y,t)/dx1dt

% 二阶导数
dy2d2x = diff(dy1dx,t)/dx1dt

% 3阶导数
dy3d3x = diff(dy2d2x,t)/dx1dt