for i = 2:length(t) % 计算温度场 T(1,:) = (2*U*d*Tout + G + q + sigma*eps*T(2,:).^4)/(2*U*d + sigma*eps*4*T(2,:).^3); T(W*10,:) = (2*U*d*Tout + G + q + sigma*eps*T((W-1)*10+1,:).^4)/(2*U*d + sigma*eps*4*T((W-1)*10+1,:).^3); T(:,1) = (2*U*d*Tout + G + q + sigma*eps*T(:,2).^4)/(2*U*d + sigma*eps*4*T(:,2).^3); T(:,H*10) = (2*U*d*Tout + G + q + sigma*eps*T(:,(H-1)*10+1).^4)/(2*U*d + sigma*eps*4*T(:,(H-1)*10+1).^3); for j = 2:W*10-1 for k = 2:H*10-1 T(j,k) = T(j,k) + (dt/(rho*cp*(1/10)^3))*((U*((Tout-T(j,k))/d) + sigma*eps*(Tout^4-T(j,k)^4) - G - q) + V*((T(j+1,k)-T(j,k))/(1/10)) - V*((T(j,k)-T(j-1,k))/(1/10)) + V*((T(j,k+1)-T(j,k))/(1/10)) - V*((T(j,k)-T(j,k-1))/(1/10))); end end end

根据代码，你尝试在第二个for循环中更新温度场，其中T的大小为50x50。然而，你在更新四个边界时，使用的是列向量T(:,1)和T(:,H*10)，这个时候T的大小是50x1，因此会导致矩阵大小不一致的错误。

为了解决这个问题，你可以将边界的赋值操作改为行向量形式，例如 T(1,:)、T(W*10,:)，这样可以确保大小一致。修改后的代码如下：

for i = 2:length(t)
    % 计算温度场
    T(1,:) = (2*U*d*Tout + G + q + sigma*eps*T(2,:).^4)/(2*U*d + sigma*eps*4*T(2,:).^3);
    T(W*10,:) = (2*U*d*Tout + G + q + sigma*eps*T((W-1)*10+1,:).^4)/(2*U*d + sigma*eps*4*T((W-1)*10+1,:).^3);
    T(:,1) = (2*U*d*Tout + G + q + sigma*eps*T(:,2).^4)/(2*U*d + sigma*eps*4*T(:,2).^3);
    T(:,H*10) = (2*U*d*Tout + G + q + sigma*eps*T(:,(H-1)*10+1).^4)/(2*U*d + sigma*eps*4*T(:,(H-1)*10+1).^3);
    for j = 2:W*10-1
        for k = 2:H*10-1
            T(j,k) = T(j,k) + (dt/(rho*cp*(1/10)^3))*((U*((Tout-T(j,k))/d) + sigma*eps*(Tout^4-T(j,k)^4) - G - q) + V*((T(j+1,k)-T(j,k))/(1/10)) - V*((T(j,k)-T(j-1,k))/(1/10)) + V*((T(j,k+1)-T(j,k))/(1/10)) - V*((T(j,k)-T(j,k-1))/(1/10)));
        end
    end
end

希望对你有所帮助！