matlab求解共线方程公式以及生成等高线

matlab课后作业2 共线方程两种格式求解（只是公式的推到） 生成等高线在TXT中，另附等高线图（读取TXT信息后绘制）

syms fa ou ka;
R1=[cos（fa） 0 -sin（fa）;0 1 0; sin（fa） 0 cos（fa）];
R2=[1 0 0;0 cos（ou） -sin（ou）;0 sin（ou） cos（ou）];
R3=[cos（ka） -sin（ka） 0;sin（ka） cos（ka） 0;0 0 1];
R=R1*R2*R3;
syms Xs Ys Zs Xa Ya Za
syms x y f
X = [Xa-Xs];
Y = [Ya-Ys];
Z = [Za-Zs];
a1 = R（11）;
a2 = R（12）;
a3 = R（13）;
b1 = R（21）;
b2 = R（22）;
b3 = R（23）;
c1 = R（31）;
c2 = R（32）;
c3 = R（33）;

x = -f*（a1*X+b1*Y+c1*Z）/（a3*X+b3*Y+c3*Z）;
y = -f*（a2*X+b2*Y+c2*Z）/（a3*X+b3*Y+c3*Z）;

dxdXs = diff（x‘Xs‘）;
dydXs = diff（y‘Xs‘）;
dxdfa = diff（x‘fa‘）;
dydfa = diff（y‘fa‘）;
dxdXs
dydXs
dxdfa
dydfa

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件     683122  2013-05-06 21:06  等高线.JPG

     文件        592  2013-05-06 23:18  共线方程.m

     文件        461  2013-05-06 23:04  共线方程2.m

     文件       1168  2013-05-06 21:10  练习2.txt

----------- ---------  ---------- -----  ----

               685343                    4