倍频程不能在循环中填充矩阵

Question

我是Octave的基础，有人能帮我写代码吗？

A = zeros(40);
format long
function xn=xnf(n,xn,c)
    if n<=40
        xn = xn^2+c;
        A(n,n)=xn;
        n++;
        xnf(n,xn,c);
    end
end
xnf(1,2,-2)

我想将xn值添加到矩阵中，但它总是只有零值-为什么？

Answer 1

对于你指定的值(n=1，xn=2，c=-2)，你将得到的矩阵是沿着主对角线的2，其他地方的0，这可以用A=diag(2*ones(40,1))来实现。

如果您想要更改参数，此代码将在Matlab中工作，因此应该在八度下工作：

N=40;
V=zeros(1,N);
V(1)=xn;
for i=2:N
    V(i)=V(i-1)^2+c;
end
A=diag(V);

Answer 2

我对Octave知之甚少，但在matlab中，如果我想要将xn的值添加到矩阵A中，并且我想要的是矩阵，那么我会这样做

function A=xnf(n,xn,c,A)
   if n<=40
       xn = xn^2+c;
       A(n,n)=xn;
       n=n+1;
       A=xnf(n,xn,c,A);
   end
end

因为函数中使用的所有变量在返回后都会被删除，所以如果需要A矩阵，请确保返回它