matlab图形实现

img


用matlab实现图形。

请把这个具体的曲线给出来，不然乱猜可不行
你说的是这种线？

% Simulation and animation of a spherical pendulum.
%
%% 
clear ; close all ; clc
%% Scenario
% Parameters
l = 3;                              % Length                [m]
% Initial conditions
theta0  = pi/6;                     % Polar angle           [rad]
phi0    = 0;                        % Azimuth               [rad]
dtheta0 = 0;                        % d/dt(theta)           [rad/s]
dphi0   = 5;                        % d/dt(phi)             [rad/s]
x0 = [theta0 phi0 dtheta0 dphi0]';
% Parameters
tf      = 300;                       % Final time            [s]
fR      = 30;                       % Frame rate            [fps]
dt      = 1/fR;                     % Time resolution       [s]
time    = linspace(0,tf,tf*fR);     % Time                  [s]
%% Simulation
[TOUT,XOUT] = ode45(@(t,x) pendulum(t,x,l),time,x0);
% Retrieving states
alpha   = XOUT(:,1);
beta    = XOUT(:,2);
dalpha  = XOUT(:,3);
dbeta   = XOUT(:,4);
% Coordinates
rx = l*sin(alpha).*cos(beta);
ry = l*sin(alpha).*sin(beta);
rz = -l*cos(alpha);
%% Animation
% Min-max axis
min_x = min(rx)-l/2;
max_x = max(rx)+l/2;
min_y = min(ry)-l/2;
max_y = max(ry)+l/2;
min_z = min(rz)-l/2;
max_z = 0;
figure
% set(gcf,'Position',[50 50 1280 720]) % YouTube: 720p
% set(gcf,'Position',[50 50 854 480]) % YouTube: 480p
set(gcf,'Position',[50 50 640 640]) % Social
% Create and open video writer object
v = VideoWriter('spherical_pendulum.mp4','MPEG-4');
v.Quality = 100;
open(v);
    
hold on ; grid on ; axis equal
set(gca,'CameraPosition',[42.0101   30.8293   16.2256])
set(gca,'XLim',[min_x max_x])
set(gca,'YLim',[min_y max_y])
set(gca,'ZLim',[min_z max_z])
set(gca,'XTickLabel',[],'YTickLabel',[],'ZTickLabel',[])
for i = 1:length(rx)
    
    cla
    % Vertical line
    plot3([0 0],[0 0],[0 -l],'k--')
    % Point fix
    p = plot3(0,0,0,'Marker','*','Color','k','MarkerSize',10);
    % Pendulum trajectory
    plot3(rx(1:i),ry(1:i),rz(1:i),'b')
    % Pendulum rod
    plot3([0 rx(i)],[0 ry(i)],[0 rz(i)],'r')
    % Pendulum sphere
    plot3(rx(i),ry(i),rz(i),'Marker','o','Color','k','MarkerFaceColor','r','MarkerSize',10);
    % Projections
    plot3(min_x*ones(1,i),ry(1:i),rz(1:i),'g')
    plot3(rx(1:i),min_y*ones(1,i),rz(1:i),'g')
    plot3(rx(1:i),ry(1:i),min_z*ones(1,i),'g')
    frame = getframe(gcf);
    writeVideo(v,frame);
    
end
close(v);
 
function dx = pendulum(~,x,l)
    
    % Parameters
    g = 9.81;                       % Gravity               [m/s2]
    % States
    x1 = x(1);
    % x2 = x(2);
    x3 = x(3);
    x4 = x(4);
    % State equations
    dx1 = x3;
    dx2 = x4;
    dx3 = (l*x4^2*sin(x1)*cos(x1) - g*sin(x1))/l;
    dx4 = -2*x3*x4/tan(x1);
    dx = [dx1 dx2 dx3 dx4]';
end

img