Una delgada placa rectangular de acero
de 9cm x 12 cm se encuentra sujeta a las
siguientes condiciones de frontera:
T (x,0)=
42C
T (0,y)=
89C
T(X,12)=
35C
T (9,y)=
180C
Calcular el perfil de Temperatura Interior si
x=y=0.01m
ECUACIONES
4T1=T12+T2+180+42
4T2=T13+T3+180+T1
4T3=T14+T4+180+T2
4T4=T15+T5+180+T3
4T5=T16+T6+180+T4
4T6=T17+T7+180+T5
4T7=T18+T8+180+T6
4T8=T19+T9+180+T7
4T9=T20+T10+180+T8
4T10=T21+T11+180+T9
4T11=T22+35+180+T10
4T12=T23+T13+T1+42
4T13=T24+T14+T2+T12
4T14=T25+T15+T3+T13
4T15=T26+T16+T4+T14
4T16=T27+T17+T5+T15
4T17=T28+T18+T6+T16
4T18=T29+T19+T7+T17
4T19=T30+T20+T8+T18
4T20=T31+T21+T9+T19
4T21=T32+T22+T10+T20
4T22=T33+35+T11+T21
4T23=T34+T24+T12+42
4T24=T35+T25+T13+T23
4T25=T36+T26+T14+T24
4T26=T37+T27+T15+T25
4T27=T38+T28+T16+T26
4T28=T39+T29+T17+T27
4T29=T40+T30+T18+T28
4T30=T41+T31+T19+T29
4T31=T42+T32+T20+T30
4T32=T43+T33+T21+T31
4T33=T44+35+T22+T32
4T34=T45+T35+T23+42
4T35=T46+T36+T24+T34
4T36=T47+T37+T25+T35
4T37=T48+T38+T26+T36
4T38=T49+T39+T27+T37
4T39=T50+T40+T28+T38
4T40=T51+T41+T29+T39
4T41=T52+T42+T30+T40
4T42=T53+T43+T31+T41
4T43=T54+T44+T32+T42
4T44=T55+35+T33+T43
4T45=T56+T46+T34+42
4T46=T57+T47+T35+T45
DESPEJANDO
4T1-T12-T2=162
4T2-T13-T3-T1=180
4T3-T14-T4-T2=180
4T4-T15-T5-T3=180
4T5-T16-T6-T4=180
4T6-T17-T7-75=180
4T7-T18-T8-T6=180
4T8-T19-T9-T7=180
4T9-T20-T10-T8=180
4T10-T21-T11-T9=180
4T11-T22-T10=215
4T12-T23-T13-T11=42
4T13-T24-T14-T2-T12=0
4T14-T25-T15-T3-T13=0
4T15-T26-T16-T4-T14=0
4T16-T27-T17-T5-T15=0
4T17-T28-T18-T6-T16=0
4T18-T29-T19-T7-T17=0
4T19-T30-T20-T8-T18=0
4T20-T31-T21-T9-T19=0
4T21-T32-T22-T10-T20=0
4T22-T33-T11-T21=35
4T23-T34-T24-T12=42
4T24-T35-T25-T13-T23=0
4T25-T36-T26-T14-T24=0
4T26-T37-T27-T15-T25=0
4T27-T38-T28-T16-T26=0
4T28-T39-T29-T17-T27=0
4T29-T40-T30-T18-T28=0
4T30-T41-T31-T19-T29=0
4T31-T42-T32-T20-T30=0
4T32-T43-T33-T21-T31=0
4T33-T44-T22-T32=35
4T34-T45-T35-T23=4235
4T35-T46-T36-T24-T34=0
4T36-T47-T37-T25-T35=0
4T37-T48-T38-T26-T36=0
4T38-T49-T39-T27-T37=0
4T39-T50-T40-T28-T38=0
4T40-T51-T41-T29-T39=0
4T41-T52-T42-T30-T40=0
4T42-T53-T43-T31-T41=0
4T43-T54-T44-T32-T42=0
4T44-T55-T33-T43=35
4T45-T56-T46-T34=42
4T46-T57-T47-T35-T45=0
4T47=T58+T48+T36+T46
4T48=T59+T49+T37+T47
4T49=T60+T50+T38+T48
4T50=T61+T51+T39+T49
4T51=T62+T52+T40+T50
4T52=T63+T53+T41+T51
4T53=T64+T54+T42+T52
4T54=T65+T55+T43+T53
4T55=T66+35+T44+T54
4T56=T67+T57+T45+42
4T57=T68+T58+T46+T56
4T58=T69+T59+T47+T57
4T59=T70+T60+T48+T58
4T60=T71+T61+T49+T59
4T61=T72+T62+T50+T60
4T62=T73+T63+T51+T61
4T63=T74+T64+T52+T62
4T64=T75+T65+T53+T63
4T65=T76+T66+T54+T64
4T66=T77+35+T55+T65
4T67=T78+T68+T56+42
4T68=T79+T69+T57+T67
4T69=T80+T70+T58+T68
4T70=T81+T71+T59+T69
4T71=T82+T72+T60+T70
4T72=T83+T73+T61+T71
4T73=T84+T74+T62+T72
4T74=T85+T75+T63+T73
4T75=T86+T76+T64+T74
4T76=T87+T77+T65+T75
4T77=T88+35+T67+T76
4T78=89+T79+T67+42
4T79=89+T80+T68+T178
4T80=89+T81T69+T179
4T81=89+T82+T70+T80
4T82=89+T83+T71+T81
4T83=89+T84+T72+T82
4T84=89+T85+T73+T83
4T85=89+T86+T74+T84
4T86=89+T87+T75+T85
4T87=89+T88+T76+T86
4T88=89+35+T77+T87
4T47-T58-T48-T36-T46=0
4T48-T59-T49-T37-T47=0
4T49-T60-T50-T38-T48=0
4T50-T61-T51-T39-T49=0
4T51-T62-T52-T40-T50=0
4T52-T63-T53-T41-T51=0
4T53-T64-T54-T42-T52=0
4T54-T65-T55-T43-T53=0
4T55-T66-T44-T54=35
4T56-T67-T57-T45=42
4T57-T68-T58-T46-T56=0
4T58-T69-T59-T47-T57=0
5T59-T70-T60-T48-T58=0
4T60-T71-T61-T49-T59=0
4T61-T72-T62-T50-T60=0
4T62-T73-T63-T51-T61=0
4T63-T74-T64-T52-T62=0
4T64-T75-T65-T53-T63=0
4T65-T76-T66-T54-T64=0
4T66-T77-T55-T65=35
4T67-T78-T68-T56=42
4T68-T79-T69-T57-T67=0
4T69-T80-T70-T58-T68=0
4T70-T81-T71-T59-T69=0
4T71-T82-T72-T60-T70=0
4T72-T83-T73-T61-T71=0
4T73-T84-T74-T62-T72=0
4T74-T85-T75-T63-T73=0
4T75-T86-T76-T64-T74=0
4T76-T87-T77-T65-T75=0
4T77-T88-T66-T76=35
4T78-T79-T67=131
4T79-T80-T68-T78=89
4T80-T81-T69-T79=89
4T81-T82-T70-T80=89
4T82-T83-T71-T81=89
4T83-T84-T72-T82=89
4T84-T85-T73-T83=89
4T85-T86-T74-T84=89
4T86-T87-T75-T85=89
4T87-T88-T76-T86=89
4T88-T77-T87=124
MATLAB
>> A=[1]
A=
1
>> B=[1]
B=
1
>> T=[A\B]
T=
84.8970
125.6662
143.6718
152.2823
156.3733
157.7956
157.1393
154.2307
147.9728
135.2519
107.1080
51.9220
94.0958
116.7388
129.0839
135.4154
137.6699
136.5309
131.8106
122.4088
105.9266
78.1803
58.5799
82.0563
100.1039
111.8989
118.5347
120.9375
119.5040
114.0719
103.9253
87.8653
64.6866
58.3413
75.4457
89.7213
99.8734
105.8868
108.0414
106.4756
101.0478
91.3552
76.9228
57.7006
57.3394
71.6640
83.4623
91.9865
97.0976
98.8657
97.3093
92.2885
83.5249
70.7699
54.1932
57.3525
70.4084
80.4775
87.5128
91.6514
93.0145
91.6073
87.2720
79.6861
68.4388
53.3024
59.6621
72.1398
80.5264
85.9357
88.9808
89.9337
88.8335
85.5059
79.5086
69.9970
55.5774
67.1561
77.9622
83.5528
86.7227
88.4023
88.9058
88.2872
86.4096
82.8453
76.4631
64.0101
>> G=[1]
G=
>> XP=0: 0.01:0.09
XP =
Columns 1 through 8
0 0.0100 0.0200 0.0300 0.0400
0.0500 0.0600 0.0700
Columns 9 through 10
0.0800 0.0900
YP =
Columns 1 through 8
0 0.0100 0.0200 0.0300 0.0400
0.0500 0.0600 0.0700
Columns 9 through 13
0.0800 0.0900 0.1000 0.1100 0.1200
>> surf(XP,YP,G)
>> hold on
>> Title('PERFIL DE TEMPERATURA
INTERIOR')
>> xlabel('X,m')
>> ylabel('Y,m')
>> zlabel('temperatura,C')
>> axis{[0.0,0.09,0.12,25,100])
shading interp
