南华大学MATLAB实验报告4 - 图文

方法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

六、实验总结