5/27/2010
5/25/2010
5/17/2010
5/14/2010
First wikipedia article
Take two people for a few days, but looks nice now.
http://en.wikipedia.org/wiki/Shear_wave_splitting
http://en.wikipedia.org/wiki/Shear_wave_splitting
5/11/2010
Try a different aproach
The way to use em and muon particles have been failed,
may be I should try the area of muon and total signal,
give it a try!
may be I should try the area of muon and total signal,
give it a try!
5/09/2010
shear wave splitting_matlab_code
% animation.m
%
%
% Author : Pieter Coppens
% Date : 8-11-2008
% Modified by:Guofeng Yuan
% Date:5-08-2010
% close all windows and clear all variables
clear;
clc;
sr=0.2;% sampling time interval
fr=1/sr; %sampling rate
max=15; %total time
split=5; %the time when the wave splits
fs = 100;
t_end =1; % end of the time interval
t = 0:1/fs:t_end;
j = length(t);
z1=zeros(1,101);
z2=z1;
%creat the vedio file
mov = avifile('shear_wave_splitting.avi','quality',50,'compression','None','fps',5)
for i=0:sr:max
if i<=split %before the splits
x0=t+i;
y0= sin(2*pi*t)*cos(pi/4);
z0=sin(2*pi*t)*sin(pi/4);
plot3(x0(1:j),y0(1:j),z0(1:j),'-r','LineWidth',1.5)
title(['Shear wave splitting']) % title
xlabel(['Time = ',num2str(i,'%.2f'),'s']) % the x-axis shows the time-instant in seconds.
hold off
axis([0 25 -1 1 -1 1 -1 1 ])
grid on
axis on
pause(0.1)
else %after the splits
x1=t+i; %the slow wave
y1=0.5*sin(2*pi*t);
x2=t+1.8*(i-split)+split; %the fast wave
y2=0.87*sin(2*pi*t);
plot3(x1(1:j),y1(1:j),z1(1:j),x2(1:j),z2(1:j),y2(1:j),'-m','LineWidth',1.2)
hold off
title(['Shear wave splitting']) % title
xlabel(['Time = ',num2str(i,'%.2f'),'s']) % the x-axis shows the time-instant in seconds.
axis([0 25 -1 1 -1 1 -1 1 ])
grid on
axis on
pause(0.1)
end
F = getframe(gcf); % save the current figure
mov = addframe(mov,F); % add it as the next frame of the movie
end
% close the open windows and clear all variables
mov = close(mov);
close all
clear all
% Open the .avi-file in the external .avi-movie application of windows.
% Comment this if you're not running this file under a Windows environment.
winopen('shear_wave_splitting.avi') %open the movie
订阅:
评论 (Atom)
