方法2： x=-3:.01:3; y=x;
ezmesh('4*x.*exp(-x.^2-y.^2)') hidden off
3.（1）t=0:.01:2*pi;
r=3*(1-cos(t)); polar(t,r,'r')
（2）t=0:.01:2*pi; r=2*(1+cos(t));
polar(t,r,'r')（3）t=0:.01:2*pi; r=2*(1+sin(t));
polar(t,r,'r')
（4）t=0:.01:2*pi; r=cos(3*t); polar(t,r,'r')
（5）t=0:.01:2*pi; r=exp(4*pi*t);
polar(t,r,'r')
4.y1=@(x)(0.2).*exp((-0.5).*x).*cos((4*pi).*x); y2=@(x)2*exp(-0.5*x)*cos(pi*x); hold on
fplot(@(x)y1(x),[0 5],'b','-') fplot(@(x)y2(x),[0 5],'r','--') X=[];
for x0=[0.47 1.51 2.49 3.53 4.45]
X(end+1)=fsolve(@(x)y1(x)-y2(x),x0); end Y=y1(X);
plot(X,Y,'.','Markersize',30); for i=1:length(X)
text(X(i),Y(i),[' (',num2str(X(i)),',',num2str(Y(i)),')'])
end
六、实验总结
