Продолжение приложения в

amax=std(Y1);

iii=alf;

end

end

set(gcf,'units','normalized','position',[0.15 0.15 .8 .75]);

print(gcf,'-dpng',[p filenameX '.png'],'-r300')

%close(gcf)

% prompt = {'alfa'};

% dlg_title = 'Чемуравноalfa?';

% num_lines = 1;

% def = {num2str(iii)};

% answer = inputdlg(prompt,dlg_title,num_lines,def);

%% вращаем, рисуем ещё раз

am=99999999999999999999999;

amax=0;

figure(2);

fori=1:18

alf=iii-18+(i-1)*2;

[ X1,Y1 ] = Func_rotate(dataXfC,dataYfC,alf );

subplot('position',[.05 .05*(i+.5) .2 .049]);

plot(1:length(X1),X1)

grid on, axis tight, ylim([-r r]);

ylabel(num2str(alf));

ifi>1, set(gca,'xticklabel',''), end

subplot('position',[.3 .05*(i+.5) .2 .049]);

plot(1:length(Y1),Y1)

Продолжение приложения в

gridon, axistight, ylim([-rr]);

ylabel(num2str(alf));

ifi>1, set(gca,'xticklabel',''), end

subplot('position',[.55 .05*(i+.5) .2 .049]);

plot(1:length(X1),X1.^3)

grid on, axis tight, ylim([-r.^3 r.^3]);

ylabel(num2str(alf));

ifi>1, set(gca,'xticklabel',''), end

subplot('position',[.8 .05*(i+.5) .2 .049]);

plot(1:length(Y1),Y1.^3)

grid on, axis tight, ylim([-r.^3 r.^3]);

ylabel(num2str(alf));

ifi>1, set(gca,'xticklabel',''), end

ifstd(Y1)<am

am=std(Y1);

ii=alf;

end

ifstd(Y1)>amax

amax=std(Y1);

iii=alf;

end

end

set(gcf,'units','normalized','position',[0.15 0.15 .8 .75]);

msgbox(num2str(iii))

print(gcf,'-dpng',[p filenameX '__.png'],'-r300')

%close(gcf)

Продолжение приложения В

function [ X1,Y1 ] = Func_rotate( X,Y,alfa )

a=alfa*pi/180;

X1=X*cos(a)-Y*sin(a);

Y1=X*sin(a)+Y*cos(a);

end

ПРИЛОЖЕНИЕ Г

(обязательное)

Программа для расчета азимута землетрясений, разработанная в системе matlab

function sp01(fname,fs,wl,overlap, meaning,zmin,zmax)

%load file

delimeter=' ';

nStrok=9;

data666 = importdata(fname, delimeter, nStrok);

data1=data666.data;

clear data666;

% computing

if meaning>1

data1=AntiTrendFast(data1,meaning);

end

[~,F,T,P]=spectrogram(data1,wl,overlap,wl,fs);

F(1:3,:)=[];

P(1:3,:)=[];

%plot 1st

figure;

subplot('position',[0.04 0.75 0.94 0.22]);

set(gca,'fontSize',9)

plot((1:length(data1))./fs,data1); axis tight; grid on;

aaa=get(gca,'ylim');

aaa(1)=aaa(1)-0.02*(aaa(2)-aaa(1));

aaa(2)=aaa(2)+0.02*(aaa(2)-aaa(1));

Ylim(aaa);

%plot 2nd

subplot('position',[0.04 0.05 0.94 0.62]);

set(gca,'fontSize',9)

surf(T,F,10*log10(P),'edgecolor','none');

aa1=0; aa2=length(data1)/fs; aa3=F(1,1); aa22=F(size(F)); aa4=aa22(1);

axis ([aa1 aa2 aa3 aa4]);

set(gca,'yscale','log');

colorbar('east');

colormap(jet(4096));

set(gca,'clim',[zminzmax]);

Продолжение приложения г

holdon;

a=[.0001 .0001 .0002 .0002 .0003 .0003 .0004 .0004 .0005 .0005 .0006 .0006 .0007 .0007 .0008 .0008 .0009 .0009 ...

.001 .001 .002 .002 .003 .003 .004 .004 .005 .005 .006 .006 .007 .007 .008 .008 .009 .009 ...

.01 .01 .02 .02 .03 .03 .04 .04 .05 .05 .06 .06 .07 .07 .08 .08 .09 .09 ...%9*2*6=108

.1 .1 .2 .2 .3 .3 .4 .4 .5 .5 .6 .6 .7 .7 .8 .8 .9 .9 ...

1 1 2 2 3 3 4 4 5 5 6 6 7 7 8 8 9 9 ...

10 10 20 20 30 30 40 40 50 50 60 60 70 70 80 80 90 90];

b=[aa1-1 aa2+1 aa2+1 aa1-1 aa1-1 aa2+1 aa2+1 aa1-1 aa1-1 aa2+1 aa2+1 aa1-1 aa1-1 aa2+1 aa2+1 aa1-1 aa1-1 aa2+1 ...

aa2+1 aa1-1 aa1-1 aa2+1 aa2+1 aa1-1 aa1-1 aa2+1 aa2+1 aa1-1 aa1-1 aa2+1 aa2+1 aa1-1 aa1-1 aa2+1 aa2+1 aa1-1 ...

aa1-1 aa2+1 aa2+1 aa1-1 aa1-1 aa2+1 aa2+1 aa1-1 aa1-1 aa2+1 aa2+1 aa1-1 aa1-1 aa2+1 aa2+1 aa1-1 aa1-1 aa2+1 ...

aa2+1 aa1-1 aa1-1 aa2+1 aa2+1 aa1-1 aa1-1 aa2+1 aa2+1 aa1-1 aa1-1 aa2+1 aa2+1 aa1-1 aa1-1 aa2+1 aa2+1 aa1-1 ...

aa1-1 aa2+1 aa2+1 aa1-1 aa1-1 aa2+1 aa2+1 aa1-1 aa1-1 aa2+1 aa2+1 aa1-1 aa1-1 aa2+1 aa2+1 aa1-1 aa1-1 aa2+1 ...

aa2+1 aa1-1 aa1-1 aa2+1 aa2+1 aa1-1 aa1-1 aa2+1 aa2+1 aa1-1 aa1-1 aa2+1 aa2+1 aa1-1 aa1-1 aa2+1 aa2+1 aa1-1];

c=3000*ones(1,108);

plot3(b,a,c,':b'); hold off;