Multiple facecolor in one fsurf plot
1 Ansicht (letzte 30 Tage)
Ältere Kommentare anzeigen
MINATI PATRA
am 23 Jul. 2023
Kommentiert: MINATI PATRA
am 25 Jul. 2023
M = 1; m = 0.1; n = 3; Pr = 3; Ec = 0.1; N = 0.5; lp = 0.01; Bi = 0.5; wf = 100; om = 0.1; %w1 = 10; w2 = 10;
syms w1 w2
for n = [3 3.7 5 5.7 8.9]
Cpf = 3594;rhof = 1063;kf = 0.492;Cps1 = 960;sigmaf = 4.3*(10^-5);rhos1 = 2810;ks1 = 173;sigmas1 = 26.77*10^6;Cps2 = 893;rhos2 = 2720;ks2 = 222;sigmas2 = 34.83*10^6;
phi1 = (w1/rhos1)./( (w1/rhos1) + (w2/rhos2) + (wf/rhof));
phi2 = (w2/rhos2)./( (w1/rhos1) + (w2/rhos2) + (wf/rhof)) ;
phi = phi1 + phi2;rhos = ((rhos1*w1) + (rhos2*w2))./(w1 + w2);Cps = ((Cps1 * w1) + (Cps2 * w2))./(w1 + w2);
A1 = (1-phi)^-2.5;A2 = ( 1- phi + phi*(rhos/rhof));A4 = ( 1- phi+phi*((rhos*Cps)/(rhof*Cpf)));
C1 = ( (ks1+(n-1)*kf-(n-1)*phi1*(kf-ks1))/(ks1+(n-1)*kf+phi1*(kf-ks1)) )*kf;
A3 = (ks2+(n-1)*C1-(n-1)*phi2*(C1-ks2))/(ks2+(n-1)*C1+phi2*(C1-ks2));
C2 = ( (sigmas1+2*sigmaf-2*phi1*(sigmaf-sigmas1))/(sigmas1+2*sigmaf+phi1*(sigmaf-sigmas1)))*sigmaf;
A5 = (sigmas2+2*C2-2*phi2*(C2-sigmas2))/(sigmas2+2*C2+phi2*(C2-sigmas2));
figure(1),fsurf(A3,[0 10 0 10],'LineWidth',2,'EdgeColor','none'), hold on,%view(30,30);
xlabel ('\bfw_1','color','blue','Fontsize',14,'FontWeight', 'bold'); ylabel ('\bfw_2','color','blue','Fontsize',10,'FontWeight', 'bold');
zlabel ('\bfk_{hnf}\k_f','color','blue','Fontsize',14,'FontWeight','bold'); axis on; box on; ax = gca; ax.BoxStyle = 'full'; ax.XAxis.LineWidth = 2; ax.YAxis.LineWidth = 2; ax.ZAxis.LineWidth = 2;
legend({'\color{green}\bfSpherical','\color{red}\bfBrick','\color{blue}\bfCylindrical','\color{black}\bfPlatelet','\color{magenta}\bfBlade'},'Box','off');
end
%%% %% I have TWO questions: (i) I want the values of w1 = w2 = 10 (fsurf(A3,[0 10 0 10]), IS THIS RIGHT?) to run the fsurf plot
%% (ii) I need 5 different facecolor in one plot as given in legend
0 Kommentare
Akzeptierte Antwort
Walter Roberson
am 23 Jul. 2023
Bearbeitet: Walter Roberson
am 23 Jul. 2023
M = 1; m = 0.1; n = 3; Pr = 3; Ec = 0.1; N = 0.5; lp = 0.01; Bi = 0.5; wf = 100; om = 0.1; %w1 = 10; w2 = 10;
syms w1 w2
facecolors = {'g', 'r', 'b', 'k', 'm'};
nvalues = [3 3.7 5 5.7 8.9];
for nidx = 1:length(nvalues)
n = nvalues(nidx);
Cpf = 3594;rhof = 1063;kf = 0.492;Cps1 = 960;sigmaf = 4.3*(10^-5);rhos1 = 2810;ks1 = 173;sigmas1 = 26.77*10^6;Cps2 = 893;rhos2 = 2720;ks2 = 222;sigmas2 = 34.83*10^6;
phi1 = (w1/rhos1)./( (w1/rhos1) + (w2/rhos2) + (wf/rhof));
phi2 = (w2/rhos2)./( (w1/rhos1) + (w2/rhos2) + (wf/rhof)) ;
phi = phi1 + phi2;rhos = ((rhos1*w1) + (rhos2*w2))./(w1 + w2);Cps = ((Cps1 * w1) + (Cps2 * w2))./(w1 + w2);
A1 = (1-phi)^-2.5;A2 = ( 1- phi + phi*(rhos/rhof));A4 = ( 1- phi+phi*((rhos*Cps)/(rhof*Cpf)));
C1 = ( (ks1+(n-1)*kf-(n-1)*phi1*(kf-ks1))/(ks1+(n-1)*kf+phi1*(kf-ks1)) )*kf;
A3 = (ks2+(n-1)*C1-(n-1)*phi2*(C1-ks2))/(ks2+(n-1)*C1+phi2*(C1-ks2));
C2 = ( (sigmas1+2*sigmaf-2*phi1*(sigmaf-sigmas1))/(sigmas1+2*sigmaf+phi1*(sigmaf-sigmas1)))*sigmaf;
A5 = (sigmas2+2*C2-2*phi2*(C2-sigmas2))/(sigmas2+2*C2+phi2*(C2-sigmas2));
fsurf(A3,[0 10 0 10],'LineWidth',2,'EdgeColor','none', 'facecolor', facecolors{nidx} );
hold on
end
view(30,30);
xlabel ('$\textbf{w}_1$','color','blue','Fontsize',14,'FontWeight', 'bold', 'interpreter', 'latex');
ylabel ('$\textbf{w}_2$','color','blue','Fontsize',10,'FontWeight', 'bold', 'interpreter', 'latex');
zlabel ('$\textbf{k}_{hnf}\textit{f}$','color','blue','Fontsize',14,'FontWeight','bold', 'interpreter', 'latex');
axis on; box on;
ax = gca;
ax.BoxStyle = 'full';
ax.XAxis.LineWidth = 2;
ax.YAxis.LineWidth = 2;
ax.ZAxis.LineWidth = 2;
legend({'\color{green}\bfSpherical','\color{red}\bfBrick','\color{blue}\bfCylindrical','\color{black}\bfPlatelet','\color{magenta}\bfBlade'},'Box','off');
I had to make some guesses about what you were trying to do with the labels, I could not figure out the \k_f -- I had to guess italics.
Weitere Antworten (0)
Siehe auch
Kategorien
Mehr zu Axis Labels finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!