引力波_LISA_GW_deformation_effect_animation_matlab

clc,clear; step = pi/15; phase = 10*pi:-step:step; H = 31; W = 31; [Y,Z] = meshgrid(1:W,1:H); Y = Y - floor(0.5*W); Z = Z - floor(0.5*H); Amplitude = 12; x_shift = -25; % Rotation matrix x axis Rm = @(Theta)[cos(Theta),-sin(Theta); sin(Theta),cos(Theta)]; smaxa = linspace(1,2,floor(0.5*length(phase))+1); % semi-major axis smaxa = [smaxa fliplr(smaxa(1,2:end-1))]; smaxa = [smaxa(floor(0.75*length(smaxa))+1:end) smaxa(1:floor(0.75*length(smaxa)))]; sminb = linspace(2,1,floor(0.5*length(phase))+1); % semi-minor axis sminb = [sminb fliplr(sminb(1,2:end-1))]; sminb = [sminb(floor(0.75*length(sminb))+1:end) sminb(1:floor(0.75*length(sminb)))]; % Display parameters time_lapse = 0.2; title_text = '引力波变对LISA航天器运动的影响'; title_on = true; filename = '地球引力波对LISA影响.gif'; polarization = '+'; az = 49; el = 27; % Display settings h = figure; set(h,'Position',get(0,'ScreenSize')); set(gcf,'Color',[0 0 0]); axis tight manual; for k = 1:length(phase) signal = Gw_model2(Amplitude,phase(k),polarization); M = sin_2D_wave(H,W,phase(1,length(phase)-k+1),smaxa(1,k),sminb(1,k),polarization)+x_shift; % space delay in x axis; default value: -25 [C,V,ny,nz,Tc,nb_pts] = LISA_triangle(M,Y,Z); mesh(M,Y,Z,M), hold on; colormap('bone'); alpha 0; % Draw wave if strcmp(polarization,'+') if k <= floor(0.5*length(phase)) line(signal(1,:),1+signal(2,:),1+signal(3,:),'Color',[0 1 1],'Linewidth',1), hold on; line(signal(1,:),1+signal(4,:),1+signal(5,:),'Color',[0 0 1],'Linewidth',2), hold on; line(signal(1,:),1-signal(4,:),1-signal(5,:),'Color',[0 0 1],'Linewidth',2), hold on; elseif k > floor(0.5*length(phase)) line(signal(1,:),1+signal(3,:),1+signal(2,:),'Color',[0 1 1],'Linewidth',1), hold on; line(signal(1,:),1+signal(5,:),1+signal(4,:),'Color',[0 0 1],'Linewidth',2), hold on; line(signal(1,:),1-signal(5,:),1-signal(4,:),'Color',[0 0 1],'Linewidth',2), hold on; end elseif strcmp(polarization,'x') if k <= floor(0.5*length(phase)) line(signal(1,:),1+signal(2,:),1+signal(3,:),'Color',[0 1 1],'Linewidth',1), hold on; line(signal(1,:),1+signal(4,:),1+signal(5,:),'Color',[0 0 1],'Linewidth',2), hold on; line(signal(1,:),1-signal(4,:),1-signal(5,:),'Color',[0 0 1],'Linewidth',2), hold on; elseif k > floor(0.5*length(phase)) line(signal(1,:),1+signal(2,:),1-signal(3,:),'Color',[0 1 1],'Linewidth',1), hold on; line(signal(1,:),1+signal(4,:),1-signal(5,:),'Color',[0 0 1],'Linewidth',2), hold on; line(signal(1,:),1-signal(4,:),1+signal(5,:),'Color',[0 0 1],'Linewidth',2), hold on; end end % Signal axis plot3([signal(1,1) signal(1,end)],[1 1],[1 1],'Color',[1 0 1]), hold on; plot3([signal(1,end) signal(1,end)],[1 -Amplitude],[1 1],'Color',[1 0 1]), hold on; plot3([signal(1,end) signal(1,end)],[1 1],[1 Amplitude],'Color',[1 0 1]), hold on; % Draw LISA triangles for m = 1:3 T = Rm(phase(1,k))*V; X = interp2(Y,Z,reshape(C(1,:),nz,ny),1+T(2,:),1+T(1,:),'cubic'); line(X,1+T(1,:),1+T(2,:),'Color',[1 0 0],'Linewidth',3), hold on; plot3(X(1,1),1+T(1,1),1+T(2,1),'o','Color',[1 0 1],'Linewidth',3), hold on; plot3(X(1,(nb_pts+2)/3),1+T(1,(nb_pts+2)/3),1+T(2,(nb_pts+2)/3),'o','Color',[0 1 1],'Linewidth',3), hold on; plot3(X(1,2*(nb_pts+2)/3-1),1+T(1,2*(nb_pts+2)/3-1),1+T(2,2*(nb_pts+2)/3-1),'o','Color',[1 1 0],'Linewidth',3), hold on; end % Features & other settings ax = gca; ax.Clipping = 'off'; set(ax,'Color',[0 0 0]); axis square, axis equal, axis off; if title_on text(-2,12,19,title_text,'Color',[1 1 1],'FontSize',16); end view([az,el]); camdolly(0,0.09,0); zoom(1.5); drawnow; frame = getframe(h); im = frame2im(frame); [imind,cm] = rgb2ind(im,256); % Write to the .gif file if k == 1 imwrite(imind,cm,filename,'gif', 'Loopcount',Inf,'DelayTime',time_lapse); else imwrite(imind,cm,filename,'gif','WriteMode','append','DelayTime',time_lapse); end clf; end