plot_traj.m

来自「一个matlab的将军模型」· M 代码 · 共 22 行

function plot_traj(file_name,SCSB_name,initial_conditions)
pwd
for k=1:length(initial_conditions)
    vert_init=vertices(initial_conditions{k});
    hold on
    place_X0=strcat(file_name,'/',SCSB_name);
    for i=1:length(vert_init)
        verpoint=vert_init(i);
        set_param(place_X0,'X0',strcat('[',num2str(verpoint'),']'));
        sim(file_name);
        if length(verpoint)==2
            plot(simula(:,1),simula(:,2));
        elseif length(verpoint)==3
            plot3(simula(:,1),simula(:,2),simula(:,3));
        else
            printf('dimension invalid\n')
            return
        end
    end
end

return