Matlab - cómo desplazar ejes de coordenadas

 
Vista:
sin imagen de perfil
Val: 2
Ha aumentado su posición en 69 puestos en Matlab (en relación al último mes)
Gráfica de Matlab

cómo desplazar ejes de coordenadas

Publicado por carmen (3 intervenciones) el 23/04/2020 18:37:42
Hola, os dejo por aquí mi duda a ver si alguien pudiera ayudarme:
Estoy haciendo un código de matlab con el cual debo representar una variable, la temperatura en funcion del ángulo theta, de manera que quede así:

1

sin embargo, el código que he hecho plotea lo siguiente:

2

Esto, posiblemente se deba a que los sistemas de coordenadas están movidos. LO QUE TENGO QUE HACER ES DESPLAZARLOS Y HACER QUE EL MÁXIMO DE TEMPERATURAS COINCIDA CON THETA=0º (puedo hacerlo jugando con la posicion donde se localice la temperatura máxima).
¿alguien sabe como se puede programar eso con matlab o alguna idea para hacerlo?
Valora esta pregunta
Me gusta: Está pregunta es útil y esta claraNo me gusta: Está pregunta no esta clara o no es útil
0
Responder
Imágen de perfil de JOSE JEREMIAS CABALLERO
Val: 8.212
Oro
Ha mantenido su posición en Matlab (en relación al último mes)
Gráfica de Matlab

cómo desplazar ejes de coordenadas

Publicado por JOSE JEREMIAS CABALLERO (5683 intervenciones) el 23/04/2020 21:01:00
Podría mostrar su código.
Valora esta respuesta
Me gusta: Está respuesta es útil y esta claraNo me gusta: Está respuesta no esta clara o no es útil
0
Comentar

cómo desplazar ejes de coordenadas

Publicado por CARMEN (3 intervenciones) el 23/04/2020 23:17:00
Si claro, Aqui le dejo el código:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
clear all
close all
 
% LEER DATOS
resultados= importdata('solucion.txt');
 
% VARIABLES ALMACENADAS EN VECTORES
NN=resultados.data(:,1);
x=resultados.data(:,2);
y=resultados.data(:,3);
z=resultados.data(:,4);
presion=resultados.data(:,5);
densidad=resultados.data(:,6);
velocidad=resultados.data(:,7);
x_vel=resultados.data(:,8);
y_vel=resultados.data(:,9);
z_vel=resultados.data(:,10);
temperatura=resultados.data(:,11);
y_star=resultados.data(:,12);
y_plus=resultados.data(:,13);
flujo_calor=resultados.data(:,14);
cell_volume=resultados.data(:,15);
 
% SISTEMA DE COORDENADAS CILÍNDRICO
[theta,rho,z]=cart2pol(x,y,z);
 
 
vector_operado=abs(unique(z)-5);
z_unicos=unique(z);
 
SC=z_unicos(find(min(vector_operado)==vector_operado));
 
a_z=find(z==SC);
 
 
 
 
x_corte=x(a_z);
y_corte=y(a_z);
z_corte=z(a_z);
presion_corte=presion(a_z);
densidad_corte=densidad(a_z);
velocidad_corte=velocidad(a_z);
x_vel_corte=x_vel(a_z);
y_vel_corte=y_vel(a_z);
z_vel_corte=z_vel(a_z);
temperatura_corte=temperatura(a_z);
y_star_corte=y_star(a_z);
y_plus_corte=y_plus(a_z);
flujo_calor_corte=flujo_calor(a_z);
cell_volume_corte=cell_volume(a_z);
theta_corte=theta(a_z);
rho_corte=rho(a_z);
 
% DIFERENCIAR ZONAS --> hay tubo, aire y sal
% DENSIDAD:
%           sal--> densidad>1000
%           tubo--> densidad>0.01 (es < que 0.01 porque aunque la densidad
%                   del haynes230 es 8970, la exporta como 0!!)
%           aire--> 0.01<densidad<1000
 
% SAL:
 
x_sal=x_corte(find(densidad_corte>1000));
y_sal=y_corte(find(densidad_corte>1000));
z_sal=z_corte(find(densidad_corte>1000));
presion_sal=presion_corte(find(densidad_corte>1000));
densidad_sal=densidad_corte(find(densidad_corte>1000));
velocidad_sal=velocidad_corte(find(densidad_corte>1000));
x_vel_sal=x_vel_corte(find(densidad_corte>1000));
y_vel_sal=y_vel_corte(find(densidad_corte>1000));
z_vel_sal=z_vel_corte(find(densidad_corte>1000));
temperatura_sal=temperatura_corte(find(densidad_corte>1000));
y_star_sal=y_star_corte(find(densidad_corte>1000));
y_plus_sal=y_plus_corte(find(densidad_corte>1000));
flujo_calor_sal=flujo_calor_corte(find(densidad_corte>1000));
cell_volume_sal=cell_volume_corte(find(densidad_corte>1000));
theta_sal=theta_corte(find(densidad_corte>1000));
rho_sal=rho_corte(find(densidad_corte>1000));
 
 
% TUBO:
 
x_tubo=x_corte(find(densidad_corte<0.01));
y_tubo=y_corte(find(densidad_corte<0.01));
z_tubo=z_corte(find(densidad_corte<0.01));
presion_tubo=presion_corte(find(densidad_corte<0.01));
densidad_tubo=densidad_corte(find(densidad_corte<0.01));
velocidad_tubo=velocidad_corte(find(densidad_corte<0.01));
x_vel_tubo=x_vel_corte(find(densidad_corte<0.01));
y_vel_tubo=y_vel_corte(find(densidad_corte<0.01));
z_vel_tubo=z_vel_corte(find(densidad_corte<0.01));
temperatura_tubo=temperatura_corte(find(densidad_corte<0.01));
y_star_tubo=y_star_corte(find(densidad_corte<0.01));
y_plus_tubo=y_plus_corte(find(densidad_corte<0.01));
flujo_calor_tubo=flujo_calor_corte(find(densidad_corte<0.01));
cell_volume_tubo=cell_volume_corte(find(densidad_corte<0.01));
theta_tubo=theta_corte(find(densidad_corte<0.01));
rho_tubo=rho_corte(find(densidad_corte<0.01));
 
 
%AIRE:
 
x_aire=x_corte(find(densidad_corte>0.01 & densidad_corte<1000));
y_aire=y_corte(find(densidad_corte>0.01 & densidad_corte<1000));
z_aire=z_corte(find(densidad_corte>0.01 & densidad_corte<1000));
presion_aire=presion_corte(find(densidad_corte>0.01 & densidad_corte<1000));
densidad_aire=densidad_corte(find(densidad_corte>0.01 & densidad_corte<1000));
velocidad_aire=velocidad_corte(find(densidad_corte>0.01 & densidad_corte<1000));
x_vel_aire=x_vel_corte(find(densidad_corte>0.01 & densidad_corte<1000));
y_vel_aire=y_vel_corte(find(densidad_corte>0.01 & densidad_corte<1000));
z_vel_aire=z_vel_corte(find(densidad_corte>0.01 & densidad_corte<1000));
temperatura_aire=temperatura_corte(find(densidad_corte>0.01 & densidad_corte<1000));
y_star_aire=y_star_corte(find(densidad_corte>0.01 & densidad_corte<1000));
y_plus_aire=y_plus_corte(find(densidad_corte>0.01 & densidad_corte<1000));
flujo_calor_aire=flujo_calor_corte(find(densidad_corte>0.01 & densidad_corte<1000));
cell_volume_aire=cell_volume_corte(find(densidad_corte>0.01 & densidad_corte<1000));
theta_aire=theta_corte(find(densidad_corte>0.01 & densidad_corte<1000));
rho_aire=rho_corte(find(densidad_corte>0.01 & densidad_corte<1000));
 
 
% LIMITAR SUPERFICIES
% interior: rho_tubo=<0.01485
% exterior: rho_tubo=>0.0165
 
 
 
% Obtengo los puntos en los bordes del tubo: TE (tubo exterior)
x_TE=x_tubo(find(flujo_calor_tubo~=0 & y_plus_tubo==0 & rho_tubo>=0.0165));
y_TE=y_tubo(find(flujo_calor_tubo~=0 & y_plus_tubo==0 & rho_tubo>=0.0165));
z_TE=z_tubo(find(flujo_calor_tubo~=0 & y_plus_tubo==0 & rho_tubo>=0.0165));
temperatura_TE=temperatura_tubo(find(flujo_calor_tubo~=0 & y_plus_tubo==0 & rho_tubo>=0.0165));
y_plus_TE=y_plus_tubo(find(flujo_calor_tubo~=0 & y_plus_tubo==0 & rho_tubo>=0.0165));
flujo_calor_TE=flujo_calor_tubo(find(flujo_calor_tubo~=0 & y_plus_tubo==0 & rho_tubo>=0.0165));
cell_volume_TE=cell_volume_tubo(find(flujo_calor_tubo~=0 & y_plus_tubo==0 & rho_tubo>=0.0165));
theta_TE=theta_tubo(find(flujo_calor_tubo~=0 & y_plus_tubo==0 & rho_tubo>=0.0165));
rho_TE=rho_tubo(find(flujo_calor_tubo~=0 & y_plus_tubo==0 & rho_tubo>=0.0165));
 
 
 
% Obtengo los puntos en los bordes del tubo: TI (tubo interior)
x_TI=x_tubo(find(flujo_calor_tubo~=0 & y_plus_tubo==0 & rho_tubo<=0.01485));
y_TI=y_tubo(find(flujo_calor_tubo~=0 & y_plus_tubo==0 & rho_tubo<=0.01485));
z_TI=z_tubo(find(flujo_calor_tubo~=0 & y_plus_tubo==0 & rho_tubo<=0.01485));
temperatura_TI=temperatura_tubo(find(flujo_calor_tubo~=0 & y_plus_tubo==0 & rho_tubo<=0.01485));
y_plus_TI=y_plus_tubo(find(flujo_calor_tubo~=0 & y_plus_tubo==0 & rho_tubo<=0.01485));
flujo_calor_TI=flujo_calor_tubo(find(flujo_calor_tubo~=0 & y_plus_tubo==0 & rho_tubo<=0.01485));
cell_volume_TI=cell_volume_tubo(find(flujo_calor_tubo~=0 & y_plus_tubo==0 & rho_tubo<=0.01485));
theta_TI=theta_tubo(find(flujo_calor_tubo~=0 & y_plus_tubo==0 & rho_tubo<=0.01485));
rho_TI=rho_tubo(find(flujo_calor_tubo~=0 & y_plus_tubo==0 & rho_tubo<=0.01485));
 
plot(temperatura_TE,theta_TE,'s');
 
posicion_max_temperatura_TE=find(temperatura_TE==max(max(temperatura_TE)));
 
posicion_0_theta_TE=find(theta_TE==min(abs(theta_TE)));
Valora esta respuesta
Me gusta: Está respuesta es útil y esta claraNo me gusta: Está respuesta no esta clara o no es útil
0
Comentar