实验一、复变函数与特殊函数图形的绘制.docVIP

实验一、复变函数与特殊函数图形的绘制 一、复变函数图形的绘制 例题：编程绘制出复变函数的图形。 解： %experiment1.m close all clear all m=30; r=(0:m)/m; theta=pi*(-m:m)/m; z=r*exp(i*theta); w=z.^3; blue=0.2; x=real(z); y=imag(z); u=real(w); v=imag(w); v=v/max(max(abs(v))); %%函数值虚部归一化 M=max(max(u)); m=min(min(u)); axis([-1 1 -1 1 m M]) caxis([-1 1]) %%指定颜色值的范围 s=ones(size(z)); subplot(131) mesh(x,y,m*s,blue*s) %%画投影图 hold on surf(x,y,u,v) %%画表面图 xlabel(x) ylabel(y) zlabel(u) title(z^3) hold off colormap(hsv(64)) %%画色轴 w=z.^(1/3); x=real(z); y=imag(z); subplot(132) for k=0:2 rho=abs(w); phi=angle(w)+k*2*pi/3; u=rho.*cos(phi); v=rho.*sin(phi); v=v/max(max(abs(v))); %%函数值虚部归一化 M=max(max(max(M,u))); m=min(min(min(m,u))); surf(x,y,u,v) %%画表面图 axis([-1 1 -1 1 m M]) hold on end s=ones(size(z)); mesh(x,y,m*s,blue*s) %%画投影图 xlabel(x) ylabel(y) zlabel(u) title(z^{1/3}) colormap(hsv(64)) %%画色轴 w=1./z; w(z==0)=NaN; x=real(z); y=imag(z); u=real(w); v=imag(w); v=v/max(max(abs(v))); %%函数值虚部归一化 M=max(max(max(M,u))); m=min(min(min(m,u))); subplot(133) surf(x,y,u,v) %%画表面图 hold on axis([-1 1 -1 1 m M]) s=ones(size(z)); mesh(x,y,m*s,blue*s) %%画投影图 xlabel(x) ylabel(y) zlabel(u) title(1/z) colormap(hsv(64)) %%画色轴 二、特殊函数图形的绘制 1、函数的绘制 % Fig1d15.m x=-3:0.01:3; y=gamma(x); plot(x,y,linewidth,4) grid on axis([-3 3 -5 5]) xlabel(x) ylabel(y) title(\Gamma函数) 2、勒让德函数的绘制 l阶勒让德多项式的定义是： 其中， 连带勒让德函数的定义是 其中，，，而是的m阶导数。 MATLAB计算连带勒让德函数的指令是 在给定N，X值以后，它将计算所有N阶连带勒让德函数在X处的函数值。如果X是矢量，所得的结果P是矩阵，而P（m+1,i）则是连带勒让德函数在X（i）处的函数值。 例如，在指令窗中输入 legendre(2,0:0.1:0.2) 产生的矩阵是： -0.5000 -0.4850 -0.4400 0 -0.2985 -0.5879 3.0000 2.9700 2.8800 它表示的结果是 例题：画出所有3阶连带勒让德函数的图形。 解： % Fig1d17.m x=0:0.01:1; y=legendre(3,x); plot(x,y(1,:),-,x,y(2,:),-.,x,y(3,:),:,x,y(4,:),--) title(勒让德多项式) legend(P_3^0,P_3^1,P_3^2,P_3^3) 运行程序，得到如下的图形： 3、贝塞尔函数的绘制 MATLAB有5种计算贝塞尔函数的指令， 计算指令 所计算的函数 J=besselj((,z)