Light_beam_deflection_Constraint_parabole_CS_Rbf_matlab

clc,clear; taille = 41; xs = 0; ys = 0; [nx,ny] = meshgrid(1:taille,1:taille); u = 1:13; u = [u fliplr(u(1,1:end-1))]; u = u(1,1:end-1); % Display parameters time_lapse = 0.25; title_text = '光束偏转与引力场的关系'; filename = '光束偏转与引力场的关系图.gif'; titleString = '引力场强度'; cmap = true; % Display settings h = figure; set(h,'Position',get(0,'ScreenSize')); set(gcf,'Color',[0 0 0]); axis tight manual; for t = 1:length(u) sigma = 0.04*u(1,t); % RBF feature parameter Phi = CS_Rbf(taille,xs,ys,-u(t),sigma); mesh(Phi), hold on; alpha 0; milieu = ceil(0.5*taille)*ones(2,1); W1 = Constraint_parabole(-u(1,t),0,sqrt(0.25*((taille-1)^2 + (taille-1)^2)),0.25*pi); W2 = Constraint_parabole(-u(1,t),0,sqrt(0.25*((taille-1)^2 + (taille-1)^2)),1.25*pi); % Compute interpolated altitude W1z = griddata(ny-milieu(1,1),nx-milieu(2,1),Phi,W1(1,:),W1(2,:),'cubic'); W2z = griddata(ny-milieu(1,1),nx-milieu(2,1),Phi,W2(1,:),W2(2,:),'cubic'); W1 = cat(1,W1,W1z); W2 = cat(1,W2,W2z); nb = 0.5*length(W1); for s = 1:nb if s > 1 plot3([W1(1,s-1),W1(1,s)]+milieu(1,1),[W1(2,s-1),W1(2,s)]+milieu(2,1),[W1(3,s-1),W1(3,s)],'Color','y','Linewidth',2.5), hold on; plot3([W1(1,nb+s-1),W1(1,nb+s)]+milieu(1,1),[W1(2,nb+s-1),W1(2,nb+s)]+milieu(2,1),[W1(3,nb+s-1),W1(3,nb+s)],'Color','y','Linewidth',2.5), hold on; plot3([W2(1,s-1),W2(1,s)]+milieu(1,1),[W2(2,s-1),W2(2,s)]+milieu(2,1),[W2(3,s-1),W2(3,s)],'Color','y','Linewidth',2.5), hold on; plot3([W2(1,nb+s-1),W2(1,nb+s)]+milieu(1,1),[W2(2,nb+s-1),W2(2,nb+s)]+milieu(2,1),[W2(3,nb+s-1),W2(3,nb+s)],'Color','y','Linewidth',2.5), hold on; end end plot3(milieu(1,1),milieu(2,1),max(max(Phi)),'o','Color',[0 0 1],'Linewidth',s), hold on; ax = gca; ax.Clipping = 'off'; set(ax,'Color',[0 0 0]); axis([1 51 1 51 -10 0]), axis off; title(title_text,'FontSize',16,'Color',[1 1 1]); if cmap colormap('hot'); hcb = colorbar; set(hcb,'YDir','reverse'); colorTitleHandle = get(hcb,'Title'); set(colorTitleHandle ,'String',titleString,'Color',[1 1 1],'Fontsize',14); else colormap([1 1 1]); end view(-22,80); drawnow; frame = getframe(h); im = frame2im(frame); [imind,cm] = rgb2ind(im,256); % Write to the .gif file if t == 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