计算二阶线性偏微分方程的变换系数.DOC

function A=tran %本程序用于计算二阶线性偏微分方程的变换系数 syms a b c x y f1 f2 disp(函数变换f1(x,y)的具体形式：) f1=input(f1=); %输入变换(=((x,y) f2=input(f2=); %输入变换(=((x,y) a11=input(a11=); a12=input(a12=); a22=input(a22=); b1=input(b1=); b2=input(b2=); % 以下计算一阶、二阶导数 d1x=diff(f1,x);d1xx=diff(f1,x,2); d1y=diff(f1,y);d1yy=diff(f1,y,2); d2x=diff(f2,x);d2xx=diff(f2,x,2); d2y=diff(f2,y);d2yy=diff(f2,y,2); d1xy=diff(diff(f1,x),y); d2xy=diff(diff(f2,x),y); % 以下计算变换后的系数 A11=a11*d1x^2+2*a12*d1x*d1y+a22*d1y^2; A22=a11*d2x^2+2*a12*d2x*d2y+a22*d2y^2; A12=a11*d1x*d2x+a12*(d1x*d2y+d1y*d2x)+a22*d1y*d2y; B1=a11*d1xx+2*a12*d1xy+a22*d1yy+b1*d1x+b2*d1y; B2=a11*d2xx+2*a12*d2xy+a22*d2yy+b1*d2x+b2*d2y; % A= simplify([A11,A12,A22,B1,B2]); 二次线性偏微分方程化简的例题 由特征方程的解可以定出变换，但抛物型需另找与 (= ((x, y) 无关的(=((x, y) 代换后的方程系数和原来系数的关系 例1 化简 因为 方程是抛物型的。特征方程 即只有一族实的特征线 y ( x = 常数 变换设为 ( = y ( x，( = x (只要( 不满足就行) A=tran 函数变换f1(x,y)的具体形式： f1=y-x f2=x a11=a a12=a a22=a b1=b b2=c A = [ 0, 0, a, -b+c, b] 化简后的方程形式为