practise.rar_从小到大排列

function [z] = practise( input_args ) %UNTITLED Summary of this function goes here % Detailed explanation goes here %load('D:\My Documents\MATLAB\totalv.txt') format long set(0,'defaultfigurecolor','w'); fni=input('wenjian','s'); fid=fopen(fni,'r'); x1=fscanf(fid,'%f',[3,inf]); x=x1'; m=length(x); y=x(:,[1,3]); for j=1:m-1 for i=1:m-j if y(i,1)>y(i+1,1) z(j,1)=y(i+1,1); z(j,2)=y(i+1,2); else n1=y(i+1,1); n2=y(i+1,2); y(i+1,1)=y(i,1); y(i+1,2)=y(i,2); y(i,1)=n1; y(i,2)=n2; z(j,1)=y(i+1,1); z(j,2)=y(i+1,2); end end end z set(gca,'XTick',-0.5:0.1:0.5); z(m,1)=y(1,1); z(m,2)=y(1,2); p=z(:,1); q=z(:,2); plot(p,q,'.'); set(gca,'LineWidth',1.5); hold on xx=-0.5:1/1000:0.5; yy=spline(z(:,1),z(:,2),xx); plot(xx,yy,'r','LineWidth',1.4) z1=-sqrt(0.5^2/2); z2=-sqrt(0.2^2/2); z3=-sqrt(0.1^2/2); z4=sqrt(0.1^2/2); z5=sqrt(0.2^2/2); z6=sqrt(0.5^2/2); hold on YL = ylim; ymin=YL(1); ymax=YL(2); N=(ymax-ymin)/1000; t=ymin:N:ymax; plot(z1,t,'k') hold on plot(z2,t,'g') hold on plot(z3,t,'y') hold on plot(z4,t,'y') hold on plot(z5,t,'g') hold on plot(z6,t,'k') hold off xlabel('x坐标') ylabel('位移') set(get(gca,'xlabel'),'fontsize',21); set(get(gca,'ylabel'),'fontsize',21); end