二元二阶微分方程组求解，并画出极坐标图（matlab）

使用matlab解决二元二阶微分方程组的求解问题，并画出包括极坐标图在内的多幅变量间的关系图

% testode45.m
%-----------------------------------------------------
% 清空所有变量
clear
% 清空屏幕
clc

% 时间跨度取0-0.5，可以空格分隔，也可以用逗号分隔
tspan = [01.5];
% 初始值
y0 = [1.732-50.2pi/2];
% 调用语句
[TY] = ode45（ @odefun tspan y0 ）;
% 绘图
figure（1）;
plot（TY（:1）‘-m‘）
xlabel（‘时间/s‘‘Fontsize‘18）;
ylabel（‘R的变化率（m/s）‘‘Fontsize‘18）;
title（‘R的变化率随时间推移的值‘‘Fontsize‘18）
legend（‘R’‘）
figure（2）;
plot（TY（:2）‘-r‘）
xlabel（‘时间/s‘‘Fontsize‘18）;
ylabel（‘θ的变化率（rad/s）‘‘Fontsize‘18）;
title（‘θ的变化率随时间推移的值‘‘Fontsize‘18）
legend（‘θ’‘）
figure（3）;
plot（TY（:3）‘-b‘）
xlabel（‘时间/s‘‘Fontsize‘18）;
ylabel（‘R（m）‘‘Fontsize‘18）;
title（‘R随时间推移的值‘‘Fontsize‘18）
legend（‘R‘）
figure（4）;
polarplot（Y（:4）Y（:3）‘-k‘）
rlim（[0 0.5]）;
thetalim（[0 360]）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2019-03-31 00:20  二元二阶微分方程组matlab\
     文件       21500  2019-03-19 19:31  二元二阶微分方程组matlab\1.jpg
     文件       21673  2019-03-19 19:36  二元二阶微分方程组matlab\2.jpg
     文件       19286  2019-03-19 19:36  二元二阶微分方程组matlab\3.jpg
     文件       17184  2019-03-19 19:36  二元二阶微分方程组matlab\4.jpg
     文件       16685  2019-03-19 20:34  二元二阶微分方程组matlab\5.jpg
     文件         863  2019-03-19 20:34  二元二阶微分方程组matlab\main.m
     文件         139  2019-03-19 18:55  二元二阶微分方程组matlab\odefun.m
     文件       54984  2019-03-31 00:20  二元二阶微分方程组matlab\题目.docx
     文件       21500  2019-03-19 19:31  动力学matlab\1.jpg
     文件       21673  2019-03-19 19:36  动力学matlab\2.jpg
     文件       19286  2019-03-19 19:36  动力学matlab\3.jpg
     文件       17184  2019-03-19 19:36  动力学matlab\4.jpg
     文件       16685  2019-03-19 20:34  动力学matlab\5.jpg
     文件         863  2019-03-19 20:34  动力学matlab\main.m
     文件         139  2019-03-19 18:55  动力学matlab\odefun.m
     目录           0  2019-03-19 22:25  动力学matlab\