用matlab写一个热管的函数，输入是冷源温度T1，热源温度T2和T3，初始气压P1,最终气压P2,求最终热源温度T22和T32

热管是一种利用毛细管效应和蒸发-冷凝循环来传递热量的装置。下面是一个简单的Matlab函数，用于计算热管的温度分布，输入为冷源温度T1、热源温度T2和T3、初始气压P1和最终气压P2，输出为最终热源温度T22和T32。

function [T22, T32] = heat_pipe(T1, T2, T3, P1, P2)
% T1: 冷源温度
% T2, T3: 热源温度
% P1: 初始气压
% P2: 最终气压
% T22, T32: 最终热源温度

% 热管长度
L = 1; % 假设热管长度为1m

% 热管内径和外径
D_i = 0.002; % 热管内径为2mm
D_o = 0.01; % 热管外径为10mm

% 热管壁厚度
delta = (D_o - D_i) / 2;

% 热管截面积
A = pi * (D_o^2 - D_i^2) / 4;

% 空气的摩尔质量和气体常数
M = 0.0289644; % 空气的摩尔质量
R = 8.314; % 气体常数

% 求解热管内部温度分布
syms T(x)
eqn = diff(T,x,2) == -2*pi*delta*A*(P2-P1)*((T(x)-T1)/(P1*T(x))+(T(x)-T2)/(P2*T2));
cond = [T(0)==T1, T(L)==T3];
T(x) = dsolve(eqn, cond);

% 计算最终热源温度
T22 = T(L);
T32 = (T2-T22)*P2/P1+T22;
end

使用方法：

T1 = 300; % 冷源温度
T2 = 400; % 热源温度1
T3 = 200; % 热源温度2
P1 = 101325; % 初始气压
P2 = 80000; % 最终气压

[T22, T32] = heat_pipe(T1, T2, T3, P1, P2)

注意，这个函数只是一个简单的模型，实际情况可能更为复杂，因此需要根据具体情况进行适当的调整。