对项目的规格感到困惑

时间:2015-10-19 21:16:58

标签: matlab gaussian

我对项目的内容感到困惑。它说如下: enter image description here

我做了#1,但在#2我需要在代码中插入Aij = cos(...)吗? 我有代码,但不知道我应该在哪里插入给定的条目Aij。任何帮助将不胜感激!谢谢!

Code

 n=8;
    for i=1:n
      for j=1:n
        a(i,j)=1/(i+j-1);
      end
    end
    c=ones(n,1);
    b=a*c;
    aorig=a;borig=b;
    %Elimination
    for i=1:n-1
      for j=i+1:n
        m=a(j,i)/a(i,i);
        for k=i:n
          a(j,k)=a(j,k)-m*a(i,k);
        end
        b(j)=b(j)-m*b(i);
      end
    end
    %Back substitution
    for i=n:-1:1
      for j=i+1:n
        b(i)=b(i)-a(i,j)*x(j);
      end
      x(i)=b(i)/a(i,i);
    end
    x'

1 个答案:

答案 0 :(得分:1)

用#2中的余弦表达式替换希尔伯特矩阵定义a(i,j)=1/(i+j-1)

相关问题