contpath.m

来自「一个基于打靶法的最优控制求解软件 求解过程中采用参数延续算法」· M 代码 · 共 31 行

% Discrete continuation path script


clear all; 
format long; 
 
filename=input('Data file prefix: ','s') 
fid=fopen([filename '.contpath']); 
% Discrete continuation path
text=fgets(fid)
% Problem size and path length
text=fgets(fid)
n = fscanf(fid,'%d',1) 
imax = fscanf(fid,'%d',1) 
% Path
blank=fgets(fid);
text=fgets(fid)
for i=1:imax
  path(i,:) = fscanf(fid,'%f',n+1);
end
fclose(fid) 

for i=1:n
figure;
hold on;
title(['PATH COMPONENT ',int2str(i)],'FontSize',12) 
plot(path(:,i),path(:,n+1),'b+','LineWidth',2);
set(gca,'FontSize',16) 
end