exn525a.m

来自「MATLAB及在电子信息课程中的应用.感觉很有用,本来也不花钱,有想用的来拿吧」· M 代码 · 共 17 行

%《MATLAB及其在理工课程中的应用指南》(十一五规划版)第五章例5-2-5程序a
% 空间任意曲面的交线
% 西安电子科技大学出版社出版  陈怀琛编著 2007年7月
%
clear,clf
s1=input('输入方程1 ','s');		% z1方程
s2=input('输入方程2 ','s');		% z2方程
xmax=input('xmax=ymax= ');dx=input('dx=dy= ');
[x,y]=meshgrid(-xmax:dx:xmax);
z1 = eval(s1);
z2 = eval(s2);
mesh(x,y,z1);hold on;mesh(x,y,z2);	%分别画出两个曲面（保持坐标比例不变）
r0=abs(z1-z2)<=.1;	%求两曲面z坐标差小于0.1的网格矩阵
zz=r0.*z1;yy=r0.*y;xx=r0.*x; 	%求这些网格上的坐标值,即交线坐标值
plot3(xx,yy,zz,'x');	%画出这些点
colormap(gray),hold off
set(gcf,'color','w')