用MATLAB解偏微分方程和系统辨识

imm

QUESTION:用MATLAB解热传导方程：T'（t）-k（T"（x)+T"(y)+T"(z))=0,   其中k=导热率/（密度乘以等体比热）
           初始条件T(0)=300K,边界条件280k，300k，边界为长方体，一角为坐标原点，在第一卦限中，长x=0.4米，
           宽y=0.2m，高z=0.2m，长方体材料为水，容器为陶瓷，热源在容器上，功率0—100W，无盖，可解否？
          另外需实验数据，输入为热源功率，输出为若干点水温，用MATLAB最小二乘法,随即逼近法，极大似然法 (ML)
     预报误差法(PE)仿真，阶次用汉克尔矩阵法，残差平方和法估计，得出差分方程。
           此时我的课题，有高手会解的话，麻烦告诉我，当感激不尽。