MATLAB时域绘图代码.docxVIP

function PlotTime% data pre-processing for JZ202 Physical model% clear all% close all% clc% sensor=mnode.txt;% sensorfilename=mnode;% datapath=signal.ASC;% datafilename=signal;% judge=0;% datastart=1500;% dataend=2000;% column=[1:50];% temp=1; global datapath; global datafilename;global datastart;global dataend;global judge;;global sensor;global sensorfilename;str1=[load sensor];eval(str1);str2=[arran= sensorfilename ;];eval(str2);str3=[load datapath];eval(str3);str4=[xx= datafilename ;];eval(str4);% judge=0;% datastart=2000;% dataend=2500;% xy=datapath;% if temp==1% column=[];% for i=1:(size(xy,2)-1)% column=[column i];% end;% end;column=arran(:,1);xy=[]; for i=1:length(column); xy=[xy xx(:,(column(i)+1))];endif judge==1 datastart=(datastart-xx(1,1))/(xx(2,1)-xx(1,1))+1; dataend=(dataend-xx(1,1))/(xx(2,1)-xx(1,1))+1;elseif judge==0 datastart=(xx(1,1)-xx(1,1))/(xx(2,1)-xx(1,1))+1; dataend=(xx(end,1)-xx(1,1))/(xx(2,1)-xx(1,1))+1;else datastart=datastart; dataend=dataend;endsensor_dir=arran(:,3)./abs(arran(:,3));sensor_dir=sensor_dir;%sign_sensor_jz202 = [1 -1 -1 1 -1 1 -1 -1 -1 1 -1 1 -1 1 -1 1 1 -1 1 -1 1 1 1 1 1 1 -1 1 -1 -1 1 1 1 -1 -1 1 1 -1 -1 1 1 1 1 1 -1 -1 -1 -1 1 -1]; % according to the node orderxy=[xx(:,1) xy];yy=xy(:,2:end);% for i=1:length(column)% yy=[yy xy(:,column(i))];% end;for i=1:size(yy,2); yy(:,i)=yy(:,i).*sensor_dir(i);end;% yy(:,nt1)=yy(:,nt1).*sin(84/180*pi);yy(:,nt2)=yy(:,nt2).*sin(86/180*pi);yy(:,nt3)=yy(:,nt3).*sin(83/180*pi);% yy(:,nt4)=yy(:,nt4).*sin(84/180*pi);yy(:,nt5)=yy(:,nt5).*sin(83/180*pi);% for i=1:length(nt1);yy(:,nt1(i))=yy(:,nt1(i)).*sin(83/180*pi);end;% y2=[yy(:,:)]./800.*9.8;y2=[y2 y(:,1)];%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% xch=1:size(yy,2); yy=detrend(yy(:,xch),constant); yy=[yy xy(:,1)];% [va,nn]=max(abs(yy(:,1))); dt=yy(2,end)-yy(1,end); % ys=1.09;% nn1=ys/dt;nn2=(ys+