洛马普里埃塔地震

1、命令行键入：load quake e n vwhos e n v

2、按“Enter”键。如图1所示。

5、命令行键入：t1 = 8*[1;1];t2 = 15*[1;1];subplot(3,1,1)hold onplot([t1 t2],yrange,'k','LineWidth',2)hold offsubplot(3,1,2)hold onplot([t1 t2],yrange,'k','LineWidth',2)hold offsubplot(3,1,3)hold onplot([t1 t2],yrange,'k','LineWidth',2)hold off

6、按“Enter”键。如图3所示。

9、命令行键入：subplot(1,1,1)k = length(t);k = round(max(1,k/2-100):min(k,k/2+100));plot(e(k),n(k),'.-')xlabel('East'), ylabel('North');title('Acceleration During a One Second Period');

10、按“Enter”键。如图5所示。

13、由位醅呓择锗置数据定义的轨迹可以用三个不同的二维投影显示出来。下面是第一个带有t注解的值的例子。命令行键入：subplot(1,1,1);cla;subplot(2,2,1)plot(npos,vpos,'b');na = max(abs(npos));na = 1.05*[-na na];ea = max(abs(epos));ea = 1.05*[-ea ea];va = max(abs(vpos));va = 1.05*[-va va];axis([na va]);xlabel('North');ylabel('Vertical');nt = ceil((max(t)-min(t))/6);k = find(fix(t/nt)==(t/nt))';for j = k text(npos(j),vpos(j),['o ' int2str(t(j))]);end

14、按“Enter”键。如图7所示。

17、第四个子情节是轨迹的三维视图。命令行键入：subplot(2,2,4)plot3(npos,epos,vpos,'k')for j = k text(npos(j),epos(j),vpos(j),['o ' int2str(t(j))]);endaxis([na ea va]);xlabel('North');ylabel('East')zlabel('Vertical');box on

18、按“Enter”键。如图9所示。