MATLAB: How to add colours on mesh

3d plotsmeshplot

A = 0.06;
k_l = 26400; %Linear stiffness

m = 483; %Mass

f = @(t,x,Om,l,k_s,d) [ x(2); ...
            -(2*k_s*(x(1)-(A*sin(Om*t))))* ...
            (sqrt((l-d)^2 + (x(1)-(A*sin(Om*t)))^2) - l)/ ...
            (m*(sqrt((l-d)^2 + (x(1)-(A*sin(Om*t)))^2))) ];
%%

Om_array = linspace(0,20,11); %in rad/s-1

l_array = linspace(0.2,1,11);
[om_array, L_array] = meshgrid(Om_array, l_array);
d = linspace(-0.005, -0.03, 5);
%d = -1;
Response_amp = zeros([size(Om_array), numel(d)]);
T = 150;
x0 = [0,0];
for k=1:numel(d)
    for i=1:numel(Om_array)
        for j=1:numel(l_array)
            Om = om_array(i,j);
            l = L_array(i,j);
            k_s = -(k_l*(l-d))/(4*d); %Spring stiffness

            
            [t, x] = ode45(@(t,x) f(t,x,Om,l,k_s,d(k)),[100,T],x0);
            Response_amp(i,j,k) = (max(x(:,1)) - min(x(:,1)))/2;
        end
    end
end
%% plot

figure(1);
ax = axes();
view(3);
hold(ax);
view([30 33]);
grid on
for i=1:size(Response_amp,3)
    mesh(om_array/(2*pi),L_array,Response_amp(:,:,i));
end
hold on
A = 0.06;
k_l = 26400; %Linear stiffness
m = 483; %Mass
f = @(t,x,Om,l,k_s,d) [ x(2); ...
            -(2*k_s*(x(1)-(A*sin(Om*t))))* ...
            (sqrt((l-d)^2 + (x(1)-(A*sin(Om*t)))^2) - l)/ ...
            (m*(sqrt((l-d)^2 + (x(1)-(A*sin(Om*t)))^2))) ];
%%
Om_array = linspace(0,20,21); %in rad/s-1
l_array = linspace(0.2,1,21);
[om_array, L_array] = meshgrid(Om_array, l_array);
d = -1;
Response_amp = zeros([size(Om_array), numel(d)]);
T = 150;
x0 = [0,0];
for k=1:numel(d)
    for i=1:numel(Om_array)
        for j=1:numel(l_array)
            Om = om_array(i,j);
            l = L_array(i,j);
            k_s = -(k_l*(l-d))/(4*d); %Spring stiffness
            
            [t, x] = ode45(@(t,x) f(t,x,Om,l,k_s,d(k)),[100,T],x0);
            Response_amp(i,j,k) = (max(x(:,1)) - min(x(:,1)))/2;
        end
    end
end
%% plot
% figure(1);
% ax = axes();
% view(3);
% hold(ax);
% view([30 33]);
% grid on
for i=1:size(Response_amp,3)
    mesh(om_array/(2*pi),L_array,Response_amp(:,:,i), 'EdgeColor', 'black');
end
xlabel('Frequency (Hz)')
ylabel('Length of the spring (m)')
zlabel('Response Amplitude (m)')
set(gca,'FontSize',15)
hold off
Hi, I managed to plot a black mesh function when d = -1 in this code. However, when d is set to be linspace(-0.005, -0.03, 5), how do I add colour for indiviual mesh (when d = -0.005, d = -0.01, … d = -0.03)?

Best Answer

Currently i changed EdgeColor for every iteration by setting it to rand rgb values. you can create an array equal to length of size(Response_amp_3) x 3 and index it in for loop
%% plot 1
figure(1);
ax = axes();
view(3);
hold(ax);
view([30 33]);
grid on
for i=1:size(Response_amp,3)
    mesh(om_array/(2*pi),L_array,Response_amp(:,:,i),'EdgeColor',rand(1,3));
end
hold on
%% plot2
% %figure(1);
% ax = axes();
view(3);
% hold(ax);
view([30 33]);
grid on
for i=1:size(Response_amp,3)
    mesh(om_array/(2*pi),L_array,Response_amp(:,:,i),'EdgeColor',rand(1,3));
end
xlabel('Frequency (Hz)')
ylabel('Length of the spring (m)')
zlabel('Response Amplitude (m)')
set(gca,'FontSize',15)
hold off
Related Question