Can't plot contour graph on a different plane
2 Ansichten (letzte 30 Tage)
Ältere Kommentare anzeigen
I plotted a contour graph on the x-y plane using this code with no problem:
figure(1)
iz = 1;
contourf(x_mesh(:,:,iz),y_mesh(:,:,iz),squeeze(g(:,:,iz)))
But once I tried to do it for the x-z plane, using this:
figure(2)
iy =1;
contourf(x_mesh(:,iy,:),y_mesh(:,iy,:),squeeze(g(:,iy,:)))
I got the following error:
"Error using contourf
Input arguments must have at most 2 dimensions."
Can anyone tell me why this is and how I can fix it? Here's the full code if it helps:
clear all
close all
clc
%Constant
rho = 4420; %kg/m^3
Cp = 550; %J/kg?K
T0 = 303.15; %K
A = 0.5; %[Absorbtivity]
k = 7.2; %W/m/K
alpha = 2.96*10^-6; %m^2/s
D = alpha;
P = 100; %W
v = 1; %m/s
u = v;
Tm = 1933; %K
d_laser = 0.0001; %m
r_laser = d_laser/2; %m
a = r_laser;
p = D/(u*a);
%Define
n = 100;
x = linspace(-0.00025,0.00125,n);
y = linspace(-0.00025,0.00025,n);
z = linspace(-0.00008,0.0005,n);
%Normalized
x_nor = x/a;
y_nor = y/a;
z_nor = z/(D*a/u).^0.5
[x_mesh,y_mesh,z_mesh] = ndgrid(x_nor,y_nor,z_nor);
fun = @(t) exp((-z_mesh.^2./(4*t))-((y_mesh.^2+(x_mesh-t).^2)./(4*p.*t+1)))./((4.*p.*t+1).*sqrt(t));
g = integral(fun,0,Inf,'ArrayValued',true);
figure(1)
hold on
iz = 1;
[C,h]= contourf(x_mesh(:,:,iz),y_mesh(:,:,iz),squeeze(g(:,:,iz)))
[C1,h1] = contour(x_mesh(:,:,iz),y_mesh(:,:,iz),squeeze(g(:,:,iz)),LevelList=0.16713,Color='red')
hold off
clabel(C,h,'manual','Color','white')
clabel(C1,h1,'manual','Color','white')
xticks(-5:5:25)
yticks(-5:2.5:5)
axis([-5 25 -5 5])
title('x\primey\prime plane')
xlabel('x\prime')
ylabel('y\prime')
colorbar
figure(2)
iy =1;
contourf(x_mesh(:,iy,:),y_mesh(:,iy,:),squeeze(g(:,iy,:)))
axis([-5 25 -10 0])
title('x\primez\prime plane')
xlabel('x\prime')
ylabel('z\prime')
colorbar
0 Kommentare
Akzeptierte Antwort
Weitere Antworten (0)
Siehe auch
Kategorien
Mehr zu Contour Plots finden Sie in Help Center und File Exchange
Produkte
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
