迭代和子图

Question

我想每次绘制 5 个具有不同 x 值的子图（1 到 5）。而不是重复代码 5 次（每次更改 x）我认为使用 for 循环来绘制每个子图会更优雅。但是，我无法使代码正常工作。

进一步说明，当我手动绘制五个图时，由于某种原因，轴标签仅出现在 Y 轴上。如何让它们出现在所有轴上？

% constants
a=30;
b=0.5;

% parameters
x=1:1:5 ;
e1 = 0:10:10000; 
e2 = 0:10:10000; 

[e1,e2]=meshgrid(e1,e2);

% plots

t(x)=(a.^(6.*x)/factorial(6.*x))*...
    (exp(-b*x.*e2));

u(x)=((a.^(a-(6.*x)))/factorial(a-(6.*x)))*...
    exp(-b*(a-(6.*x)).*e1);

p(x)=t(x)./(t(x)+u(x));

%FIGURES:

for i=length(x)

    figure
    subplot(3,2,i);
    mesh(e1,e2,p(i));
    title('X=i');
    axis([1 10000 1 1000 10^-300 1])
    xlabel('e1');
    ylabel('e2');
    zlabel('p');

end

Answer 1

好的，您的代码有几个错误。这是一个有效的方法：

% constants
a=30;
b=0.5;

% parameters
x=1:1:5 ;
e1 = 0:10:10000; 
e2 = 0:10:10000; 

[e1,e2]=meshgrid(e1,e2);

% plots

t=@(x)(a.^(6.*x)./factorial(6.*x))*...
    (exp(-b*x.*e2));

u=@(x)((a.^(a-(6.*x)))/factorial(a-(6.*x)))*...
    exp(-b*(a-(6.*x)).*e1);

p=@(x)t(x)./(t(x)+u(x));

%FIGURES:
figure
for i=x


    subplot(3,2,i);
    mesh(e1,e2,p(i));
    title(['X=',int2str(i)]);
    axis([1 10000 1 1000 10^-300 1])
    xlabel('e1');
    ylabel('e2');
    zlabel('p');

end

问题是：

1. 匿名函数由t=@(x)([function here])而不是t(x)=[function here]定义
2. 您需要遍历 x 的所有元素，因此 for 循环需要为 for i=x
3. figure 语句需要在 for 循环之外，否则每次迭代都会有一个新的数字。
4. title() 不解析字符串中的变量。

希望有所帮助。