exm0563_1.m

matlab程序

%exm0563_1.m
clear
%（0）从三维坐标初步观察两函数图形相交情况
x=-2:0.05:2;y=x;[X,Y]=meshgrid(x,y);	%产生x-y平面上网点坐标
F1=sin(X-Y);F2=cos(X+Y);
F0=zeros(size(X));
figure(1)
surf(X,Y,F1),
xlabel('x'),ylabel('y'),
view([-31,62]),hold on,
surf(X,Y,F2),surf(X,Y,F0),
shading interp,
hold off  


%（1）在某区域观察两函数0等位线的交点情况
clear;
x=-2:0.5:2;y=x;[X,Y]=meshgrid(x,y);	%产生x-y平面上网点坐标
F1=sin(X-Y);F2=cos(X+Y);
v=[-0.2, 0, 0.2];	%指定三个等位值，是为了更可靠地判断0等位线的存在。
figure(2)
contour(X,Y,F1,v)						%画F1的三条等位线。
hold on,contour(X,Y,F2,v),hold off	%画F2的三条等位线。  


%（2）从图形获取零点的初始近似值
clc
disp('用鼠标在图上获取2个零点初始数据')
[x0,y0]=ginput(2);  					%在图上取两个点的坐标
disp([x0,y0])  


%（3）利用fsolve求精确解。以求（0.7926,7843）附近的解为例。
fun='[sin(x(1)-x(2)),cos(x(1)+x(2))]';								%<12>
xy=fsolve(fun,[x0(2),y0(2)])											%<13>  


%（4）检验
fxy1=sin(xy(1)-xy(2));
fxy2=cos(xy(1)+xy(2));
disp([fxy1,fxy2])