嘿，这个问题我刚好在做材料相场模拟（巧了，你说的代理问题我也接触过）的时候琢磨过，咱们拆解开来聊：

首先说你纠结的正弦/余弦基能不能处理非齐次Dirichlet边界的问题：

• 正弦基本身是天生对应齐次Dirichlet边界（T(0)=T(L)=0）的，因为sin(nπx/L)在两个边界点上的值都是0，直接用的话没法直接适配非齐次的情况；
• 余弦基则是对应齐次Neumann边界（边界处导数为0），或者说它的边界值是1和(-1)^n，也没法直接匹配非齐次Dirichlet的给定值。

但这不代表三角基就用不了，有个很实用的小技巧：变量替换消去非齐次边界。假设你的边界条件是T(0,t)=a(t)，T(L,t)=b(t)，你可以构造一个辅助函数：
T(x,t) = U(x,t) + (b(t)-a(t))x/L + a(t)
这样U(x,t)就自动满足齐次Dirichlet边界U(0,t)=U(L,t)=0了。把这个代入原热方程，你会得到一个关于U的PDE，这时候就可以用你熟悉的正弦伪谱法求解U，最后再通过上面的式子还原出T就行。这个方法不用换基，成本很低，很适合你这种用热方程做代理的情况。

再说说切比雪夫基的选项：
切比雪夫配置法（尤其是用包含边界点的Gauss-Lobatto配置点）是可以直接处理非齐次Dirichlet边界的，不需要做变量替换。因为配置点本身就包含区间的两个端点（标准区间[-1,1]的话就是x=-1和x=1），你可以直接把边界条件作为约束嵌入到离散后的方程组里。举个例子，在构造离散的拉普拉斯算子矩阵时，把对应边界点的行直接替换成边界条件对应的方程（比如x=0对应的行，只有该配置点的位置系数为1，其余为0，方程右边就是给定的边界值a(t)）。切比雪夫基的优势是能处理更一般的边界条件，而且用FFT加速导数计算的效率也很高，不过如果你已经熟悉三角基的话，上面的变量替换方法可能更顺手。

最后给你几本靠谱的参考资料，都是谱方法领域的经典：

• Spectral Methods in Fluid Dynamics（Canuto等人著）：里面详细覆盖了三角基和切比雪夫基的伪谱法实现，包括非齐次边界条件的各种处理技巧；
• Numerical Solution of Partial Differential Equations: An Introduction（Morton & Mayers著）：有专门章节讲谱方法的边界条件处理，变量替换的方法讲得通俗易懂；
• Spectral Methods: Fundamentals in Single Domains（Canuto等人著）：偏入门级的经典，伪谱配置法的细节和边界约束的实现都有清晰的推导。

备注：内容来源于stack exchange，提问作者Doug