How to make animtion of 3D surface plot?

I'm having my code here where I want to put animation. Can we show this animation wrt to any parameter?

[Constantino Carlos Reyes-Aldasoro]

To create an animation you will need to plot continuously different plots. Say you want to have 10 or 20 frames in your animation you could do this with a for loop:

for k=1:10
    % use your code here, change the variables and surf each time
end

Now, you will need to grab each plot, so after each surf add the following:

drawnow
F(k) = getframe;

The variable F will contain your animation and then you can save as a video like this

 v = VideoWriter('video_name'), 'MPEG-4');
            open(v);
            writeVideo(v,F);
            close(v);

Hope this helps. If this does not solve your problem, please let me know. If it does, please accept the answer.

[priya anjali]

Thank you, but it's showing some error about the function.

Answer 1

[Ameer Hamza]

See this example using variable 't'. Adapt it according to your requirement

%suface plot;
x=-20:1:20;
y=-20:1:20;
[x,y]=meshgrid(x,y);
c1=0.1576;
A1=0.9706;
A2=0.9572;
A3=0.4854;
A4=2*A3;
T = linspace(0, 1, 100);

fig = figure;
ax = axes();
hold(ax);
view(3);
grid on
colormap jet;
xlabel('\bf{x}','fontsize',18);
ylabel('\bf{y}','fontsize',18);
zlabel('\bf{u (x,y,t)}','fontsize',18);
s = surf(x,y,zeros(size(x)),'Facealpha','1'); % temporary surface to get the handle

v = VideoWriter('test.mp4', 'MPEG-4');
open(v);
for i=1:numel(T)
    t=T(i);
    b=cos(t);
    A=log(b);
    s.ZData = c1./((x./(A1+3*t).^1./3+A-(A2+y)./(A1+3*t).^1./3-A4)*(A1+3*t).^1./3);
    drawnow;
    
    frame = getframe(fig);
    writeVideo(v, frame)
end
close(v);

[priya anjali]

Thankyou very much. Is this possible in 2D?

[Ameer Hamza]

Try this

t=-20:1:20;
A4=0.0967;
A5=0.818;
c6=0.8175;
a=0.765;
b=0.795;
k=0.186;
a1=0.489;
x=1;
Y=linspace(0,2,100);
C=0.6324;
B1= sign(t+A5).*log(1+abs(t+A5)/10^C);

fig = figure;
ax = axes();
hold(ax);
grid on;
xlabel('\bf{t}','fontsize',26);
ylabel('\bf{u(x,y,t)}','fontsize',26);
p = plot(t,zeros(size(t)));

v = VideoWriter('test.mp4', 'MPEG-4');
open(v);
for i=1:numel(Y)
    y = Y(i);
    p.YData = (k*a1.^2+1)./(2*a)+(-x./(2*a)+(a1*y)./(2*a)+(A4*B1)./((2*a)*(t+A5))+c6)./(t+A5);
    title(sprintf('\bf{x=1,y=%.2f}',y),'fontsize',26);
    drawnow;
    
    frame = getframe(fig);
    writeVideo(v, frame)
end
close(v);

[priya anjali]

And how to save it, so that it's visible in pdf file also?

[Ameer Hamza]

You want to add images to pdf file? See export graphics() to save images at few values of 'y'.

[priya anjali]

Yes, in pdf file which i generated from latex file. Here i can put a command on .mp4 file to attach. So, is it possible to save it as a vedio ?