www.gusucode.com > elfun18工具箱matlab源码程序 > elfun18/elfun18v1_3/examples/plot/plotLambdaZeta.m

    fa = 0.75;

figure(1)
set(gcf, 'Position', [488 298 560 560]);
clf
hold on
x = -1:0.005:1;
k = [0, 0.25, 0.5, 0.75, 1];
for n = 1:length(k)
    clg{n} = num2str(k(n));
    plot(x,pHeumanLambda0(pi*x,k(n)),'LineWidth',2)
end
hlg = legend(clg,'Location','southeast');
title(hlg, 'k','FontWeight','normal')
title('Heuman \Lambda_0')
xlabel('\phi/\pi')
ylabel('\Lambda_0(\phi,k)');
%ylim([-0.5,1])
grid on
hold off

figure(2)
set(gcf, 'Position', [488 298 560 560]);
clf
hold on
x = -1:0.005:1;
k = [0, 0.25, 0.5, 0.75, 1];
for n = 1:length(k)
    clg{n} = num2str(k(n));
    plot(x,pJacobiZeta(pi*x,k(n)),'LineWidth',2)
end
hlg = legend(clg,'Location','southeast');
title(hlg, 'k','FontWeight','normal')
title('Jacobi Z')
xlabel('\phi/\pi')
ylabel('Z(\phi,k)');
%ylim([-0.5,1])
grid on
hold off

figure(3)
set(gcf, 'Position', [488 298 560 560]);
clf
hold on
title('Jacobi Z(\phi,k)')
xlabel('\phi/\pi')
ylabel('k')
view(3)
x = -1:0.01:1;
k = -1:0.01:1;
[X,K] = meshgrid(x,k);
hs = surfc(X,K,pJacobiZeta(pi*X,K),'EdgeColor','none','FaceAlpha',fa);
hc = hs(2);
hc.ContourZLevel = -1;
%hc.LineWidth = 1;
hc.LevelList = -1:0.1:1;
caxis([-1 1])
%zlim([-1.1,1.1])
grid on
hold off

figure(4)
set(gcf, 'Position', [488 298 560 560]);
clf
hold on
title('Jacobi Z(x,k)')
xlabel('x')
ylabel('k')
view(3)
x = -2:0.01:2;
k = -1:0.01:1;
[X,K] = meshgrid(x,k);
hs = surfc(X,K,JacobiZeta(X,K),'EdgeColor','none','FaceAlpha',fa);
hc = hs(2);
hc.ContourZLevel = -1;
%hc.LineWidth = 1;
hc.LevelList = -1:0.1:1;
caxis([-1 1])
%zlim([-1.1,1.1])
grid on
hold off